xiaoxue 2020-05-09
??在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(‘原始数据‘,‘拟合直线‘);
结果示例: