通信应用系统实验报告.docx

通信应用系统实验报告学院：计算机科学与工程专业：通信工程组员：1241903105 靳晓春1241903106 李学红 1241903107 钮佳超1241903108 王 娟题目一AM、DSB、SSB调制题目1. 设信源为 m(t)=2 cos 2?t ，载波为 s (t ) ? A cos 20? tA ? 2 ，AM 调制信号A ? 1 ，DSB 调制信号A ? 1 ，SSB 调制信号信道中各自加入进过带通滤波器后的窄带高斯白噪声，功率为 0.1，解调 每个信号，并画出滤除前后，双边带信号和单边带信号的频谱。设计流程图源程序代码ts = 0.01;t0 = 2;t =-t0:ts:t0; A = 2;fa = 1;N = length(t); w = 0:2*pi/(N-1):2*pi; %载波和调制信号波形%mt = sqrt(2)*cos(2*pi*fa.*t);st = A*cos(20*pi*t); %AM信号调制程序%figure(1);amt = (mt).*st;subplot(321);plot(t,amt);xlabel(AM时域波形图);subplot(322);plot(w,fft(amt));axis([0,pi,0,500]);xlabel(AM频域波形图); %DSB信号调制程序%dsbt = mt.*st;subplot(323);plot(t,dsbt);xlabel(DSB时域波形图);subplot(324);plot(w,fft(dsbt));axis([0,pi,0,500]);xlabel(DSB频域波形图); %SSB信号调制程序%ssbt = real(hilbert(mt).*exp(j*20*pi*t));subplot(313);plot(t,ssbt);hold on;plot(t,mt,r--);xlabel(t);axis([0 5 -2 2]);subplot(325);plot(t,ssbt);xlabel(SSB时域波形图);subplot(326);plot(w,fft(ssbt));axis([0 pi 0 200]);xlabel(SSB频域波形图); %在信道中加入窄带高斯白噪声%figure(2);fangcha = 0.1;[h,l] = size(mt);nit = sqrt(fangcha).*randn(h,l);amtn = amt + nit;subplot(321);plot(t,amtn);xlabel(白噪声AM时域波形图);subplot(322);plot(w,fft(amtn));axis([0,pi,0,500]);xlabel(白噪声AM频域波形图);dsbtn = dsbt + nit;subplot(323);plot(t,dsbtn);xlabel(白噪声DSB时域波形图);subplot(324);plot(w,fft(dsbtn));axis([0,pi,0,500]);xlabel(白噪声DSB频域波形图);ssbtn = ssbt + nit;subplot(325);plot(t,ssbtn);xlabel(白噪声SSB时域波形图);subplot(326);plot(w,fft(ssbtn));axis([0,pi,0,500]);xlabel(白噪声SSB频域波形图); %已调信号的解调%figure(3);amtm = amtn.*st;amtm = (amtm-mean(amtm))*2;subplot(311);plot(t,amtm);xlabel(AM解调);dsbtm = (dsbtn.*st)*2;subplot(312);plot(t,dsbtm);xlabel(DSB解调);ssbtm = (ssbtn.*st)*2;subplot(313);plot(t,ssbtm);xlabel(SSB解调); %解调信号通过低通滤波器%Rp = 0.1;Rs = 50;Wp = 5/100;Ws = 10/100; %Rp、Rs分别为通带最大波纹和阻带最小衰减% %Wp、Ws分别为通带边界频率和阻带边界频率% [n,Wn] = ellipord(Wp,Ws,Rp,Rs); [b,a] = ellip(n,Rp,Rs,Wn);figure(4);x1 = filter(b,a,amtm);subplot(311);plot(t,x1);xlabel(AM通过低通滤波器);x2 = filter(b,a,dsbtm);subplot(312);plot(t,x2);xlabel(DSB通过低通滤波器);x3 = filter(b,a,ssbtm); subplot(313);plot(t,x3)