close all;
%t<-[0,2pi]
t=linspace(0,2*pi,256);
%信源頻率為1Hz的余弦
m=cos(2*pi*t);
%載波10Hz，調(diào)制
sm=m.*cos(2*pi*10*t);
%同步解調(diào)，與載波相乘
mo=sm.*cos(2*pi*10*t);
%切比雪夫?yàn)V波器，我不會(huì)做，下面一行從別人的程序中拷來(lái)的
[b,a]=cheby1(12,0.5,100/500);  
%濾波
mo1=filter(b,a,mo);   
figure(1);
subplot(4,2,1);
%畫出信源的圖形
plot(m);
subplot(4,2,2);
%畫出信源的頻譜幅度圖形
plot(abs(fft(m,256)));
subplot(4,2,3);
%畫出與載波相乘后的信號(hào)圖形
plot(sm);
subplot(4,2,4);
%畫出sm的頻譜幅度
plot(abs(fft(sm,256)));
subplot(4,2,5);
%畫出解調(diào)時(shí)與載波相乘后的圖形
plot(mo);
subplot(4,2,6);
%畫出mo的頻譜幅度
plot(abs(fft(mo,256)));
subplot(4,2,7);
%畫出輸出
plot(mo1);
subplot(4,2,8);
%畫出mo1的幅度頻譜
plot(abs(fft(mo1,256)));