MATLAB的離散傅里葉變換的仿真

1、應用MATLAB對信號進行頻譜分析及濾波設計目的要求學生會用MATLA語言進行編程，繪出所求波形，并且運用FFT求對連續(xù)信號進行分析。一、設計要求1、用Matlab產(chǎn)生正弦波，矩形波，并顯示各自的時域波形圖；2、進行FFT變換，顯示各自頻譜圖，其中采樣率、頻率、數(shù)據(jù)長度自選， 要求注明；3、繪制三種信號的均方根圖譜；4、用IFFT回復信號，并顯示恢復的正弦信號時域波形圖。二、系統(tǒng)原理用FFT對信號作頻譜分析是學習數(shù)字信號處理的重要內(nèi)容。經(jīng)常需要進行頻譜分析的信號是模擬信號和時域離散信號。頻譜分辨率直接和FFT的變換區(qū)間N有關，因為FFT能夠?qū)崿F(xiàn)頻率分辨率是2n /N。x(n)是一個長度為M的有

2、限長序列，貝U x(n)的N點離散傅立葉變換為：NX(k)=DFTx (n )=n1x( n)W0,k=0,1,.,N-1.2Wne J N1 N1knX逆變換：x(n) =IDFTX(k)=N n0,k=0,1,. ,N-1但FFT是一種比DFT更加快速的一種算法，提高了 DFT的運算速率,為數(shù)字 信號處理技術(shù)應用于各種信號處理創(chuàng)造了條件，大大提高了數(shù)字信號處理技術(shù)的開展。本實驗就是采用FFT, IFFT對信號進行譜分析。三、程序設計fs=i nput('please in put the fs:');%設定采樣頻率N=i nput('please in put th

3、e N:');%設定數(shù)據(jù)長度t=0:0.001:1;f=100;%設定正弦信號頻率%生成正弦信號x=sin(2*pi*f*t);figure(1);subplot(211);plot(t,x);% 作正弦信號的時域波形axis(0,0.1,-1,1);title(' 正弦信號時域波形 ');z=square(50*t);subplot(212)plot(t,z) axis(0,1,-2,2);title(' 方波信號時域波形 ');grid;泌行FFT變換并做頻譜圖y=fft(x,N);% 進行 fft 變換mag=abs(y);% 求幅值f=(0:N-

4、1)*fs/N;% 橫坐標頻率的表達式為 f=(0:M-1)*Fs/M; figure(2);subplot(211);plot(f,mag);% 做頻譜圖axis(0,1000,0,200);title(' 正弦信號幅頻譜圖 ');y仁fft(z,N);% 進行fft變換mag=abs(y1);% 求幅值f=(0:N-1)*fs/N;% 橫坐標頻率的表達式為 f=(0:M-1)*Fs/M;subplot(212);plot(f,mag);% 做頻譜圖axis(0,1000,0,200);title(' 方波信號幅頻譜圖 ');grid;%求功率譜sq=abs(

5、y);power=sq42;figure(3)subplot(211);plot(f,power);title(' 正弦信號功率譜 ');grid; sq1=abs(y1);power1=sq1.A2;subplot(212);plot(f,power1);title(' 方波信號功率譜 ');grid;%用 IFFT恢復原始信號 xifft=ifft(y);magx=real(xifft); ti=0:length(xifft)-1/fs;figure(4);subplot(211);plot(ti,magx);axis(0,0.1,-1,1);title(&

6、#39;通過IFFT轉(zhuǎn)換的正弦信號波形');zifft=ifft(y1);magz=real(zifft); ti1=0:length(zifft)-1/fs;subplot(212); plot(ti1,magz);title(' 通過IFFT轉(zhuǎn)換的方波信號波形);grid;四、仿真結(jié)果及分析由圖可以看出正弦波周期T=0.01，采樣點N=1024程序為:x=s in (2*pi*f*t);figure(1);subplot(211);plot(t,x);%作正弦信號的時域波形axis(0,0.1,-1,1);title('正弦信號時域波形');z=square

7、(50*t);subplot(212)plot(t,z)axis(0,1,-2,2);title'方波信號時域波形'grid;2、對正弦波、方波信號進行FFT變換程序：y=fft(x,N);% 進行fft 變換mag=abs(y);% 求幅值f=(0:N-1)*fs/N;%橫坐標頻率的表達式為 f=(0:M-1)*Fs/M;figure (2);subplot(211);plot(f,mag);% 做頻譜圖axis(0,1000,0,200);title('正弦信號幅頻譜圖');y1=fft(z,N);% 進行fft 變換mag=abs(y1);% 求幅值f=(

8、0:N-1)*fs/N;%橫坐標頻率的表達式為 f=(0:M-1)*Fs/M;subplot(212);plot(f,mag);% 做頻譜圖axis(0,1000,0,200);title('方波信號幅頻譜圖');grid;正弦信號、方波信號功率譜程序:sq=abs(y);power=sq.A2;figure©)subplot(211);plot(f,power);title('正弦信號功率譜');grid;sq 1= abs(y1);power仁 sq1.A2;subplot(212);plot(f,power1);title('方波信號功率譜');grid;對兩個信號進行恢復程序:xifft=ifft(y);magx=real(xifft);ti=O:le ngth(xifft)-1/fs;figure(4);subplot(211);plot(ti,magx);axis(0,0.1,-1,1);title('通過IFFT轉(zhuǎn)換的正弦信號波形');zifft=ifft(y1);ma