关于fsolve求解非线性方程的问题

caiwhph

本人采用以下程序运算非线性方程组：
定义函数：
function F=myfun(x)
r0=4.8839
l1=9.5006
v1=0.0436
l2=9.5167
v2=0.0785
F(1)=x(1)*(tan(v1)*(x(1)^2+l1^2-r0^2)+sqrt(4*x(1)^2*l1^2-(x(1)^2+l1^2-r0^2)^2))-x(3)*(x(1)^2+l1^2-r0^2-sqrt(4*x(1)^2*l1^2-(x(1)^2+l1^2-r0^2)^2)*tan(v1));

F(2)=x(2)*(tan(v2)*(x(2)^2+l2^2-r0^2)+sqrt(4*x(2)^2*l2^2-(x(2)^2+l2^2-r0^2)^2))-x(3)*(x(2)^2+l2^2-r0^2-sqrt(4*x(2)^2*l2^2-(x(2)^2+l2^2-r0^2)^2)*tan(v2));

F(3)=x(2)-x(1)-0.381;
调用函数运算：
clc
clear all
close all
options=optimset('Display','on');
[r1,val1]=fsolve('myfun',[5,5,0.5])
运行结果：
r1 =
   4.6157 - 0.0008i   4.6306 + 0.0000i   0.2168 - 0.0456i

val1 =
   1.0282 - 2.3754i  12.9889 - 5.0570i  -0.3662 + 0.0008i
求出来为虚数，而且不满足方程组F(3).
求教问题出在哪里。

shamohu

用1stOpt得到下面两组，不知对否？
第一组：
x1        13.0745742771892
x2        13.455574277189
x3        4.67463223720975

第二组：
x1        -7.52671059593533
x2        -7.14571059593532
x3        -4.91214349097977

caiwhph

2# shamohu
x1和x2取值应该是在4附近

messenger

1# caiwhph

命令没有问题，只是x1和x2在4附近找不到解，如果你的方程有实际意义要求x1和x2在4附近取值，检查一下你的方程是否正确。