版權(quán)說(shuō)明:本文檔由用戶(hù)提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、<p> 課 程 設(shè) 計(jì) 報(bào) 告</p><p> 課程名稱(chēng) 信號(hào)與系統(tǒng)課程設(shè)計(jì) </p><p> 指導(dǎo)教師 </p><p><b> 設(shè)計(jì)起止日期 </b></p><p> 學(xué) 院 信息與通信工程 </p>
2、;<p> 專(zhuān) 業(yè) 電子信息工程 </p><p> 學(xué)生姓名 </p><p> 班級(jí)/學(xué)號(hào) </p><p> 成 績(jī) </p><p> 指導(dǎo)老師簽字
3、 </p><p><b> 目 錄</b></p><p> 1、課程設(shè)計(jì)目的1</p><p> 2、課程設(shè)計(jì)要求1</p><p> 3、課程設(shè)計(jì)任務(wù)1</p><p> 4、課程設(shè)計(jì)內(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è)計(jì)目的和要求</p><p><b> 目的:</b></p><p> 掌握MATLAB
5、的使用方法?!靶盘?hào)與系統(tǒng)”是一門(mén)重要的專(zhuān)業(yè)基礎(chǔ)課,MATLAB作為信號(hào)處理強(qiáng)有力的計(jì)算和分析工具是電子信息工程技術(shù)人員常用的重要工具之一。本課程設(shè)計(jì)基于MATLAB完成信號(hào)與系統(tǒng)綜合設(shè)計(jì)實(shí)驗(yàn),以提高學(xué)生的綜合應(yīng)用知識(shí)能力為目標(biāo),是“信號(hào)與系統(tǒng)”課程在實(shí)踐教學(xué)環(huán)節(jié)上的必要補(bǔ)充。通過(guò)課設(shè)綜合設(shè)計(jì)實(shí)驗(yàn),激發(fā)學(xué)生理論課程學(xué)習(xí)興趣,提高分析問(wèn)題和解決問(wèn)題的能力。</p><p><b> 要求: </
6、b></p><p> ?。?)運(yùn)用MATLAB編程得到簡(jiǎn)單信號(hào)、簡(jiǎn)單信號(hào)運(yùn)算、復(fù)雜信號(hào)的頻域響應(yīng)圖;</p><p> (2)通過(guò)對(duì)線(xiàn)性時(shí)不變系統(tǒng)的輸入、輸出信號(hào)的時(shí)域和頻域的分析,了解線(xiàn)性時(shí)不變系統(tǒng)的特性,同時(shí)加深對(duì)信號(hào)頻譜的理解。</p><p> 2、課程設(shè)計(jì)原理</p><p> 1離散時(shí)間信號(hào)卷積使用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)對(duì)輸入x進(jìn)行濾波,系統(tǒng)的輸出為y。 </p><p> 3.在MATLAB中,zplane函數(shù)——繪制系統(tǒng)零極點(diǎn)圖 </p>
8、;<p> 4.在MATLAB中,有專(zhuān)門(mén)的frepz()函數(shù)來(lái)求解。調(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的冪次方,以便提高計(jì)算速度。 </p><p> 5.MATLAB
9、提供了如abs(),angle()等基本函數(shù)用來(lái)計(jì)算DTFT幅度,相位。</p><p> 3、課程設(shè)計(jì)所用設(shè)備</p><p><b> ?。?)個(gè)人電腦一臺(tái)</b></p><p> (2)MATLAB2012設(shè)計(jì)程序</p><p> 4、課程設(shè)計(jì)內(nèi)容和步驟</p><p> ?。?)
10、根據(jù)設(shè)計(jì)題目的要求,熟悉相關(guān)內(nèi)容的理論基礎(chǔ),理清程序設(shè)計(jì)的措施和步驟;</p><p> ?。?)根據(jù)設(shè)計(jì)題目的要求,提出各目標(biāo)的實(shí)施思路、方法和步驟;</p><p> ?。?)根據(jù)相關(guān)步驟完成MATLAB程序設(shè)計(jì),所編程序應(yīng)能完整實(shí)現(xiàn)設(shè)計(jì)題目的要求;</p><p> (4)調(diào)試程序,分析相關(guān)理論;</p><p> (5)編寫(xiě)設(shè)計(jì)報(bào)
11、告。</p><p> 5、設(shè)計(jì)過(guò)程、實(shí)現(xiàn)過(guò)程、實(shí)現(xiàn)結(jié)果</p><p><b> ?。ㄒ唬┗静糠?lt;/b></p><p> ?。?)信號(hào)的時(shí)頻分析</p><p> 任意給定單頻周期信號(hào)的振幅、頻率和初相,要求準(zhǔn)確計(jì)算出其幅度譜,并準(zhǔn)確畫(huà)出時(shí)域和頻域波形,正確顯示時(shí)間和頻率。</p><p&g
12、t;<b> 設(shè)計(jì)思路:</b></p><p> 首先給出橫坐標(biāo),即時(shí)間,根據(jù)設(shè)定的信號(hào)的振幅、頻率和初相,寫(xiě)出時(shí)域波形的表達(dá)式;然后對(duì)時(shí)域波形信號(hào)進(jìn)行傅里葉變化,得到頻域波形;最后使用plot函數(shù)繪制各個(gè)響應(yīng)圖。</p><p><b> 源程序:</b></p><p> %基本題1,時(shí)頻分析,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('時(shí)域波形圖')</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('信號(hào)幅度譜')</p><p> subplot(224),plot(f,phaseX)</p><p> xlabel('f(Hz)'),ylabel('angle[X(jw)]'),title('信號(hào)相位譜')</p>&
17、lt;p><b> 結(jié)果圖:</b></p><p><b> 結(jié)果分析:</b></p><p> cos函數(shù)波形為周期信號(hào),其頻域響應(yīng)為兩個(gè)沖擊函數(shù),并且符合對(duì)偶性。</p><p> ?。?)傅里葉級(jí)數(shù)分析</p><p> 分析周期鋸齒波的傅里葉級(jí)數(shù)系數(shù),用正弦信號(hào)的線(xiàn)性組合構(gòu)
18、成鋸齒波,要求諧波次數(shù)可以任意輸入,分析不同諧波次數(shù)所構(gòu)成的鋸齒波,解釋是否存在吉伯斯現(xiàn)象。</p><p><b> 設(shè)計(jì)思路:</b></p><p> 根據(jù)原理——任意周期信號(hào)都可以表達(dá)成傅里葉級(jí)數(shù)的形式,對(duì)周期舉行鋸齒波進(jìn)行傅里葉級(jí)數(shù)分解,利用for循環(huán)完成級(jí)數(shù)求和運(yùn)算。當(dāng)N值較小時(shí),傅里葉級(jí)數(shù)的逼近效果不是很理想,隨著N值變大,傅里葉級(jí)數(shù)越來(lái)越接近理想值
19、。</p><p><b> 源程序:</b></p><p> %基本題2,傅里葉級(jí)數(shù)分析,2011010512</p><p> clc; clear; close all; %初始化工作空間</p><p> t=-6:6/1000:6; %時(shí)間</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值的增大,傅里葉級(jí)數(shù)與理想的周期鋸齒波越來(lái)越接近,符合預(yù)期效果。</p><p><b> ?。?)系統(tǒng)分析</b></p><p> 任意給定微分方程或差分方程描述的系統(tǒng),畫(huà)出系統(tǒng)的幅頻響應(yīng)和相頻響應(yīng)。</p><p><b> 設(shè)計(jì)思路:</b>&l
24、t;/p><p> 根據(jù)微分方程、差分方程與系統(tǒng)函數(shù)的對(duì)應(yīng)關(guān)系,結(jié)合matlab自帶的freqs和freqz兩個(gè)函數(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ù)時(shí)間
27、系統(tǒng)頻率響應(yīng)</p><p><b> figure(1)</b></p><p> subplot(2,1,1);zplane(b,a);title('連續(xù)時(shí)間系統(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ù)時(shí)間系統(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ù)時(shí)間系統(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); %求離散時(shí)間系統(tǒng)頻率響應(yīng)</p><p><b> figure(2)</b></p><p> subplot(2,1,1);zplane(b,a);title('離散時(shí)間系統(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('離散時(shí)間系統(tǒng)幅頻響應(yīng)');</p><p> subplot(2,2,4);plot(w/pi,angle(Hw))</p><p> xlabel('w
32、/pi');ylabel('Phase');title('離散時(shí)間系統(tǒng)相頻響應(yīng)');</p><p><b> 結(jié)果圖;</b></p><p> 結(jié)果分析:所顯示零極點(diǎn)等圖符合驗(yàn)算結(jié)果</p><p> (4)音樂(lè)合成程序設(shè)計(jì)</p><p> 對(duì)于任意一小段音樂(lè),利用“
33、十二平均律”計(jì)算該音樂(lè)中各個(gè)樂(lè)音的頻率,產(chǎn)生并播放這些樂(lè)音。分析音樂(lè)的頻譜,從中識(shí)別出不同的樂(lè)音。</p><p><b> 設(shè)計(jì)思路:</b></p><p> 根據(jù)時(shí)間長(zhǎng)短來(lái)區(qū)別各個(gè)音符拍長(zhǎng)短,根據(jù)頻率高低來(lái)區(qū)別各個(gè)音符的音調(diào),具體細(xì)節(jié)參考“十二平均律”,將每段樂(lè)音連接起來(lái),使用sound函數(shù)播放樂(lè)音。對(duì)每段樂(lè)音進(jìn)行傅里葉變換,分析其頻譜。</p>
34、<p><b> 源程序:</b></p><p> %基本題4,音樂(lè)合成程序設(shè)計(jì),一分錢(qián),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é)果分心:所聽(tīng)音樂(lè)符合預(yù)期。</p><p><b> ?。?)調(diào)制</b></p><p> 分析單位沖激響應(yīng)為的系統(tǒng)的濾波特性,畫(huà)出其幅頻響應(yīng)
39、曲線(xiàn)。</p><p><b> 設(shè)計(jì)思路:</b></p><p> 根據(jù)傅里葉變換特性,時(shí)域相乘對(duì)應(yīng)頻域卷積,將題目中的時(shí)域波形進(jìn)行分解,并將兩個(gè)函數(shù)分別對(duì)應(yīng)的傅里葉變換進(jìn)行卷積。參考常用的傅里葉變換對(duì),所求的頻域特性為方波和沖擊函數(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é)果符合演算過(guò)程,與
44、預(yù)期一致。</p><p><b> (二)提高部分</b></p><p> 題目1、音樂(lè)合成 2011010512</p><p> 設(shè)計(jì)思路:通過(guò)定義每個(gè)音來(lái)寫(xiě)樂(lè)譜,在加包絡(luò)諧波,在播放音樂(lè),后畫(huà)出波形。 </p><p><b> 源程序</b></p><p&
45、gt;<b> %提高題 音樂(lè)合成</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> %% 簡(jiǎn)譜 歡樂(lè)頌</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個(gè)音符(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('時(shí)間/s');ylabel('幅值');title('普通 無(wú)包絡(luò) 無(wú)諧波 波形圖');</p
69、><p> subplot(312);</p><p> plot((1:length_num_song)/8000,wav_envelope_only(1:length_num_song));</p><p> xlabel('時(shí)間/s');ylabel('幅值');title('加包絡(luò) 無(wú)諧波 波形圖');<
70、;/p><p> subplot(313);</p><p> plot((1:length_num_song)/8000,wav_envelope_harmonic(1:length_num_song));</p><p> xlabel('時(shí)間/s');ylabel('幅值');title('加包絡(luò) 加諧波 波形圖
71、9;);</p><p> %% 抽取前13個(gè)音符(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('普通音無(wú)諧波頻譜圖');</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個(gè)音符(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、致 聽(tīng)到音樂(lè)也和預(yù)期一致。</p><p><b> 6、設(shè)計(jì)總結(jié)</b></p><p> 通過(guò)這次試驗(yàn),我明白了學(xué)到的知識(shí)只有會(huì)用才有價(jià)值,只有在應(yīng)用中才能找到自己的不足并加以彌補(bǔ),如果我們不去將知識(shí)應(yīng)用,那我們學(xué)的東西將失去價(jià)值。</p><p> 在設(shè)計(jì)過(guò)程中,自己學(xué)會(huì)了Matlab軟件的部分功能,在今后的生活中我希望能夠?qū)W到
79、更多!</p><p> 在設(shè)計(jì)工程中我遇到了很多困難,但是我堅(jiān)持作了下來(lái),這對(duì)我來(lái)說(shuō)是一次磨練,令我更加堅(jiān)強(qiáng)。在今后的生活中,當(dāng)我們遇到困難我們不應(yīng)該輕言放棄而是應(yīng)該堅(jiān)持,努力實(shí)現(xiàn)自己的目標(biāo)。</p><p><b> 參考文獻(xiàn)</b></p><p> [1] 鄭君里、應(yīng)啟珩 、楊為理,信號(hào)與系統(tǒng)引論,北京:高等教育出版社,2009年
80、3月.</p><p> [2] 谷源濤、應(yīng)啟珩、鄭君里,信號(hào)與系統(tǒng)——MATLAB綜合實(shí)驗(yàn),北京:高等教育出版社,2008年1月.</p><p> [3] 梁虹等,信號(hào)與系統(tǒng)分析及Matlab實(shí)現(xiàn),北京:電子工業(yè)出版社,2002年2月.</p><p> [4]《信號(hào)分析和處理MATLAB語(yǔ)言及應(yīng)用》國(guó)防科技大學(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&
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 眾賞文庫(kù)僅提供信息存儲(chǔ)空間,僅對(duì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 信號(hào)與系統(tǒng)課程設(shè)計(jì)報(bào)告
- 信號(hào)與系統(tǒng)課程設(shè)計(jì)報(bào)告
- 信號(hào)與線(xiàn)性系統(tǒng)課程設(shè)計(jì)報(bào)告--周期信號(hào)的分解與合成
- 信號(hào)與線(xiàn)性系統(tǒng)課程設(shè)計(jì)報(bào)告
- 信號(hào)與系統(tǒng)課程設(shè)計(jì)
- 信號(hào)與系統(tǒng)課程設(shè)計(jì)
- 信號(hào)與系統(tǒng)課程設(shè)計(jì)
- 信號(hào)與系統(tǒng)課程設(shè)計(jì)報(bào)告---matlab在信號(hào)與系統(tǒng)中的使用
- 電路、信號(hào)與系統(tǒng)課程設(shè)計(jì)
- 課程設(shè)計(jì)----音樂(lè)信號(hào)濾波去噪
- 信號(hào)與線(xiàn)性系統(tǒng)課程設(shè)計(jì)
- 信號(hào)與系統(tǒng)課程設(shè)計(jì)--信號(hào)的采樣與恢復(fù)
- 信號(hào)與系統(tǒng)課程設(shè)計(jì)大綱doc
- 信號(hào)取樣與恢復(fù)課程設(shè)計(jì)報(bào)告
- 音樂(lè)播放系統(tǒng)設(shè)計(jì)--單片機(jī)課程設(shè)計(jì)報(bào)告
- 數(shù)字信號(hào)課程設(shè)計(jì)報(bào)告
- 課程設(shè)計(jì)--駝峰信號(hào)自動(dòng)控制課程設(shè)計(jì)報(bào)告
- 信號(hào)處理分析課程設(shè)計(jì)報(bào)告
- 信號(hào)與系統(tǒng)課程設(shè)計(jì)--信號(hào)發(fā)生器的設(shè)計(jì)與實(shí)現(xiàn)
- 濾波器設(shè)計(jì)與信號(hào)濾波課程設(shè)計(jì)報(bào)告
評(píng)論
0/150
提交評(píng)論