歡迎對人類大腦與意識有興趣的同學，加入腦與意識實驗室！
【徵才】國立臺灣大學心理學系謝伯讓老師誠徵學士級或碩士級專任研究助理一名
【徵求條件】
1. 對意識及大腦神經科學研究有高度興趣。
2. 具備 Python 及 Matlab 等電腦程式語言編撰經驗。
【工作內容】設計實驗及數據分析相關程式撰寫。
【工作時間】2020年9月起。
【工作地點】台大心理學系南館。
【工作待遇】週休二日，勞保、健保。薪資依台灣大學學士級/碩士級專任研究助理相關規定支給，相關年資可累計。年終獎金 1.5 個月，依當年度到職月份比例折算。
【應徵方式】請將成績單及個人履歷自傳（說明學經歷、研究興趣、相關軟體使用經驗、及未來生涯規劃）以 e-mail 寄至謝伯讓老師：[email protected]。
【實驗室網頁】https://sites.google.com/site/brainconsciousnesslab/
1.0 e 03 matlab 在 Re: [問題]解非線性方程式.....三個函數- 看板MATLAB 的推薦與評價
此題可以用fsolve來解喔~~
先設函數:
function y=prac1(x)
y=zeros(1,3);
y(1) =
x(1)*(1-x(2)*exp(-2600/x(3))-exp(-2600/x(3))...
+x(2)*exp(-3000/x(3))*exp(137.6/x(3)))-873;
y(2) =
x(1)*(1-x(2)*exp(-1600/x(3))-exp(-1600/x(3))...
+x(2)*exp(-3000/x(3))*exp(137.6/x(3)))-555;
y(3) =
x(1)*(1-x(2)*exp(-800/x(3))-exp(-800/x(3))...
+x(2)*exp(-3000/x(3))*exp(137.6/x(3)))-91;
接著來解~
clear
clc
x01=[1,1,1];
a1=fsolve('prac1',x01)
y1=prac1(a1)
x02=[100,100,100];
a2=fsolve('prac1',x02)
y2=prac1(a2)
x03=[10000,10000,10000];
a3=fsolve('prac1',x03)
y3=prac1(a3)
% fsolve的用法是 fsolve('待解函數',初始值)
(*初始值可任意代 但最好盡量靠近答案(mo,f,T1) 就是要你先猜大概的(mo,f,T1)
若與真實答案的級數差太多
則可能找不到你要的解)
x01,x02,x03就是我設的三種初始值
而a1,a2,a3就是各解出的(mo,f,T1)值
y1,y2,y3為各代回function 看這三個方程式的值是否都等於0
來判斷解出的答案是否是我們要的 若與0差距太多 則可能是初始值設的離實際值太遠
讓我們來看看結果如何吧~~
a1 =
506.3333 1.0000 1.0000
y1 =
-366.6667 -48.6667 415.3333
______________________________________
a2 =
1.0e+003 *
1.0833 0.0008 1.4808
y2 =
1.0e-010 *
-0.1489 0.4025 0.6632
_____________________________________
a3 =
1.0e+003 *
1.0833 0.0008 1.4808
y3 =
1.0e-010 *
0.5537 0.5082 -0.5024
______________________________________
可看出a2和a3一樣 而y2和y3也確實趨近於0
至於a1的值和a2,a3差蠻多的 但也有可能是另一組解
不過y1就不趨近於0了 可見a1這組不是我們要的答案!
所以 此題之 mo=1083.3 f=0.8 T1=1480.8
呼~~~ 解完囉~~~~
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 118.169.14.166
... <看更多>