• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

MATLAB最小二乘法

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

 MATLAB最小二乘法

作者:凯鲁嘎吉 - 博客园
http://www.cnblogs.com/kailugaji/

三、实验程序

四、实验内容

设有如下数据:

3次多项式拟合这组数据.

五、解答(按如下顺序提交电子版)

1.(程序)

LSM1.m:

function p=LSM1(x,y,m)  %x,y为序列长度相等的数据向量,m为拟合多项式次数
format short;
A=zeros(m+1,m+1);
for i=0:m
    for j=0:m
        A(i+1,j+1)=sum(x.^(i+j));
    end
    b(i+1)=sum(x.^i.*y);
end
a=A\b\';
p=fliplr(a\');

2.(运算结果)

>> x=[-3,-2,-1,0,1,2,3];
>> y=[-1.76,0.42,1.2,1.34,1.43,2.25,4.38];
>> p=LSM1(x,y,3)

p =

    0.1133   -0.0018    0.0035    1.3300

3.(拓展(方法改进、体会等))

MATLAB中有关于最小二乘法的现成的函数,如下编写程序:

ploy1.m:

function [p]=ploy1(x,y)  
P=polyfit(x,y,3);  
xi=-4:.2:4;  
p=polyfit(x,y,3);
yi=polyval(P,xi);  
plot(xi,yi,x,y,\'r*\'); 

结果:

>> x=[-3,-2,-1,0,1,2,3];
y=[-1.76,0.42,1.2,1.34,1.43,2.25,4.38]; [p]=ploy1(x,y)

p =

    0.1133   -0.0018    0.0035    1.3300

则y=0.1133*x^3-0.0018*x^2+0.0035*x+1.33

 注:非线性曲线拟合见:MATLAB实例:非线性曲线拟合


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
微软正在开发基于 Rust 的安全编程语言发布时间:2022-07-18
下一篇:
matlab *与.*的区别发布时间:2022-07-18
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap