圖片1 圖片3 答案1：： function [t,omg,FT,IFT] = prefourier(Trg,N,OMGrg,K) % 輸入參數： % Trg : 二維矢量，兩個元素分別表示時域信號的起止時間； % N : 時域抽樣數量； % OMGrg: 二維矢量，兩個元素分別表示頻譜的起止頻率； % K : 頻域抽樣數量。 % 輸出參數： % t : 抽樣時間； % omg : 抽樣頻率； % FT : 實現傅裏葉變換的矩陣~U~及系數； % IFT : 實現傅裏葉逆變換的矩陣~V~及系數。 T = Trg(2)-Trg(1); t = linspace(Trg(1),Trg(2)-T/N,N)'; OMG = OMGrg(2)-OMGrg(1); omg = linspace(OMGrg(1),OMGrg(2)-OMG/K,K)'; FT = T/N*exp(-j*kron(omg,t.')); IFT = OMG/2/pi/K*exp(j*kron(t,omg.')); end 在另壹個腳本文件中： clc;clear ;close all; N=1024*8; K=500; OMGrg=[0,100]; Trg=[0,1]; [t,omg,FT,IFT] = prefourier(Trg,N,OMGrg,K); % f0=10; % f=sin(2*pi*f0*t); f=randn(N,1); F=FT*f; figure; plot(t,f); figure; plot(omg/2/pi,abs(F).^2); 高斯白噪聲的功率譜理論上為壹直線，妳給出的功率譜圖不對，除非它是在某些特定情況下成立，比如經過了濾波器。功率譜密度為傅裏葉變換的平方，在我的程序中，先求出傅裏葉變換 F，然後求模的平方，得到功率譜。結果的圖像看起來是雜亂的，圍繞 :::::::::::::::::::請參考以下相關問題:::::::::::::::::::: MATLAB畫功率譜密度的圖形...求高手 :::::::::::::::::::請參考以下相關問題:::::::::::::::::::: Matlab畫圖問題，自相關函數，功率譜密度，求具體代碼。 :::::::::::::::::::請參考以下相關問題:::::::::::::::::::: 為什麽matlab畫出的 高斯白噪聲功率譜密度不是直線 :::::::::::::::::::請參考以下相關問題:::::::::::::::::::: 為什麽用matlab畫出的高斯白噪聲功率譜密度不是直線 :::::::::::::::::::請參考以下相關問題:::::::::::::::::::: 誰會matlab編程顯示功率譜密度呢，跪求... 答案1：： function [t,omg,FT,IFT] = prefourier(Trg,N,OMGrg,K) % 輸入參數： % Trg : 二維矢量，兩個元素分別表示時域信號的起止時間； % N : 時域抽樣數量； % OMGrg: 二維矢量，兩個元素分別表示頻譜的起止頻率； % K : 頻域抽樣數量。 % 輸出參數： % t : 抽樣時間； % omg : 抽樣頻率； % FT : 實現傅裏葉變換的矩陣~U~及系數； % IFT : 實現傅裏葉逆變換的矩陣~V~及系數。 T = Trg(2)-Trg(1); t = linspace(Trg(1),Trg(2)-T/N,N)'; OMG = OMGrg(2)-OMGrg(1); omg = linspace(OMGrg(1),OMGrg(2)-OMG/K,K)'; FT = T/N*exp(-j*kron(omg,t.')); IFT = OMG/2/pi/K*exp(j*kron(t,omg.')); end 在另壹個腳本文件中： clc;clear ;close all; N=1024*8; K=500; OMGrg=[0,100]; Trg=[0,1]; [t,omg,FT,IFT] = prefourier(Trg,N,OMGrg,K); % f0=10; % f=sin(2*pi*f0*t); f=randn(N,1); F=FT*f; figure; plot(t,f); figure; plot(omg/2/pi,abs(F).^2); 高斯白噪聲的功率譜理論上為壹直線，妳給出的功率譜圖不對，除非它是在某些特定情況下成立，比如經過了濾波器。功率譜密度為傅裏葉變換的平方，在我的程序中，先求出傅裏葉變換 F，然後求模的平方，得到功率譜。結果的圖像看起來是雜亂的，圍繞
上一篇:關於ul和dl的區別下一篇:Windows Phone 7程序設計的譯者序