|
信號與系統(tǒng)大作業(yè)大作業(yè)要求 課程大作業(yè)采用小組形式進(jìn)行，每4~6名學(xué)生一組，每組推選組長1名，由組長負(fù)責(zé)組織小組成員分工合作，完成語音信號時頻認(rèn)知。具體要求如下，①錄制一段語音信號，對錄制的信號進(jìn)行采樣，并畫出采樣后語音信號的時域波形和頻譜圖；②對語音信號進(jìn)行加噪，畫出加噪后時域波形和頻譜；③回放語音信號；實現(xiàn)快錄慢放；④撰寫設(shè)計報告及完成仿真，并進(jìn)行驗收答辯。 錄制一段語音信號，對錄制的信號進(jìn)行采樣，并畫出采樣后語音信號的時域波形和頻譜圖；
1、用window的附件中的“錄音機(jī)”錄制一段語音信號并存盤。（例如：d:\1\2.wav）
2、用MATLAB軟件進(jìn)行該原始語音信號的時頻分析。
2.1代碼
%對語言信號做原始的時域波形分析和頻譜分析
[y,fs]=audioread('d:\1\2.wav');
% sound(y,fs) % 回放語音信號
n=length(y) %選取變換的點數(shù)
y_p=fft(y,n); %對n 點進(jìn)行傅里葉變換到頻域
f=fs*(0:n/2-1)/n; % 對應(yīng)點的頻率
figure(1)
subplot(2,1,1);
plot(y); %語音信號的時域波形圖
title('原始語音信號采樣后時域波形');
xlabel('時間軸')
ylabel('幅值 A')
subplot(2,1,2);
plot(f,abs(y_p(1:n/2))); %語音信號的頻譜圖
title('原始語音信號采樣后頻譜圖');
xlabel('頻率Hz');
ylabel('頻率幅值');
2.2 圖像
對語音信號進(jìn)行加噪，畫出加噪后時域波形和頻譜；
2.1代碼
%對音頻信號產(chǎn)生噪聲
L=length(y) %計算音頻信號的長度
noise=0.1*randn(L,2); %產(chǎn)生等長度的隨機(jī)噪聲信號(這里的噪聲的大小取決于隨機(jī)函數(shù)的幅度倍數(shù)）
y_z=y+noise; %將兩個信號疊加成一個新的信號——加噪聲處理
%sound(y_z,fs)
%對加噪后的語音信號進(jìn)行分析
n=length(y); %選取變換的點數(shù)
y_zp=fft(y_z,n); %對n 點進(jìn)行傅里葉變換到頻域
f=fs*(0:n/2-1)/n; % 對應(yīng)點的頻率
figure(2)
subplot(2,1,1);
plot(y_z); %加噪語音信號的時域波形圖
title('加噪語音信號時域波形');
xlabel('時間軸')
ylabel('幅值 A')
subplot(2,1,2);
plot(f,abs(y_zp(1:n/2))); %加噪語音信號的頻譜圖
title('加噪語音信號頻譜圖');
xlabel('頻率Hz');
ylabel('頻率幅值');
回放語音信號；實現(xiàn)快錄慢放；
在MATLAB寫好的代碼末尾輸入sound(y,fs);語句 可實現(xiàn)語音正常放；
若輸入sound(y,1.5*fs); 可實1.5倍快播放
若輸入sound(y,0.5*fs); 可實現(xiàn)0.5倍慢播放
完整的Word格式文檔51黑下載地址：
信號與系統(tǒng)大作業(yè).docx
(114.58 KB, 下載次數(shù): 24)
2017-12-4 22:55 上傳
點擊文件名下載附件
下載積分: 黑幣 -5
| 
QQ好友和群
QQ空間
騰訊微博
騰訊朋友
收藏
淘帖
頂
踩
