没有合适的资源?快使用搜索试试~ 我知道了~
最小二乘法曲线拟合_原理及matlab实现.docx
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 189 浏览量
2022-11-03
09:02:40
上传
评论
收藏 353KB DOCX 举报
温馨提示
试读
12页
...
资源推荐
资源详情
资源评论
曲线拟合(curve-fitting):
工程实践中,用测量到的一些离散的数据
{(x , y ),i 0,1, 2,...m}求一个近似的函数 ( ) 来拟合这组数据,要求所得的拟合曲
x
i
i
线能最好的反映数据的基本趋势(即使 ( ) 最好地逼近 f x ,而不必满足插值
x
(x )
= y ,只要使 (x ) y 尽可能地小)。
原则。因此没必要取
i
i
i
i
i
原理:
给定数据点
{(x , y ),i
0,1,
2,...
m}
。求近似曲线 x 。并且使得近似曲线与 f x 的偏差最
( )
i
i
( )
x y
小。近似曲线在该点处的偏差
,i=1,2,...,m。
i
i
i
常见的曲线拟合方法:
1.使偏差绝对值之和最小
2.使偏差绝对值最大的最小
3.使偏差平方和最小
最小二乘法:
按偏差平方和最小的原则选取拟合曲线,并且采取二项式方程为拟合曲线的方法,称为
最小二乘法。
推导过程:
1. 设拟合多项式为:
(x) a a x ... a x
k
0
1
k
2. 各点到这条曲线的距离之和,即偏差平方和如下:
3. 问题转化为求待定系数a ... a 对等式右边求a 偏导数,因而我们得到了:
0
k
i
.......
4、 把这些等式化简并表示成矩阵的形式,就可以得到下面的矩阵:
5. 将这个范德蒙得矩阵化简后可得到:
6. 也就是说 X*A=Y,那么 A = (X'*X)-1*X'*Y,便得到了系数矩阵A,同时,我们也就得
到了拟合曲线。
MATLAB 实现:
MATLAB 提供了 polyfit()函数命令进行最小二乘曲线拟合。
调用格式:p=polyfit(x,y,n)
[p,s]= polyfit(x,y,n)
[p,s,mu]=polyfit(x,y,n)
x,y 为数据点,n 为多项式阶数,返回 p 为幂次从高到低的多项式系数向量 p。x
必须是单调的。矩阵 s 包括 R(对 x 进行 QR 分解的三角元素)、df(自由度)、
normr(残差)用于生成预测值的误差估计。
[p,s,mu]=polyfit(x,y,n)在拟合过程中,首先对 x 进行数据标准化处理,以在拟合
中消除量纲等影响,mu 包含标准化处理过程中使用的 x 的均值和标准差。
polyval( )为多项式曲线求值函数,调用格式: y=polyval(p,x)
[y,DELTA]=polyval(p,x,s)
y=polyval(p,x)为返回对应自变量 x 在给定系数 P 的多项式的值。
[y,DELTA]=polyval(p,x,s) 使用 polyfit 函数的选项输出 s 得出误差估计 Y DELTA。它
假设 polyfit 函数数据输入的误差是独立正态的,并且方差为常数。则Y DELTA将
至少包含 50%的预测值。
如下给定数据的拟合曲线:
x=[,,,,,],
y=[,,,,,]。
解:MATLAB 程序如下:
剩余11页未读,继续阅读
资源评论
不吃鸳鸯锅
- 粉丝: 8292
- 资源: 2万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- Picasso_v3.1 2.ipa
- chromedriver-mac-arm64.zip
- 蓝zapro.apk
- chromedriver-linux64.zip
- UCAS研一深度学习实验-MNIST手写数字识别python源码+详细注释(高分项目)
- 基于Python和PyTorch框架完成的一个手写数字识别实验源码(带MINIST手写数字数据集)+详细注释(高分项目)
- 基于Matlab在MNIST数据集上利用CNN完成手写体数字识别任务,并实现单层CNN反向传播算法+源代码+文档说明(高分项目)
- NVIDIA驱动、CUDA和Pytorch及其依赖
- 基于SVM多特征融合的微表情识别python源码+项目说明+详细注释(高分课程设计)
- html动态爱心代码一(附源码)
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功