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