%設(shè)計(jì)一個(gè)通帶為0.35π到0.65π，阻帶下截止為0.2π，上截止為0.8π
N=40;
alfa=（40-1）/2;
k=0:N-1;
w1=（2*pi/N）*k;
T1=0.109021;?T2=0.59417456;
hrs=[zeros（15）T1T2ones（17）T2T1zeros（19）T1T2ones（17）T2T1zeros（14）];
hdr=[001100];?wd1=[00.20.350.650.81];
k1=0:floor（（N-1）/2）;??????k2=floor（（N-1）/2）+1:N-1;
angH=[-alfa*（2*pi）/N*k1alfa*（2*pi/N*（N-k2））];
H=hrs.*exp（j*angH）;
h=real（ifft（H））;
[dbmagphagrdw]?=?freqz_m（h1）;
[HrwwaL]?=Hr_Type2（h）;
figure（1）
subplot（221）
plot（w1（1:21）/pihrs（1:21）‘o‘wd1hdr）
axis（[01-0.11.1]）;
title（‘帶通：N=40，T1=0.109021，?T2=0.59417456‘）
ylabel（‘Hr（k）‘）;
set（gca‘XTickMode‘‘manual‘‘XTick‘[00.20.350.650.81]）
set（gca‘YTickMode‘‘manual‘‘YTick‘[00.0590.1091]）;
grid??%繪制帶網(wǎng)格的圖像
subplot（222）;?
stem（kh）;
axis（[-1N-0.40.4]）
title（‘脈沖響應(yīng)‘）;??ylabel（‘h（n）‘）;?text（N+1-0.4‘n‘）
subplot（223）;?plot（ww/piHrw1（1:21）/pihrs（1:21）‘o‘）;
axis（[01-0.11.1]）;title（‘振幅響應(yīng)‘）
xlabel（‘頻率?（單位：pi）‘）;?ylabel（‘Hr（w）‘）
set（gca‘XTickMode‘‘manual‘‘XTick‘[00.20.350.650.81]）;
set（gca‘YTickMode‘‘manual‘‘YTick‘[00.0590.1091]）;
grid
subplot（224）;?plot（w/pidb）;?axis（[01-10010]）;
grid
title（‘幅度響應(yīng)‘）;xlabel（‘頻率?（單位：pi）‘）;?ylabel（‘分貝‘）
set（gca‘XTickMode‘‘manual‘‘XTick‘[00.20.350.650.81]）
set（gca‘YTickMode‘‘manual‘‘YTick‘[-60;0]）;
set（gca‘YTickLabelMode‘‘manual‘‘YTickLabels‘[60;0]）;

%信號取樣頻率為8000HZ
[sfsnbits]=wavread（‘111.wav‘）;
x=s（:1）;
sound（x8000）;
L=length（x）;
f=fs*（0:L-1）/L;
t=0:1/fs:（L-1）/fs;?%將所加噪聲信號的點(diǎn)數(shù)調(diào)整到與原始信號相同
%Au=1　　　　　
d=0.02*abs（max（x））*[cos（2*pi*500*t）+cos（2*pi*3300*t）]‘;?%噪聲為500和3300Hz的余弦信號
%dz=cos（0.5*pi*fs*t）‘;%載波
xz=x+d;
%pause（10）;
sound（xzfs）;?%播放加噪聲后的語音信號
X=fft（x）;?%求信號的頻譜
XZ=fft（xz）;??%信號加噪后的頻譜

figure（2）
subplot（211）;plot（tx）
title（‘未加噪的信號‘）;?xlabel（‘time?s‘）;ylabel（‘幅度?‘）;
subplot（212）;plot（txz）
title（‘加噪后的信號‘）;?xlabel（‘time?s‘）;ylabel（‘幅度?‘）;
%subplot（313）;plot（tx2）
%title（‘調(diào)制后的信號‘）;?xlabel（‘time?n‘）;ylabel（‘fuzhi?n‘）;

figure（3）
subplot（211）;
plot（fabs（X））;
title（‘原始語音信號頻譜‘）;
xlabel（‘頻率（單位：Hz）‘）;
ylabel（‘幅度‘）;
subplot（212）;
plot（fabs（XZ））;
title（‘加噪后的信號頻譜‘）;
xlabel（‘頻率（單位：Hz）‘）;
ylabel（‘幅度‘）;

y?=?fftfilt（hxz）;
Y=fft（y）;
sound（3*yfs）;


figure（4）
subplot（211）;plot（t3*y）
title（‘濾波后的信號‘）;?xlabel（‘time?s‘）;ylabel（‘幅度‘）;
subplot（212）;plot（fabs（Y））;
title（‘濾波后的信號頻譜‘）;?xlabel（‘頻率（單位：Hz）‘）;ylabel（‘幅度‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????427??2011-06-16?21:10??yuyincaiyang.m

?????文件??????66452??2011-06-16?22:04??111.wav

?????文件?????307198??2011-06-16?22:04??112.wav

?????文件??????55308??2011-06-16?22:04??113.wav

?????文件??????67839??2011-06-16?20:19??12312.wma

?????文件???????2641??2011-06-16?22:34??BPF.m

?????文件???????1387??2011-06-16?22:42??freFIR.m

?????文件????????769??2011-06-12?01:02??freqz_m.m

?????文件????????138??2011-06-12?00:09??Hr_Type2.m

-----------?---------??----------?-----??----

???????????????502159????????????????????9