SSB-SC通信系统仿真..doc

SSB-SC通信系统MATLAB仿真 SSB（单边带）信号的产生： 设单频调制信号为：m（t）=AmcosΩmt 载波为：c（t）=cosΩct 则DSB信号的时域表示式为：SDSB= AmcosΩmt cosΩct= 1/2Amcos (Ωc+Ωm)t+ 1/2Amcos (Ωc-Ωm)t 则SSB信号的时域表示式为：SSSB=1/2AmcosΩmt cosΩct- 1/2AmsinΩmtsinΩct SSSB=1/2AmcosΩmt cosΩct+1/2AmsinΩmtsinΩct 其中：“-”表示上边带信号；“+”表示下边带信号。 SSB（单边带）信号的解调 因SSB信号是抑制载波的已调信号，其包络不直接表示调制信号，因而在这里我们采用相干解调（同步检波）。相干解调时，为了无失真地恢复原基带信号，接收端必须提供一个与接收的已调载波严格同步（同频通相）的本地载波（相干波），它与接收的已调信号相乘后，经低通滤波器取出低频分量，即可得到原始的基带信号。 送入解调器的已调信号的一般表达式为： Sm（t）=SI（t）cosΩct+SQ（t）sinΩct 与同频同相的相干载波c（t）相乘后，得 SP（t）= Sm（t）cosΩct= 1/2 SI（t）+1/2 SI（t）cos2Ωct+1/2 SQ（t）sin2Ωct 经低通滤波器（LPF）后，得 Sd（t）=1/2 SI（t）∝m（t） 即可得到基带信号：m（t） MATLAB程序仿真 ts=1.e-4;%时间采样间隔 t=-0.04:ts:0.04; Ta=0.01; fc=300; m_sig=triangl((t+0.01)/0.01)-triangl((t-0.01)/0.01);%产生基带信号 Lm_sig=length(m_sig);%采样点数 Lfft=length(t);Lfft=2^ceil(log2(Lfft));% M_fre=fftshift(fft(m_sig,Lfft));%对基带信号进行傅里叶变换由时域变到频域，并让正半轴部分和负半轴部分的图像分别关于各自的中心对称，因为直接用fft得出的数据与频率不是对应的，fftshift可以纠正过来 freqm=(-Lfft/2:Lfft/2-1)/(Lfft*ts);%由时域到频域时每个点所对应的频率 B_m=150;%基带信号的带宽 h=fir1(40,[B_m*ts])%返回6dB截止频率为B_m*ts的40阶（单位脉冲响应h长度N=40+1）FIR低通（wc为标量）滤波器系数向量h，默认选用哈明窗。; s_dsb=m_sig.*cos(2*pi*fc*t);%双边带信号的产生 Lfft=length(t);Lfft=2^ceil(log2(Lfft)+1); S_dsb=fftshift(fft(s_dsb,Lfft));%对双边带信号进行傅里叶变换 L_lsb=floor(fc*ts*Lfft);%以载波频率为基准向下取整 SSBfilt=ones(1,Lfft); SSBfilt(Lfft/2-L_lsb+1:Lfft/2+L_lsb)=zeros(1,2*L_lsb); %理想低通滤波器 S_ssb=S_dsb.*SSBfilt;%滤波法得到下边带 freqs=(-Lfft/2:Lfft/2-1)/(Lfft*ts); %由时域到频域时每个点所对应的频率 s_ssb=real(ifft(fftshift(S_ssb)));%取下边带信号实部 s_ssb=s_ssb(1:Lm_sig); s_dem=s_ssb.*cos(2*pi*fc*t)*2;%相干解调 S_dem=fftshift(fft(s_dem,Lfft));%解调后的信号由时域变换到频域 s_rec=filter(h,1,s_dem);%输入s_dem为滤波前序列，s_rec为滤波结果序列，h/1 提供滤波器系数 S_rec=fftshift(fft(s_rec,Lfft));%将接收到的信号由时域转换到频域 Trange=[-0.025 0.025 -1 1];%坐标轴的范围 figure(1) subplot(221);td1=plot(t,m_sig); axis(Trange);set(td1,Linewidth,1.5); xlabel({\it t} (sec));ylabel({\it m}({\it t})); title(message signal); subplot(222);td2=plot(t,s_ssb); axis(Trange);set(td2,Linewidth,1.5); xlabel({\it t}(sec)