[gqfs]=audioread（‘piano.wav‘）;????????????????????????????????????????????%讀取音頻文件
n=length（gq）;??????????????????????????????????????????????????????????????%計(jì)算信號(hào)長(zhǎng)度
gq1=fft（gqn）;?????????????????????????????????????????????????????????????%對(duì)音頻數(shù)據(jù)進(jìn)行快速傅里葉變化?????
len?=?round（（length（gq1））/2）;??????????????????????????????????????????????%取一半數(shù)據(jù)

figure（1）;?????????????????????????????????????????????????????????????????%建立一個(gè)畫(huà)圖界面
subplot（221）;????????????????????????????????????????????????????????????%該界面中圖片為兩行兩列，這里畫(huà)第一個(gè)
plot（gq）;??????????????????????????????????????????????????????????????????%畫(huà)出音頻時(shí)域圖像
title（‘初始信號(hào)波形‘）;??????????????????????????????????????????????????????%設(shè)置標(biāo)題
xlabel（‘時(shí)間‘）;?????????????????????????????????????????????????????????????%設(shè)置X軸單位
ylabel（‘幅度‘）;?????????????????????????????????????????????????????????????%設(shè)置Y軸單位
grid?on????????????????????????????????????????????????????????????????????%添加網(wǎng)格線

subplot（222）;????????????????????????????????????????????????????????????%畫(huà)第二個(gè)圖
plot（abs（gq1（1:len）））;?????????????????????????????????????????????????????%畫(huà)出音頻的頻譜圖
title（‘初始信號(hào)頻譜‘）;??????????????????????????????????????????????????????%設(shè)置標(biāo)題
xlabel（‘頻率‘）;?????????????????????????????????????????????????????????????%設(shè)置X軸單位
ylabel（‘幅度‘）;?????????????????????????????????????????????????????????????%設(shè)置Y軸單位
grid?on????????????????????????????????????????????????????????????????????%添加網(wǎng)格線

load（‘E:\lowpass.mat‘）;?????????????????????????????????????????????????????%導(dǎo)入濾波器
[ba]=tf（Hd）;??????????????????????????????????????????????????????????????%得出傳遞系數(shù)

x=filter（bagq）;??????????????????????????????????????????????????????????%濾波器進(jìn)行濾波
n=length（x）;???????????????????????????????????????????????????????????????%計(jì)算濾波后音頻數(shù)據(jù)長(zhǎng)度
X1=fft（xn）;???????????????????????????????????????????????????????????????%對(duì)濾波后音頻進(jìn)行傅里葉變化
len?=?round（（length（X1））/2）;???????????????????????????????????????????????%選取一半長(zhǎng)度

subplot（223）;????????????????????????????????????????????????????????????%畫(huà)第三個(gè)圖
plot（x）;???????????????????????????????????????????????????????????????????%畫(huà)出濾波后音頻的時(shí)域圖
title（‘濾波后信號(hào)的波形‘）;??????????????????????????????????????????????????%設(shè)置標(biāo)題
xlabel（‘時(shí)間‘）;?????????????????????????????????????????????????????????????%設(shè)置X軸單位
ylabel（‘幅度‘）;?????????????????????????????????????????????????????????????%設(shè)置Y軸單位
grid?on????????????????????????????????????????????????????????????????????%添加網(wǎng)格線

subplot（224）;????????????????????????????????????????????????????????????%畫(huà)第四個(gè)圖
plot（abs（X1（1:len）））;??????????????????????????????????????????????????????%畫(huà)出濾波后音頻的頻譜圖
title（‘濾波后信號(hào)的頻譜‘）;???????????????????????????????????????????????????%設(shè)置標(biāo)題
xlabel（‘頻率‘）;?????????????????????????????????????????????????????????????%設(shè)置X軸
ylabel（‘幅度‘）;?????????????????????????????????????????????????????????????%設(shè)置Y軸
grid?on?????????????????????????????????????????????????????????????????????%添加網(wǎng)格線

fvtool（Hd‘Analysis‘‘magestimate‘）;??????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????3513??2020-05-19?12:10??gangqin.m