用Matlab模拟SSB调制波形及解调过程

用Matlab产生一个频率为1Hz、功率为1的余弦信源m(t)，设载波频率为10Hz，试画出（1）SSB调制信号；（2）该调制信号的功率谱密度；（3）相干解调后的信号波形。Matlab代码：SSB.m%显示模拟调制的波形及解调方法SSB%信源closeall;clearall;dt=0.001;%采样时间间隔fm=1;%信源最高频率fc=10;%载波中心频率N=4096;T=(N-1)*dt;%信号时长t=0:dt:T;mt=sqrt(2)*cos(2*pi*fm*t);%信源s_ssb=real(hilbert(mt).*exp(j*2*pi*fc*t));%SSB单边带抑制载波调幅B=2*fm;figure(1);subplot(311);plot(t,s_ssb,'k-');holdon;%画出SSB信号波形plot(t,mt,'g--');%画出m(t)信号波形title('SSB调制信号');xlabel('t');gridon;legend('ssb信号波形','m(t)信号波形');%SSBdemodulationrt=s_ssb.*cos(2*pi*fc*t);rt=rt-mean(rt);[f,rf]=T2F(t,rt);[t,rt]=lpf(f,rf,B);subplot(312);plot(t,rt,'b-');holdon;plot(t,mt/2,'r--');title('相干解调后的信号波形与输入信号的比较');xlabel('t');gridon;subplot(313);[f,sf]=T2F(t,s_ssb);%单边带信号频谱psf=(abs(sf).^2)/T;%单边带信号功率谱plot(f,psf);axis([-2*fc2*fc0max(psf)]);title('SSB信号功率谱');xlabel('f');gridon;F2T.mfunction[t,st]=F2T(f,sf)df=f(2)-f(1);Fmx=(f(end)-f(1)+df);dt=1/Fmx;N=length(sf);T=dt*N;%t=-T/2:dt:T/2-dt;t=0:dt:T-dt;sff=fftshift(sf);st=Fmx*ifft(sff);Lpf.mfunction[t,st]=lpf(f,sf,B)df=f(2)-f(1);T=1/df;hf=zeros(1,length(f));bf=[-floor(B/df/2):floor(B/df/2)]+floor(length(f)/2);hf(bf)=1;yf=hf.*sf;[t,st]=F2T(f,yf);st=real(st);T2F.mfunction[f,sf]=T2F(t,st)dt=t(2)-t(1);T=t(end);df=1/T;N=length(st);f=-N/2*df:df:N/2*df-df;sf=fft(st);%sf=T/N*sf;sf=T/N*fftshift(sf);SSB.m、F2T.m、Lpf.m、T2F.m四个文件要在一个文件夹里才能出来正确的运行结果运行结果：00.511.522.533.544.5-202SSB调制信号tssb信号波形m(t)信号波形00.511.522.533.544.5-101相干解调后的信号波形与输入信号的比较t相干解调后的信号输入信号-20-15-10-505101520012SSB信号功率谱f