线性最小二乘法拟合 matlab程序,曲线拟合的线性最小二乘法及其MATLAB程序
1 曲线拟合的线性最小二乘法及其MATLAB 程序
例7.2.1 给出一组数据点),(i i y x 列入表7–2中,试用线性最小二乘法求拟合曲线,并用(7.2),(7.3)和(7.4)式估计其误差,作出拟合曲线.
表7–2 例7.2.1的一组数据),(y x
解 (1)在MATLAB 工作窗口输入程序
>> x=[-2.5 -1.7 -1.1 -0.8 0 0.1 1.5 2.7 3.6];
y=[-192.9 -85.50 -36.15 -26.52 -9.10 -8.43 -13.12 6.50 68.04];
plot(x,y,'r*'),
legend('实验数据(xi,yi)')
xlabel('x'), ylabel('y'),
title('例7.2.1的数据点(xi,yi)的散点图')
运行后屏幕显示数据的散点图(略).
(3)编写下列MATLAB 程序计算)(x f 在),(i i y x 处的函数值,即输入程序
>> syms a1 a2 a3 a4
x=[-2.5 -1.7 -1.1 -0.8 0 0.1 1.5 2.7 3.6];
fi=a1.*x.^3+ a2.*x.^2+ a3.*x+ a4
运行后屏幕显示关于a 1,a 2, a 3和a 4的线性方程组
fi =[ -125/8*a1+25/4*a2-5/2*a3+a4,
-4913/1000*a1+289/100*a2-17/10*a3+a4,
-1331/1000*a1+121/100*a2-11/10*a3+a4,
-64/125*a1+16/25*a2-4/5*a3+a4,
a4, 1/1000*a1+1/100*a2+1/10*a3+a4,
27/8*a1+9/4*a2+3/2*a3+a4, 19683/1000*a1+729/100*a2+27/10*a3+a4, 5832/125*a1+324/25*a2+18/5*a3+a4]
编写构造误差平方和的MATLAB 程序
>> y=[-192.9 -85.50 -36.15 -26.52 -9.10 -8.43 -13.12 6.50 68.04];
fi=[-125/8*a1+25/4*a2-5/2*a3+a4,
-4913/1000*a1+289/100*a2-17/10*a3+a4,
-1331/1000*a1+121/100*a2-11/10*a3+a4,
-64/125*a1+16/25*a2-4/5*a3+a4, a4, 1/1000*a1+1/100*a2+1/10*a3+a4,
27/8*a1+9/4*a2+3/2*a3+a4,
19683/1000*a1+729/100*a2+27/10*a3+a4,
5832/125*a1+324/25*a2+18/5*a3+a4];
fy=fi-y; fy2=fy.^2; J=sum(fy.^2)
运行后屏幕显示误差平方和如下
J=
(-125/8*a1+25/4*a2-5/2*a3+a4+1929/10)^2+(-4913/1000*a1+2
89/100*a2-17/10*a3+a4+171/2)^2+(-1331/1000*a1+121/100*a2-11/10*a3+a4+723/20)^2+(-64/125*a1+16/25*a2-4/5*a3+a4+663/25)^2+(a4+91/10)^2+(1/1000*a1+1/100*a2+1/10*a3+a4+843/100)^2+(27/8*a1+9/4*a 2+3/2*a3+a4+328/25)^2+(19683/1000*a1+729/100*a2+27/10*a3+a4-13/
2)^2+(5832/125*a1+324/25*a2+18/5*a3+a4-1701/25)^2
为求4321,,,a a a a 使J 达到最小,只需利用极值的必要条件0=??k
a J )4,3,2,1(=k ,
线性最小二乘法拟合 matlab程序,曲线拟合的线性最小二乘法及其MATLAB程序相关推荐
- c++椭圆最小二乘法原理_利用最小二乘法拟合椭圆方程的理论推导,附有matlab代码...
为了很好的进行椭圆方程拟合,本文先对椭圆基本知识进行复习,后进行非标准椭圆方程拟合公式推导,最后有matlab代码的实现. 1. 用最小二乘法做椭圆拟合 1.1. 椭圆标准方程 对椭圆印象最深的就是高 ...
- 插值与拟合 (二) : 曲线拟合的线性最小二乘法
目录 1 线性最小二乘法 最小二乘准则 系数 的确定 1.2 函数 的选取 常用的拟合曲线:直线.多项式曲线.双曲线.指数曲线 2 最小二 ...
- Matlab 隐函数方程求解最小二乘法拟合一阶线性拟合二阶拟合传感器实验
九层妖塔 起于垒土 Matlab 最小二乘法拟合一阶线性拟合&传感器实验 一.代码 二.数据处理结果 三.Notes 一.代码 %电容传感器位移实验数据 最小二乘法一阶线性拟合 x = ...
- matlab最小二乘法拟合图旋转,【Matlab】—{最小二乘法拟合一阶线性拟合传感器实验}...
[Matlab]-{最小二乘法拟合一阶线性拟合传感器实验} [Matlab]-{最小二乘法拟合一阶线性拟合&传感器实验} ???九层妖塔?起于垒土 [Matlab]-{最小二乘法拟合一阶线性拟 ...
- matlab 最小二乘法拟合_高中数学:最小二乘法与线性回归方程
1.怎样的拟合直线最好?--与所有点都近,即与所有点的距离之和最小. 最小二乘法可以帮助我们在进行线性拟合时,如何选择"最好"的直线.要注意的是,利用实验数据进行拟合时,所用数据的 ...
- 最小二乘法线性拟合和2次曲线拟合算法
最近由于项目要求,应用了最小二乘法线性拟合和2次曲线拟合算法,现总结如下: 最小二乘法线性拟合应用已有的采样时间点,再现这些点所描述的线性变化,即求出一个线性方程y=ax+b(这个算法的主要问题也就是 ...
- matlab直线拟合的程序,MATLAB最小二乘法拟合直线的程序
最小二乘法拟合直线 程序: function linear_fit %最小二乘法拟合直线clear; clc; prompt={'Name of data file'}; title='Linear_ ...
- 曲线拟合的线性最小二乘法
最小二乘法拟合 最小二乘法拟合解方程组方法多项式拟合 解方程组方法 栗子:最小二乘法求一个形如: 的经验公式. x = [19 25 31 38 44]'; y = [19.0 32.3 49.0 ...
- 【Matlab】最小二乘法拟合多项式
前言 在最近的电机项目中,有遇到有传感器数据并不线性的问题,然后想要用最小二乘法做个曲线拟合,反过来去校准不线性的传感器的数据,因此记录一下使用最小二乘法来拟合多项式的曲线的步骤.本篇从最小二乘法的原 ...
最新文章
- SQLplus 和mysql区别_mysql和oracle的区别有哪些
- 【easy】349. Intersection of Two Arrays
- 洛谷1008 三连击
- CSMA/CA与CSMA/CD
- Intel IPP密码库 IPPCP 2018 开发笔记与总结(全)
- credential provider filter注意
- 空手套白狼之——兼职中介
- 在elementUI中sort-orders排序,默认为三种,怎么改成两种
- Solidworks二次开发平台 --- RyS.SwWorks [2015-09-18更新]
- R语言和Julia以及Python使用Feather包共享数据
- 微信公众号开发系统入门教程之公众号与小程序、APP的区别
- 小精灵无尽的长廊_神奇宝贝之无尽系统
- 引领企业电销革新,外呼系统是不可缺的电销工具
- python 变量聚类 proc varclus_使用SAS进行简单的聚类分析讲解(转载)
- 一头扎进Mysql视频教程 + 源码
- 【干货】关于电商售后(退货、退款、换货)的相关流程
- windows中Linux ubuntu子系统 安装talib 编译talib python 安装 talib 报错 - 简书
- 广度优先遍历类似于二叉树的_深度优先遍历类似于二叉树的()
- C语言 一条语句max = x>y? (x>z?x:z) : (y>z?y:z)实现求 x, y, z 最大值
- excel两个表格数据对比_快速对比excel表格两个sheet中不同的数据,极大的降低重复工作量...