程序

zhuchang1121

function dz=dafen(t,z)
m1=4;
m2=32.1;
R=0.00025;
L=0.00012;
c1=1050;
c2=2100;
k=2.5*10^7;
kc=2*10^6;
R0=0.00016;
f=0.1;
u=0.018;
c=0.00002;
w=900;
g=9.8;
e=0.000005;
%油膜力
a=(w*z(3)+2*z(2))/2*(sqrt((c^2-z(1)^2-z(3)^2)^3));
b=3*z(1)*(z(1)*z(2)+z(3)*z(4))/sqrt((c^2-z(1)^2-z(3)^2)^5);
p=(2*z(4)-w*z(1))/2*(sqrt((c^2-z(1)^2-z(3)^2)^3));
d=3*z(3)*(z(1)*z(2)+z(3)*z(4))/sqrt((c^2-z(1)^2-z(3)^2)^5);
Fx=-u*pi*R*L^3*(a+b);
Fy=-u*pi*R*L^3*(p+d);
%碰摩力
T=sqrt(z(5)^2+z(7)^2);
if T>=R0
    N=kc*(T-R0);
    Px=-N*(z(5)-f*z(7))/T;
    Py=-N*(f*z(5)+z(7))/T;
else
    N=0;
    Px=0;
    Py=0;
end
dz=zeros(8,1);
dz(1)=z(2);
dz(2)=(Fx-k*(z(1)-z(5))-c1*z(2))/m1;
dz(3)=z(4);
dz(4)=(Fy-k*(z(3)-z(7))-c1*z(4)-m1*g)/m1;
dz(5)=z(6);
dz(6)=(Px+m2*e*w^2*cos(w*t)-2*k*(z(5)-z(1))-c2*z(6))/m2;
dz(7)=z(8);
dz(6)=(Py+m2*e*w^2*sin(w*t)-m2*g-2*k*(z(7)-z(3))-c2*z(6))/m2;
自己编写的一个程序,但是运行的时候出现这个提示Warning: Imaginary parts of complex X and/or Y arguments ignored..到底是怎么回师啊!

sigma665

开根号的时候有负的

sogooda

请楼主发帖时注意标题应正确表达帖子主题，参见会员守则及本版置顶帖。