[問題] 解八元非線性聯立方程組~>"<(已爬文)已回收

看板MATLAB作者 (Don't tell)時間16年前 (2009/08/19 17:26), 編輯推噓4(403)
留言7則, 3人參與, 最新討論串1/1
我是要進行最小平方法求解未知變數 可是執行完後卻出現以下訊息: 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
dsolve有解出收斂的答案?
08/21 23:55, 7F
文章代碼(AID): #1AYyL6rp (MATLAB)
文章代碼(AID): #1AYyL6rp (MATLAB)