FindRoot,起始值

时间:2015-05-09 23:27:23

标签: numerical-methods

我有一个简单的方程组,我试图用数字解决,

    ClearAll["Global`*"]
     b1 = 1
     T = 2
     t = .5

FindRoot[{(1 - xbar1) 2 ra1 == n1 - (1 - n1)/T,
ra1 == (pa1^2) (b1^2)/pw,
pw == (1 - xbar1) ra1 - 1,
xbar1 == 1/t  Log[(t n1 + ra1)/ra1],
t[1 - xbar1] + Log[T] - Log[ra1] + Log[t (1 - n1)/(Exp[t] - 1)] == 
 0},
{{n1, .5}, {ra1, .5}, {pa1, .5}, {pw, 1.5}, {xbar1, .5}}, 
AccuracyGoal -> 11]

这是我得到的输出,

FindRoot::nlnum: "The function value {0.25,0.333333,2.25,-0.31093,0.432752 
+0.5[0.5]} is not a list of numbers with dimensions {5} at 
{n1,ra1,pa1,pw,xbar1} = {0.5,0.5,0.5,1.5,0.5}."de here`{n1,ra1,pa1,pw,xbar1} =                    `   {0.5,0.5,0.5,1.5,0.5}."``

我想知道我是不是正确指定了初始起始值。我在使用NSolve的类似问题上做了这个,我没有问题得到结果,所以我很难解决问题所在。

0 个答案:

没有答案