典型序列頻譜分析

1、第1章 設(shè)計(jì)任務(wù)與要求211對(duì)于三種典型序列-單位采樣序列、實(shí)指數(shù)序列、矩形序列21.2自行設(shè)計(jì)一個(gè)周期序列2第2 章 原理及過(guò)程31設(shè)計(jì)原理3第3 章 設(shè)計(jì)內(nèi)容51.1單位采樣序列51.1.1時(shí)域波形51.2傅里葉變換51.3幅度譜及相位譜61.4頻移71.5時(shí)移82.1時(shí)域圖形82.2傅里葉變換92.3幅度譜與相位譜102.4頻移112.5時(shí)移113.1時(shí)域圖形123.2傅里葉變換123.3幅度譜與相位譜133.4時(shí)移143.5頻移154.1幅度特性曲線154.4周期序列的DFS164.5傅里葉變換17第4章 心得與體會(huì)18參考文獻(xiàn)1914第1章 設(shè)計(jì)任務(wù)與要求11對(duì)于三種典型序列-單位采

2、樣序列、實(shí)指數(shù)序列、矩形序列要求： （1）畫(huà)出以上序列的時(shí)域波形圖； （2）求出以上序列的傅里葉變換； （3）畫(huà)出以上序列的幅度譜及相位譜，并對(duì)相關(guān)結(jié)果予以理論分析； （4）對(duì)以上序列分別進(jìn)行時(shí)移，畫(huà)出時(shí)移后序列的頻譜圖，驗(yàn)證傅里葉變換的時(shí)移性質(zhì)； （5）對(duì)以上序列的頻譜分別進(jìn)行頻移，求出頻移后頻譜所對(duì)應(yīng)的序列，并畫(huà)出序列的時(shí)域波形圖，驗(yàn)證傅里葉變換的頻移性質(zhì)。1.2自行設(shè)計(jì)一個(gè)周期序列要求： （1）畫(huà)出周期序列的時(shí)域波形圖； （2）求周期序列的DFS，并畫(huà)出幅度特性曲線； （3）求周期序列的FT，并畫(huà)出幅頻特性曲線； （4）比較DFS和FT的結(jié)果，從中可以得出什么結(jié)論。 第2 章 原理及過(guò)程

3、1設(shè)計(jì)原理運(yùn)用Matlab及相關(guān)軟件，對(duì)三種典型序列進(jìn)行頻譜分析，得到并驗(yàn)證相應(yīng)的結(jié)果。MATLAB是由美國(guó)mathworks公司發(fā)布的主要面對(duì)科學(xué)計(jì)算、可視化以及交互式程序設(shè)計(jì)的高科技計(jì)算環(huán)境。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以及非線性動(dòng)態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在一個(gè)易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計(jì)語(yǔ)言（如C、Fortran）的編輯模式，代表了當(dāng)今國(guó)際科學(xué)計(jì)算軟件的先進(jìn)水平。1 MATLAB和Mathematica、Maple并稱為三大數(shù)學(xué)軟件。它在數(shù)學(xué)類科技應(yīng)

4、用軟件中在數(shù)值計(jì)算方面首屈一指。MATLAB可以進(jìn)行矩陣運(yùn)算、繪制函數(shù)和數(shù)據(jù)、實(shí)現(xiàn)算法、創(chuàng)建用戶界面、連 接其他編程語(yǔ)言的程序等，主要應(yīng)用于工程計(jì)算、控制設(shè)計(jì)、信號(hào)處理與通訊、圖像處理、信號(hào)檢測(cè)、金融建模設(shè)計(jì)與分析等領(lǐng)域。 MATLAB的基本數(shù)據(jù)單位是矩陣，它的指令表達(dá)式與數(shù)學(xué)、工程中常用的形式十分相似，故用MATLAB來(lái)解算問(wèn)題要比用C，F(xiàn)ORTRAN等語(yǔ)言完成相同的事情簡(jiǎn)捷得多，并且MATLAB也吸收了像Maple等軟件的優(yōu)點(diǎn)，使MATLAB成為一個(gè)強(qiáng)大的數(shù)學(xué)軟件。在新的版本中也加入了對(duì)C，F(xiàn)ORTRAN，C+，JAVA的支持?？梢灾苯诱{(diào)用,用戶也可以將自己編寫(xiě)的實(shí)用程序?qū)氲組ATLAB

5、函數(shù)庫(kù)中方便自己以后調(diào)用，此外許多的MATLAB愛(ài)好者都編寫(xiě)了一些經(jīng)典的程序，用戶可以直接進(jìn)行下載就可以用。利用傅里葉變換的方法對(duì)振動(dòng)的信號(hào)進(jìn)行分解，并按頻率順序展開(kāi)，使其成為頻率的函數(shù)，進(jìn)而在頻率域中對(duì)信號(hào)進(jìn)行研究和處理的一種過(guò)程，稱為頻譜分析。將信號(hào)在時(shí)間域中的波形轉(zhuǎn)變?yōu)轭l率域的頻譜，進(jìn)而可以對(duì)信號(hào)的信息作定量解釋。對(duì)信號(hào)進(jìn)行頻譜分析，是對(duì)其進(jìn)行傅里葉變換，得到其振幅譜與相位譜。分析軟件主要為Matlab。 對(duì)于信號(hào)來(lái)說(shuō)，分模擬信號(hào)與數(shù)字信號(hào)。進(jìn)行頻譜分析時(shí)，對(duì)于模擬信號(hào)來(lái)說(shuō)，首先對(duì)其進(jìn)行抽樣，使其離散化，然后利用離散傅里葉變換（DFT）或者快速傅里葉變換（FFT），然后對(duì)其幅度（ABS）

6、和相位（ANGLE）的圖像進(jìn)行分析，而對(duì)于數(shù)字信號(hào)來(lái)說(shuō)，則可直接進(jìn)行離散傅里葉變換或快速傅里葉變換第3 章 設(shè)計(jì)內(nèi)容1.1單位采樣序列1.1.1時(shí)域波形 程序：n=-5:10;y=zeros(1,5),1,zeros(1,10);stem(n,y)axis(-5,10,0,2);title('單位取樣序列') 圖11.2傅里葉變換 程序：clearclc% 輸入信號(hào)n=-5:10;%時(shí)域信號(hào)的時(shí)間范圍x=zeros(1,5),1,zeros(1,10);%時(shí)域信號(hào)xw=0:1e-2:2;%想要觀察的頻率范圍% 預(yù)定義y=w;a=w;j=sqrt(-1);%先定義變量維度，提高運(yùn)

7、算速度% 計(jì)算頻點(diǎn)for i=1:length(w)f=trapz(n,x.*exp(-j*w(i)*n);y(i)=abs(f);a(i)=angle(f);end% 輸出subplot(3,1,1),plot(n,x)subplot(3,1,2),plot(w,y)subplot(3,1,3),plot(w,a)圖2 1.3幅度譜及相位譜n=1:50; % 定義序列的長(zhǎng)度是50 x=zeros(1,50); % 注意：MATLAB 中數(shù)組下標(biāo)從1 開(kāi)始 x(1)=1;close all; subplot(3,1,1);stem(x);title('單位沖擊信號(hào)序列'); k

8、=-25:25; X=x*(exp(-j*pi/12.5).(n'*k); magX=abs(X); % 繪制x(n) 的幅度譜 subplot(3,1,2);stem(magX);title('單位沖擊信號(hào)的幅度譜'); angX=angle(X); % 繪制x(n) 的相位譜 subplot(3,1,3);stem(angX) ; title ('單位沖擊信號(hào)的相位譜') ;圖31.4頻移程序：fs=256;N=256; %采樣頻率和數(shù)據(jù)點(diǎn)數(shù)n=0:N-1;t=n/fs; %時(shí)間序列f=20; %定義信號(hào)中心頻率(hz)x=zeros(1,5),1,

9、zeros(1,10); %信號(hào)%程序區(qū)y1=fft(x,N); %對(duì)信號(hào)進(jìn)行快速Fourier變換y2=fftshift(y1);mag2=abs(y2); f1=n*fs/N; %頻率序列f2=n*fs/N-fs/2; %shift后頻率序列subplot(4,1,1),plot(f2,mag2,'r'); %繪出隨頻率變化的振幅圖4 1.5時(shí)移程序N=50;n=0:N-1;M=20;xn=1 zeros(1,N); nm=mod(n-M),N);xm=xn(nm+1);subplot(3,1,1),stem(xm);title('|X(k)|');k=-2

10、5:25;X=x*(exp(-j*pi/25).(nm)'*k); magX=abs(X); % 繪制x(m) 的幅度譜 subplot(3,1,2);stem(magX);title('矩形序列的幅度譜'); angX=angle(X); % 繪制x(m) 的相位譜 subplot(3,1,3);stem(angX) ; title ('矩形序列相位譜') ;圖52.1時(shí)域圖形程序：n=-5:10;y=zeros(1,5),ones(1,5),zeros(1,6);plot(n,y)stem(n,y)axis(-5,10,0,2);title('

11、;矩形序列');圖62.2傅里葉變換程序：clearclc% 輸入信號(hào)n=-5:10;%時(shí)域信號(hào)的時(shí)間范圍x=zeros(1,5),1,zeros(1,10);%時(shí)域信號(hào)xw=0:1e-2:2;%想要觀察的頻率范圍% 預(yù)定義y=w;a=w;j=sqrt(-1);%先定義變量維度，提高運(yùn)算速度% 計(jì)算頻點(diǎn)for i=1:length(w)f=trapz(n,x.*exp(-j*w(i)*n);y(i)=abs(f);a(i)=angle(f);end% 輸出subplot(3,1,1),plot(n,x)subplot(3,1,2),plot(w,y)subplot(3,1,3),plo

12、t(w,a) 圖72.3幅度譜與相位譜程序：n=1:50;x=sign(sign(10-n)+1); close all;subplot(3,1,1); stem(x);title('矩形信號(hào)序列'); k=-25:25; X=x*(exp(-j*pi/25).(n'*k); magX=abs(X); % 繪制x(n) 的幅度譜 subplot(3,1,2);stem(magX);title('矩形序列的幅度譜'); angX=angle(X); % 繪制x(n) 的相位譜 subplot(3,1,3);stem(angX) ; title ('

13、矩形序列相位譜') ;圖82.4頻移程序：fs=256;N=256; %采樣頻率和數(shù)據(jù)點(diǎn)數(shù)n=0:N-1;t=n/fs; %時(shí)間序列f=20; %定義信號(hào)中心頻率(hz)x=zeros(1,5),ones(1,5),zeros(1,6); %信號(hào)%程序區(qū)y1=fft(x,N); %對(duì)信號(hào)進(jìn)行快速Fourier變換y2=fftshift(y1);mag2=abs(y2); f1=n*fs/N; %頻率序列f2=n*fs/N-fs/2; %shift后頻率序列subplot(4,1,1),plot(f2,mag2,'r'); %繪出隨頻率變化的振幅圖92.5時(shí)移程序：N=5

14、0;n=0:N-1;M=20;xn=sign(sign(10-n)+1); nm=mod(n-M),N);xm=xn(nm+1);subplot(3,1,1),stem(xm);title('|X(k)|');k=-25:25;X=x*(exp(-j*pi/25).(nm)'*k); magX=abs(X); % 繪制x(m) 的幅度譜 subplot(3,1,2);stem(magX);title('矩形序列的幅度譜'); angX=angle(X); % 繪制x(m) 的相位譜 subplot(3,1,3);stem(angX) ; title (&

15、#39;矩形序列相位譜') ;圖104.1幅度特性曲線程序：n=1:50; % 定義序列的長(zhǎng)度是50 x=sin(pi*n); % 注意：MATLAB 中數(shù)組下標(biāo)從1 開(kāi)始 close all; subplot(3,1,1);stem(x);title('單位沖擊信號(hào)序列'); k=-25:25; X=x*(exp(-j*pi/12.5).(n'*k); magX=abs(X); % 繪制x(n) 的幅度譜 subplot(3,1,2);stem(magX);title('單位沖擊信號(hào)的幅度譜'); angX=angle(X); % 繪制x(n)

16、 的相位譜 subplot(3,1,3);stem(angX) ; title ('單位沖擊信號(hào)的相位譜') ;圖154.4周期序列的DFS程序：n=0:9;k=0:1:9;x=sin(5*pi*n);WN=exp(-j*2*pi/9);nk=n'*k;WNnk=WN.nk;Xk=x*WNnk; xlabel('Xk')4.5傅里葉變換程序：fs=100;%設(shè)定采樣頻率 N=128; n=0:N-1; t=n/fs; %設(shè)定正弦信號(hào)頻率 %生成正弦信號(hào) x=sin(pi *t); figure(1); subplot(231); plot(t,x);%作

17、正弦信號(hào)的時(shí)域波形 xlabel('t'); ylabel('y'); title('正弦信號(hào)y=2*pi*10t時(shí)域波形'); grid; %進(jìn)行FFT變換并做頻譜圖 y=fft(x,N);%進(jìn)行fft變換 mag=abs(y);%求幅值 f=(0:length(y)-1)'*fs/length(y);%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換 figure(1); subplot(232); plot(f,mag);%做頻譜圖 axis(0,100,0,80); xlabel('頻率(Hz)'); ylabel('幅值'); title('正弦信號(hào)y=2*pi*10t幅頻譜圖N=128');圖16 第4章 心得與體會(huì)這次的課設(shè)設(shè)計(jì)對(duì)我來(lái)說(shuō)時(shí)間是比較緊促的，但是門(mén)課程是一門(mén)動(dòng)手能力要求很高的課程，需要很多實(shí)踐操作才能真正的去掌握它，這次的課程設(shè)計(jì)為我提供了很好的機(jī)會(huì)。通過(guò)這次的課程