代码:
syms r
r1=34.5;
a=asin(8.2385e-6*r^3 + 0.000084205*r^2 - 0.0041311*r + 0.38649);
f=@(r)int(-tan(a)*(sin(a)+r*cos(a))*diff(a),-6.62,r);
v=@(r)(f-r*cos(a))/r1;
fun1=@(r)r*cos(a-v)+r1*sin(v);
fun2=@(r)r*sin(a-v)+r1*cos(v);
b=linspace(-6.62,6.76,15);
fn1=arrayfun(fun1,b);
fn2=arrayfun(fun2,b);
plot(fn1,fn2)
运行结果图片: