多项式拟合:
polyfit(x,y,n):对散点列(x,y)做n次多项式拟合,返回多项式系数矩阵
polyval(b, x):返回x关于系数矩阵为b的多项式对应的y
来个例子:
一次拟合:
x = [1 2 3 4 5];
y = [1.1 1.9 3.1 3.99 5.01];
plot(x, y, 'o');
hold on;
b = polyfit(x, y, 1);
ny = polyval(b, x);
plot(x, ny);
求相关系数:
corrcoef(x, y) :返回向量x,y的相关系数矩阵
R = corrcoef(x, y);
插值:
interp1(x, y, nx, 'cubic'):在点列(x,y)中插值,返回nx所对应的ny
来个例子:
x = 0:2:64;
y = [35.2 34.4 28.8 18.4 6.4 8.8 20.0 28.0 33.0 34.4 29.2 ...
19.2 7.2 8.8 19.2 28.0 33.6 33.6 28.0 19.2 8.0 6.4 ...
17.6 27.2 32.8 32.8 27.2 19.2 7.2 7.2 18.4 26.4 32.0];
plot(x, y, 'o');
hold on;
nx = 0:0.1:64;
ny = interp1(x, y, nx, 'cubic');
plot(nx, ny);