检波波形
2、编写matlab程序实现DSB振幅调制与解调制的设计与仿真；
1）设计DSB振幅调制与解调制仿真电路，要求调制信号的幅度A和频率F可变；载波信号的幅度A和频率F可变；
2）绘制调制信号u21、载波信号uc2和已调波信号udsb1的时域波形图和频谱图（要求谱线清晰）；
3）要求调制信号为三个以上正弦波信号的合成，幅度和频率均可变，绘制调制信号u22、载波信号uc2和已调波信号udsb2的时域波形图和频谱图（要求谱线清晰）；观察相位突变点处的波形；
4）用同步检波对已调波信号udsb1进行解调制，在同一幅图中绘制原调制信号u21和解调后的信号y21，并绘制这两个信号的频谱图；5）用同步检波对已调波信号udsb2进行解调制，在同一幅图中绘制原调制信号u22和解调后的信号y22，并绘制这两个信号的频谱图。
1.单调制信号的DSB波
function[u21,uc2,udsb1] =dantizozhiDSBbo(a1,a2,f1,f2) ? 调制信号的幅度;
? 调制信号的频率; ￠ 载波信号的幅度; ò 载波信号的频率; fs=20*max(f1,f2);
t=0:1/fs:2/min(f1,f2);
u21=a1*cos(2*pi*f1*t); %设调制信号频率 uc2=a2*cos(2*pi*f2*t); %设载波频率 udsb1= u21.* uc2; %调制信号与载波相乘 subplot(3,1,1);plot(t, u21);grid on subplot(3,1,2);plot(t, uc2);grid on subplot(3,1,3);plot(t, udsb1);grid on X=fftshift(fft(u21)); Y=fftshift(fft(uc2)); Z=fftshift(fft(udsb1)); figure;
F=linspace(-fs/2,fs/2,length(t)); subplot(3,1,1);plot(F,abs(X));grid on subplot(3,1,2);plot(F,abs(Y));grid on subplot(3,1,3);plot(F,abs(Z));grid on s1=udsb1.*uc2; Rp=3;
Rs=60;%信号衰减幅度 Wp=40/500;%通带截止频率 Ws=200/500;%阻带截止频率
[n,Wn]=ellipord(Wp,Ws,Rp,Rs);%阶数n
[b,a]=ellip(n,Rp,Rs,Wn);%传递函数分子分母b，a y21=5*filter(b,a,s1); A=fftshift(fft(u21)); B=fftshift(fft(y21));
F=linspace(-fs/2,fs/2,length(t)); figure;
subplot(4,1,1);plot(t,u21);grid on subplot(4,1,2);plot(t,y21);grid on
subplot(4,1,3);plot(F(1000:1400),abs(A(1000:1400)));grid on subplot(4,1,4);plot(F(1000:1400),abs(B(1000:1400)));grid on 在命令窗口出入：
[u21,uc2,udsb1] =dantizozhiDSBbo(1,1,500,30000) 得到：
三个调制信号的调制波形
频谱图
检波波形
2.三个输入的DSB波
function[u21,uc2,udsb1] =duogexinhaodsbbo (a11,a12,a13,a2,f11,f12,f13,f2)
?1 调制信号1的幅度; ?2 调制信号2的幅度; ?3 调制信号3的幅度; ?1 调制信号1的频率; ?2 调制信号2的频率; ?3 调制信号3的频率; ￠ 载波信号的幅度; ò 载波信号的频率; fs=20*f2; k=1;
t=0:1/fs:2/f11;
u21=a11*cos(2*pi*f11*t)+ a12*cos(2*pi*f12*t)+ a13*cos(2*pi*f13*t); %设调制信号频率
uc2=a2*cos(2*pi*f2*t); %设载波频率
udsb1=k* u21.* uc2; %调制信号与载波相乘 subplot(3,1,1);plot(t, u21);grid on subplot(3,1,2);plot(t, uc2);grid on subplot(3,1,3);plot(t, udsb1);grid on X=fftshift(fft(u21)); Y=fftshift(fft(uc2)); Z=fftshift(fft(udsb1));
figure;
F=linspace(-fs/2,fs/2,length(t));
subplot(3,1,1);plot(F(1000:1400),abs(X(1000:1400)));grid on subplot(3,1,2);plot(F(1000:1400),abs(Y(1000:1400)));grid on subplot(3,1,3);plot(F(1000:1400),abs(Z(1000:1400)));grid on s1=udsb1.*uc2; Rp=3;
Rs=60;%信号衰减幅度 Wp=40/500;%通带截止频率 Ws=200/500;%阻带截止频率
[n,Wn]=ellipord(Wp,Ws,Rp,Rs);%阶数n
[b,a]=ellip(n,Rp,Rs,Wn);%传递函数分子分母b，a y21=5*filter(b,a,s1); A=fftshift(fft(u21)); B=fftshift(fft(y21));
F=linspace(-fs/2,fs/2,length(t)); figure;
subplot(4,1,1);plot(t,u21);grid on subplot(4,1,2);plot(t,y21);grid on
subplot(4,1,3);plot(F(1000:1400),abs(A(1000:1400)));grid on subplot(4,1,4);plot(F(1000:1400),abs(B(1000:1400)));grid on 在命令窗口出入：
[u21,uc2,udsb1] =duogexinhaodsbbo (1, 2, 3,1,500,700,1000,30000) 得到：
三个调制信号的调制波形