我有一个函数g=m*(m1/m)*(m2/m)*(r^2)*rm*Wc1*cos(2*u),我需要一个程序求这个函数在在(-2*pi,2*pi)上所有的根。
其中u是自变量,其他所有参数都是已知。已知参数如下:
kx=10037;ky=10037;kP=1500;kD=1500;kE=1500;
f1=0;f2=0;fx=125;fy=125;fP=15;fD=15;fE=15;
m=100;m1=4;m2=4;M=m+m1+m2;
r=0.05;
JP=10;JD=10;JE=10;
rm=0.04;
lz=0.3;ly=0.3;
Omiga=157;
kt=0;
Wc1=abs(((Omiga^4)*(M^2))/((((Omiga^2)*M)^2)+(Omiga^2)*(fx^2)))...
*cos(atan(fx*Omiga/((Omiga^2)*M-kx)))+...
abs(((Omiga^4)*(M^2))/((((Omiga^2)*M)^2)+(Omiga^2)*(fy^2)))*cos(atan(fy*Omiga/((Omiga^2)*M-ky)))-...
abs(((Omiga^4)*(JP^2))/((((Omiga^2)*JP)^2)+(Omiga^2)*(fP^2)))*((lz^2)*M/JP)*cos(atan(fP*Omiga/((Omiga^2)*JP-kP)))-...
abs(((Omiga^4)*(JD^2))/((((Omiga^2)*JD)^2)+(Omiga^2)*(fD^2)))*((lz^2)*M/JD)*cos(atan(fD*Omiga/((Omiga^2)*JD-kD)))+...
abs(((Omiga^4)*(JE^2))/((((Omiga^2)*JE)^2)+(Omiga^2)*(fE^2)))*((ly^2)*M/JE)*cos(atan(fE*Omiga/((Omiga^2)*JE-kE)));
求大神帮我写一个程序求这个函数在在(-2*pi,2*pi)区间内所有的根