MATLAB的数值计算(三)曲线拟合与插值

一、曲线拟合

在许多实际应用中,常常只能测得一些分散的数据点。为了能从这些分散的数据点中找出它们内在的变化规律,人们就依据这些分散的数据点,运用最小二乘法,生成一个函数来逼近这些已知的数据点,称为曲线拟合。曲线拟合对应的命令是polyfit,调用格式为:

polyfit(x,y,n)

其中,x,y为样本数据;n为拟合多项式的阶数;返回值是一个多项式系数的行向量。

在使用polyfit函数的时候,x,y为已知样本数据,n需要设计者自己确定。如果选择n=1,即把样本数据作为直线对待,称为线性回归。

例:已知热敏电阻实验数据如下表所示,利用polyfit函数求其最小二乘法拟合曲线,并计算其误差程度。

|-------------------------------------------------------|------|------|------|------|------|
| 温度t/() | 20.5 | 32.7 | 51.0 | 73.0 | 95.7 |
| 电阻 | 765 | 826 | 873 | 942 | 1032 |

为了对比拟合的多项式与原来样本数据的拟合程度,通过绘图进行观察。

误差计算有很多种方法,在这里使用均方差(Mean Square Error,MSE)的方法来计算,具体计算公式为

其中,n为样本数量,为拟合值,y为样本值。

误差的计算代码如下:

方法1:直接使用一次曲线进行计算。

方法2:使用二次曲线进行计算。

二、插值

插值是在已知数据点之间计算估计值。通过这些插值,可以将无规律、没有联系的数据点,连接得到一个近似的曲线。对应一、二维插值的调用命令如下:

yi = interp1(x,y,xi,'method')

zi =interp2 (x,y,z,xi,yi,'method')

其中,x,y为样本数据;xi是计算的x的位置;method为所用的插值方法,无论是一维插值还是二维插值,method有以下4种选择方法。

|----|---------|------------|--------|-------|
| 方法 | Nearest | Linear | Spline | Cubic |
| 说明 | 邻近点插值 | 线性插值(默认方法) | 3次样条插值 | 立方插值 |

例:通过下表中的已知数据求自变量x在[0:0.25:10]等位置的插值,并且画出相应的曲线。

|---|---|--------|--------|--------|---------|---------|---------|--------|--------|--------|---------|
| x | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
| y | 0 | 0.8415 | 0.9039 | 0.1411 | -0.7568 | -0.9589 | -0.2794 | 0.6570 | 0.9894 | 0.4121 | -0.5440 |

相关推荐
我是一只小青蛙88824 分钟前
JavaScript DOM操作全解析
开发语言·javascript·ecmascript
软件算法开发8 小时前
基于改进麻雀优化的LSTM深度学习网络模型(ASFSSA-LSTM)的一维时间序列预测算法matlab仿真
深度学习·matlab·lstm·一维时间序列预测·改进麻雀优化·asfssa-lstm
你怎么知道我是队长8 小时前
C语言---头文件
c语言·开发语言
期待のcode8 小时前
Java虚拟机的运行模式
java·开发语言·jvm
hqwest9 小时前
码上通QT实战25--报警页面01-报警布局设计
开发语言·qt·qwidget·ui设计·qt布局控件
a程序小傲9 小时前
京东Java面试被问:动态规划的状态压缩和优化技巧
java·开发语言·mysql·算法·adb·postgresql·深度优先
HellowAmy9 小时前
我的C++规范 - 玩一个小游戏
开发语言·c++·代码规范
徐先生 @_@|||10 小时前
Palantir Foundry 五层架构模型详解
开发语言·python·深度学习·算法·机器学习·架构
tang7778910 小时前
爬虫如何绕过绕过“5秒盾”Cloudflare:从浏览器指纹模拟到Rust求解之不完全指南
开发语言·爬虫·rust·cloudflare
Yuer202510 小时前
什么是 Rust 语境下的“量化算子”——一个工程对象的最小定义
开发语言·后端·rust·edca os·可控ai