. AM、DSB调制及解调用matlab产生一个频率为1Hz，功率为1的余弦信源，设载波频率，,试画出：l AM及DSB调制信号的时域波形；l 采用相干解调后的AM及DSB信号波形；l AM及DSB已调信号的功率谱；l 调整载波频率及m0，观察分的AM的过调与DSB反相点现象l 在接收端带通后加上窄带高斯噪声，单边功率谱密度，重新解调一 图表输出及结果分析1 AM及DSB调制信号的时域波形AM信号的包络波形与其调制信号的波形完全一样DSB调制信号在载波处反向2 采用相干解调后的AM及DSB信号波形AM信号可用包络检波器检波，滤除直流后可恢复出原信号与AM相比，DSB信号由于不存在载波分量，全部功率都用来传输信息3 AM及DSB已调信号的功率谱AM信号所需的传输带宽为调制信号带宽的两倍DSB节省了载波功率，但所需的传输带宽仍为调制信号带宽的两倍4 调整载波频率及m0，观察分的AM的过调与DSB反相点现象AM信号用于包络检波器检波时，当调制信号的最大值大于直流信号时，会出现“过调制”现象，此时会发生失真。
DSB相位在0.2π角频率处频率处明显存在反相点现象5 在接收端带通后加上窄带高斯噪声，单边功率谱密度，重新解调加入相同功率的高斯白噪声信号，DSB信号解调较AM信号解调误差较大二 代码clc;close all;clear all;fm=1;fc=10;T=5;t=linspace(-5,5,1024);%采样频率a=sqrt(2);mt=a*cos(2*pi*fm*t);% 带通滤波器A=2;s_db=mt.*cos(2*pi*fc*t);s_am=(A+mt).*cos(2*pi*fc*t);figure(1);subplot(2,1,1);plot(t,s_am);hold on; plot(t,A+mt,r--); title(AM调制信号);xlabel(t/s);ylabel(幅度);subplot(2,1,2);plot(t,s_db);title(DSB调制信号);xlabel(t/s);ylabel(幅度);hold on; plot(t,mt,r--); sp=s_am.*cos(2*pi*fc*t);sp2=s_db.*cos(2*pi*fc*t);fs1=0.3;fp1=0.5;fpu=1.8;fsu=2.0;Fs=50;Rp=1;Rs=42;ws1=2*pi*fs1/Fs;wp1=2*pi*fp1/Fs;wpu=2*pi*fpu/Fs;wsu=2*pi*fsu/Fs;wd1=0.5*(wp1+ws1);wd2=0.5*(wpu+wsu);bt=min(abs(ws1-wp1),abs(wsu-wpu));n=ceil(6.6*pi/bt);wn=[wd1/pi,wd2/pi];hn=fir1(n-1,wn,bandpass,hamming(n));y=conv(sp,hn);y2=conv(sp2,hn);y=2*y;y2=2*y2;figure(2);subplot(2,1,1);plot(t,y(1:1024));hold on;plot(t,mt,r--);axis([0 5 -1.5 1.5]);xlabel(t/s),ylabel(m(t)),title(AM相干解调);subplot(2,1,2);plot(t,y2(1:1024));hold on;plot(t,mt,r--);axis([0 5 -1.5 1.5]);xlabel(t/s),ylabel(m(t)),title(DSB相干解调);S=fft(s_am,1024);p=S.*conj(S)/1024; S2=fft(s_db,1024);p2=S2.*conj(S2)/1024; f=0:0.1:102.3;figure(3);subplot(2,1,1);plot(f,p);title(AM已调信号功率谱);xlabel(f/hz);ylabel(w);subplot(2,1,2);plot(f,p2);title(DSB已调信号功率谱);xlabel(f/hz);ylabel(w);s_am1=(A-1+mt).*cos(2*pi*fc*t);figure(4);subplot(2,1,1);plot(t,s_am);hold on; plot(t,A+mt,r--); title(AM不过调);xlabel(t/s);ylabel(幅度);subplot(2,1,2);plot(t,s_am1);title(AM过调);xlabel(t/s);ylabel(幅度);f=(0:1023)/1024;S1=fft(s_db,1024);xp=abs(S1);gl=angle(S1);figure(5);subplot(2,1,1);plot(f,xp);title(DSB调幅波);xlabel(w/(2*pi));ylabel(幅度);subplot(2,1,2);plot(f,gl);xlabel(w/(2*pi));ylabel(相位);hold on; f=0:0.1:102.3;s_am=s_am+wgn(1,length(t),0.1);s_db=s_db+wgn(1,length(t),0.1);sp1=s_am.*cos(2*pi*fc*t);sp2=s_db.*cos(2*pi*fc*t);fs1=0.3;fp1=0.5;fpu=1.8;fsu=2.0;Fs=50;Rp=1;Rs=42;ws1=2*pi*fs1/Fs;wp1=2*pi*fp1/Fs;wpu=2*pi*fpu/Fs;wsu=2*pi*fsu/Fs;wd1=0.5*(wp1+ws1);wd2=0.5*(wpu+wsu);bt=min(abs(ws1-wp1),abs(wsu-wpu));n=ceil(6.6*pi/bt);wn=[wd1/pi,wd2/pi];hn=fir1(n-1,wn,bandpass,hamming(n));y=conv(sp1,hn);y2=conv(sp2,hn);y=2*y;y2=2*y2;figure(6);subplot(2,1,1);plot(t,y(1:1024));hold on;plot(t,mt,r--);axis([0 5 -1.5 1.5]);xlabel(t/s),ylabel(m(t)),title(AM相干解调);subplot(2,1,2);plot(t,y2(1:1024));hold on;plot(t,mt,r--);axis([0 5 -1.5 1.5]);xlabel(t/s),ylabel(m(t)),title(DSB相干解调); . . . 。