数学模型实验课（二）.pptVIP

下载本文档

1
0
约1.18千字
约 10页
2017-08-30 发布于上海
举报
版权申诉

数学模型实验课（二）.ppt

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

数学模型实验课（二）最小二乘法与直线拟合康尔乃奶粉 32.4元 400g; 67.1元 900g. 4 k1 + 42/3 k2 = 32.4 9 k1 + 92/3 k2 = 67.1 解得: k1 = 5.3791, k2 = 4.3192 模型: C(W)=5.3791 W + 4.3192 W2/3. 预测: W=1800, C(W) = 126.49. W=2500, C(W) = 154.36 实际: W=1800, C(W)=115.9 W=2500, C(W)=146.85 x=[4,9,18,25];y=[32.4,67.1,115.9,146.85]; x1=[x(1),x(2)]’,x2=x1.^(2/3); y1=[y(1),y(2)]’;X=[x1,x2]; k=X\y1 x1=[x(1),x(2),x(3)];x2=x1.^(2/3); y1=[y(1),y(2),y(3)];X=[x1,x2]; k=X\y1 y=[17,16,17,23,26,27,41,49]; x=[1953,2015,2015,2821,3049,3049,5133,5592]; X=[ones(8,1),x’];b=X\y’; y1=b(1)+b(2).*x; plot(x,y1,’r’,x,y,’*r’) 2. 直线拟合的效果（相关系数） X1=[x’,y’];m=size(X1);n=m(1);p=m(2); X=[ones(n,1),x’];B=X\y’;a=B(1);b=B(2:p); y1=a+b*x;plot(x,y1,’r’,x,y,’*r’); Cov=X1’*(eye(n)-ones(n)./n)*X1; Lyy=Cov(p,p);Lxx=Cov(1:(p-1),1:(p-1)); R=(b’*Lxx*b)/Lyy;r=sqrt(R); 3. 多项式拟合 b=polyfit(x,y,1); 4.线性回归 [b,r,j,k,l]=regress(y’,X); * * . 求解线性方程组：A x = y 如果设：Ax+ε= y，求 x 使得误差 ||ε|| 最小。 A’A x = A’y x = (A’A)-1A’y 二. 直线拟合 1. 矩阵左除数据（xi, yi）,拟合直线 y = a + bx, 误差最小模型：yi = a + bxi +ε X b = y