[問題] 解八元非線性聯立方程組~>"<(已爬文)已回收
我是要進行最小平方法求解未知變數
可是執行完後卻出現以下訊息:
Error in ==>C:\MATLAB6p5\toolbox\symbolic\maplemex.dll
Error in ==> C:\MATLAB6p5\toolbox\symbolic\maple.m
On line 116 ==> [result,status] = maplemex(statement);
Error in ==> C:\MATLAB6p5\toolbox\symbolic\solve.m
On line 115 ==> [R,stat] = maple('solve',eqns,vars);
Error in ==> C:\MATLAB6p5\toolbox\symbolic\@sym\solve.m
On line 49 ==> [varargout{1:max(1,nargout)}] = solve(S{:});
實在不清楚哪裡有問題了
身為matlab新手的我 爬文後還是沒頭緒>"<
請各位高手指點小妹一下
以下為我的程式碼
====================================================
syms Xb Yb Xc Yc Xd Yd Xe Ye;
Vab=sym('((Xb-1000)^2+(Yb-1000)^2)^0.5-319.236');
Vac=sym('((Xc-1000)^2+(Yc-1000)^2)^0.5-550.933');
Vad=sym('((Xd-1000)^2+(Yd-1000)^2)^0.5-610.974');
Vae=sym('((Xe-1000)^2+(Ye-1000)^2)^0.5-407.268');
Vbc=sym('((Xc-Xb)^2+(Yc-Yb)^2)^0.5-276.962');
Vbd=sym('((Xd-Xb)^2+(Yd-Yb)^2)^0.5-466.565');
Vbe=sym('((Xe-Xb)^2+(Ye-Yb)^2)^0.5-404.487');
Vcd=sym('((Xd-Xc)^2+(Yd-Yc)^2)^0.5-285.577');
Vce=sym('((Xe-Xc)^2+(Ye-Yc)^2)^0.5-393.802');
Vde=sym('((Xe-Xd)^2+(Ye-Yd)^2)^0.5-255.056');
S=[0.026 0.024 0.025 0.030 0.022 0.015 0.036 0.016 0.025 0.022];
Q=S.*S;
W=1./Q;%權矩陣
PHI=W(1)*Vab^2+W(2)*Vac^2+W(3)*Vad^2+W(4)*Vae^2+W(5)*Vbc^2+W(6)*Vbd^2
+W(7)*Vbe^2+W(8)*Vcd^2+W(9)*Vce^2+W(1)*Vde^2;
dPHI_dXb=diff(PHI,Xb);
dPHI_dYb=diff(PHI,Yb);
dPHI_dXc=diff(PHI,Xc);
dPHI_dYc=diff(PHI,Yc);
dPHI_dXd=diff(PHI,Xd);
dPHI_dYd=diff(PHI,Yd);
dPHI_dXe=diff(PHI,Xe);
dPHI_dYe=diff(PHI,Ye);
[Xb,Yb,Xc,Yc,Xd,Yd,Xe,Ye]=solve(dPHI_dXb,dPHI_dYb,dPHI_dXc,dPHI_dYc,
dPHI_dXd,dPHI_dYd,dPHI_dXe,dPHI_dYe)
--
為你的難過而快樂的,是敵人,為你的快樂而快樂的,是朋友,
而為你的難過而難過的,就是那些,該放進心裡的人。
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.66.130.37
※ 編輯: tingyili 來自: 61.66.130.37 (08/19 17:27)
推
08/19 20:15, , 1F
08/19 20:15, 1F
→
08/19 20:23, , 2F
08/19 20:23, 2F
→
08/19 20:24, , 3F
08/19 20:24, 3F
推
08/20 12:32, , 4F
08/20 12:32, 4F
→
08/21 21:51, , 5F
08/21 21:51, 5F
推
08/21 23:53, , 6F
08/21 23:53, 6F
推
08/21 23:55, , 7F
08/21 23:55, 7F
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章