prony1.txt

PRONG的应用程序,一般都可以用的吧.啊哦

x=0:0.05:10;%原理想信号
y=220*exp(-0.25*x).cos(2*pi*1.5*x+pi)+110*exp(0.25*x).cos(2*pi*10*x+pi/2);
p=6% p为模型阶数
sf=1/0.001 %sf为采样频率
[F,D,A,theta]=exprony_m a(y,p,sf)  %调用prony算法函数
function [F,D,A,theta]=exprony_m a(y,p,sf);
nm=2*p;
n=fix(length(y)/2);
dt=1/sf; %由采样频率计算时间间隔
L=length(h);
M=L/2;
for k=1:nm
x1(:,k)=h(k:M-1-K);
end
for k=1:M
x2(k,:)=-h(nm+k);
end
B=x1\x2 %最小二乘求解prony多项式系数
B(nm+1)=1；
B1=B(nm+1:-1:1);
V=roots(B1);%求出特征方程的根
F1=abs(log(v)))/(2*pi*dt); %初步计算频率
for k=0:(2*n-1)
Va(k+1,:)=conj(V %生成范德梦矩阵
end
S1=2*(inv(conj(Va)*Va)conj(Va)*h);
[F2,I]=sort(F1);将平率F1从小到大排序
m=0;
for k=1:nm-1
if F2(K)=F2(K+1) %消除其中的共轭相
continue；
end
n=m+1；
ii=I(k);%提取有效的元素序数
F(m)=F1(ii);
D(m)=D1(ii);
A(m)=abs(S1(ii));
theta(m)=angle(S1(ii));
end