功率谱估计 matlab arma,MATLAB/SIMULINK 功率谱估计

淩亂°似流年 2023-01-15 08:06 153阅读 0赞

丁老师，您好！我测量了脉动速度，采样频率为4000Hz，测量时间为4.09575s，测量点数16384，想用FFT对其进行功率谱分析，我参考同学的程序(其采样频率10Hz，测量点数2000)进行了编程，但功率谱这项的程序感觉有点问题，与一些书上的程序不同，不太确定，请丁老师帮忙！还有我的频率是不是太大，测量时间是不是太小了？

程序如下：

%FFT变换,获得采样数据时域和频域图

%这里的向量都用行向量，被测变量是速度，单位为m/s

clear;

close all;

load data.txt                                 %仪器测量的原始数据，存储为data.txt中

A=data;                                          %将测量数据赋给A，A为N×2的数组

TimeSequence=A(:,1);

TimeSequence=TimeSequence';

FanData=A(:,2);

FanData=FanData';

%显示原始数据曲线图(时域)

subplot(2,1,1);

plot(TimeSequence,FanData)

axis(\[min(TimeSequence) max(TimeSequence) 1.1\*floor(min(FanData)) 1.1\*ceil(max(FanData))\])

xlabel('时间 (s)');

ylabel('脉动速度 (m/s)');

title('原始信号(时域)');

grid on;

SampleCapacity=16384

DeltaT=0.00025

AverageSpeed=mean(FanData)

PulsatileVector=FanData-AverageSpeed;

%频谱分析

WaveCounting=0

![19b0fe3299812d8066e982a21c3cbca5.gif][]SampleCapacity-1);

FrequencyVector=WaveCounting/(SampleCapacity\*DeltaT);    %频率向量

fftPulsatileVector=fft(PulsatileVector);

PowerSpectrumVector=2\*DeltaT\*((abs(fftPulsatileVector)).^2)/SampleCapacity;     % 功率谱向量

LogFrequencyVector=log10(FrequencyVector);

LogPowerSpectrumVector=log10(PowerSpectrumVector);

%显示频谱图(频域)

subplot(2,1,2);

plot(LogFrequencyVector(1:SampleCapacity/2),LogPowerSpectrumVector(1:SampleCapacity/2),'r')

axis(\[min(LogFrequencyVector(1:SampleCapacity/2)) max(LogFrequencyVector(1:SampleCapacity/2)) 1.1\*floor(min(LogPowerSpectrumVector(1:SampleCapacity/2))) 1.1\*ceil(max(LogPowerSpectrumVector(1:SampleCapacity/2)))\])

xlabel('lg f')

ylabel('lg E(f)')

title('频谱图(频域)')

grid on

![a70cbf5f56cb187f20fb09bae08ed3de.gif][]

2011-3-22 16:34 上传

点击文件名下载附件

258.22 KB, 下载次数: 7774

[19b0fe3299812d8066e982a21c3cbca5.gif]: /images/20221022/a6b67d9d67414e9bb793115d7e946103.png
[a70cbf5f56cb187f20fb09bae08ed3de.gif]: /images/20221022/fe09ba03a24548619bf154d8c7b0620a.png

发表评论取消回复

表情：

评论列表（有 0 条评论，153人围观）

还没有评论，来说两句吧...

相关阅读

相关 Python 图像频谱：探究图像的功率谱密度

Python 图像频谱：探究图像的功率谱密度分析图像频谱是计算机视觉和图像处理领域中的重要话题。在本文中，我们将深入探讨如何使用 Python 对图像执行傅立叶变换得到其频

刺骨的言语ヽ痛彻心扉/ 2024年03月26日 20:49/ 0 赞/ 52 阅读

相关基于Blackman方法的正弦信号功率谱估计Matlab仿真

基于Blackman方法的正弦信号功率谱估计Matlab仿真在信号处理中，功率谱估计是一项重要的任务。其中，Blackman方法是广泛使用的一种窗函数技术，它可以用于对信号

痛定思痛。/ 2024年03月24日 12:52/ 0 赞/ 91 阅读

相关 LS估计的AR参数估计：Cadzow谱估计与信号功率谱密度（Matlab实现）

LS估计的AR参数估计：Cadzow谱估计与信号功率谱密度（Matlab实现） AR模型是一种常用的时序数据建模方法，它可以用于信号分析、预测和滤波等领域。在AR模型中，参数

ゝ一世哀愁。/ 2024年03月04日 02:11/ 0 赞/ 82 阅读

相关 OFDM和F-OFDM的功率谱与峰均功率比仿真

OFDM和F-OFDM的功率谱与峰均功率比仿真摘要： OFDM(Orthogonal Frequency Division Multiplexing)是第四代通信技术

向右看齐/ 2024年02月19日 19:13/ 0 赞/ 55 阅读

相关无线发射功率

无线电发射功率无线电发射机输出的射频信号，通过馈线（电缆）输送到天线，由天线以电磁波形式辐射出去。电磁波到达接收地点后，由天线接收下来（仅仅接收很小很小一部分功率），

朴灿烈づ我的快乐病毒、/ 2023年10月12日 20:04/ 0 赞/ 81 阅读

相关 matlab 倒谱ceptrum,倒谱法求共振峰

我用自己的语音，读出来之后就出现 Index exceeds matrix dimensions. 求解为什么。我的Fs=22050，我想请问共振峰的个数与什么有关 >>

末蓝、/ 2023年01月16日 06:21/ 0 赞/ 140 阅读

相关 matlab功率谱补零,matlab如何在功率谱密度图中添加一条信度值为0.9的置信曲线

下面的数据是从1770年至1869年这100年间每年12个月所记录到的太阳黑子出现次数的平均值。 a. 用周期图法做该数据的功率谱，从该曲线您是否可以看出太阳黑子活动的周期？

╰半橙微兮°/ 2023年01月15日 09:29/ 0 赞/ 140 阅读

相关功率谱估计 matlab arma,MATLAB/SIMULINK 功率谱估计

丁老师，您好！我测量了脉动速度，采样频率为4000Hz，测量时间为4.09575s，测量点数16384，想用FFT对其进行功率谱分析，我参考同学的程序(其采样频率10Hz，测量

淩亂°似流年/ 2023年01月15日 08:06/ 0 赞/ 154 阅读

相关 matlab高斯白噪声功率谱密度,matlab-正弦波-高斯白噪声-均匀白噪声-功率谱密度-自相关函数.doc...

现代通信原理作业一利用matlab完成：产生正弦波信号、均匀白噪声以及高斯白噪声并分别将两种噪声叠加到正弦波信号上，绘出波形。分别求取均匀白噪声序列和高斯白噪声序

痛定思痛。/ 2023年01月13日 13:54/ 0 赞/ 173 阅读

相关算法--傅里叶变换，非周期信号频谱分析,周期信号的功率谱

算法--傅里叶变换，非周期信号频谱分析,周期信号的功率谱非周期信号进行分析的思路：从周期信号的傅里叶级数出发。时域：当周期T1--无穷大，周期信号成为非周期信号；频域

曾经终败给现在/ 2022年09月30日 07:02/ 0 赞/ 191 阅读