f=100;     %取樣頻率
N=1024;             %總點數
Ts=1/f;               %取樣間隔
T=Ts*N;             %取樣總時間
n=0:N-1;
t=n*Ts;
x1=cos(pi*t);
x2=cos(9*pi*t);
x3=x1.*x2;
subplot(511),plot(x1);% 繪出調制信號
grid;title('調制信號');
subplot(512),plot(x3);% 繪出己調信號
grid;title('己調信號');
Y=abs(fft(x3,N));      %已調信號頻譜
subplot(513)
plot(n,Y);title('己調信號頻譜')
x4=x3.*cos(9*pi*t);     %解調信號
Wp=1;Ws=3;Rp=2;Rs=30;
[N,Wn]=buttord(Wp/50,Ws/50,Rp,Rs);
[b,a]=butter(N,Wn);
y=filter(b,a,x4);         %濾波
Y2=abs(fft(y))
Y3=fft(y);              %濾波后信號頻譜
subplot(514)
plot(Y2);title('濾波后信號頻譜');
Y4=ifft(Y3);
subplot(515)
plot(t,Y4);             %濾波后得到原信號
title('濾波后得到原信號');