SSB调制解调(matlab)

dt=0.001; %时间采样间隔 fm=1; %信源最高频率 fc=10; %载波中心频率 T=5; %信号时长 t=0:dt:T; mt=cos(2*pi*fm*t); %信源 ct=cos(2*pi*fc*t); %载波信号 figure(1) plot(t,mt); title('调制信号时域波形') axis([0 5 -1.01 1.01]); line([0,5],[0,0],'color','k'); figure(2) [fmt,Fmt]=T2F(t,mt); plot(fmt,Fmt); title('调制信号频域波形') axis([-4 4 0 2.6]); figure(3) plot(t,ct); title('载波信号时域波形'); axis([0 5 -1.01 1.01]); line([0,5],[0,0],'color','k'); figure(4) [fct,Fct]=T2F(t,ct); plot(fct,Fct); title('调制信号频域波形') axis([-12 12 0 2.8]); %SSB调制 s=modulate(mt,fc,1000,'amssb'); %已调信号 figure(5) plot(t,s); %画出DSB信号波形 title('SSB调制信号'); axis([0 5 -1.2 1.2]); line([0,5],[0,0],'color','k'); figure(6) [frt,Frt]=T2F(t,s); plot(frt,Frt); title('SSB调制信号频谱'); axis([-14 14 0 2.8]); figure(7) psf=(abs(Frt).^2)/T; plot(frt,psf); axis([-20 20 0 max(psf)]); title('SSB信号功率谱'); sssb1=awgn(s,5); sssb2=awgn(s,10); sssb3=awgn(s,15); figure(8) subplot(311) plot(t,sssb1); %画出DSB信号波形 title('SSB调制信号 s/N=5dB'); axis([0 5 -2.4 2.4]); line([0,5],[0,0],'color','k'); subplot(312) plot(t,sssb2); %画出DSB信号波形 title('SSB调制信号 s/N=10dB'); axis([0 5 -2.4 2.4]); line([0,5],[0,0],'color','k'); subplot(313) plot(t,sssb3); %画出DSB信号波形 title('SSB调制信号 s/N=15dB'); axis([0 5 -2.4 2.4]); line([0,5],[0,0],'color','k'); %DSB相干解调 rt1=sssb1.*ct; rt2=sssb2.*ct; rt3=sssb3.*ct; figure(9); subplot(311) plot(t,rt1); hold on plot (t,mt,'r-'); %画出原信号 title('SSB调制信号相干解调与调制信号 S/N=5dB') line([0,5],[0,0],'color','k'); subplot(312) plot(t,rt2); hold on plot (t,mt,'r-'); %画出原信号 title('SSB调制信号相干解调与调制信号 S/N=10dB') line([0,5],[0,0],'color','k'); subplot(313) plot(t,rt3); hold on plot (t,mt,'r-'); %画出原信号 title('SSB调制信号相干解调与调制信号 S/N=15dB') line([0,5],[0,0],'color','k'); [f1,rf1]=T2F(t,rt1); %傅里叶变换 [t,rt1]=lpf(f1,rf1,fm); %低通滤波 [f2,rf2]=T2F(t,rt2); %傅里叶变换 [t,rt2]=lpf(f2,rf2,fm); %低通滤波 [f3,rf3]=T2F(t,rt3); %傅里叶变换 [t,rt3]=lpf(f3,rf3,fm); %低通滤波 figure(10) subplot(311) plot(t,rt1); line([0,5],[0,0],'color','k'); title('经过低通滤波的相干解调信号波形 S/N=5dB'); subplot(312) plot(t,rt2); line([0,5],[0,0],'color','k'); title('经过低通滤波的相干解调信号波形 S/N=10dB'); subplot(313) plot(t,rt3); line([0,5],[0,0],'color','k'); title('经过低通滤波的相干解调信号波形 S/N=15dB');