[問題]+[討論]用fsolve找到不同參數的解 已回收

看板MATLAB作者 (^^)時間17年前 (2009/01/09 04:31), 編輯推噓0(000)
留言0則, 0人參與, 最新討論串1/1
m-file名稱fun function f=fun(X) r=1; f=zeros(3,1); f=[0.5*0.1-X(1)*0.1-1000*2.1*0.1*0.1*X(1)*X(2)-1000*0.00012*X(1); 0.5*0.1*r-X(2)*0.1-1000*2.1*0.1*0.1*X(1)*X(2)-1000*1.5*0.1*0.1*X(2)*X(3); -X(3)*0.1+1000*2.1*0.1*0.1*X(1)*X(2)-1000*1.5*0.1*0.1*X(2)*X(3)]; _______________在command window___________________ >> xzero=fsolve('fun',[0 0 0]) Optimization terminated: first-order optimality is less than options.TolFun. xzero = 0.0960 0.0143 0.0917 __________________________________________________ [問題]+[討論] m-file內的意思就是把fun的函數寫成一個A=F 左A為3*1矩陣 右F為三條方程式3*1的排列方式 即[0 [F1(Xi) 0 = F2(Xi) 0] F3(Xi)] 可是此為r=1的情況下得到的答案 如今要寫r=2,r=3,r=4.... 甚至r=linspace(1,10,100)所對應的不同答案 我想可能要寫 r=1 [0 [F1(Xi) 0 = F2(Xi) 0] F3(Xi)] 得其解 再 r=2 [0 [F1(Xi) 0 = F2(Xi) 0] F3(Xi)] 得其解 再 r=3 [0 [F1(Xi) 0 = F2(Xi) 0] F3(Xi)] 得其解 ... 可是問題是用fsolve得其解 只能一次控制一個所以不能"全解" 該怎嚜寫比較好咧@@" -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 119.14.154.20
文章代碼(AID): #19PcAVgp (MATLAB)
文章代碼(AID): #19PcAVgp (MATLAB)