| 查看: 511 | 回复: 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似乎是颠倒的。到底什么原因啊，折磨了我一下午一晚上也没找出原因，难道是我人品不好啊，  ，求高人解答。 |
回复此楼» 猜你喜欢
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有191人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有1人回复
-
投稿chemical physical letters不送审？
已经有2人回复
-
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这个积分函数能用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