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

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

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

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

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

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

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

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

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