1200字范文,内容丰富有趣,写作的好帮手!
1200字范文 > MATLAB —— polyfit()多项式曲线拟合(线性拟合/线性回归)

MATLAB —— polyfit()多项式曲线拟合(线性拟合/线性回归)

时间:2023-06-10 05:39:22

相关推荐

MATLAB —— polyfit()多项式曲线拟合(线性拟合/线性回归)

文章目录

(一)多项式曲线拟合函数(二)多项式曲线拟合函数(计算标准误差估计值)

打开工具 - 方法1:MATLAB - APP - Curve Fitting

打开工具 - 方法2:命令行窗口:cftool(Curve Fitting Tool)

多项式曲线拟合公式:p(x) = p(1)*x^n + p(2)*x^(n-1) + ... + p(n)*x + p(n+1).

定义1:用连续曲线近似地刻画或比拟平面上离散点组所表示的坐标之间的函数关系。

定义2:在数值分析中,曲线拟合就是用解析表达式逼近离散数据,即离散数据的公式化。

最小二乘法:通过" 最小化误差的平方和 "寻找与数据匹配的最佳函数,使得求得的函数(解析表达式)与实际数据之间误差的平方和为最小。

MathWorks官网:多项式曲线拟合(polyfit)

(一)多项式曲线拟合函数

p = polyfit(x, y, n);返回阶数为 n 的多项式系数,p为多项式系数向量,向量中的系数按降幂排列。该阶数是 y 中数据的最佳拟合(在最小二乘方式中)。

y = polyval(p, x);计算多项式(p为系数向量)在(x)点出得函数值(y)

clear;% 删除工作区中项目,释放系统内存clc;% 清空命令行窗口close all;% 关闭句柄可见的所有图窗%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% (1)构建原始数据x = 1:20;y = -1.2*x + 2*randn(1, length(x));%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% (2)多项式曲线拟合p = polyfit(x, y, 1);% 多项式曲线拟合f = polyval(p, x);% 多项式计算% p = polyfit(log(x), y, 2);% 多项式(对数)曲线拟合%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% (3)开始画图figure; % 创建图窗窗口scatter(x, y);% 原始数据。scatter:绘制散点图。hold on;% 添加新绘图时,保留当前绘图plot(x, f, 'Color', 'r');% 拟合曲线。plot:画图xlabel('x');% 添加 x 轴标签ylabel('y');% 添加 y 轴标签title('Polynomial curve fitting');% 添加标题legend('data', 'linear fitting');% 添加图例

(二)多项式曲线拟合函数(计算标准误差估计值)

[p, S] = polyfit(x, y, n)将返回一个结构体 S,可用作 polyval 的输入来获取误差估计值。

[y, delta] = polyval(p, x, S)使用 polyfit 生成的可选输出结构体 S 来生成误差估计值。delta 是使用 p(x) 预测 x 处的未来观测值时的标准误差估计值。

clear;% 删除工作区中项目,释放系统内存clc;% 清空命令行窗口close all;% 关闭句柄可见的所有图窗%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%x = 1:100; y = -0.3*x + 2*randn(1,100); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%[p, S] = polyfit(x, y, 1); [f, delta] = polyval(p, x, S);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 绘制原始数据、线性拟合、95%预测区间(y±2Δ)。plot(x, y, 'bo')hold onplot(x, f, 'r-')plot(x, f+2*delta, 'm--', x, f-2*delta, 'm--')title('Linear Fit of Data with 95% Prediction Interval')legend('Data','Linear Fit','95% Prediction Interval')

本内容不代表本网观点和政治立场,如有侵犯你的权益请联系我们处理。
网友评论
网友评论仅供其表达个人看法,并不表明网站立场。