用matlab求解方程

来源：百度知道 编辑：UC知道 时间：2024/06/15 13:00:57
解如下方程
1.01*t_5+(1.84*t_5+2500)*0.622*602.4*10^(7.45*t_5/(235+t_5))/(518400-602.4*10^(7.45*t_5/(235+t_5))=23.7
一看就想说用solve的大哥们好好看一下再说
谢谢

这个方程用sovle是求不出解析解的，求它的数值解：

>> t_5=10:0.1:20;
y=1.01*t_5+(1.84*t_5+2500)*0.622*602.4.*10.^(7.45*t_5./(235+t_5))./(518400-602.4.*10.^(7.45*t_5./(235+t_5)))-23.7;
plot(t_5,y,t_5,0,'r')

看出解大概在17到18之间。

>> f=inline('1.01*t_5+(1.84*t_5+2500)*0.622*602.4*10^(7.45*t_5/(235+t_5))/(518400-602.4.*10^(7.45*t_5/(235+t_5)))-23.7','t_5')

f =

Inline function:
f(t_5) = 1.01*t_5+(1.84*t_5+2500)*0.622*602.4*10^(7.45*t_5/(235+t_5))/(518400-602.4.*10^(7.45*t_5/(235+t_5)))-23.7

>> t_5=fsolve(f,17)
Optimization terminated: first-order optimality is less than options.TolFun.

t_5 =

17.494553672455257

解得的结果如上，不过和楼上用的图像法得到的结果很近似，哈哈，不过我的作法是比较正规的。

其实遇到非线性度很高的方程有一个很简单的方法，或许很多人都不屑于这么做，但是确实很有效，你用matlab做个图找它和x轴的交点，就算是用牛顿法数值解，其原理也大致相同。
解出来t_5=17.49455364409869