求解T的值, 这个算法该如何解?

[zta8507]

表达式中, d,c1,c2,c3,c4,c5,c6,c7都是已知数,请教各位大虾,该如何计算T的值?
其中, ln(T) 可通过数学函数 log(T) 计算.

[枫叶饭团]

数学家在哪里。。。
我数学不好= =:(

[meng0987]

牛顿切线法？

[tangboyun]

要算法还是软件还是库？
算法可以参考http://en.wikipedia.org/wiki/List_of_algorithms 中2.5.8节的Root finding 一节相关介绍。（要实际应用不要自己去写这类算法。）
库的话可以读下 GNU Scientific Library (or GSL )的manual，至于软件，几乎所有数学软件都有相关部分，搜下manual

[cuihao]

……用万能的牛顿法。
对于这道具体的题目，导函数算是已知的，好写多了。

[zta8507]

谢谢大家的帮助,这个公式已经通过牛顿迭代法解决了,下面是C语言代码片段:

double f(double c1, double c2, double c3, double c4, double c5, double c6, double c7, double d, double T)
{
return (c1 + c2*T + c3*pow(T,2) + c4*pow(T,3) + c5*pow(T,4) + c6*pow(T,5) + c7*log(T)*T - d*T);
}

double f1(double c2, double c3, double c4, double c5, double c6, double c7, double d, double T)
{
return (c2 + 2*c3*T + 3*c4*pow(T,2) + 4*c5*pow(T,3) + 5*c6*pow(T,4) + c7*log(T) - d);
}

double newton_ddf()
{
double T;
doubel eps;

T = (double)473.15; // 初值设定很重要
eps = 1E-5;
//初始化c1,c2,c3,c4,c5,c6,c7,d

// 原始表达式 : d= c1/T + c2 + c3*T + c4*T^2 + c5*T^3 + c6*T^4 + c7*ln(T)
// 转换表达式1 : 0 = c1/T + c2 + c3*T + c4*T^2 + c5*T^3 + c6*T^4 + c7*ln(T) - d
// 转换表达式2 : 0 = c1 + c2*T + c3*T^2 + c4*T^3 + c5*T^4 + c6*T^5 + c7*ln(T)*T - d*T
do
{
T1 = T;
T = T1 - f(c1, c2, c3, c4, c5, c6, c7, d, T1) / f1(c2, c3, c4, c5, c6, c7, d, T1);
} while (fabs(T-T1) >= eps);


printf("T: %lf\n", T);

[zhugejunqing]

感觉是数值计算方法里的题目。。。一般用matlab编程，自动求解。。。

[b33e]

嗯，是我就用matlab

[月下叹逍遥]

Fortran和mathematica...