马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
%计算Y序列的自相关系数,n为序列长度
sum1=1/n*Y*Y';
for k=1:n
sum3(k)=0;
for m=1:n-k
sum3(k)=sum3(k)+1/(n-k)*Y(m)*Y(m+k);
end
r0(k)=sum3(k)/sum1; %最后求得的自相关系数存在r0向量中。
end
%计算Y序列的偏相关系数
f0(1,1)=r0(1);
pr0(1)=f0(1,1);
for k=2:n-1
sum4=0;
sum5=0;
for j=1:k-1
sum4=sum4+f0(k-1,j)*r0(k-j);
sum5=sum5+f0(k-1,j)*r0(j);
end
f0(k,k)=(r0(k)-sum4)/(1-sum5);
for j=1:k-1
f0(k,j)=f0(k-1,j)-f0(k,k)*f0(k-1,k-j);
end
end
for i=2:n-1
pr0(i)=f0(i,i); %最后的偏相关系数存在向量pro中。
end
大侠帮我看看,程序有没有问题,谢谢了 。
[ 本帖最后由 xuebx 于 2006-9-1 17:01 编辑 ] |