VSB调制仿真及其抗噪声性能分析

VSB调制及解调仿真（包含三角函数及语音信号测试） 包括各模块代码讲解，时域波形，频谱分析以及信噪比增益曲线

function AM_simulation

%AM调制
fc=10;
f0=1;
A=2;
t0=5;%如果频率为1，那至少要信号时长5个周期
snr=15;
dt=0.003;
fs=1/dt;
B=2;
df=0.0003;
t=0:dt:t0;
Lt=length（t）;
snr_lin=10^（snr/10）;
figure（1）;
subplot（321）%调制信号波形
m1=cos（2*pi*f0*t）;
plot（tm1（1:length（t）））
R=2*max（m1）;
axis（[0 t0 -R/2 R/2]）
ylabel（‘基波信号‘）;
subplot（322）%调制信号频谱
[M1m1df1f]=T2F（m1dtdffs）;
plot（ffftshift（abs（M1）））;
axis（[-2 2 0 3]）;
xlabel（‘f‘）;
ylabel（‘调制信号频谱‘）;
subplot（323）%载波波形
m2=cos（2*pi*fc*t）;
plot（tm2（1:length（t）））;
ylabel（‘载波信号波形‘）;
axis（[0 0.4 -R/2 R/2]）
subplot（324）;%载波频谱
[M2m2df1f]=T2F（m2dtdffs）;
plot（ffftshift（abs（M2）））;
xlabel（‘f‘）;
ylabel（‘载波信号频谱‘）;
axis（[-20 20 0 3]）
subplot（325）%已调波波形
m=（m1（1:length（t））+A）.*m2（1:length（t））;
plot（tm（1:length（t）））;
R=2*max（m）+A;
axis（[0 3 -R/2 R/2]）
ylabel（‘已调信号波形‘）;
subplot（326）
[Mmdf1f]=T2F（mdtdffs）;
plot（ffftshift（abs（M）））;
ylabel（‘已调信号频谱‘）;
xlabel（‘f‘）;
axis（[-15 15 0 3]）
%将已调信号送入信道
signal_power=power_x（m（1:Lt））;
noise_power=（signal_power*fs）/（snr_lin*2*B）;
noise_std=sqrt（noise_power）;
noise=noise_std*randn（1Lt）;
figure（2）;
subplot（321）%噪声波形
plot（tnoise）
axis（[0 5 -10 10]）;
xlabel（‘t‘）;
ylabel（‘噪声信号‘）;
subplot（322）%噪声信号频谱
[noisefnoisedf1f]=T2F（noisedtdffs）;
plot（ffftshift（abs（noisef）））;
axis（[-50 50 0 1]）
xlabel（‘f‘）;
ylabel（‘噪声信号频谱‘）;
subplot（323）%信道中的信号
sam=m（1:Lt）+noise（1:Lt）;
plot（tsam）;
axis（[0 1 -R R]）
xlabel（‘t‘）;
ylabel（‘信道中的信号‘）;
subplot（324）%信道中信号的频谱
[samfsamdf1f]=T2F（samdtdffs）;
plot（ffftshift（abs（samf）））;
axis（[-17 17 0 6]）;
ylabel（‘信道中信号的频谱‘）;
xlabel（‘f‘）;
subplot（326）%带通滤波器
f_start=fc-B;
f_cutoff=fc+B;
[H f]=bp_f（length（sam）f_startf_cutoffdf1fsA）;
plot（ffftshift（abs（H）））;
axis（[-15 15 0 2.5 ]）
xlabel（‘f‘）;
ylabel（‘带通滤波器‘）;
%经过带通滤波器
DEM=H.*samf;
[dem]=F2T（DEMfs）;
figure（3）
subplot（321）
plot（tdem（1:Lt））
axis（[0 3 -5 5]）
xlabel（‘t‘）
ylabel（‘理想BPF输出信号‘）;
[demfdemdf1f]=T2F（dem（1:Lt）dtdffs）;
subplot（322）%经过理想带通滤波器后信号频谱
plot（ffftshift（abs（demf）））;
xlabel（‘f‘）;
ylabel（‘理想BPF输出信号频谱‘）;
axis（[-15 15 0 5]）;
% %与本地载波相乘，混频
% subplot（323）%本地载波
% plot（tm2（1:Lt））;
% axis（[0 0.4 -2 2]）
% xlabel（‘t‘）;
% ylabel（‘本地载波‘）;
% subplot（324）%本地载波频谱
% plot（ffftshift（abs（M2）））;
% xlabel（‘f‘）;
% ylabel（‘载波信号频谱‘）;
% axis（[-20 20 0 3]）
% der=dem（1:Lt）.*m2（1:Lt）;%混频
% subplot（325）;
% plot（tder）;
% xlabel（‘t‘）;
% ylabel（‘混频后的信号‘）;
% axis（[0 3 -1 R]）;
% subplot（326）%混频后信号频谱
% [derfderdf1f]=T2F（derdtdffs）;
% plot（ffftshift（abs（derf）））
% xlabel（‘f‘）;
% ylabel（‘混频后信号频谱‘）;
% axis（[-30 30 0 7]）
% %经过低通滤波器
% figure（4）
% [LPF f]=lp_f（length（der）Bdf1fs2）;
% subplot（322）
% plot（ffftshift（abs（LPF）））
% xlabel（‘f‘）
% ylabel（‘理想LPF‘）
% axis（[-3 3 0 3]）
% DM=LPF.*derf;
% [dm]=F2T（DMfs）;
% subplot（323）
% plot（tdm（1:Lt））;
% xlabel（‘t‘）
% ylabel（‘LPF输出信号‘）
% subplot（324）
% [dmfdmdf1f]=T2F（dm（1:Lt）dtdffs）;
% plot（ffftshift（dmf））;
% xlabel（‘f‘）
% ylabel（‘LPF输出信号频谱‘）
% ax

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3669  2017-11-09 13:20  VSB_simulation\AM_simulation.m

     文件        292  2017-11-13 14:12  VSB_simulation\bp_f.m

     文件        665  2017-11-20 18:26  VSB_simulation\coherent_demodulation_power_transf.m

     文件         42  2017-11-05 15:30  VSB_simulation\F2T.m

     文件        185  2017-11-05 15:29  VSB_simulation\fftseq.m

     文件        219  2017-11-13 14:12  VSB_simulation\lp_f.m

     文件       6493  2017-11-21 22:21  VSB_simulation\modulation_performance_comparing.m

     文件     128044  2017-11-18 09:48  VSB_simulation\output.wav

     文件         48  2017-11-09 08:58  VSB_simulation\power_x.m

     文件         96  2017-11-05 15:30  VSB_simulation\signalband.m

     文件     128078  2017-11-15 15:16  VSB_simulation\sound.wav

     文件        107  2017-11-05 17:08  VSB_simulation\T2F.m

     文件       4166  2017-11-17 12:44  VSB_simulation\VSB.m

     文件        340  2017-11-17 12:36  VSB_simulation\vsb_lpf.m

     文件       5713  2017-11-21 22:21  VSB_simulation\VSB_simulation.m

     目录          0  2017-11-20 19:52  VSB_simulation

----------- ---------  ---------- -----  ----

               278157                    16