图片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，然后求模的平方，得到功率谱。结果的图像看起来是杂乱的，围绕
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024