通过已知的点求拟合曲线及方程

1、编辑器编译：

p1 = polyfit(xr,yr,5); % 使用 polyfit函数 拟合横、纵轴数据得到拟合多项式，储存在p中
y5 = polyval(p1,xr); % polyval用于计算出每个横轴坐标(x)在拟合多项式(p)中对应的函数值
p1=vpa(poly2sym(p1),5);%显示5阶多项式

2、在命令行窗口输入p5，按回车，得到5阶多项式

3、要使横纵坐标的单位刻度是相同的，则在plot函数下一行输入axis equal；设置坐标轴范围axis([xmin,xmax,ymin,ymax]);

figure(2);%创建图形窗
plot(xr,yr,'o',xr,y5,'r:');%o是标记符号圆圈，r是颜色符号
axis equal;%使横纵坐标的单位刻度相同
axis([0,100,0,80]);%设置坐标轴范围

4、求曲率

h1 = abs(diff(xr)) ;%对xr求一阶导数，绝对值
h = [h1 h1(end)];
ht = h;
yapp1 = gradient(yr)./ht;
yapp2 = del2(yr)./ht;           
%.^:按元素求幂,即计算 A 中每个元素在 B 中对应指数的幂。A 和 B 的大小必须相同或兼容
kappar = abs(yapp2)./(1+yapp1.^2).^(3/2);%前面加.表示"按元素"