mengxiao6861om 发表于 2011-4-5 14:20

matlab 方程组求解

本帖最后由 mengxiao6861om 于 2011-4-5 14:59 编辑


syms U V W U1 V1 W1 bb1 bb2 bb3 B1 B2 h
eq1=h*(-V1*W-V*U*W1)*B2+h*V*U1*B1-bb1;
eq2=h*U1*W1*B2+h*U*B1-bb2;
eq3=h*(V*W-V1*U*W1)*B2+h*V1*U1*B1-bb3;
eq1=subs(eq1,{h,B1,B2},{0.75,sin(55.8/180*pi),cos(55.8/180*pi)});
eq2=subs(eq2,{h,B1,B2},{0.75,sin(55.8/180*pi),cos(55.8/180*pi)});
eq3=subs(eq3,{h,B1,B2},{0.75,sin(55.8/180*pi),cos(55.8/180*pi)});
eq1=subs(eq1,{U1,V1,W1},{(1-U^2)^(1/2),(1-V^2)^(1/2),(1-W^2)^(1/2)})
eq2=subs(eq2,{U1,V1,W1},{(1-U^2)^(1/2),(1-V^2)^(1/2),(1-W^2)^(1/2)})
eq3=subs(eq3,{U1,V1,W1},{(1-U^2)^(1/2),(1-V^2)^(1/2),(1-W^2)^(1/2)})
eq1=vpa(eq1,4)
eq2=vpa(eq2,4)
eq3=vpa(eq3,4)
=solve(eq1,eq2,eq3,'U,V,W')


eq1 =

-.4216*(1.-1.*V^2)^(1/2)*W-.4216*V*U*(1.-1.*W^2)^(1/2)+.6203*V*(1.-1.*U^2)^(1/2)-1.*bb1



eq2 =

.4216*(1.-1.*U^2)^(1/2)*(1.-1.*W^2)^(1/2)+.6203*U-1.*bb2



eq3 =

.4216*V*W-.4216*(1.-1.*V^2)^(1/2)*U*(1.-1.*W^2)^(1/2)+.6203*(1.-1.*V^2)^(1/2)*(1.-1.*U^2)^(1/2)-1.*bb3



bb1 bb2 bb3是符号常量,请问怎么求出三个未知数(U V W )的解啊?
我用solve求出来是Warning: Explicit solution could not be found.
U =

[ empty sym ]



V =

   []


W =

   []
希望高手给解决下,我都困扰了好久了!

ChaChing 发表于 2011-4-7 00:30

不是說了"Warning: Explicit solution could not be found.":@)
页: [1]
查看完整版本: matlab 方程组求解