用matlab求出最小曲率半径曲线方程实例:
y=(-5.324e-015*x.^7)+(-8.452e-012*x^6)+(-5.204e-009*x^5)+(-1.558 e-006*x^4)+(-0.0002253*x^3)+(-0.01515*x^2)+(-1.464*x)+222.9
求在区间(0,2)之间的最小曲率半径。
Clear all
clc
syms x;
f=(-5.324e-15*x^7)+(-8.452e-12*x^6)+(-5.204e-9*x^5)+(-1.558e-6*x^4)+(-0.0002253* x^3)+(-0.01515*x^2)+(-1.464*x)+222.9;
f1=diff(f);
f2=diff(f1);
g=(1+f1^2)^(3/2)/abs(f2);
F=@(x)eval(g);%将括号内的字符视为语句运行
[X,fval]=fminbnd(F,0,2)%X为最小曲率半径事的x的值,fval为最小曲率半径