LS算法最小二乘法 算法
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
clear all;
close all;
N=128;
n=1:128;
a1=sqrt(20);
a2=sqrt(2);
x=a1*sin(2*pi*0.2*n)+a2*sin(2*pi*0.213*n); figure(2);
R=xcorr(x,'unbiased');
Rx=0;
Rx=hankel(R); %生成自相关函数矩阵
p=4;
A=Rx(1:p,1:p); %构造p*p汉克矩阵
b=-Rx(2:p+1,p);
Ls=inv(A'*A)*(A'*b);%最小二乘法解出来参数Ls=[1;Ls];
freqz(1,Ls,128,1);
title('Ls');