网易首页 > 网易号 > 正文 申请入驻

最小二乘法:多项式拟合实例(matlab实现)

0
分享至

一、什么是最小二乘法

最小二乘法在生活中的例子:

拿测一个人的身高为例,现实生活中我们无法用测量仪器完全准确地测得一个人的身高,相应地,用任何仪器测量都会存在误差。当然,我们求某人的身高时也不需要完全准确,我们只要使误差在一定范围之内就行了。比如我们生活中并不认为1.831米和1.832米身高的人有什么区别。但是我们会觉得1.83米的人和1.85米的人有所差别,比如很多女生会更喜欢身高为后者的男生,因而我们要把测量身高的误差控制在0.01米以内。这里我们用v代表误差。

那么这就引出了一个问题,如何来降低误差。

很直观的想法就是去求样本的平均数,比如测五次同一个人的身高,求五个样本的平均数,然后就认为这个平均数就是这个人的身高。那么问题来了,为什么要用算数平均数,而不是中位数、几何平均数还有调和平均数之类的。

要解决这个问题,首先我们得假设一个值就是这个人的准确身高。比如我们假设某个值y就是身高的准确值,那么误差v就是 |y-yi| ,(yi代表样本,就是每次的身高测量值)(如下图)

黑色直线长度即代表误差大小,既然是长度,那就是一个绝对值,但由于|y-yi|不方便计算,故而直接用平方来代表这个误差,误差为(y-yi)^2,总误差v即所有样本误差之和,因此最小二乘法的目标就是求出使总误差v尽量小的y。

要求总误差v的最小值,我们对v求导,使导数为0,此时对应的y恰好就是样本yi的算数平均数。

二、最小二乘法的应用

上面的测身高只展示了最小二乘法的一个简单的应用,事实上,最小二乘法为我们提供了一个思路:我们可以先假设一个“准确值y”,然后用这个准确值和其它样本,将误差表示出来,再对误差求导,以此求出误差最小时的y,这个y画成线,就可以代表样本中变量的关系。(如下图)

这里我们就用蓝线拟合了红色的点,我们也就可以用蓝线的方程来表示横轴的变量和纵轴的变量之间的关系(本例子为y=2*x)。

下面我们再举一个例子:matlab里面的polyfit函数可以用最小二乘法来对给定的数据进行拟合,其有三个参数ployfit(x,y,n)表示用最高次数为n次幂的函数来拟合y与x的关系。(如上图就是用了一次函数y=2*x来进行拟合)

我们随便输入一些数据,如代码所示:

clear all

clc

x=[0.3 0.35 0.4 0.7 0.92 1.2 1.92 2.4 2.8 3.6 4.5];

y=[1 3 2 4 5 7 6 7 6 9 10];

plot(x,y,'bo');

axis([-1 6 0 11]);

(x和y的散点图如下图)

我们对其进行拟合:

先用最高次为一次幂的方法。

一次幂就是y=a*x+b,polyfit函数求的就是a和b的值,更高次幂依次类推。代码如下:

clear all

clc

x=[0.3 0.35 0.4 0.7 0.92 1.2 1.92 2.4 2.8 3.6 4.5];

y=[1 3 2 4 5 7 6 7 6 9 10];

p1=polyfit(x,y,1) %求拟合曲线的参数。这里的1代表用1次函数拟合

y2=polyval(p1,x) %求拟合曲线下,y的值

figure;

plot(x,y,'ro');

hold on;

plot(x,y2,'b-');

axis([0 5 0 12]);

这其中,p1=polyfit(x,y,1) 是核心语句,用于求拟合曲线的参数。这里的1代表用1次函数拟合(曲线见上图)。当我们把1变为2或者更大

时。

值得注意的是,选取次数的时候,次数长度不得超过length(x)-1。

特别声明:以上内容(如有图片或视频亦包括在内)为自媒体平台“网易号”用户上传并发布,本平台仅提供信息存储服务。

Notice: The content above (including the pictures and videos if any) is uploaded and posted by a user of NetEase Hao, which is a social media platform and only provides information storage services.

/阅读下一篇/

创新杯大数据挑战赛被国内两大权威媒体齐报道!含金量直线飙升!

返回网易首页 下载网易新闻客户端
相关推荐
热点推荐
上海交大解剖500名脑梗死者,惊讶发现患脑梗的人,有6大特点

上海交大解剖500名脑梗死者,惊讶发现患脑梗的人,有6大特点

白宸侃片
2025-06-21 12:40:49
一景区突发垮塌!亲历游客:离死亡最近的一次……最新消息:部分区域暂时关闭

一景区突发垮塌!亲历游客:离死亡最近的一次……最新消息:部分区域暂时关闭

极目新闻
2025-06-22 18:49:11
一投资专家发动态挑衅特朗普:如果你敢打伊朗,我开直播表演切掉小弟弟

一投资专家发动态挑衅特朗普:如果你敢打伊朗,我开直播表演切掉小弟弟

不掉线电波
2025-06-22 14:39:37
伊外长质问:伊朗何以“重返”谈判桌?

伊外长质问:伊朗何以“重返”谈判桌?

新京报
2025-06-22 17:19:09
小伙帮邻居死亡后续:被曝内脏蒸熟 死前画面流出 求救无门很绝望

小伙帮邻居死亡后续:被曝内脏蒸熟 死前画面流出 求救无门很绝望

鋭娱之乐
2025-06-21 19:15:32
特朗普电视讲话:伊朗核设施已被“彻底摧毁”,若反击将面临更猛烈袭击

特朗普电视讲话:伊朗核设施已被“彻底摧毁”,若反击将面临更猛烈袭击

澎湃新闻
2025-06-22 11:16:29
常州一大排档因赞助苏超爆火,大量网友前来打卡“世界品牌”,老板娘忙得求饶:过几天再来吧

常州一大排档因赞助苏超爆火,大量网友前来打卡“世界品牌”,老板娘忙得求饶:过几天再来吧

极目新闻
2025-06-22 17:55:59
匿名捐款1.8亿?演戏的味道扑鼻而来

匿名捐款1.8亿?演戏的味道扑鼻而来

不主流讲话
2025-06-21 22:41:18
利好来了!A股重磅新规,明日实施!

利好来了!A股重磅新规,明日实施!

证券时报
2025-06-22 14:57:22
伊朗外长:核设施摧毁一两个也能重建,技术知识炸不垮!如果特朗普决定让美国打击伊朗,伊朗保留报复的权利

伊朗外长:核设施摧毁一两个也能重建,技术知识炸不垮!如果特朗普决定让美国打击伊朗,伊朗保留报复的权利

和讯网
2025-06-21 10:20:08
被央媒怒批,目不识丁,脑袋空空,这5位“绝望的文盲”凭啥走红

被央媒怒批,目不识丁,脑袋空空,这5位“绝望的文盲”凭啥走红

卷史
2025-06-21 19:57:06
常州五连败被“剃光头”,一景区推出“光头免票”:没头发就行,性别无限制

常州五连败被“剃光头”,一景区推出“光头免票”:没头发就行,性别无限制

极目新闻
2025-06-22 18:26:28
伊朗人民苦神棍久矣!

伊朗人民苦神棍久矣!

廖保平
2025-06-22 08:56:17
国家公布违规吃喝认定标准,从此以后官场吃喝风将彻底一去不复返

国家公布违规吃喝认定标准,从此以后官场吃喝风将彻底一去不复返

深析古今
2025-06-22 21:14:03
联合国安理会将就伊朗局势举行紧急会议

联合国安理会将就伊朗局势举行紧急会议

界面新闻
2025-06-22 21:58:50
伊朗导弹撕裂以军防空,特朗普彻底明白:解放军已经惹不起!

伊朗导弹撕裂以军防空,特朗普彻底明白:解放军已经惹不起!

文雅笔墨
2025-06-21 00:18:38
北京酒店HIV试纸后续:全是人血两条杠!客人紧急就医,细节曝光

北京酒店HIV试纸后续:全是人血两条杠!客人紧急就医,细节曝光

哄动一时啊
2025-06-22 07:55:03
快评丨不扶倒于积水老人被处理,民警的“民”字究竟该如何写?

快评丨不扶倒于积水老人被处理,民警的“民”字究竟该如何写?

红星新闻
2025-06-22 16:02:35
社科院研究员和他的“旷世奇文”

社科院研究员和他的“旷世奇文”

地球公民金建国
2025-06-21 19:42:33
外交部发布视频:回家

外交部发布视频:回家

新华社
2025-06-22 18:24:30
2025-06-23 00:07:00
数模乐园官方
数模乐园官方
专注于数学建模,分享干货知识
1214文章数 708关注度
往期回顾 全部

教育要闻

高考后留学申请哪所世界名校?独家汇总2022-2026QS全球前200名校

头条要闻

伊朗投入5万亿美元的核项目遭摧毁 哈梅内伊面临抉择

体育要闻

中国女排0-3意大利 香港站3胜1负收官

娱乐要闻

离婚四年!赵丽颖被曝新恋情惹争议

财经要闻

苹果后院起火

科技要闻

"你应该靠嘴吃饭",罗永浩自述被梁文锋劝退

汽车要闻

首搭华为双王牌/6月24日预售 岚图FREE+正式下线

态度原创

艺术
本地
手机
时尚
公开课

艺术要闻

故宫珍藏的墨迹《十七帖》,比拓本更精良,这才是地道的魏晋写法

本地新闻

被贵妃带火的“唐代顶流”,如今怎么不火了

手机要闻

消息称小米 MIX Flip 2 小折叠手机配备 5100mAh 电池

伊姐周日热推:电视剧《完美的救赎》;电视剧《书卷一梦》......

公开课

李玫瑾:为什么性格比能力更重要?