资源描述:
《0144366陈官亮实验3》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、实验课程名称:通信系统实验与设计学院:软件与通信工程学院实验班级专业:通信工程指导教师:奎刚学生姓名:陈官亮学号:01443662017年4月13日实验报告正文学生实验报告正文一、实验内容:消息信号为t,0<1=<—t+2,Kf<20,othert载波信号为频率为25Hz,初始相位为0,幅度为1的余弦信号。1)、用上述信号分别进行①DSB②LSSB③常规AM(调制指数a=0.2)④FM(kf=50)调制。用Matlab编程求出对应的消息信号功率,己调信号功率,并绘制消息信号及已调信号的时域波形和频谱。2)、设上述各已调信号经过AWGN信道,解调器
2、输入信噪比为100,用Matlab对上述各已调信号进行解调,绘制加噪声后的接收信号时间波形、频谱,解调信号的时域波形并和m⑴对比。二、实验仪器及设备:电脑以及matlab软件三、实验结果如下:程序1,DSB%MatlabdemonstrationscriptforDSB-AMmodulation.Themessagesignalt=[0:ts:tO];snr_lin=l(T(snr/10);*(t>=O&t=l&t<2):c=cos(2本pi本fc.*t):u=m•氺c;[M,m,dfl]=fftseq(m,ts,df);
3、M=M/fs;%timevector%linearSNR%messagesignal%carriersignal%modulatedsignal%Fouriertransform%scaling[U,u,dfl]=fftseq(u,ts,df)%Fouriertransformechoont0=2;%signaldurationts=0.001:%samplingintervalfc=25;%carrierfrequencysnr=20;%SNRindB(logarithmic)fs=l/ts;%samplingfrequencydf=O.3;%de
4、siredfreq.resolutioU=U/fs;[C,c,dfl]=fftseq(c,ts,df);f=[0:df1:dfl*(length(m)-l)]-fs/2;message_power=spower(m(l:length(t)));1signalpower^spower(u(1:1ength⑴));alnoise_power=signal_power/snr_lin;noisc_std=sqrt(noisc_powcr):deviationnoise=noisestd^randn(1,length(u)):r=u+noise;tedsig
5、nal[R,r,dfl]=fftseq(r,ts,df):noise%scaling%Fouriertransform%freq.vector%powerinmessagesigna%powerinmodulatedsign%computenoisepower%computenoisestandard%generatenoise%addnoisetothemodula%spectrumofthesignal+R=R/fs;%scalingpause%Pressakeytoshowthemodulatedmessagepowermessage_powe
6、rpause%Pressakeytoshowthemodulatedsignalpowersignal_powerpause%Pressanykeytoseeaplotofthemessageelfplot(t,m(l:length(t)))xlabel(’Time’)title(’Themessagesignal’)pause%Pressanykeytoseeaplotofthecarriersubplot(2,2,2)plot(t,c(1:length(t)))xlabelCTime’)title(’Thecarrier’)pause%Press
7、anykeytoseeaplotofthemodulatedsignalsubplot(2,2,3)plot(t,u(l:length(t)))xlabelCTime’)title(’Themodulatedsignal’)pause%Pressanykeytoseeaplotsofthemagnitudeofthemessageandthe%modulatedsignalinthefrequencydomain,subplot(2,1,1)plot(f,abs(fftshift(M)))xlabelCFrequency’)title(’Spectr
8、umofthemessagesignal’)axis([-5050,01]);subplot(2,1,2)p