北邮通信原理软件实验报告北邮通信原理软件实验报告通信原理软实验实验报告学院:信息与通信工程学院班级:@@@@@@@@@@姓名:execf2t.sci;fs=800;//采样速率T=200;//截短时间N=T*fs;//采样点数dt=1/fs;//时域采样间隔t=[-T/2:dt:T/2-dt];//时域采样点df=1/T;//频域采样间隔f=[-fs/2:df:fs/2-df];//频域采样点数fm1=1;//待观测正弦波频率,单位KHz,下同fm2=0.5; //待观测余弦波频率fc=20;//载波频率//以上为初始化参数设置m1=sin((2*%pi)*fm1*t);//待观测正弦波部分M1=t2f(m1,fs);//傅里叶变换MH1=-%i*sign(f).*M1;//希尔伯特变换mh1=real(f2t(MH1,fs));//希尔伯特反变换m2=2*cos((2*%pi)*fm2*t);//待观测余弦波部分M2=t2f(m2,fs);//傅里叶变换MH2=-%i*sign(f).*M2;//希尔伯特变换mh2=real(f2t(MH2,fs));//希尔伯特反变换s1=(1+(m1+m2)/abs(max(m1+m2))).*cos((2*%pi)*fc*t); //AM信号时域表达式S1=t2f(s1,fs);//AM信号频域表达式s2=(m1+m2).*cos((2*%pi)*fc*t);//DSB-SC信号时域表达式S2=t2f(s2,fs);//DSB-SC信号频域表达式s3=(m1+m2).*cos((2*%pi)*fc*t)-(mh1+mh2).*sin((2*%pi)*fc*t);//SSB信号时域表达式,以上边带为例S3=t2f(s3,fs);//SSB信号上边带频域表达式//以上是仿真计算部分//以下为绘图部分//AM信号xset(window,1)plot(f,abs(S1))title(AM信号频谱)xlabel(f)ylabel(S(f))mtlb_axis([-25,25,0,max(abs(S1))]); xset( window,2)plot(t,s1)title(AM信号波形)xlabel(t)ylabel(s(t))mtlb_axis([-3,3,-3,3]);//DSB-SC信号window,3)plot(f,abs(S2))title(DSB-SC信号频谱)xlabel(f)ylabel(S(f))mtlb_axis([-25,25,0,max(abs(S2))]); xset( window,4)plot(t,s2)title(DSB-SC信号波形)xlabel(t)s(t))mtlb_axis([-1,4,-3,3]);//SSB信号(以上边带为例)xset(window,5)plot(f,abs(S3))title(SSB信号频谱)xlabel(f)ylabel(S(f))mtlb_axis([-25,25,0,max(abs(S3))]) xset(window,6)plot(t,s3)title(SSB信号波形)xlabel(t)ylabel(s(t))mtlb_axis([0,6,-3,3]) 2、产生的波形图:1)AM信号:2)AM信号频谱:3)DSB-SC信号:4)DSB-SC信号的频谱:5)SSB信号波形:6)SSB信号频谱:3、实验心得:做这个实验时的我的理论知识已经准备得比较充分了,所以难点主要在编程方面。