作者kabii460 (卡比獸有460公斤哦)
看板MATLAB
標題[問題]急問!解聯立方程式
時間Mon Apr 15 13:41:30 2019
>> syms nx ny nz
>> eq1=-7.26*nx+5*ny+2*nz==0;
>> eq2 eq1=5*nx-14.26*ny+3*nz==0;
>> eq3=2*nx+3*ny+11.26*nz==0;
>> eq4=nx^2+ny^2+nz^2==1;
>> lmnx=solve(eq1,eq2,eq3,eq4,nx,ny,nz)
lmnx =
struct with fields:
nx: [0x1 sym]
ny: [0x1 sym]
nz: [0x1 sym]
這是我打的code
想解三個未知數 nx ny nz
因為是初學者,還很生疏,雖然試了很多方法但還是解不出來,拜託請大家的賜教了>_<
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 223.137.186.26
※ 文章網址: https://webptt.com/m.aspx?n=bbs/MATLAB/M.1555306892.A.82D.html
1F:→ tenpoinyuki: 用矩陣寫(?) 04/16 16:10
2F:→ kabii460: 目前只想到將前面三式做矩陣列運算找出nx ny nz 的比例 04/16 16:31
3F:→ kabii460: 關係,然後再帶回第四式求解。想知道有沒有更簡單的方法 04/16 16:31
4F:→ kabii460: 可以使用>_< 04/16 16:31
5F:→ kevio: 前三個方程式不是就只有一個唯一解嗎? 不一定能滿足eq4 04/17 13:22
6F:推 erre: 這不是手算就出來了嗎? 04/18 02:24
7F:推 aQQQQQQQQQQ: overdetermined system在matlab有close form 解 04/18 22:04
8F:→ aQQQQQQQQQQ: : 將係數列成行列式 A 和各方程式的和b 04/18 22:09
9F:→ aQQQQQQQQQQ: close form solution x = b\A 04/18 22:10
10F:→ aQQQQQQQQQQ: *close form solution x =A\b 04/18 22:13
11F:→ aQQQQQQQQQQ: 發現有非線性的方程式... 可能還是要用fsolve 04/18 22:34