%生成調制信號，并加入加性噪聲和窄帶干擾

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）;