最小二乘法matlab程序.docVIP

下载本文档

170
0
约1.81千字
约 5页
2019-08-07 发布于广东
举报
版权申诉

最小二乘法matlab程序.doc

1、本文档共5页，可阅读全部内容。
2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

PAGE PAGE 1 理解最小二乘法的基本原理，通过计算机解决相关问题；用matlab程序解决问题。实验目的： (1)理解最小二乘法的基本原理，通过计算机解决相关问题； (2)用matlab解决问题。实验内容：给定一组实验数据（Xi,Yi）(i=1,2,…n)编写一个通用程序，计算线性拟合和二次拟合（x,y）的值。 3. 程序源代码及运行结果：（一）线性拟合 x=[1 2 3 4 5 6 7 8] y=[15.3 20.5 27.4 36.6 49.1 65.6 87.8 117.6] x1=0; x12=0; y1=0; xy=0; n=8; for i=1:n x1=x1+x(i); x12=x12+x(i)^2 y1=y1+y(i); xy=xy+x(i)*y(i); end x1 y1 x12 xy a0=(y1*x12-x1*xy)/(n*x12-x1*x1); a1=(n*xy-x1*y1)/(n*x12-x1*x1); a0 a1 yy =a0+a1*x plot(x,y,o) plot(yy) 运行结果： x = 1 2 3 4 5 6 7 8 y = 15.3000 20.5000 27.4000 36.6000 49.1000 65.6000 87.8000 117.6000 x12 = 1 5 14 30 55 91 140 204 x1 = 36 y1 = 419.9000 x12 = 204 xy = 2.4794e+003 a0 = -10.7107 a1 = 14.0440 yy = Columns 1 through 7 3.3333 17.3774 31.4214 45.4655 59.5095 73.5536 87.5976 Column 8 101.6417 （二）二次拟合 xx=0; yy=0; x2=0; x3=0; x4=0; xy=0; y2=0; n=8; x1=1:8; y1=[15.3 20.5 27.4 36.6 49.1 65.6 87.8 117.6]; for i=1:n xx=xx+x1(i); x2=x2+x1(i)^2; x3=x3+x1(i)^3; x4=x4+x1(i)^4; yy=yy+y1(i); xy=xy+x1(i)*y1(i); y2=y2+x1(i)*x1(i)*y1(i); end xx x2 x3 x4 yy xy y2 A=[n xx x2;xx x2 x3;x2 x3 x4]; B=[yy;xy;y2]; C=A\B; C 运行结果： xx = 36 x2 = 204 x3 = 1296 x4 = 8772 yy = 419.9000 xy = 2.4794e+003 y2 = 1.6347e+004 C = 18.8518 -3.6935 1.9708 （三）线性拟合和二次拟合： function ZXE(x,y,m) S=zeros(1,2*m+1); T=zeros(m+1,1); for k=1:2*m+1 S(k)=sum(x.^(k-1)); end for k=1:m+1 T(k)=sum(x.^(k-1).*y); end A=zeros(m+1,m+1); a=zeros(m+1,1); for i=1:m+1 for j=1:m+1 A(i,j)=S(i+j-1); end end a=A\T; for k=1:m+1 fprintf(a[%d]=%f\n,k,a(k)); end 运行结果： ZXE(x,y,1); a[1]=-10.710714 a[2]=14.044048 ZXE(x,y,2); a[1]=18.851786 a[2]=-3.693452 a[3]=1.970833 图像：