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

╰半橙微兮° 2023-01-15 09:29 145阅读 0赞

下面的数据是从1770年至1869年这100年间每年12个月所记录到的太阳黑子出现次数的平均值。

a. 用周期图法做该数据的功率谱，从该曲线您是否可以看出太阳黑子活动的周期？

b. 对该数据建立一个三阶的AR模型，分别用Levinson-Durbin法和Burg法求解AR模型的参数，并得到功率谱曲线，从该曲线来发现太阳黑子的活动周期。改变AR模型的阶次，看不同阶次下所估计出的周期的差别。

x = \[101 82 66 35 31 7 20 92 154 125 85 68 38 23 10 24 83 132 131 118 90 67 60 47 41 21 16 6 4 7 14 34 45 43 48 42 28 10 8 2 0 1 5 12 14 35 46 41 30 24 16 7 4 2 8 17 36 50 62 67 71 48 28 8 13 57 122 138 103 86 63 37 24 11 15 40 62 98 124 96 66 64 54 39 21 7 4 23 55 94 96 77 59 44 47 30 16 7 37 74\];

下面是matlab程序：

clear;

clc;

x = \[101 82 66 35 31 7 20 92 154 125 85 68 38 23 10 24 83 132 131 118 90 67 60 47 41 21 16 6 4 7 14 34 45 43 48 42 28 10 8 2 0 1 5 12 14 35 46 41 30 24 16 7 4 2 8 17 36 50 62 67 71 48 28 8 13 57 122 138 103 86 63 37 24 11 15 40 62 98 124 96 66 64 54 39 21 7 4 23 55 94 96 77 59 44 47 30 16 7 37 74\];

N=length(x);

R=ZXG(x);

%han=zeros(1,1024);

%定义一个窗长为N的汉明窗

for i=1:N

han(i)=100\*(0.54-0.46\*cos(2\*pi\*i/N));

end

%用周期图法求数据的功率谱估计P

X=abs(fft(x,1024));

for i=1:1024

P(i)=X(i)\*X(i)/N;

end

plot(P);

%改进之后的周期图法,u为汉明窗

for i=1:99

u(i)=R(i)\*han(i);

end

UU=abs(fft(u,1024));

figure(2);

plot(UU);

axis(\[0 100 0 1000000\]);

grid on

%三阶AR模型参数得其功率谱

a\_11=-R(2)/R(1);

D\_1=(1-abs(a\_11)\*abs(a\_11))\*R(1);

a\_22=-(R(3)+a\_11\*R(2))/D\_1;

a\_21=a\_11+a\_22\*a\_11;

D\_2=(1-abs(a\_22)\*abs(a\_22))\*D\_1;

a\_33=-(R(4)+a\_21\*R(3)+a\_22\*R(2))/D\_2;

a\_32=a\_22+a\_33\*a\_21;

a\_31=a\_21+a\_33\*a\_22;

D\_3=(1-abs(a\_33)\*abs(a\_33))\*D\_2;

A=\[a\_31 a\_32 a\_33\];

D=abs(fft(A,256));

for i=1:256

PP(i)=D\_3/((1+D(i))\*(1+D(i)));

end

figure(3);

plot(PP);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%因为AR谱估计高阶算法与三阶算法类似， %

%%%而这种方法与线性预测谱估计等效，故 %

%%%可以直接直接调用求线性预测的系统函数LPC(X,N)

L=LPC(x,70);

LL=abs(fft(L,256));

for i=1:256

PPP(i)=D\_3/((1+LL(i))\*(1+LL(i)));

end

figure(4);

plot(PPP);

%axis(\[0 128 0 5\]);

此程序主要用于输入数据的功率谱估计，该程序由三部分组成：burg\_Kp.m，Mainfun.m，ZXG.m.其中：Mainfun.m为主函数，burg\_Kp.m为Burg算法函数，ZXG.m为求输入信号自相关的函数(估计)。

问题在于，如何在功率谱密度图中添加一条信度值为0.9的置信曲线？

\[本帖最后由 vincent92 于 2009-3-10 13:42 编辑\]

发表评论取消回复

表情：

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

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

相关阅读

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

相关 matlab求雷克子波相位谱,求雷克子波的振幅谱和相位谱的MATLAB程序，谢谢！！！！！...

%% 自己写了如下程序，望采纳！ clc;clear;close all; %Ricker子波 f=40;dt=0.004; nw=6./f/dt; nw=2\flo

布满荆棘的人生/ 2023年01月15日 08:23/ 0 赞/ 182 阅读

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

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

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

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

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

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

相关 matlab求方差、均方值和功率

% 求方差、均方值、功率等等，今天要弄明白了。 clc;clear all;close all; %实数域 x = [1;3

绝地灬酷狼/ 2022年09月10日 06:08/ 0 赞/ 205 阅读