声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 942|回复: 1

[编程技巧] matlab 方程组求解

[复制链接]
发表于 2011-4-5 14:20 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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

NOHF}SV)[J`KSFJ52E380.jpg
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)
[U,V,W]=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 =

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

回复
分享到:

使用道具 举报

发表于 2011-4-7 00:30 | 显示全部楼层
不是說了"Warning: Explicit solution could not be found.":@)
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-16 09:47 , Processed in 0.062990 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表