- 积分
- 82
- 注册时间
- 2004-5-11
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2004-7-22 10:12:00
|
显示全部楼层
来自 江西景德镇
回复: 【讨论】求曲线最小曲率半径的坐标
x =[ 12.5 25 50 100 200 400 800 1698.243652 3200 ];
y =[ 1.132 1.128 1.110 1.074 1.013 0.940 0.843 0.734 0.634 ];
p=polyfit(x,y,2);
vpa(p)
plot(x,y)
hold on
x=12.5:10:3200;
y=0.666010537473e-7*x.^2-0.359694992821e-3*x+1.11222063984;
plot(x,y,'r')
syms x
y=0.666010537473e-7*x.^2-0.359694992821e-3*x+1.11222063984;
k2=diff(y,2);
k1=diff(y,1);
k=k2./((1+k1.^2).^(3/2));
x =[ 12.5 25 50 100 200 400 800 1698.243652 3200 ];
kk=subs(k);
kk=double(kk);
mk=max(kk)
以上仅作参考,不是你需要的。你的数据间隔太大,并且数据太少,单凭这么几个数据很难有好的结果,必须拟合,但按二乘法拟合的精度不高,结果不会很满意。希望你对自己的表达式有个大概的了解,是不是指数或其它。 |
|