參與評論

評論

登錄后參與評論

全部評論(3條)

• 

  卡米拉tks 2008-05-11 00:00:00

  matlab 基本工具箱中的data analysis中的fourier analysis中就有fft函數(shù)包括快速傅立葉變換還有反變換什么的。

  贊(11)

  回復(fù)(0)

  評論

  評論

  登錄后參與評論

• 

  北暖29 2008-05-10 00:00:00

  Image Process 里好像有，不過help fft2試試

  贊(16)

  回復(fù)(0)

  評論

  評論

  登錄后參與評論

• 

  夢色未央 2018-05-06 14:47:25

  　　1、采樣數(shù)據(jù)導(dǎo)入Matlab 。 　　采樣數(shù)據(jù)的導(dǎo)入至少有三種方法。 　　diyi就是手動將數(shù)據(jù)整理成Matlab支持的格式，這種方法僅適用于數(shù)據(jù)量比較小的采樣。 　　第二種方法是使用Matlab的可視化交互操作，具體操作步驟為：File --> Import Data，然后在彈出的對話框中找到保存采樣數(shù)據(jù)的文件，根據(jù)提示一步一步即可將數(shù)據(jù)導(dǎo)入。這種方法適合于數(shù)據(jù)量較大，但又不是太大的數(shù)據(jù)。 　　第三種方法，使用文件讀入命令。數(shù)據(jù)文件讀入命令有textread、fscanf、load等，如采樣數(shù)據(jù)保存在txt文件中，則推薦使用 textread命令。如[a，b]=textread('data.txt'，'%f%*f%f'); 這條命令將data.txt中保存的數(shù)據(jù)三個三個分組，將每組的diyi個數(shù)據(jù)送給列向量a，第三個數(shù)送給列向量b，第二個數(shù)據(jù)丟棄。命令類似于C語言，詳細可查看其幫助文件。文件讀入命令錄入采樣數(shù)據(jù)可以處理任意大小的數(shù)據(jù)量，且錄入速度相當(dāng)快，一百多萬的數(shù)據(jù)不到20秒即可錄入。 　　2、對采樣數(shù)據(jù)進行頻譜分析 。 　　頻譜分析自然要使用快速傅里葉變換FFT了，對應(yīng)的命令即 fft ，簡單使用方法為：Y=fft(b，N)，其中b即是采樣數(shù)據(jù)，N為fft數(shù)據(jù)采樣個數(shù)。一般不指定N，即簡化為Y=fft(b)。Y即為FFT變換后得到的結(jié)果，與b的元素數(shù)相等，為復(fù)數(shù)。以頻率為橫坐標(biāo)，Y數(shù)組每個元素的幅值為縱坐標(biāo)，畫圖即得數(shù)據(jù)b的幅頻特性；以頻率為橫坐標(biāo)，Y數(shù)組每個元素的角度為縱坐標(biāo)，畫圖即得數(shù)據(jù)b的相頻特性。典型頻譜分析M程序舉例如下： clc fs=100; 　　t=[0:1/fs:100]; 　　N=length(t)-1;%減1使N為偶數(shù) %頻率分辨率F=1/t=fs/N 　　p=1.3*sin(0.48*2*pi*t)+2.1*sin(0.52*2*pi*t)+1.1*sin(0.53*2*pi*t)... +0.5*sin(1.8*2*pi*t)+0.9*sin(2.2*2*pi*t); 　　%上面模擬對信號進行采樣，得到采樣數(shù)據(jù)p，下面對p進行頻譜分析 　　figure(1) plot(t，p); grid on 　　title('信號 p(t)'); xlabel('t') ylabel('p') Y=fft(p); 　　magY=abs(Y(1:1:N/2))*2/N; f=(0:N/2-1)'*fs/N; figure(2) 　　%plot(f，magY); 　　h=stem(f，magY，'fill'，'--'); 　　set(h，'MarkerEdgeColor'，'red'，'Marker'，'*') grid on 　　title('頻譜圖 （理想值：[0.48Hz，1.3]、[0.52Hz，2.1]、[0.53Hz，1.1]、[1.8Hz，0.5]、[2.2Hz，0.9]） '); xlabel('f (Hz)') ylabel('幅值') 　　對于現(xiàn)實中的情況，采樣頻率fs一般都是由采樣儀器決定的，即fs為一個給定的常數(shù)；另一方面，為了獲得一定精度的頻譜，對頻率分辨率F有一個人為的規(guī)定，一般要求F<0.01，即采樣時間ts>100秒；由采樣時間ts和采樣頻率fs即可決定采樣數(shù)據(jù)量，即采樣總點數(shù)N=fs*ts。這就從理論上對采樣時間ts和采樣總點數(shù)N提出了要求，以保證頻譜分析的jing準(zhǔn)度。

  贊(1)

  回復(fù)(0)

  評論

  評論

  登錄后參與評論