Matlab曲线polyfit拟合求参数

在已知输入和输出的数据基础上，获取输入输出具体的函数形式，可采用曲线拟合：

一：已知函数阶数

polyfit（x,y,N）

x:输入的数组 y：输出的数组  N：函数的阶数

构造X，Y。这里其实就是一阶函数。故意试下3阶：

clc;
x=1:1:10;
y=2*x;
coefficient=polyfit(x,y,3)

结果：

从右向左阶数升高。结果正确。

二：未知函数阶数

未知就求解，

构造一个复杂点的函数，然后利用误差平方向大小判断最佳的阶次。函数阶次一般都在5次以下

clc;
x=1:1:10;
y=x.^4+x.^3+2*x+6;
for i=1:5
    Coef=polyfit(x,y,i);
    y2=polyval(Coef,x);
    if sum((y2-y).^2)<0.1
        N=i;
        break
    end
end
Coef=polyfit(x,y,N)

结果如下：

收工，把x和y换成自己的数据就行了。

一般数据多的话，都是从Excel中导入，使用函数

xlsread（path，sheet，zone）；

path：路径/文件名

sheet：位于excel的第几张表

zone：区域使用A2：A5形式

clc;
set(0,'defaultAxesFontName', '<宋体>');	
X1=xlsread('data.xls','sheet1','A2:A44');
Y1=xlsread('data.xls','sheet1','B2:B44');
plot(X1,Y1);
%legend('左前轮','右前轮','左后轮','右后轮');   %右上角标注
xlabel('记数')  %x轴坐标描述
ylabel('脉冲/秒') %y轴坐标描述

第一行set，是为了防止matlab画图时出现中文乱码。