prony.m

来自「基于prony算法的低频振荡模式识别」· M 代码 · 共 11 行

%基于prony算法的低频振荡模式识别
clc
clear
x=0:199;
%低频振荡信号
y=220*exp(-0.25*x*0.05).*cos(2*pi*1.5*x*0.05+pi)+110*exp(0.25*x*0.05).*cos(2*pi*1.0*x*0.05+pi/2);
plot(y)
p=6 %p为模型阶数
sf=1/0.05%sf为采样频率
[F,D,A,theta]=exprony_ma(y,p,sf)