我有一个简单的方程组,我试图用数字解决,
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的类似问题上做了这个,我没有问题得到结果,所以我很难解决问题所在。