最小二乘多项式拟合程序matlab,最小二乘法的多项式拟合(matlab实现)
1、用最小二乘法进行多项式拟合(matlab实现)西安交通大学 徐彬华算法分析:对给定数据 (i=0 ,1,2,3,.,m),一共m+1个数据点,取多项式P(x),使函数P(x)称为拟合函数或最小二乘解,令似的 使得 其中,a0,a1,a2,an为待求未知数,n为多项式的最高次幂,由此,该问题化为求的极值问题。由多元函数求极值的必要条件:j=0,1,n得到:j=0,1,n这是一个关于a0,a1,a2,an的线性方程组,用矩阵表示如下:因此,只要给出数据点 及其个数m,再给出所要拟合的参数n,则即可求出未知数矩阵(a0,a1,a2,an)试验题1编制以函数 为基的多项式最小二乘拟合程序,并用于对下列。
2、数据作三次多项式最小二乘拟合(取权函数wi1)xi-1.0 -0.50.00.51.01.52.0yi-4.447-0.4520.5510.048-0.4470.5494.552总共有7个数据点,令m=6第一步:画出已知数据的的散点图,确定拟合参数n;x=-1.0:0.5:2.0;y=-4.447,-0.452,0.551,0.048,-0.447,0.549,4.552;plot(x,y,*)xlabel x轴ylabel y轴title 散点图hold on因此将拟合参数n设为3.第二步:计算矩阵A= 注意到该矩阵为(n+1)*(n+1)矩阵,多项式的幂跟行、列坐标(i,j)的关系为i+j。
3、-2,由此可建立循环来求矩阵的各个元素,程序如下:m=6;n=3;A=zeros(n+1);for j=1:n+1for i=1:n+1for k=1:m+1A(j,i)=A(j,i)+x(k)(j+i-2)endendend;再来求矩阵B=B=0 0 0 0;for j=1:n+1for i=1:m+1B(j)=B(j)+y(i)*x(i)(j-1)endend第三步:写出正规方程,求出a0,a1,an.B=B;a=inv(A)*B;第四步:画出拟合曲线x=-1.0:0.0001:2.0;z=a(1)+a(2)*x+a(3)*x.2+a(4)*x.3;plot(x,z) legend(离散点。
4、,y=a(1)+a(2)*x+a(3)*x.2+a(4)*x.3)title(拟合图)总程序附下:x=-1.0:0.5:2.0;y=-4.447,-0.452,0.551,0.048,-0.447,0.549,4.552;plot(x,y,*)xlabel x轴ylabel y轴title 散点图hold onm=6;n=3;A=zeros(n+1);for j=1:n+1for i=1:n+1for k=1:m+1A(j,i)=A(j,i)+x(k)(j+i-2)endendend;B=0 0 0 0;for j=1:n+1for i=1:m+1B(j)=B(j)+y(i)*x(i)(j-1)endendB=B;a=inv(A)*B;x=-1.0:0.0001:2.0;z=a(1)+a(2)*x+a(3)*x.2+a(4)*x.3;plot(x,z) legend(离散点,y=a(1)+a(2)*x+a(3)*x.2+a(4)*x.3)title(拟合图。
最小二乘多项式拟合程序matlab,最小二乘法的多项式拟合(matlab实现)相关推荐
- matlab最小二乘法拟合二次多项式,最小二乘法的多项式拟合(matlab实现)
用最小二乘法进行多项式拟合(matlab 实现) 西安交通大学 徐彬华 算法分析: 对给定数据|(斗』i=0 ,1,2,3,..,m), -共m+1个数据点,取多项式P(x), 使 战 m 刃:ITb ...
- matlab最小二乘法_基于最小二乘法的线性回归拟合
阅读本文需要的知识储备: 高等数学 概率论与数理统计 Python基础 线性回归,其实生活中有很多这样的例子,比如:票价与行车距离.服务质量之间的关系,买房时房价与面积.地域等的关系.给我们一组这样的 ...
- matlab 赋权法,基于MATLAB的非线性曲线拟合赋权法
成都航空职业技术学院学报 Journal of Chengdu Aeronautic Polytechnic 2018 年 12 月第 4 期(总第 117 期) Vol.34 No.4(Serial ...
- 多项式在matlab中的应用,matlab的应用-多项式函数及多项式拟合
matlab的应用-多项式函数及多项式拟合 Matlab 的应用- 多项式函数及多项式拟合 本节将向大家简要介绍 matlab 在多项式处理方面的应用. 多项式函数主要有: roots 求多项式的根 ...
- MatLab中的多项式拟合(lsqcurvefit函数的非线性曲线拟合)
MatLab中的多项式拟合(lsqcurvefit函数的非线性曲线拟合) lsqcurvefit函数 函数原型 x=lsqcurvefit(fun, x0, xdata, ydata) == 其中的参 ...
- MATLAB中的多项式拟合
下面来给大家介绍一下如何在MATLAB中进行多项式拟合,多项式拟合是用一个多项式展开去拟合包含数个分析格点的一小块分析区域中的所有观测点,得到观测数据的客观分析场.展开系数用最小二乘拟合确定.但此方法 ...
- 最佳二次逼近多项式MATLAB代码,数学实验“Chebshev多项式最佳一致逼近,最佳平方逼近”实验报告(内含matlab程序).doc...
数学实验"Chebshev多项式最佳一致逼近,最佳平方逼近"实验报告(内含matlab程序) 西京学院数学软件实验任务书 课程名称数学软件实验班级数0901学号0912020107 ...
- matlab求勒让德多项式零点,有没有勒让德多项式导数 零点程序
求N+1次勒让德多项式的m(m=0,1,...)阶导数零点Matlab程序 子程序: function x=jp(N,alpha,beta) n=1:N; a(1)=(alpha+beta+2)/2; ...
- 基于MATLAB 2014b的多项式除法程序
多项式运算所涉及的函数见下表 conv 实现卷积和多项式乘法 deconv 实现去卷积和多项式除法 poly 求具有指定根的多项式 polyder 多项式求导 polyeig 求多项式的特征值 pol ...
最新文章
- 伪语言PHP,PHP语言之伪静态几种做法
- 利润中心, 成本中心, 工作中心
- 知乎 开源机器学习_使用开源数据和机器学习预测海洋温度
- 【Xamarin挖墙脚系列:Xamarin.IOS机制原理剖析】
- stream().map().collect()用法
- LINUX C#开发教程
- Java 重载和重写
- dedecms php用不了,织梦DEDECMS安装360漏洞补丁之后不能够运行PHP代码的问题
- 小白学 Python 爬虫(18):Requests 进阶操作
- 水晶易表 Xcelsius 2008 安装指南 完美支持office2010
- RS485以及MODBUS知识积累
- 微信小程序PNG图片去白底
- Studio 3T for MongoDB 破解无限使用脚本
- ubuntu 18.04 配置深度学习环境
- The Sultan's Successors UVA - 167
- TCP/IP 协议详解内容总结
- 如何从Matlab中导出清晰的结果图片
- 2022-2028全球与中国智能家居安全市场现状及未来发展趋势
- 数据处理之连续数据离散化
- [Nowcoder] [Wannafly26B] 冥土追魂 [贪心]