[xFs]?=?audioread（‘D:\北科大課程\信號分析與處理\信號分析處理大作業\程序\002.mp3‘）;
%sound（xFs）;???????????%播放采集到的語音信號
N?=?length（x）;
t0?=?linspace（0N/FsN）;
n?=?0?:?N?-?1;

figure（1）;
subplot（211）;?
plot（t0x）;?????????????%畫出原始語音信號的時域波形圖
title（‘原始音頻信號時域波形圖‘）;
xlabel（‘時?間?t/s‘）;
ylabel（‘幅?度?n‘）;
subplot（212）;
f?=?fft（xN）;????????????%對原始信號進行快速傅立葉變換
plot（abs（n）abs（f））;?????%畫出原始語音信號經過傅立葉變換后的頻譜
title（‘原始音頻信號幅頻圖‘）;
xlabel（‘頻?率?f/Hz‘）;
ylabel（‘幅?度?n‘）;

%FIR低通濾波器程序
fb?=?290;fc?=?300;As?=?120;??%確定技術指標
wp?=?2?*?pi?*?fb?/?Fs;
ws?=?2?*?pi?*?fc?/?Fs;????????%求歸一化頻率，其中Fs由語音信號提取函數audioread返回值得到
wo?=?ws?-?wp;???????????%求過渡帶
beta?=?0.1102?*?（As?-?0.87）;??%求凱澤窗函數參數
N?=?ceil（（As?-?8）?/?2.285?/?wo）;??%由過渡帶決定N的值
wc?=?（wp?+?ws）?/?2?/?pi;
b?=?fir1（Nwckaiser（N?+?1beta））;
[Hww]?=?freqz（b1）;
figure（2）;
plot（w?*?Fs?/?2?/?pi?20?*?log10（abs（Hw）））;??%繪制頻率響應曲線
grid;
xlabel（‘頻?率f/Hz‘）;
ylabel（‘H（f）‘）;
title（‘FIR低通濾波器的幅頻特性‘）;

%FIR高通濾波器程序
fp?=?6000;?fs?=?5900;?As?=?135;%確定技術指標
wp?=?2?*?pi?*?fp?/?Fs;
ws?=?2?*?pi?*?fs?/?Fs;????????%求歸一化頻率，其中Fs由語音信號提取函數audioread返回值得到
wo?=?wp?-?ws;????????????%求過渡帶
beta?=?0.1102?*?（As?-?0.87）;??%求凱澤窗函數參數
N?=?ceil（（As?-?8）?/?2.285?/?wo）;??%由過渡帶決定N的值
wc?=?（wp?+?ws）?/?2?/?pi;
b1?=?fir1（Nwc‘high‘kaiser（N+1beta））;
[Hww]?=?freqz（b11）;
figure（3）;
plot（w?*?Fs?/?2?/?pi?20?*?log10（abs（Hw）））;??%繪制頻率響應曲線
grid;
xlabel（‘頻?率f/Hz‘）;
ylabel（‘H（f）‘）;
title（‘FIR高通濾波器的幅頻特性‘）;

%FIR帶通濾波器程序
fp1?=?980;?fc1?=?1000;?fp2?=?4200;?fc2?=?4300;?rs?=?120;????%確定技術指標
wp1?=?2?*?pi?*?fc1?/?Fs;
ws1?=?2?*?pi?*?fp1?/?Fs;
wp2?=?2?*?pi?*?fp2?/?Fs;
ws2?=?2?*?pi?*?fc2?/?Fs;??????%求歸一化頻率，其中Fs由語音信號提取函數audioread返回值得到
wo1?=?wp1?-?ws1;
wo2?=?wp2?-?ws2;
wo?=?max（wo1?wo2）;????????????????????????%求過渡帶
beta?=?0.1102?*?（rs?-?0.87）;??????????????????????%求凱澤窗函數參數
N?=?ceil（（rs?-?8）?/?2.285?/?wo）;?????????????????????%由過渡帶決定N的值
wn1?=?（fp1?+?fc1）?/?Fs;
wn2?=?（fp2?+?fc2）?/?Fs;
wn?=?[wn1wn2];
b2?=?fir1（Nwn‘band‘kaiser（N?+?1beta））;
[Hww]?=?freqz（b21512Fs）;
figure（4）;
plot（w20?*?log10（abs（Hw）））;??????????????????%繪制頻率響應曲線
grid;
xlabel（‘頻?率f/Hz‘）;
ylabel（‘H（f）‘）;
title（‘FIR帶通濾波器的幅頻特性‘）;

figure（4）;
y?=?filter（b1x）;
subplot（211）;
plot（t0x）;
title（‘原始音頻信號時域波形圖‘）;
xlabel（‘時?間?t/s‘）;
ylabel（‘幅?度?n‘）;
subplot（212）;
plot（t0y）;
title（‘低通濾波后音頻信號時域波形圖‘）;
xlabel（‘時?間?t/s‘）;
ylabel（‘幅?度?n‘）;

figure（5）;
y1?=?filter（b11x）;
y2?=?filter（b21x）;
subplot（211）;
plot（t0y1）;
title（‘高通濾波后音頻信號時域波形圖‘）;
xlabel（‘時?間?t/s‘）;
ylabel（‘幅?度?n‘）;
subplot（212）;
plot（t0y2）;
title（‘帶通濾波后音頻信號時域波形圖‘）;
xlabel（‘時?間?t/s‘）;
ylabel（‘幅?度?n‘）;

figure（6）;
N1?=?length（x）;
t1?=?linspace（0N1/fsN1）;
n1?=?0?:?N1?-?1;
X?=?fft（xN1）;????????????%對原始信號進行快速傅立葉變換
Y?=?fft（yN1）;
subplot（211）;
plot（abs（n1）abs（X））;???%畫出原始語音信號經過傅立葉變換后的頻譜
title（‘原始音頻信號幅頻圖‘）;
xlabel（‘頻?率?f/Hz‘）;
ylabel（‘幅?度?n‘）;
subplot（212）;
plot（abs（n）abs（Y））;???%畫出原始語音信號經過傅立葉變換后的頻譜
title（‘低通濾波后音頻信號幅頻圖‘）;
xlabel（‘頻?率?f/Hz‘）;
ylabel（‘幅?度?n‘

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????193051??2019-11-06?22:03??程序\002.mp3

?????文件???????4093??2019-11-07?08:19??程序\Xinhao_1.m

?????文件??????36763??2019-11-07?21:15??程序\低通.jpg

?????文件??????58018??2019-11-07?08:10??程序\低通濾波時域.jpg

?????文件??????34802??2019-11-07?08:11??程序\低通濾波頻域對比.jpg

?????文件?????709748??2019-11-07?21:50??程序\信號大作業.pdf

?????文件??????47940??2019-11-07?08:30??程序\原始.jpg

?????文件??????53992??2019-11-07?08:19??程序\高通_帶通.jpg

?????文件??????37704??2019-11-07?08:20??程序\高通_帶通_頻域.jpg

?????目錄??????????0??2020-11-02?19:09??程序

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

??????????????1176111????????????????????10