matlab fsolve用法 在 コバにゃんチャンネル Youtube 的精選貼文
「matlab fsolve用法」的推薦目錄：
- 關於matlab fsolve用法 在 コバにゃんチャンネル Youtube 的精選貼文
- 關於matlab fsolve用法 在 大象中醫 Youtube 的最佳貼文
- 關於matlab fsolve用法 在 大象中醫 Youtube 的最佳解答
- 關於matlab fsolve用法 在 Re: [問題]解非線性方程式.....三個函數- 看板MATLAB 的評價
- 關於matlab fsolve用法 在 MATLAB教學- 11方程式求根 的評價
- 關於matlab fsolve用法 在 請教matlab 非線性聯立方程式的問題 - Mobile01 的評價
- 關於matlab fsolve用法 在 Re: [問題]解非線性方程式.....三個函數- 看板MATLAB - PTT數位 ... 的評價
matlab fsolve用法 在 請教matlab 非線性聯立方程式的問題 - Mobile01 的推薦與評價
請熟悉matlab的高手大大看一下那兒出錯指點嗎?謝謝。 ... 你的用法似乎比較接近solve, 我不確定在fsolve能不能這樣用. 或許先修正kc,kd, ... <看更多>
matlab fsolve用法 在 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
... <看更多>