[問題]解非線性方程式.....三個函數

看板MATLAB作者 (我要戒煙!!!)時間14年前 (2011/12/07 16:15), 編輯推噓0(0推 0噓 1→)
留言1則, 1人參與, 最新討論串1/4 (看更多)
之前有解過2個未知數的非線性方程式 y=solve('873=M0*(1-2*exp(-TI3/T1)+exp(-TR/T1)*exp(TA/T1))', '555=M0*(1-2*exp(-TI2/T1)+exp(-TR/T1)*exp(TA/T1))'); 大概需要將近十分鐘才能解出值 現在想增加一個函數來解方程式 如下 y=solve('873 = M0*(1-f*exp(-2600/T1)-exp(-2600/T1)+f*exp(-3000/T1)* exp(137.6/T1))', '555 = M0*(1-f*exp(-1600/T1)-exp(-1600/T1)+f*exp(-3000/T1)* exp(137.6/T1))', '91 = M0*(1-f*exp(-800/T1)-exp(-800/T1)+f*exp(-3000/T1)* exp(137.6/T1))' ); 以上為我三個方程式的聯立 大概解了一天了還解不出答案~~電腦一直呈現在BUSY的狀態 有人建議我用lsqcurvefit來解 爬了文也沒看到多少相關的文章 873,555,91.....總共有六組不同的值 -2600,-1600,-800 也有六組不同的值 請問我該如何去設定輸入跟輸出 因為help的example我看不太懂....所以想請問一下板上的高手 謝謝 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.117.157.80
→
12/07 18:56, , 1F
我想解f,M0,T1
12/07 18:56, 1F
‣ 更多 marco110232 的文章
文章代碼(AID): #1Eto2JdO (MATLAB)
‣ 更多 marco110232 的文章
文章代碼(AID): #1Eto2JdO (MATLAB)