python polyfit函数怎么使用
Posted
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了python polyfit函数怎么使用相关的知识,希望对你有一定的参考价值。
参考技术A 用polyfit(X,Y,1)得到的拟合函数只能得到a,b,但不能得到线性相关系数R^2。如想要得到其线性相关系数,可以用regress(y,X),其使用格式[b,bint,r,rint,stats]
=
regress(y,X);
b——拟合系数
bint——b的置信区间
r——残差值
rint——r的置信区间
stats——检验统计量,第一个就是相关系数
例如:
x=[。。。];y=[。。。]
X=[x
ones(n,1)];
%x的行数(列数)
[b,bint,r,rint,stats]
=
regress(y,X);
MATLAB中的polyfit函数的使用方法
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(‘原始数据‘,‘拟合直线‘);
结果示例:
以上是关于python polyfit函数怎么使用的主要内容,如果未能解决你的问题,请参考以下文章
Python曲线拟合(polyfit , curve_fit, interp1d插值)