matlab求解平面方程的原理

2021-05-27 19:30:07 字數 469 閱讀 8981

已知三點p1,p2,p3,求其平面方程

syms x y z

%p1,p2,p3的座標由自己定義。

p1=[x1,y1,z1];

p2=[x2,y2,z2];

p3=[x3,y3,z3];

%%那麼求解下面矩陣q行列式就是了

q=[ones(4,1),[x,y,z];p1;p2;p3]; %%

detb=det(q);

%%最後令

q=0;

這裡的求解方法是四階矩陣的行列式求解,可以從三階矩陣的行列式推導出來,那就是常規的求解方法:方程組寫成[p1;p2;p3]=d,其中d=[d,d,d]'表示平面方程的常數項,先求q=[p1;p2;p3]的行列式,再分別求d與三個座標兩個分量合成的三階矩陣的行列式,設為q1,q2,q3;那麼所求的平面方程的四個引數就是a=q1/q;b=q2/q;c=q3/q;d=d;同時乘以q/d所得結果本質上就是上面四階矩陣行列式的求解了

matlab求解振動方程

看了一篇柱塞幫浦離散化動力學建模的文章,感覺還挺有意思,於是嘗試做一下 二 matlab下的動力學方程總結 斜盤式軸向柱塞幫浦是一類常見的柱塞幫浦,本文以 型斜盤式軸向柱塞幫浦為研究物件,研究幫浦內機械振動的傳遞問題。由於該幫浦傳動軸與缸體之間為過盈配合,且柱塞滑靴元件位於缸體的柱塞腔內,因此,將傳...

利用Matlab求解Laplace方程

問題描述 求解泊松方程 u 1 求解區域為單位圓盤,邊界條件為在圓盤邊界上u 0 下面求它的數值解,編寫程式如下 1 問題定義 g circleg 單位圓 b circleb1 邊界上為零條件 c 1 a 0 f 1 2 產生初始的三角形網格 p,e,t initmesh g 3 迭代直至得到誤差允...

MATLAB 求解方程(組)

eg.解方程x 2 x 2 0 1.roots p 函式 此 matlab 函式 以列向量的形式返回 p 表示的多項式的根。輸入 p 是乙個包含 n 1 多項式係數的向量,以 xn 係數開頭。0係數表示方程中不存在的中間冪。p 1 1,2 x roots p 2.solve函式 利用solve函式求...