| 查看: 559 | 回复: 1 | ||
[求助]
薅不住了，matlab算出来的结果怎么这么悲剧啊
|
|
我使用matlab求一个方程组的解， 给出来的结果明显不对，求解如下： >> clear syms a b c A B C lx ly lz xy yz xz; a=6.5209; b=10.761; c=7.3062; A=90; B=102.058; C=90; [lx,ly,lz,xy,yz,xz]=solve('lx=a','ly^2+xy^2=b^2','lz^2+xz^2+yz^2=c^2','xy*xz+ly*yz =cosd(A)*b*c', 'xz= cosd(B)*c', 'xy=cosd(C)*b', 'lx', 'ly', 'lz', 'xy', 'yz', 'xz') eval(lx) eval(ly) eval(lz) eval(xy) eval(yz) eval(xz) lx = a a a a ly = (1-cosd(C)^2)^(1/2)*b (1-cosd(C)^2)^(1/2)*b -(1-cosd(C)^2)^(1/2)*b -(1-cosd(C)^2)^(1/2)*b lz = 1/(cosd(C)^2-1)*((cosd(C)^2-1)*(-2*cosd(C)*cosd(B)*cosd(A)+cosd(A)^2-1+cosd(C)^2+cosd(B)^2))^(1/2)*c -1/(cosd(C)^2-1)*((cosd(C)^2-1)*(-2*cosd(C)*cosd(B)*cosd(A)+cosd(A)^2-1+cosd(C)^2+cosd(B)^2))^(1/2)*c 1/(cosd(C)^2-1)*((cosd(C)^2-1)*(-2*cosd(C)*cosd(B)*cosd(A)+cosd(A)^2-1+cosd(C)^2+cosd(B)^2))^(1/2)*c -1/(cosd(C)^2-1)*((cosd(C)^2-1)*(-2*cosd(C)*cosd(B)*cosd(A)+cosd(A)^2-1+cosd(C)^2+cosd(B)^2))^(1/2)*c xy = cosd(C)*b cosd(C)*b cosd(C)*b cosd(C)*b yz = cosd(B)*c cosd(B)*c cosd(B)*c cosd(B)*c xz = c*(-cosd(C)*cosd(B)+cosd(A))/(1-cosd(C)^2)^(1/2) c*(-cosd(C)*cosd(B)+cosd(A))/(1-cosd(C)^2)^(1/2) -c*(-cosd(C)*cosd(B)+cosd(A))/(1-cosd(C)^2)^(1/2) -c*(-cosd(C)*cosd(B)+cosd(A))/(1-cosd(C)^2)^(1/2) ans = 6.5209 6.5209 6.5209 6.5209 ans = 10.7610 10.7610 -10.7610 -10.7610 ans = -7.1450 7.1450 -7.1450 7.1450 ans = 0 0 0 0 ans = -1.5263 -1.5263 -1.5263 -1.5263 ans = 0 0 0 0 matlab给出的解带入方程式明显不对。 这个方程式手算也很容易，和我手算得出的结果比较，matlab算出的xz、yz似乎是颠倒的。到底什么原因啊，折磨了我一下午一晚上也没找出原因，难道是我人品不好啊，  ，求高人解答。 |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有282人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这个积分函数能用matlab画出来吗？谢谢
已经有12人回复
- 怎样用matlab或者是mathcad求解含有变参数的四次方程？
已经有7人回复
- 如何改变matlab默认的变量最大数值
已经有3人回复
- 三元二次的三角方程如何用MATLAB求解？？
已经有10人回复
- matlab函数调用与数组的问题
已经有19人回复
- matlab拟合求值
已经有17人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 如何用Matlab画z关于x和y的函数图？谢谢大家！
已经有8人回复
- Matlab两个小问题
已经有10人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 计算强帖: 1
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
【答案】应助回帖
yangmou(金币+2): 谢谢解答，好人啊 2011-10-01 13:59:32
|
[lx,ly,lz,xy,yz,xz]=solve('lx=a','ly^2+xy^2=b^2','lz^2+xz^2+yz^2=c^2','xy*xz+ly*yz =cosd(A)*b*c', 'xz= cosd(B)*c', 'xy=cosd(C)*b', 'lx', 'ly', 'lz', 'xy', 'yz', 'xz') 你的这句有误； 注意未知数的顺序，应该是 'lx', 'ly', 'lz', 'xy', 'xz', 'yz' solve解出来的是按照字母顺序排列出来的，而你自己列出来的顺序 'yz', 'xz'颠倒了，所以显示也就颠倒了 |
2楼2011-10-01 09:53:25
