matlab信號處理實習(xí)報告

上傳人：奔*** IP屬地：河北更新時間：2024-03-02 格式：doc 頁數(shù)：23 大?。?64.50KB 人氣指數(shù)：12 舉報 版權(quán)申訴

已閱讀1頁，還剩22頁未讀，繼續(xù)免費(fèi)閱讀

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、　　MATLAB信號處理實習(xí)報告　　實習(xí)名稱 MATLAB信號處理 　　實習(xí)題目數(shù)字信號處理 　　專業(yè)班級 　　姓名 </p&g

2、t;　　學(xué) 號 　　成績評定　　電氣與信息學(xué)院　　二0一一年六月　　實習(xí)考核和成績評定辦法&

3、lt;/p>　　1、實習(xí)的考核由指導(dǎo)教師根據(jù)設(shè)計表現(xiàn)（出勤、遵守紀(jì)律情況等）、設(shè)計報告、設(shè)計成果、答辯等幾個方面，給出各項成績或權(quán)重，綜合后給出實習(xí)總成績。該設(shè)計考核須經(jīng)教研室主任審核，主管院長審批備案。　　2、成績評定采用五級分制，即優(yōu)、良、中、及格和不及格。　　3、參加本次設(shè)計時間不足三分之二或曠課三天以上者，不得參加本次

4、考核，按不及格處理。　　4、實習(xí)結(jié)束一周內(nèi)，指導(dǎo)教師提交成績和設(shè)計總結(jié)。　　5、設(shè)計過程考核和成績在教師手冊中有記載。　　注意： 　　1、實習(xí)任務(wù)書和指導(dǎo)書在實習(xí)前發(fā)給學(xué)生，設(shè)計任務(wù)書放置在設(shè)計報告封面后和正文目錄前。</p&g

5、t;　　2、為了節(jié)省紙張，保護(hù)環(huán)境，便于保管設(shè)計報告，統(tǒng)一采用A4紙，實習(xí)報告建議雙面打印（正文采用宋體五號字）或手寫，左側(cè)裝訂。　　目錄　　一、實習(xí)目的- 4 -　　二、實習(xí)任務(wù)及功能概述- 4 -　　三、

6、相關(guān)設(shè)計內(nèi)容原理及其實現(xiàn)之的MATLAB函數(shù)- 5 -　　四、調(diào)試運(yùn)行代碼，得到仿真圖- 8 -　　五、實習(xí)總結(jié)- 22 -　　六、參考文獻(xiàn)- 22 -　　MATALB信號處理實習(xí)報告　　一、實習(xí)目的&l

7、t;/b>　　“數(shù)字信號處理”是電子信息工程專業(yè)的主干課程，其理論性較強(qiáng)，學(xué)生通過理論課的系統(tǒng)學(xué)習(xí)后，應(yīng)通過MATLAB語言對其所涉及的算法進(jìn)行仿真，這不僅能幫助學(xué)生理解其抽象的物理概念，工程概念和復(fù)雜算法，加強(qiáng)感性認(rèn)識，而且能激發(fā)學(xué)生更進(jìn)一步地在該領(lǐng)域的學(xué)習(xí)和探索熱情。　　二、實習(xí)任務(wù)及功能概述　　1.

8、FT與DFT關(guān)系的MATLAB仿真　　任務(wù)要求：　　任設(shè)一個長度為N的實序列；　　編寫MATLAB代碼，實現(xiàn)該序列的傅立葉變換和離散傅立葉變換；　　調(diào)試運(yùn)行代碼，得到仿真圖形；　　（4）通過仿

9、真圖形觀察二者的關(guān)系，并得出結(jié)論。　　DFT的快速算法FFT的仿真　　任務(wù)要求：　　在熟悉FFT算法原理的基礎(chǔ)上，對所給的時域信號進(jìn)行分析；　　對信號做時域分析時要注意混疊現(xiàn)象的發(fā)生；　　對混

10、合信號做快速傅立葉變換FFT運(yùn)算，其中自設(shè)，變換點數(shù)根據(jù)所設(shè)的情況自定；　?。?）對變換結(jié)果做分析。　　3. 利用MATLAB信號處理工具箱設(shè)計FIR濾波器　　要求如下： 　　(1) 熟悉FIR濾波器的工作原理；

11、　　(2) 設(shè)計一個N階高通數(shù)字濾波器，階數(shù)和截止頻率自定；　　(3) 設(shè)計一個N階帶通數(shù)字濾波器，階數(shù)和截止頻率自定； 　　設(shè)計一個N階低通數(shù)字濾波器，階數(shù)和截止頻率自定；　　有完整的MATLAB代碼。 　　4. 利用MATLAB信號處理工具箱設(shè)計

12、IIR濾波器　　要求如下：　　熟悉IIR數(shù)字濾波器的工作原理；　　設(shè)計一個N階高通數(shù)字濾波器，階數(shù)和截止頻率自定；　　設(shè)計一個N階帶通數(shù)字濾波器，階數(shù)和截止頻率自定；　　設(shè)計一個N階高通數(shù)字濾

13、波器，階數(shù)和截止頻率自定；　　有完整的MATLAB代碼。　　5. 利用所設(shè)計的濾波器對混合信號進(jìn)行濾波　　要求如下:　　（1）設(shè)計一個IIR或FIR數(shù)字濾波器，將混合信號中的、或選擇出來，、、自定，濾波器類型自定

14、；　?。?）有完整的MATALB代碼；　　（3）仿真結(jié)果；　　（4）結(jié)果分析。　　三、相關(guān)設(shè)計內(nèi)容原理及其實現(xiàn)之的MATLAB函數(shù)　　(1) 序列的傅立葉

15、變換FT、離散傅立葉變換DFT和快速傅立葉變換FFT原理　　FFT算法基本上分為兩大類：時域抽取法FFT(Decimation In Time FFT,簡稱DIT-FFT)和頻域抽取法FFT(Decimation In Frequency FFT,簡稱DIF―FFT)。　　在MATLAB中，DFT的快速算法FFT的實現(xiàn)函數(shù)是fft(x,n)，這里x是離散

16、序列，n是變換點數(shù)，n的選擇由信號的最高頻率和頻率分辨率決定，即n>2fh/F。　　FIR濾波器的工作原理　　設(shè)希望設(shè)計的濾波器傳輸函數(shù)為Hd(ejω),hd(n)是與其對應(yīng)的單位脈沖響應(yīng)，因此 　　為了構(gòu)造一個長度為N的線性相位濾波器，只有將hd(n)截取一段，并保證截取的一段對(N-1)/2對稱。設(shè)截取

17、的一段用h(n)表示，即　　h(n)=hd(n)RN(n) 　　我們實際實現(xiàn)的濾波器的單位取樣響應(yīng)為h(n),長度為N，其系統(tǒng)函數(shù)為H(z)，　　在MATLAB中，實現(xiàn)FIR濾波的函數(shù)是　　-----濾波器的階

18、數(shù);　　----- 濾波器的截止頻率，若是帶通或帶阻濾波器，和分別代表通帶（阻帶）上、下限截止頻率，且這里的和均為歸一化截止頻率。　　提示：Matlab信號處理工具箱濾波器設(shè)計中的截止頻率均采用歸一化頻率（對采樣頻率的1/2做歸一化處理），歸一化頻率轉(zhuǎn)換為角頻率，則將歸一化頻率乘以。　　----- 濾波器的類型，

19、即低通（），高通（），帶阻（），帶通（），　　默認(rèn)值為低通或帶阻。 　　-----窗函數(shù)類型，即矩形窗（）、漢寧窗（）、巴特利窗（）、布萊克曼窗（）、海明窗（）等。默認(rèn)值為海明窗（）。　　-----返回濾波器的系數(shù)。　　濾波器的系數(shù)得到后，就可以做它的頻率響應(yīng)圖了。M

20、atlab提供了數(shù)字濾波器的頻率響應(yīng)函數(shù)，其基本調(diào)用格式為：　　-----代表濾波器傳遞函數(shù)的分子多項式系數(shù)。　　-----代表濾波器傳遞函數(shù)的分目多項式系數(shù)。　　-----代表復(fù)頻率響應(yīng)的采樣點數(shù)。　　-----濾波器的復(fù)頻率響應(yīng)。<p&g

21、t;　　-----頻率點，該是將歸一化后的頻率乘以的角頻率值。　　IIR濾波器的工作原理　　IIR濾波器設(shè)計方法有兩類，經(jīng)常用的一類設(shè)計方法是借助于模擬濾波器的設(shè)計方法進(jìn)行的。其設(shè)計步驟是：先設(shè)計模擬濾波器得到傳輸函數(shù)Ha(s)，然后將Ha(s)按某種方法轉(zhuǎn)換成數(shù)字濾波器的系統(tǒng)函數(shù)H(z)。模擬濾波器的理論和設(shè)計方法已發(fā)展得相當(dāng)成熟，且有若干典型的模擬濾

22、波器供我們選擇，如巴特沃斯(Butterworth)濾波器、切比雪夫(Chebyshev)濾波器、橢圓(Cauer)濾波器、貝塞爾(Bessel)濾波器等，這些濾波器都有嚴(yán)格的設(shè)計公式、現(xiàn)成的曲線和圖表供設(shè)計人員使用。 　　模擬低通濾波器的設(shè)計指標(biāo)及逼近方法　　模擬低通濾波器的設(shè)計指標(biāo)有αp, Ωp,αs和Ωs。其中Ωp和Ωs分別稱為通帶截止頻率和阻帶截

23、止頻率，αp是通帶Ω(=0~Ωp)中的最大衰減系數(shù)，αs是阻帶Ω≥Ωs的最小衰減系數(shù)，αp和αs一般用dB數(shù)表示。對于單調(diào)下降的幅度特性，可表示成：　　如果Ω=0處幅度已歸一化到1，即|Ha(j0)|=1,αp和αs表示為　　Ωc稱為3dB截止頻率　　濾波器的技術(shù)指標(biāo)給定后，需要設(shè)計一個傳輸函數(shù)Ha(s)，希望其

24、幅度平方函數(shù)滿足給定的指標(biāo)αp和αs，一般濾波器的單位沖激響應(yīng)為實數(shù)，因此　　在MATALB中，實現(xiàn)IIR濾波的函數(shù)有切比雪夫濾波器cheby和橢圓濾波器ellip。這里我們只選用橢圓濾波器，橢圓濾波器的調(diào)用格式為：；　　濾波器的分子多項式系數(shù)　　濾波器的分母多項式系數(shù)

25、　　濾波器的階數(shù)　　通帶最大衰減　　阻帶最小衰減　　濾波器的類型，其取值為代表高通濾波器，代表低通濾波器，　　代表帶通濾波器，代表帶阻濾波器。</p&g

26、t;　　MATLAB的濾波函數(shù)　　在MATLAB中，實現(xiàn)濾波的函數(shù)是filter(b,a,x)。這里，b和a分別是濾波器傳遞函數(shù)分子和分母多項式的系數(shù),x是被濾波的信號。濾波函數(shù)的調(diào)用格式為：　　-----代表濾波器傳遞函數(shù)的分子多項式系數(shù)。　　-----代表濾波器傳遞函數(shù)的分目多項式系數(shù)。

27、　　-----代表被濾波的信號。　　四、調(diào)試運(yùn)行代碼，得到仿真圖　　1. FT與DFT關(guān)系的MATLAB仿真　　任務(wù)要求：　　任設(shè)一個長度為N的實序列；

28、;　　編寫MATLAB代碼，實現(xiàn)該序列的傅立葉變換和離散傅立葉變換；　　調(diào)試運(yùn)行代碼，得到仿真圖形；　　通過仿真圖形觀察二者的關(guān)系，并得出結(jié)論。　　function [X,Xk]=mydft(x,N)　　w=2*pi*(0:127)/128;<

29、;p>　　X=x*exp(-j*[1:length(x)]'*w) 　　subplot(211)　　plot(w,abs(X))　　title('電子0841 張憲')　　k=

30、0:N-1;　　XK=x*exp(-j*[1:length(x)]'*(2*pi*k)/N) 　　subplot(212)　　stem(k,abs(XK))　　title('FT與DFT關(guān)系的MATLA

31、B仿真')　　end　　圖1-1　　DFT的快速算法FFT的仿真　　任務(wù)要求：　　在熟悉FFT算法原理的基礎(chǔ)

32、上，對所給的時域信號進(jìn)行分析；　　對信號做時域分析時要注意混疊現(xiàn)象的發(fā)生；　　對混合信號做快速傅立葉變換FFT運(yùn)算，其中、、自設(shè)，變換點數(shù)根據(jù)所設(shè)的情況自定；　?。?）對變換結(jié)果做分析。　　function [x,xfft]=myfft(f1,f2,f3,fs

33、,n)　　t=0:1/fs:1;　　x=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);　　subplot(211)　　plot(t,x)　　titl

34、e('DFT的快速算法FFT的仿真')　　subplot(212)　　xfft=fft(x,n)　　f=fs*(0:n/2-1)/n;　　plot(f,abs(xfft(1:n/2)))　　title('電

35、子0841 張憲')　　end　　f1=20　　f2=40　　f3=60

36、;　　fs=120　　n=64　　3. 利用MATLAB信號處理工具箱設(shè)計FIR濾波器　　要求如下： 　　(1) 熟悉FIR濾波器的工作原理；</p&g

37、t;　　(2) 設(shè)計一個N階高通數(shù)字濾波器，階數(shù)和截止頻率自定；　　(3) 設(shè)計一個N階帶通數(shù)字濾波器，階數(shù)和截止頻率自定； 　　設(shè)計一個N階低通數(shù)字濾波器，階數(shù)和截止頻率自定；　　有完整的MATLAB代碼　　帶通：</b&g

38、t;　　function [h,s,sf]=myFIRbandpassfilter(f1,f2,f3,fs)　　m=(0.3*f1)/(fs/2); 　　M=round(8/m); 　　N=M-1;

39、 　　b=fir1(N,[0.8*f2/(fs/2) 1.2*f2/(fs/2)]); 　　[h,f]=freqz(b,1,512);　　axes('position',[0.1 0.76 0.33 .2])　　plot(f*fs/(2*pi),20*l

40、og10(abs(h)))　　xlabel('frequency/Hz','fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,'fontweight','bold')

41、;　　t=0:1/fs:1;　　s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);　　axes('position',[0.1 0.43 0.33 .2])　　plot(t,s);　　xla

42、bel('time/Second','fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,'fontweight','bold')　　sf=filter(b,1,s);

43、;　　axes('position',[0.1 0.1 0.33 .2])　　plot(t,sf)　　axis([0.2 .5 -2 2])　　xlabel('time/Second','fontsize',8,'fontweight&

44、#39;,'bold')　　set(gca,'fontsize',8,'fontweight','bold')　　end　　f1=10

45、;　　f2=20　　f3=30　　fs=60　　低通：　　function [h,s,sf]=myFIRlowpassf

46、ilter(f1,f2,fs)　　m=(0.3*f1)/(fs/2); 　　M=round(8/m); 　　N=M-1; 　　b=fir1(N,0.9*f2/(fs/2));

47、　　[h,f]=freqz(b,1,512); 　　plot(f*fs/(2*pi),20*log10(abs(h))) 　　xlabel('frequency/Hz','fontsize',8,'fontweight','bold')&l

48、t;p>　　set(gca,'fontsize',8,'fontweight','bold')　　figure　　subplot(211) 　　t=0:1/fs:1;　　s=si

49、n(2*pi*f1*t)+sin(2*pi*f2*t); 　　plot(t,s); 　　xlabel('time/Second','fontsize',8,'fontweight','bold')

50、　set(gca,'fontsize',8,'fontweight','bold')　　subplot(212)　　SK=fft(s,512); 　　SXK=abs(SK);<

51、p>　　f=(0:255)*fs/512;　　plot(f,SXK(1:256)) 　　xlabel('Frequency/Hz','fontsize',8,'fontweight','bold')　　set

52、(gca,'fontsize',8,'fontweight','bold')　　figure 　　sf=filter(b,1,s); % filter s　　subplot(211)</

53、p>　　plot(t,sf) % draw filtered signal sf　　xlabel('time/Second','fontsize',8,'fontweight','bold')　　set(gca,&#

54、39;fontsize',8,'fontweight','bold')　　axis([0.2 .5 -2 2])　　subplot(212)　　SFK=fft(sf,512);　　SFXK=abs(SFK);

55、　　f=(0:255)*fs/512;　　plot(f,SFXK(1:256)) 　　xlabel('Frequency/Hz','fontsize',8,'fontweight','bold')

56、;　　set(gca,'fontsize',8,'fontweight','bold')　　end　　高通：　　function [h,s,sf]=myFIRbandpassfilter(f

57、1,f2,f3,fs)　　m=(0.3*f1)/(fs/2); 　　M=round(8/m); 　　N=M-1; 　　b=fir1(N,[0.8*f2/(fs/2) 1.2*f2/(fs/2)]),

58、9;high'; 　　[h,f]=freqz(b,1,512);　　axes('position',[0.1 0.76 0.33 .2])　　plot(f*fs/(2*pi),20*log10(abs(h)))weight','bold')&

59、lt;p>　　set(gca,'fontsize',8,'fontweight','bold')　　t=0:1/fs:1;　　s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);　　axes('position

60、',[0.1 0.43 0.33 .2])　　plot(t,s);　　xlabel('time/Second','fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,&

61、#39;fontweight','bold')　　sf=filter(b,1,s);　　axes('position',[0.1 0.1 0.33 .2])　　plot(t,sf)　　axis([0.2 .5 -2 2])</

62、p>　　xlabel('time/Second','fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,'fontweight','bold')<b&g

63、t;　　end　　xlabel('frequency/Hz','fontsize',8,'font　　4. 利用MATLAB信號處理工具箱設(shè)計IIR濾波器　　要求如下：　　熟

64、悉IIR數(shù)字濾波器的工作原理；　　設(shè)計一個N階高通數(shù)字濾波器，階數(shù)和截止頻率自定；　　設(shè)計一個N階帶通數(shù)字濾波器，階數(shù)和截止頻率自定；　　設(shè)計一個N階高通數(shù)字濾波器，階數(shù)和截止頻率自定；　　有完整的MATLAB代碼。<b&g

65、t;　　低通：　　function [x,m,xf,mf]=mylowIIRfilter(f1,f2,fs,n)　　subplot(411)　　t=0:1/fs:.5;　　x=sin(2*pi*f1*t)+sin(2*pi*f2*t);</p&g

66、t;　　plot(t,x)　　subplot(412)　　y=fft(x,n);　　m=abs(y)/max(abs(y));　　f=fs*(0:n/2-1)/n;　　p

67、lot(f,m(1:n/2));　　[b,a]=ellip(6,.1,40,.5*(f1+f2)/(fs/2));　　xf=filter(b,a,x);　　subplot(413)　　plot(t,xf);　　subplot(41

68、4)　　yf=fft(xf,n);　　mf=abs(yf)/max(abs(yf));　　f=fs*(0:n/2-1)/n;　　plot(f,mf(1:n/2));　　end

69、;　　f1=100　　f2=200　　fs=400　　n=32　　高通：<

70、/b>　　function [x,m,xf,mf]=mylowIIRfilter(f1,f2,fs,n)　　subplot(411)　　t=0:1/fs:.5;　　x=sin(2*pi*f1*t)+sin(2*pi*f2*t);<p&g

71、t;　　plot(t,x)　　subplot(412)　　y=fft(x,n);　　m=abs(y)/max(abs(y));　　f=fs*(0:n/2-1)/n;　　plot(f,m(1:n

72、/2));　　[b,a]=ellip(6,.1,40,.5*(f1+f2)/(fs/2),'high');　　xf=filter(b,a,x);　　subplot(413)　　plot(t,xf);　　subplo

73、t(414)　　yf=fft(xf,n);　　mf=abs(yf)/max(abs(yf));　　f=fs*(0:n/2-1)/n;　　plot(f,mf(1:n/2));　　end</

74、p>　　f1=100,　　f2=200,　　fs=400,　　n=32

75、帶通　　function [x,m,xf,mf]=mypassIIRfilter(f1,f2,fs,n)　　subplot(411)　　t=0:1/fs:.5;　　x=sin(2*pi*f1*t)+sin(2*pi*f2*t);

76、　　plot(t,x)　　subplot(412)　　y=fft(x,n);　　m=abs(y)/max(abs(y));　　f=fs*(0:n/2-1)/n;　　plot(

77、f,m(1:n/2));　　[b,a]=ellip(6,.1,40,[0.8*f2/(fs/2) 1.2*f2/(fs/2)],'high');　　xf=filter(b,a,x);　　subplot(413)　　plot(t,xf);

78、　　subplot(414)　　yf=fft(xf,n);　　mf=abs(yf)/max(abs(yf));　　f=fs*(0:n/2-1)/n;　　plot(f,mf(1:n/2));

79、　end　　圖1-2　　1.5使用function定義一個帶通FIR濾波器ftxuanpinfilter　　s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);通過帶通分別選出三個頻率f1,f2,f3.</p

80、>　　function [h,s,sf]=ftxuanpinfilter(f1,f2,f3,fs)　　m=(0.3*f1)/(fs/2); % define transition bandwidth 　　M=round(8/m); % define length of hamming

81、window　　N=M-1; % define order of FIR filter　　b=fir1(N,[0.8*f1/(fs/2) 1.2*f1/(fs/2)]); % 上下限截止頻率分別取f1的0.8和1.2　　[h,f]=freqz(b,1,512);</

82、p>　　axes('position',[0.1 0.73 0.8 .2])　　plot(f*fs/(2*pi),20*log10(abs(h)))　　title('濾波器頻譜圖')　　xlabel('frequency/Hz','

83、;fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,'fontweight','bold')　　t=0:1/fs:1;　　s=sin(2*pi*f1*t)+sin(2

84、*pi*f2*t)+sin(2*pi*f3*t);　　axes('position',[0.1 0.4 0.8 .2])　　plot(t,s);　　title('原始信號f(t)的波形圖')　　xlabel('time/Second

85、','fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,'fontweight','bold')　　sf=filter(b,1,s);　　axes(&

86、#39;position',[0.1 0.08 0.8 .2])　　plot(t,sf)　　title('選出f(t)中的f1的波形圖')　　axis([0.2 .5 -2 2])　　xlabel('time/Second',

87、9;fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,'fontweight','bold')　　figure　　b1=fir1

88、(N,[0.8*f2/(fs/2) 1.2*f2/(fs/2)]); % 上下限截止頻率分別取f2的0.8和1.2　　[h1,f1]=freqz(b1,1,512);　　axes('position',[0.1 0.73 0.8 .2])　　plot(f1*fs/(2*pi),20*log10(

89、abs(h)))　　title('濾波器頻譜圖')　　xlabel('frequency/Hz','fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,'

90、;fontweight','bold')　　axes('position',[0.1 0.4 0.8 .2])　　plot(t,s);　　title('原始信號f(t)的波形圖')　　xlabel('time

91、/Second','fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,'fontweight','bold')　　sf=filter(b1,1,s);

92、　　axes('position',[0.1 0.08 0.8 .2])　　plot(t,sf)　　title('選出f(t)中的f2的波形圖')　　axis([0.2 .5 -2 2])　　xlabel('time/Second&

93、#39;,'fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,'fontweight','bold')　　figure

94、　b2=fir1(N,[0.8*f3/(fs/2) 1.2*f3/(fs/2)]); % 上下限截止頻率分別取f3的0.8和1.2　　[h2,f2]=freqz(b2,1,512);　　axes('position',[0.1 0.73 0.8 .2])　　plot(f2*fs/(2*pi),2

95、0*log10(abs(h)))　　title('濾波器頻譜圖')　　xlabel('frequency/Hz','fontsize',8,'fontweight','bold')　　set(gca,'fontsize'

96、;,8,'fontweight','bold')　　axes('position',[0.1 0.4 0.8 .2])　　plot(t,s);　　title('原始信號f(t)的波形圖')　　xlabel(&

97、#39;time/Second','fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,'fontweight','bold')　　sf=filter(b2,1,s);&

98、lt;p>　　axes('position',[0.1 0.08 0.8 .2])　　plot(t,sf)　　title('選出f(t)中的f3的波形圖')　　axis([0.2 .5 -2 2])　　xlabel('time

99、/Second','fontsize',8,'fontweight','bold')　　set(gca,'fontsize',8,'fontweight','bold')　　end<p

100、>　　再在命令窗口調(diào)用此函數(shù)結(jié)果如下　　f1=30　　f2=60　　f3=130　　fs=500

101、;　　ftzhangxianpin(f1,f2,f3,fs);　　濾出f1：　　五、實習(xí)總結(jié)　　本次實習(xí)的內(nèi)容為matlab與數(shù)字信號處理，說實話，matlab和數(shù)字信號處理都是好久前學(xué)的了，多少有些生疏，我把老師教

102、過的課本認(rèn)認(rèn)真真的反復(fù)研讀，理解內(nèi)容為主，主要以課本為主，因為貪多嚼不爛。到實習(xí)的最后三天，雖然高燒，喉嚨疼痛，但是我還是沒有放棄matlab與數(shù)字信號處理的學(xué)習(xí)，我希望能好好掌握matlab這門實用課程，因為它不僅在我們的電子領(lǐng)域還是其他方面都有著廣泛的應(yīng)用。最后，感謝老師的指導(dǎo)和同學(xué)的幫助，使我能在本周的實習(xí)中圓滿完成任務(wù)。　　六、參考文獻(xiàn)</p&

眾賞文庫> 全部分類> 畢業(yè)設(shè)計

溫馨提示

1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
5. 眾賞文庫僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

matlab信號處理實習(xí)報告

文檔簡介

溫馨提示

最新文檔

評論

matlab信號處理實習(xí)報告

文檔簡介

溫馨提示

最新文檔

評論

免費(fèi)下載