提高題音樂合成信號與系統(tǒng)課程設(shè)計報告

1、<p>　　課 程 設(shè) 計 報 告</p><p>　　課程名稱 信號與系統(tǒng)課程設(shè)計 </p><p>　　指導(dǎo)教師 </p><p><b>　　設(shè)計起止日期 </b></p><p>　　學(xué) 院 信息與通信工程 </p>

2、;<p>　　專 業(yè) 電子信息工程 </p><p>　　學(xué)生姓名 </p><p>　　班級/學(xué)號 </p><p>　　成 績 </p><p>　　指導(dǎo)老師簽字

3、 </p><p><b>　　目 錄</b></p><p>　　1、課程設(shè)計目的1</p><p>　　2、課程設(shè)計要求1</p><p>　　3、課程設(shè)計任務(wù)1</p><p>　　4、課程設(shè)計內(nèi)容1</p><p><b>　　5、總

4、結(jié)11</b></p><p><b>　　參考文獻(xiàn)12</b></p><p><b>　　附錄12</b></p><p>　　1、課程設(shè)計目的和要求</p><p><b>　　目的：</b></p><p>　　掌握MATLAB

5、的使用方法?！靶盘柵c系統(tǒng)”是一門重要的專業(yè)基礎(chǔ)課，MATLAB作為信號處理強(qiáng)有力的計算和分析工具是電子信息工程技術(shù)人員常用的重要工具之一。本課程設(shè)計基于MATLAB完成信號與系統(tǒng)綜合設(shè)計實(shí)驗(yàn)，以提高學(xué)生的綜合應(yīng)用知識能力為目標(biāo)，是“信號與系統(tǒng)”課程在實(shí)踐教學(xué)環(huán)節(jié)上的必要補(bǔ)充。通過課設(shè)綜合設(shè)計實(shí)驗(yàn)，激發(fā)學(xué)生理論課程學(xué)習(xí)興趣，提高分析問題和解決問題的能力。</p><p><b>　　要求： </

6、b></p><p>　?。?）運(yùn)用MATLAB編程得到簡單信號、簡單信號運(yùn)算、復(fù)雜信號的頻域響應(yīng)圖；</p><p>　　（2）通過對線性時不變系統(tǒng)的輸入、輸出信號的時域和頻域的分析，了解線性時不變系統(tǒng)的特性，同時加深對信號頻譜的理解。</p><p>　　2、課程設(shè)計原理</p><p>　　1離散時間信號卷積使用MATLAB中

7、的conv函數(shù) </p><p>　　2.在MATLAB中，filter函數(shù)的作用是進(jìn)行數(shù)字濾波，其實(shí)也就是系統(tǒng)的響應(yīng)序列求解；y=filter(b,a,x); y=filter(b,a,x)表示向量b和a組成的系統(tǒng)對輸入x進(jìn)行濾波，系統(tǒng)的輸出為y。 </p><p>　　3.在MATLAB中，zplane函數(shù)——繪制系統(tǒng)零極點(diǎn)圖 </p>

8、;<p>　　4.在MATLAB中，有專門的frepz()函數(shù)來求解。調(diào)用形式為 [h w]=freqz(b,a,,n,fs)其中向量b和a為離散系統(tǒng)的系數(shù)，fs為采樣頻率，n為在區(qū)間[0 fs/2]頻率范圍內(nèi)選取的頻率點(diǎn)數(shù)，f記錄頻率點(diǎn)數(shù)。由于freqz函數(shù)是采用基2的FFT的算法，n常取2的冪次方，以便提高計算速度。 </p><p>　　5.MATLAB

9、提供了如abs(),angle()等基本函數(shù)用來計算DTFT幅度，相位。</p><p>　　3、課程設(shè)計所用設(shè)備</p><p><b>　?。?）個人電腦一臺</b></p><p>　?。?）MATLAB2012設(shè)計程序</p><p>　　4、課程設(shè)計內(nèi)容和步驟</p><p>　?。?）

10、根據(jù)設(shè)計題目的要求，熟悉相關(guān)內(nèi)容的理論基礎(chǔ)，理清程序設(shè)計的措施和步驟；</p><p>　?。?）根據(jù)設(shè)計題目的要求，提出各目標(biāo)的實(shí)施思路、方法和步驟；</p><p>　?。?）根據(jù)相關(guān)步驟完成MATLAB程序設(shè)計，所編程序應(yīng)能完整實(shí)現(xiàn)設(shè)計題目的要求；</p><p>　　（4）調(diào)試程序，分析相關(guān)理論；</p><p>　　（5）編寫設(shè)計報

11、告。</p><p>　　5、設(shè)計過程、實(shí)現(xiàn)過程、實(shí)現(xiàn)結(jié)果</p><p><b>　?。ㄒ唬┗静糠?lt;/b></p><p>　　（1）信號的時頻分析</p><p>　　任意給定單頻周期信號的振幅、頻率和初相，要求準(zhǔn)確計算出其幅度譜，并準(zhǔn)確畫出時域和頻域波形，正確顯示時間和頻率。</p><p&g

12、t;<b>　　設(shè)計思路：</b></p><p>　　首先給出橫坐標(biāo)，即時間，根據(jù)設(shè)定的信號的振幅、頻率和初相，寫出時域波形的表達(dá)式；然后對時域波形信號進(jìn)行傅里葉變化，得到頻域波形；最后使用plot函數(shù)繪制各個響應(yīng)圖。</p><p><b>　　源程序：</b></p><p>　　%基本題1,時頻分析，2011010

13、512</p><p>　　clc,clear,close all</p><p>　　A=4;fc=200;a=0*pi/3;</p><p><b>　　T=1/fc;</b></p><p>　　fs=15*fc;Ts=1/fs;</p><p>　　t=0:Ts:3*T-Ts;</p

14、><p>　　N=length(t);</p><p>　　x=A*sin(2*pi*fc*t+a);</p><p>　　subplot(211),plot(t,x)</p><p>　　xlabel('t(s)'),ylabel('x(t)'),title('時域波形圖')</p>

15、<p>　　X=fft(x);XX=fftshift(X);</p><p>　　magX=abs(XX);phaseX=angle(XX);</p><p>　　f=(0:N-1)*(fs/N)-fs/2;</p><p>　　subplot(223),plot(f,magX)</p><p>　　xlabel('f(H

16、z)'),ylabel('|X(jw)|'),title('信號幅度譜')</p><p>　　subplot(224),plot(f,phaseX)</p><p>　　xlabel('f(Hz)'),ylabel('angle[X(jw)]'),title('信號相位譜')</p>&

17、lt;p><b>　　結(jié)果圖：</b></p><p><b>　　結(jié)果分析：</b></p><p>　　cos函數(shù)波形為周期信號，其頻域響應(yīng)為兩個沖擊函數(shù)，并且符合對偶性。</p><p>　?。?）傅里葉級數(shù)分析</p><p>　　分析周期鋸齒波的傅里葉級數(shù)系數(shù)，用正弦信號的線性組合構(gòu)

18、成鋸齒波，要求諧波次數(shù)可以任意輸入，分析不同諧波次數(shù)所構(gòu)成的鋸齒波，解釋是否存在吉伯斯現(xiàn)象。</p><p><b>　　設(shè)計思路：</b></p><p>　　根據(jù)原理——任意周期信號都可以表達(dá)成傅里葉級數(shù)的形式，對周期舉行鋸齒波進(jìn)行傅里葉級數(shù)分解，利用for循環(huán)完成級數(shù)求和運(yùn)算。當(dāng)N值較小時，傅里葉級數(shù)的逼近效果不是很理想，隨著N值變大，傅里葉級數(shù)越來越接近理想值

19、。</p><p><b>　　源程序：</b></p><p>　　%基本題2，傅里葉級數(shù)分析，2011010512</p><p>　　clc; clear; close all; %初始化工作空間</p><p>　　t=-6:6/1000:6; %時間</p><p>　　N = inp

20、ut('please input the nunber:'); </p><p>　　w0=pi; %角頻率</p><p><b>　　XN=0; </b></p><p>　　for n=1:N; </p><p>　　fn=(-1)^(n+1)/(n*pi); </p>

21、<p>　　XN=XN+fn*sin(n*w0*t); </p><p><b>　　end </b></p><p>　　XN=XN+0.5; %抬升0.5</p><p>　　figure; plot(t,XN); </p><p>　　title(['Gibbs,N=',nu

22、m2str(N)]); </p><p>　　xlabel('Time(sec)'); </p><p>　　ylabel(['X',num2str(N),('t')]);</p><p><b>　　結(jié)果圖；</b></p><p><b>　　結(jié)果分析：<

23、;/b></p><p>　　隨著N值的增大，傅里葉級數(shù)與理想的周期鋸齒波越來越接近，符合預(yù)期效果。</p><p><b>　　（3）系統(tǒng)分析</b></p><p>　　任意給定微分方程或差分方程描述的系統(tǒng)，畫出系統(tǒng)的幅頻響應(yīng)和相頻響應(yīng)。</p><p><b>　　設(shè)計思路：</b>&l

24、t;/p><p>　　根據(jù)微分方程、差分方程與系統(tǒng)函數(shù)的對應(yīng)關(guān)系，結(jié)合matlab自帶的freqs和freqz兩個函數(shù)，分析系統(tǒng)的幅頻特性和相頻特性。</p><p><b>　　源程序：</b></p><p>　　%基本題3 系統(tǒng)分析 2011010512</p><p>　　clc,clear,close all,&

25、lt;/p><p>　　%微分方程為 r''(t)+6r'(t)+8r（t）=2*e(t)

26、 </p><p>　　b=[0,0,2]; %系統(tǒng)函數(shù)H(s)分子多項(xiàng)式系數(shù)</p><p>　　a=[1,6,8]; %系統(tǒng)函數(shù)H(s)分母多項(xiàng)式系數(shù)</p><p>　　[Hw1,w1]=freqs(b,a); %求連續(xù)時間

27、系統(tǒng)頻率響應(yīng)</p><p><b>　　figure(1)</b></p><p>　　subplot(2,1,1);zplane(b,a);title('連續(xù)時間系統(tǒng)零極點(diǎn)圖');</p><p>　　subplot(2,2,3);plot(w1/pi,abs(Hw1))</p><p>　　xlab

28、el('rad/s');ylabel('|H(jw)|');title('連續(xù)時間系統(tǒng)幅頻響應(yīng)');</p><p>　　subplot(2,2,4);</p><p>　　plot(w1/pi,angle(Hw1))</p><p>　　xlabel('rad/s');ylabel('Phas

29、e');title('連續(xù)時間系統(tǒng)相頻響應(yīng)');</p><p>　　%差分方程為 y（n）-0.5*y(n-1)+4/5*y(n-2)=3*x(n)-0.6*x(n-1);</p><p>　　b=[3,-0.6]; %系統(tǒng)函數(shù)H(z)分子多項(xiàng)式系數(shù)</p><p>　　a=[1,-0.5,4/5]; %系統(tǒng)函數(shù)H(z)分母多項(xiàng)

30、式系數(shù)</p><p>　　[Hw,w]=freqz(b,a); %求離散時間系統(tǒng)頻率響應(yīng)</p><p><b>　　figure(2)</b></p><p>　　subplot(2,1,1);zplane(b,a);title('離散時間系統(tǒng)零極點(diǎn)圖');</p><p>　　subplot(2,2

31、,3);plot(w/pi,abs(Hw))</p><p>　　xlabel('w/pi)');ylabel('|H(e^j^w)|');title('離散時間系統(tǒng)幅頻響應(yīng)');</p><p>　　subplot(2,2,4);plot(w/pi,angle(Hw))</p><p>　　xlabel('w

32、/pi');ylabel('Phase');title('離散時間系統(tǒng)相頻響應(yīng)');</p><p><b>　　結(jié)果圖；</b></p><p>　　結(jié)果分析：所顯示零極點(diǎn)等圖符合驗(yàn)算結(jié)果</p><p>　?。?）音樂合成程序設(shè)計</p><p>　　對于任意一小段音樂，利用“

33、十二平均律”計算該音樂中各個樂音的頻率，產(chǎn)生并播放這些樂音。分析音樂的頻譜，從中識別出不同的樂音。</p><p><b>　　設(shè)計思路：</b></p><p>　　根據(jù)時間長短來區(qū)別各個音符拍長短，根據(jù)頻率高低來區(qū)別各個音符的音調(diào)，具體細(xì)節(jié)參考“十二平均律”，將每段樂音連接起來，使用sound函數(shù)播放樂音。對每段樂音進(jìn)行傅里葉變換，分析其頻譜。</p>

34、<p><b>　　源程序：</b></p><p>　　%基本題4，音樂合成程序設(shè)計，一分錢，2011010512</p><p>　　clc,clear,close all,</p><p><b>　　fs=8e3;</b></p><p>　　freq=[384,512,426

35、.66,512,426.66, 320,384,288,320,384, 320,384,426.66,512, 384,426.66,384,320, 384,256,288, 320,288,256,288,320, 436.66,384,320,384,436.66 384,512,426.66,384,320,384,288, 384,288,320,256];</p><p&

36、gt;　　N=[1/2,1/2,1/4,1/4,1/2, 1/4,1/4,1/4,1/4,1, 1/4,1/4,1/4,1/4, 1/4,1/4,1/4,1/4, 1/4,1/2,1, 1/4,1/4,1/4,1/4,1, 1/4,1/4,1/4,1/4,1, 1/4,1/4,1/4,1/4, 1/4,1/4,1/2, 1/4,1/4,1/2,1/2]*fs;</p><p>

37、　　D=zeros(1,sum(N));</p><p><b>　　X=1;</b></p><p>　　for k=1:42</p><p>　　D(X:X+N(k)-1)=(sin(2*pi*freq(k)*(1/fs:1/fs:N(k)/fs)));</p><p><b>　　X=X+N(k);<

38、;/b></p><p><b>　　end;</b></p><p>　　sound(D,fs);</p><p>　　結(jié)果分心：所聽音樂符合預(yù)期。</p><p><b>　　（5）調(diào)制</b></p><p>　　分析單位沖激響應(yīng)為的系統(tǒng)的濾波特性，畫出其幅頻響應(yīng)

39、曲線。</p><p><b>　　設(shè)計思路：</b></p><p>　　根據(jù)傅里葉變換特性，時域相乘對應(yīng)頻域卷積，將題目中的時域波形進(jìn)行分解，并將兩個函數(shù)分別對應(yīng)的傅里葉變換進(jìn)行卷積。參考常用的傅里葉變換對，所求的頻域特性為方波和沖擊函數(shù)的卷積。</p><p><b>　　源程序：</b></p>&l

40、t;p>　　%基本題5，調(diào)制，2011010512</p><p>　　clc,clear,close all,</p><p><b>　　fs=300;</b></p><p>　　t=-10:1/fs:10;</p><p>　　N=length(t);</p><p>　　n=[-(

41、N-1)/2:(N-1)/2];</p><p>　　h=sin(50*pi*t).*sin(100*pi*t)./(pi*t+0.00001);</p><p>　　f1=fft(h);</p><p>　　Fn=fftshift(f1);</p><p><b>　　f=n/N*fs;</b></p>

42、<p><b>　　figure(1)</b></p><p>　　plot(t,h),</p><p>　　xlabel('t'),ylabel('Amplitude'),title('H(t)');</p><p><b>　　figure(2)</b><

43、/p><p>　　plot(f,abs(Fn)),</p><p>　　xlabel('w'),ylabel('Amplitude'),title('Modulation');</p><p><b>　　結(jié)果圖；</b></p><p>　　結(jié)果分析：所得結(jié)果符合演算過程，與

44、預(yù)期一致。</p><p><b>　　（二）提高部分</b></p><p>　　題目1、音樂合成 2011010512</p><p>　　設(shè)計思路：通過定義每個音來寫樂譜，在加包絡(luò)諧波，在播放音樂，后畫出波形。 </p><p><b>　　源程序</b></p><p&

45、gt;<b>　　%提高題 音樂合成</b></p><p>　　clear;clc;close all;</p><p>　　TUNE_BASED=220*[2^(3/12) 2^(5/12) 2^(7/12) 2^(8/12) 2^(10/12) 2^(12/12) 2^(14/12) ]; %C調(diào)</p><p>　　H_TUNE=

46、2*TUNE_BASED;</p><p>　　L_TUNE=0.5*TUNE_BASED;</p><p>　　l_dao=L_TUNE(1); %將“l(fā)_dao”宏定義為低音“1”的頻率262Hz</p><p>　　l_re =L_TUNE(2); %將“l(fā)_re”宏定義為低音“2”的頻率286Hz</p><p

47、>　　l_mi =L_TUNE(3); %將“l(fā)_mi”宏定義為低音“3”的頻率311Hz</p><p>　　l_fa =L_TUNE(4); %將“l(fā)_fa”宏定義為低音“4”的頻率349Hz</p><p>　　l_sao= L_TUNE(5); %將“l(fā)_sao”宏定義為低音“5”的頻率392Hz</p><p&g

48、t;　　l_la =L_TUNE(6); %將“l(fā)_a”宏定義為低音“6”的頻率440Hz</p><p>　　l_xi =L_TUNE(7); %將“l(fā)_xi”宏定義為低音“7”的頻率494Hz</p><p>　　%以下是C調(diào)中音的音頻宏定義</p><p>　　dao= TUNE_BASED(1); %將“dao”宏定

49、義為中音“1”的頻率523Hz</p><p>　　re= TUNE_BASED(2); %將“re”宏定義為中音“2”的頻率587Hz</p><p>　　mi= TUNE_BASED(3); %將“mi”宏定義為中音“3”的頻率659Hz</p><p>　　fa= TUNE_BASED(4); %將“fa”宏定義為中音

50、“4”的頻率698Hz</p><p>　　sao= TUNE_BASED(5); %將“sao”宏定義為中音“5”的頻率784Hz</p><p>　　la= TUNE_BASED(6); %將“l(fā)a”宏定義為中音“6”的頻率880Hz</p><p>　　xi= TUNE_BASED(7); %將“xi”宏定義為中音“7

51、”的頻率523H </p><p>　　%以下是C調(diào)高音的音頻宏定義</p><p>　　h_dao= H_TUNE(1); %將“h_dao”宏定義為高音“1”的頻率1046Hz</p><p>　　h_re= H_TUNE(2); %將“h_re”宏定義為高音“2”的頻率1174Hz</p><p>

52、　　h_mi= H_TUNE(3); %將“h_mi”宏定義為高音“3”的頻率1318Hz</p><p>　　h_fa= H_TUNE(4); %將“h_fa”宏定義為高音“4”的頻率1396Hz</p><p>　　h_sao= H_TUNE(5); %將“h_sao”宏定義為高音“5”的頻率1567Hz</p><p>

53、;　　h_la= H_TUNE(6); %將“h_la”宏定義為高音“6”的頻率1760Hz</p><p>　　h_xi= H_TUNE(7); %將“h_xi”宏定義為高音“7”的頻率1975Hz</p><p><b>　　%% 簡譜 歡樂頌</b></p><p>　　song=[mi mi fa sao..

54、.</p><p>　　sao fa mi re...</p><p>　　dao dao re mi...</p><p>　　mi re re ,...</p><p>　　mi mi fa sao...</p><p>　　sao fa mi re...</p><p>　　dao dao

55、 re mi...</p><p>　　re dao dao ,...</p><p>　　re re mi dao...</p><p>　　re mi fa mi dao...</p><p>　　re mi fa mi re...</p><p>　　dao re l_sao mi...</p>&

56、lt;p>　　mi mi fa sao...</p><p>　　sao fa mi fa re... </p><p>　　dao dao re mi...</p><p>　　re dao dao 0];</p><p>　　delay=0.5*[1 1 1 1,...</p><p>　　1 1

57、1 1,...</p><p>　　1 1 1 1,...</p><p>　　1.5 0.5 1 ,...</p><p>　　1 1 1 1,...</p><p>　　1 1 1 1,...</p><p>　　1 1 1 1,...</p><p>　　1.5 0.5 1 ,...<

58、;/p><p>　　1 1 1 1,...</p><p>　　1 0.5 0.5 1 1,...</p><p>　　1 0.5 0.5 1 1,...</p><p>　　1 1 1 1,...</p><p>　　1 1 1 1,...</p><p>　　1 1 1 0.5 0.5,...&

59、lt;/p><p>　　1 1 1 1,...</p><p>　　1.5 0.5 1 1];</p><p><b>　　%% 設(shè)置采樣精度</b></p><p><b>　　fs=8000;</b></p><p><b>　　%% 生成波形</b>&

60、lt;/p><p>　　wav_normal=0; %普通波形</p><p>　　wav_low=0; %低八度波形</p><p>　　wav_hight=0; %高八度的波形</p><p>　　wav_envelope_only=0; %包絡(luò)的波形</p><p>　　wav_envelope_harmonic=0

61、; %加包絡(luò)加諧波的波形</p><p><b>　　m=1;</b></p><p>　　for i=1:length(song) </p><p>　　n=1:fs*delay(i);</p><p>　　temp=sin(2*pi*song(i)*n/fs); %普通</p>

62、;<p>　　temp_low=sin(pi*song(i)*n/fs); %低八度</p><p>　　temp_highet=sin(4*pi*song(i)*n/fs); %高八度</p><p>　　temp_envelope=exp((-1/4000)*n); %指數(shù)衰減包絡(luò)</p>&

63、lt;p>　　temp_harmonic=0.2*sin(4*pi*song(i)*n/fs); %2倍頻諧波</p><p>　　wav_normal(m:m+fs*delay(i)-1)=temp; %正常波形</p><p>　　wav_low(m:m+fs*delay(i)-1)=temp_low; %低八度波形</p><

64、;p>　　wav_hight(m:m+fs*delay(i)-1)=temp_highet; %高八度波形</p><p>　　wav_envelope_only(m:m+fs*delay(i)-1)=temp.*temp_envelope; %只加包絡(luò)波形</p><p>　　wav_harmonic_only(m:m+fs*delay(i)-1)=temp+temp_ha

65、rmonic; %只加諧波波形</p><p>　　wav_envelope_harmonic(m:m+fs*delay(i)-1)=(temp+temp_harmonic).*temp_envelope; %包絡(luò)諧波波形</p><p>　　m=m+fs*delay(i);</p><p><b>　　end</b></p>

66、;<p><b>　　%% 播放</b></p><p>　　sound(wav_normal,fs); %播放原始聲音</p><p>　　sound(wav_envelope_only,fs); %播放加包絡(luò)聲音</p><p>　　sound(wav_envelope_harmoni

67、c,fs); %播放加包絡(luò)加諧波聲音</p><p>　　%% 抽取前13個音符(8秒)分析包絡(luò)與諧波</p><p>　　num_song=13;</p><p>　　length_num_song=fs*sum(delay(1:num_song));</p><p><b>　　figure;</b></

68、p><p>　　subplot(311);</p><p>　　plot((1:length_num_song)/8000,wav_normal(1:length_num_song));</p><p>　　xlabel('時間/s');ylabel('幅值');title('普通 無包絡(luò) 無諧波 波形圖');</p

69、><p>　　subplot(312);</p><p>　　plot((1:length_num_song)/8000,wav_envelope_only(1:length_num_song));</p><p>　　xlabel('時間/s');ylabel('幅值');title('加包絡(luò) 無諧波 波形圖');<

70、;/p><p>　　subplot(313);</p><p>　　plot((1:length_num_song)/8000,wav_envelope_harmonic(1:length_num_song));</p><p>　　xlabel('時間/s');ylabel('幅值');title('加包絡(luò) 加諧波 波形圖

71、9;);</p><p>　　%% 抽取前13個音符(8秒)分析頻域諧波效果</p><p>　　normal=abs(fftshift(fft(wav_normal(1:length_num_song)))); %普通頻譜圖</p><p>　　harmonic=abs(fftshift(fft(wav_harmonic_only(1:length_num_so

72、ng))));%加諧波頻譜圖</p><p>　　fft_rang=(0:length_num_song-1)/length_num_song*fs-fs/2; %頻域范圍</p><p><b>　　figure;</b></p><p>　　subplot(211);</p><p>　　plot(fft_rang,

73、normal);</p><p>　　xlabel('頻率/Hz');ylabel('幅值');title('普通音無諧波頻譜圖');</p><p>　　subplot(212);</p><p>　　plot(fft_rang,harmonic);</p><p>　　xlabel('

74、;頻率/Hz');ylabel('幅值');title('普通音加諧波頻譜圖');</p><p><b>　　figure;</b></p><p>　　%% 抽取前13個音符(8秒)分析頻域的高八度低八度效果</p><p>　　hight=abs(fftshift(fft(wav_hight(1:l

75、ength_num_song)))); %高八度頻譜圖</p><p>　　low=abs(fftshift(fft(wav_low(1:length_num_song)))); %低八度頻譜圖</p><p>　　subplot(311);</p><p>　　plot(fft_rang,normal);</p><p>　　xlabe

76、l('頻率/Hz');ylabel('幅值');title('普通音頻譜圖');</p><p>　　subplot(312);</p><p>　　plot(fft_rang,low);</p><p>　　xlabel('頻率/Hz');ylabel('幅值');title('

77、;低八度音頻譜圖');</p><p>　　subplot(313);</p><p>　　plot(fft_rang,hight);</p><p>　　xlabel('頻率/Hz');ylabel('幅值');title('高八度音頻譜圖');</p><p>　　結(jié)果分析：與預(yù)期一

78、致 聽到音樂也和預(yù)期一致。</p><p><b>　　6、設(shè)計總結(jié)</b></p><p>　　通過這次試驗(yàn)，我明白了學(xué)到的知識只有會用才有價值，只有在應(yīng)用中才能找到自己的不足并加以彌補(bǔ)，如果我們不去將知識應(yīng)用，那我們學(xué)的東西將失去價值。</p><p>　　在設(shè)計過程中，自己學(xué)會了Matlab軟件的部分功能，在今后的生活中我希望能夠?qū)W到

79、更多！</p><p>　　在設(shè)計工程中我遇到了很多困難，但是我堅持作了下來，這對我來說是一次磨練，令我更加堅強(qiáng)。在今后的生活中，當(dāng)我們遇到困難我們不應(yīng)該輕言放棄而是應(yīng)該堅持，努力實(shí)現(xiàn)自己的目標(biāo)。</p><p><b>　　參考文獻(xiàn)</b></p><p>　　[1] 鄭君里、應(yīng)啟珩 、楊為理，信號與系統(tǒng)引論，北京：高等教育出版社，2009年

80、3月.</p><p>　　[2] 谷源濤、應(yīng)啟珩、鄭君里，信號與系統(tǒng)——MATLAB綜合實(shí)驗(yàn)，北京：高等教育出版社，2008年1月.</p><p>　　[3] 梁虹等，信號與系統(tǒng)分析及Matlab實(shí)現(xiàn)，北京：電子工業(yè)出版社，2002年2月.</p><p>　　[4]《信號分析和處理MATLAB語言及應(yīng)用》國防科技大學(xué)出版社黃文梅楊勇2000年</p>

81、;<p><b>　　附錄</b></p><p>　　MATLAB程序清單：</p><p><b>　　basic1.m</b></p><p><b>　　basic2.m</b></p><p><b>　　basic3.m</b>&

82、lt;/p><p><b>　　basic4.m</b></p><p><b>　　basic5.m</b></p><p>　　enhanced7.m</p><p><b>　　輔助文件：</b></p><p><b>　　msg.wav&