码迷,mamicode.com
首页 > 其他好文 > 详细

MATLAB中的polyfit函数的使用方法

时间:2020-05-09 21:29:30      阅读:182      评论:0      收藏:0      [点我收藏+]

标签:描述   数学   proc   atl   返回   数据集   blog   大小   遇到   

MATLAB中的polyfit函数的使用方法

??在MATLAB中polyfit函数是用来进行多项式拟合的。其数学原理是基于最小二乘法进行拟合的。具体使用语法是:

p = polyfit(x,y,n);
% 其中x,y表示需要拟合的坐标点,大小需要一样; n表示多项式拟合的次数。
% 返回值p表示多项式拟合的系数,系数从高到低排列

具体用法示例:

1、使用polyfit函数拟合一次多项式,示例:

%% polyfit函数的使用
clear; clc; close all;
% 原始数据
x = 1:20;
y = [1,12,23,46,78,98,100,123,160,210,230,270,292,350,346,386,438,349,460,512];

p = polyfit(x,y,1);  % 进行拟合
y1 = x*p(1)+p(2);  % 生成表达式,计算y的值

figure;  % 绘图
scatter(x,y,‘.‘);
hold on;
plot(x,y1,‘Color‘,‘r‘);
% 添加说明
xlabel(‘x‘);ylabel(‘y‘);title(‘自己构建表达式‘);
legend(‘原始数据‘,‘拟合直线‘);

结果示例:

技术图片

2、在一般情况下使用polyfit函数的同时也会使用polyval函数,其中polyval函数是用来计算多项式的结果的。这样就不需要自己手动来计算多项式的结果了。其中polyval函数的使用方法如下:

p = polyfit(x,y,1);
y1 = polyval(p,x);
% 其中p就是使用polyfit函数拟合数据集x,y之后的结果,p是一个向量。
% 结果y1表示使用拟合多项式系数p来求出拟合结果y1。
具体示例:

clear; clc; close all;
% 原始数据
x = 1:20;
y = [1,12,23,46,78,98,100,123,160,210,230,270,292,350,346,386,438,349,460,512];

p = polyfit(x,y,1);  % 进行拟合
y1 = polyval(p,x);  % 生成表达式,计算y的值

figure;  % 绘图
scatter(x,y,‘.‘);
hold on;
plot(x,y1,‘Color‘,‘r‘);
% 添加说明
xlabel(‘x‘);ylabel(‘y‘);title(‘使用polyval函数‘);
legend(‘原始数据‘,‘拟合直线‘);

结果示例:

技术图片

3、在一些特殊的情况下,当x的坐标值不是从小到大进行排列的时候,只要是x,y的坐标点一一对应的,都可以使用polyfit函数进行拟合处理,示例:

%% 当拟合的数据x的坐标不是排列好顺序的时候同样可以进行拟合
clear; clc; close all;
% 遇到的数据的x轴的坐标是乱序的,不是按照从大到小进行排序,每个点是一一对应的
x = [1,2,5,7,3,6,9,10,4,8,15,16,13,11,17,12,19,18,20];
y = [1,12,78,100,23,98,160,210,46,123,346,386,292,230,438,270,460,349,512];

% 当遇到x的点是乱序的时候,同样的可以使用polyfit函数直接进行拟合
p = polyfit(x,y,1);
y1 = polyval(p,x);

figure;  % 绘图
scatter(x,y,‘.‘);
hold on;
plot(x,y1,‘Color‘,‘r‘);
% 添加说明
xlabel(‘x‘);ylabel(‘y‘);title(‘x乱序‘);
legend(‘原始数据‘,‘拟合直线‘);

结果示例:

技术图片

MATLAB中的polyfit函数的使用方法

标签:描述   数学   proc   atl   返回   数据集   blog   大小   遇到   

原文地址:https://www.cnblogs.com/zhicungaoyuan-mingzhi/p/12859769.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!