資源簡介
gangqin.zip

代碼片段和文件信息
[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
評論
共有 條評論