此题可以用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