生成调制信号，并加入加性噪声和窄带干扰

%生成调制信号，并加入加性噪声和窄带干扰

k=0:1023;

x1=20*sin（70*2*pi*k/40）;

x2=20*cos（70*2*pi*k/40）;

x=x1.*randn（11024）+x2.*randn（11024）+15*sin（70*2*pi*k/40）+15*sin（75*2*pi*k/42）+15*cos（80*2*pi*k/35）;

w=blackman（1024）;                  %通过blackman窗，减小频谱泄漏 

Y=fft（w‘.*x1024）;                  %对截取的数据进行快速傅里叶变换 

R=abs（real（Y））; 

I=abs（imag（Y））; 

Ampl=abs（Y）;                          %计算频谱幅度

subplot（121）;

plot（Ampl）;

title（‘采用窄带抑制算法前幅频响应‘）;xlabel（‘频率Hz‘）

axis（[0 1000 0 2800]）;

u=0;

for i=1:1024

    u=u+Ampl（i）;

end

mean1=u./1024;                     %计算频谱幅度的均值，便于分析频谱特性

Ampldb=10.*log10（Ampl）;            %为防止溢出将幅度取分贝

u=0;

m=0;

for i=1:1024

    u=u+Ampldb（i）;