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

淩亂°似流年 2023-01-15 08:06 145阅读 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 条评论，145人围观）

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

相关阅读

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

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

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

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

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

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

相关无线发射功率

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

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

相关 matlab全搜索运动估计,三步搜索运动估计

种运动估计的快速预测搜索算法 (S) 该算法首先用当前块的三个 /PA.... x? 2 (x ? di ) ?di x?B ( xi ) 使用基于一阶梯度下降的方法

怼烎@/ 2023年01月17日 04:48/ 0 赞/ 121 阅读

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

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

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

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

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

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

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

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

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

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

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

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

相关有功功率计算公式和无功功率计算公式

无功功率(reactive power)在具有电感和的电路里,这些储能元件在半周期的时间里把能量变成磁场(或电场)的能量存起来,在另半周期的时间里对已存的磁场(或电场)能量送还

素颜马尾好姑娘i/ 2022年11月10日 05:38/ 0 赞/ 226 阅读

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

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

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