fz.m

来自「对超声电机摩擦驱动的仿真」· M 代码 · 共 14 行

function f=f(x0)
Fn=135;
br=0.002;
hr=0.0005;
Er=5*10^8;
kr=Er*br/hr;
n=9;
e0=0.8;
Ap=n*e0;
R=0.03;
Rav=0.95*R;
rumda=2*pi*Rav/n;
k=2*pi/rumda;
f=Fn-2*Rav*kr*e0*(sin(k*x0)-k*x0*cos(k*x0));