clc;

 fm=;     %信号频率

 fc=;     %载波频率

 fs=;    %抽样频率

 Am=;

 A=;

 N=;

 K=N-;

 n=:N-;

 t=(:/fs:K/fs);

 yt=Am*cos(*pi*fm*t);

 figure()

 subplot(,,);plot(t,yt);title('频率为100HZ信号时域波形');

 y0=A+yt;

 y2=y0.*cos(*pi*fc*n/fs);

 y3=fft(y2,N);

 q1=(:N/-)*fs/N;

 mx1=abs(y3(:N/));

 figure()

 subplot(,,);

 plot(t,y2);

 title('已调信号时域波形');

 subplot(,,);

 plot(q1,mx1);

 title('已调信号频谱');

 yc=cos(*pi*fc*t);

 figure();

 subplot(,,),plot(t,yc),title('载波fc时域')

 n=:N-;

 yc1 = Am*cos(*pi*fc*n/fs);

 y3=fft(yc1,N)

 q=(:N/-)*fs/N

 mx=abs(y3(:N/));

 figure()

 subplot(,,),plot(q,mx),title('载波fc频谱')

 N=;

 n=:N-;

 y4=0.01*randn(,length(t)); %产生高斯噪声

 w=y4.^;                    %噪声功率

 figure();

 subplot(,,);

 plot(t,y4);

 title('高斯白噪声时域波形');

 y5=fft(y4,N);

 q2=(:N/-)*fs/N;

 mx2=abs(y5(:N/));

 figure();

 subplot(,,);

 plot(q2,mx2);

 title('高斯白噪声频域波形');

 y6=y2+y4;

 figure();

 subplot(,,);

 plot(t,y6);

 title('加噪声后时域信号')

 q3=q1;

 mx3=mx1+mx2;

 subplot(,,);

 plot(q3,mx3);

 title('加噪声后频谱')

 yv=y6.*yc;   %乘以载波想干解调

 Ws=yv.^;

 p1=fc-fm;

 [k,Wn,beta,ftype]=kaiserord([p1 fc],[ ],[0.05 0.01],fs);%数字低通过滤波器

 window=kaiser(k+,beta);%使用kaiser窗函数

 b=fir1(k,Wn,ftype,window,'noscale');%

 yt=filter(b,,yv);

 yssdb=yt.*-;

 figure()

 subplot(,,);

 plot(t,yssdb);

 title('经过低通已调信号时域波形')    ;%解调

 y9=fft(yssdb,N);

 mx=abs(y9(:N/));

 subplot(,,);

 plot(q,mx);

 title('已调信号频域波形')