Digital 发表于 2006-7-5 12:15

DUNLSF的收敛问题

用DUNLSF求解一个2维非线性方程组的解。参数数组Iparam和Rparam的初始值如下:

IPARAM(1) = 1
IPARAM(2) = 15
IPARAM(3) = 1E9
IPARAM(4) = 1E9
IPARAM(5) = 1E9
IPARAM(6) = 1

RPARAM(1) = 0.1D-100
RPARAM(2) = 0.1D-100
RPARAM(3) = 0.1D-100
RPARAM(4) = 0.1D-100
RPARAM(5) = 0.1D-100
RPARAM(6) = 0.1D+100
RPARAM(7) = 0.1D+100

但程序停止计算时,并没有满足收敛准则,比如结束时,残差向量FVEC的值为:
FVEC(1)= -0.57126370E+06 FVEC(2)= -0.22382156E+06
大大超过收敛准则。
数组 IPARAM 和 RPARAM 的输出结果为:

IPARAM(1) =1
IPARAM(1)=15
IPARAM(1)=1
IPARAM(1)=345
IPARAM(1)=2
IPARAM(1)=1

RPARAM(1)= 0.10000000-100
RPARAM(2) =0.10000000-100
RPARAM(3) =0.10000000-100
RPARAM(4) =0.10000000-100
RPARAM(5) =0.10000000-100
RPARAM(6) =0.10000000+100
RPARAM(7) =0.10000000+100


有没有什么办法,强制程序停止计算时,残差很小,以至于可以忽略。
thanks.
页: [1]
查看完整版本: DUNLSF的收敛问题