如何用MATLAB绘制功率谱密度图形

如何用MATLAB绘制功率谱密度图形
如何用MATLAB绘制功率谱密度图形

如何用MATLAB绘制功率谱密度图形
图片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,然后求模的平方,得到功率谱.结果的图像看起来是杂乱的,围绕