7楼: Originally posted by walk1997 at 2013-12-24 09:03:42
又试了下
不用分数代替也可以 不过是复数解(精度也只能是机器精度之内左边==0)
{t1 -> 1174.72 - 207.014 I, t2 -> 772.469 + 196.068 I}
也许还有其他复数解...

8楼: Originally posted by dingd at 2013-12-24 10:03:03
1stOpt试试，很简单，数值解的话比mathematica强的多：

Function
1*(1273 - t1) + 1.76*10^(-4)*5.67*10^(-8)*(1273^4 - (t1)^4) - ((1.5 + 0.6 + (
2*1.38*10^(-23)*t1)/(1.6*10^(-19)))*1.2*10^6*(t1)^2*
E ...

12楼: Originally posted by walk1997 at 2013-12-24 22:17:28
谢谢提醒 LZ的方程是有实数解
NSolve应该还是不够强大 对于无穷多解方程 在给定解参数范围下
没给出解 这点是不够的(对有些方程能给出给定范围的全部解)
不过FindRoot能给出上面这解，另外单独这个的t1 t2带入表 ...

8楼: Originally posted by dingd at 2013-12-24 10:03:03
1stOpt试试，很简单，数值解的话比mathematica强的多：

Function
1*(1273 - t1) + 1.76*10^(-4)*5.67*10^(-8)*(1273^4 - (t1)^4) - ((1.5 + 0.6 + (
2*1.38*10^(-23)*t1)/(1.6*10^(-19)))*1.2*10^6*(t1)^2*
E ...

13楼: Originally posted by chyanog at 2013-12-24 23:28:09
超越方程的全局求解本身就不容易，Matlab里的fsolve、1stOpt也都不是全局求解的...

14楼: Originally posted by Simewe at 2013-12-24 23:34:04
1stOpt算数值解强在哪里，是求解速度还是精度？...