| 查看: 548 | 回复: 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人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有214人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这个积分函数能用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
