[討論] solve方程組

看板MATLAB作者saniao (saniao)時間10年前 (2015/05/05 21:14)推噓0(0推 0噓 7→)

留言7則, 4人參與討論串1/1

請問一下我今天有3個未知數a,b,c待求定亦有三個方程式eq1,eq2,eq3 我使用symbolic進行解方程組為什麼每個解都是8x1 syms呢小弟數學不好請見諒 code如下 syms a b c eq1=(1+a^2-b^2-c^2)/(1+a^2+b^2+c^2)-0.99768250; eq2=(1-a^2+b^2-c^2)/(1+a^2+b^2+c^2)-0.75500072; eq3=(1-a^2-b^2+c^2)/(1+a^2+b^2+c^2)-0.75731587; [a,b,c]=solve(eq1,eq2,eq3,'a','b','c') ------ 結果如下 a = -8638450155872395594504574183690^(1/2)/7903815296897390 8638450155872395594504574183690^(1/2)/7903815296897390 -8638450155872395594504574183690^(1/2)/7903815296897390 8638450155872395594504574183690^(1/2)/7903815296897390 -8638450155872395594504574183690^(1/2)/7903815296897390 8638450155872395594504574183690^(1/2)/7903815296897390 -8638450155872395594504574183690^(1/2)/7903815296897390 8638450155872395594504574183690^(1/2)/7903815296897390 b = -10456213264794931385910795^(1/2)/3951907648448695 -10456213264794931385910795^(1/2)/3951907648448695 10456213264794931385910795^(1/2)/3951907648448695 10456213264794931385910795^(1/2)/3951907648448695 -10456213264794931385910795^(1/2)/3951907648448695 -10456213264794931385910795^(1/2)/3951907648448695 10456213264794931385910795^(1/2)/3951907648448695 10456213264794931385910795^(1/2)/395190764844869 c = -82451023629991619035344171910^(1/2)/7903815296897390 -82451023629991619035344171910^(1/2)/7903815296897390 -82451023629991619035344171910^(1/2)/7903815296897390 -82451023629991619035344171910^(1/2)/7903815296897390 82451023629991619035344171910^(1/2)/7903815296897390 82451023629991619035344171910^(1/2)/7903815296897390 82451023629991619035344171910^(1/2)/7903815296897390 82451023629991619035344171910^(1/2)/7903815296897390 還請板上高手指點迷津是我用的function有錯嗎還是正解就是這樣= _= -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.116.47.32 ※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1430831654.A.978.html