1、Matlab編程實現FFT變換及頻譜分析的程序代碼(喜歡進行電磁兼容仿真分析的朋友可以借用2007-10-08 sysop 點擊: 436Matlab編程實現FFT變換及頻譜分析的程序代碼(喜歡進行電磁兼容仿真分析的朋友可以借用Matlab編程實現FFT變換及頻譜分析的程序代碼(喜歡進行電磁兼容仿真分析的朋友可以借用內容1.用Matlab產生正弦波,矩形波,以及白噪聲信號,并顯示各自時域波形圖2.進行FFT變換,顯示各自頻譜圖,其中采樣率,頻率、數據長度自選3.做出上述三種信號的均方根圖譜,功率圖譜,以及對數均方根圖譜4.用IFFT傅立葉反變換恢復信號,并顯示恢復的正弦信號時域波形圖源程序%*

2、% FFT實踐及頻譜分析 %*%*%*1.正弦波*%fs=100;%設定采樣頻率N=128;n=0:N-1;t=n/fs;f0=10;%設定正弦信號頻率%生成正弦信號x=sin(2*pi*f0*t;figure(1;subplot(231;plot(t,x;%作正弦信號的時域波形xlabel('t'ylabel('y'title('正弦信號y=2*pi*10t時域波形'grid;%進行FFT變換并做頻譜圖y=fft(x,N;%進行fft變換mag=abs(y;%求幅值f=(0:length(y-1'*fs/length(y;%進行對應的頻

3、率轉換figure(1;subplot(232;plot(f,mag;%做頻譜圖axis(0,100,0,80;xlabel('頻率(Hz'ylabel('幅值'title('正弦信號y=2*pi*10t幅頻譜圖N=128' grid;%求均方根譜sq=abs(y;figure(1;subplot(233;plot(f,sq;xlabel('頻率(Hz'ylabel('均方根譜'title('正弦信號y=2*pi*10t均方根譜' grid;%求功率譜power=sq.2;figure(1;subp

4、lot(234;plot(f,power;xlabel('頻率(Hz'ylabel('功率譜'title('正弦信號y=2*pi*10t功率譜'grid;%求對數譜ln=log(sq;figure(1;subplot(235;plot(f,ln;xlabel('頻率(Hz'ylabel('對數譜'title('正弦信號y=2*pi*10t對數譜'grid;%用IFFT恢復原始信號xifft=ifft(y;magx=real(xifft;ti=0:length(xifft-1/fs;figure(1;

5、subplot(236;plot(ti,magx;xlabel('t'ylabel('y'title('通過IFFT轉換的正弦信號波形' grid;%*2.矩形波*% fs=10;%設定采樣頻率t=-5:0.1:5;x=rectpuls(t,2;x=x(1:99;figure(2;subplot(231;plot(t(1:99,x;%作矩形波的時域波形xlabel('t'ylabel('y'title('矩形波時域波形'grid;%進行FFT變換并做頻譜圖y=fft(x;%進行fft變換mag=ab

6、s(y;%求幅值f=(0:length(y-1'*fs/length(y;%進行對應的頻率轉換figure(2;subplot(232;plot(f,mag;%做頻譜圖xlabel('頻率(Hz'ylabel('幅值'title('矩形波幅頻譜圖'grid;%求均方根譜sq=abs(y;figure(2;subplot(233;plot(f,sq;xlabel('頻率(Hz'ylabel('均方根譜'title('矩形波均方根譜'grid;%求功率譜power=sq.2;figure(2;s

7、ubplot(234;plot(f,power;xlabel('頻率(Hz'ylabel('功率譜'title('矩形波功率譜'grid;%求對數譜ln=log(sq;figure(2;subplot(235;plot(f,ln;xlabel('頻率(Hz'ylabel('對數譜'title('矩形波對數譜'grid;%用IFFT恢復原始信號xifft=ifft(y;magx=real(xifft;ti=0:length(xifft-1/fs;figure(2;subplot(236;plot(ti

8、,magx;xlabel('t'ylabel('y'title('通過IFFT轉換的矩形波波形'grid;%*3.白噪聲*%fs=10;%設定采樣頻率t=-5:0.1:5;x=zeros(1,100;x(50=100000;figure(3;subplot(231;plot(t(1:100,x;%作白噪聲的時域波形xlabel('t'ylabel('y'title('白噪聲時域波形'grid;%進行FFT變換并做頻譜圖y=fft(x;%進行fft變換mag=abs(y;%求幅值f=(0:length

9、(y-1'*fs/length(y;%進行對應的頻率轉換figure(3;subplot(232;plot(f,mag;%做頻譜圖xlabel('頻率(Hz'ylabel('幅值'title('白噪聲幅頻譜圖'grid;%求均方根譜sq=abs(y;figure(3;subplot(233;plot(f,sq;xlabel('頻率(Hz'ylabel('均方根譜'title('白噪聲均方根譜'grid;%求功率譜power=sq.2;figure(3;subplot(234;plot(f,p

10、ower;xlabel('頻率(Hz'ylabel('功率譜'title('白噪聲功率譜'grid;%求對數譜ln=log(sq;figure(3;subplot(235;plot(f,ln;xlabel('頻率(Hz'ylabel('對數譜'title('白噪聲對數譜'grid;%用IFFT恢復原始信號xifft=ifft(y;magx=real(xifft;ti=0:length(xifft-1/fs;figure(3;subplot(236;plot(ti,magx;xlabel('t'ylabel('y'title('通過IFFT轉換的白噪聲波形'grid;本文責任編輯: sy