杏林学院通信原理课程设计(注释版)

《杏林学院通信原理课程设计(注释版)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、课程设计报告1、常规AM调制基带信号1，0≤t

2、已调信号加噪声（u(t)+n(t)）的频谱图；（7）比较示范课题中的m(t)和u(t)的波形及其频谱，写出它们的不同点。说明：（1）注意到，因此有（2）本题求已调信号功率时不能用下面关系因为本题中m(t)不是零均值的信号。程序如下：%MatlabdemonstrationscriptforDSBmodulation.%Themessagesignalis+1for0

3、0.001;时间间隔fc=250;载波频率snr=20;信噪比为20DBfs=1/ts;抽样频率df=0.3;采样间隔t=(0:ts:t0);时间范围snr_lin=10^(snr/20);信噪比单位转换%messagesignalm=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)];m的表达式；ones是全1矩阵；zeros是全0矩阵c=cos(2*pi*fc.*t);载波表达式u=(2+0.85*m).*c;已调信号表达式%To

4、specifythepowerofXtobe0dBWandaddnoisetoproduce%anSNRof10dB,use:%X=sqrt(2)*sin(0:pi/8:6*pi);awgn噪声函数调用n=awgn(u,10,0)-u;噪声表达式a=u+n;加噪声后已调信号表达式[M,m,df1]=fftseq(m,ts,df);调用F变换求m的频谱M=M/fs;按fs比例缩放[U,u,df1]=fftseq(u,ts,df);调用F变换求u的频谱U=U/fs;按fs比例缩放[C,c,df1]=fftse

5、q(c,ts,df);调用F变换求c的频谱f=[0:df1:df1*(length(m)-1)]-fs/2;[N,n,df1]=fftseq(n,ts,df);调用F变换求n的频谱N=N/fs;按fs比例缩放[A,a,df1]=fftseq(a,ts,df);调用F变换求a的频谱A=A/fs;按fs比例缩放signal_power=spower(u(1:length(t)));信号功率计算noise_power=signal_power/snr_lin;噪声功率计算'显示信号功率和噪声功率'pause按一

6、个键显示signal_power信号功率noise_power噪声功率'显示基带信号波形图'pause按一个键看图clf清空当前figurefigure(1)plot(t,m(1:length(t)))以t为横坐标，m为纵坐标画图xlabel('Time')横坐标为'Time'title('Themessagesignal')标题为'Themessagesignal''显示载波信号和调制信号波形图'pause按一个键看图figure(2)subplot(2,1,1)把多个图画到同一个平面上，位置为1行1列

7、plot(t,c(1:length(t)))以t为横坐标，c为纵坐标画图xlabel('Time')横坐标为'Time'title('Thecarrier')标题为'Thecarrier'subplot(2,1,2)把多个图画到同一个平面上，位置为1行2列plot(t,u(1:length(t)))以t为横坐标，u为纵坐标画图xlabel('Time')横坐标为'Time'title('Themodulatedsignal')标题为'Themodulatedsignal''显示基带信号和调制信号的频谱图'

8、pausefigure(3)subplot(2,1,1)把多个图画到同一个平面上，位置为1行1列plot(f,abs(fftshift(M)))以f为横坐标，M为纵坐标画图xlabel('frequency')横坐标为'frequency'title('spectrumofthemessagesignal')标题为'spectrumofthemessagesignal'subplot(2,1,2)plot(f,abs(ffts