houshang 发表于 2006-3-5 11:01

非线性方程组的解法

请教有那位高手能帮忙解下下面这个方程:<BR>y(1)=0.1472*960-2.13e-7*((t(1)+t(2))/2)^3*(t(2)-t(1))-2.13e-7*((t(2)+t(3))/2)^3*(t(2)-t(3));<BR>y(2)=0.11*960-19.71*(t(1)-t(3))-2.13e-7*((t(1)+t(3))/2)^3*(t(1)-t(3))-2.13e-7*((t(1)+t(2))/2)^3*(t(1)-t(2));<BR>y(3)=0.022*860-2.13e-7*((t(2)+t(3))/2)^3*(t(3)-t(2))-2.04e-8*((t(4)+t(3))/2)^3*(t(3)-t(4))-1.38*(t(3)-t(4));<BR>y(4)=1.38*(t(3)-t(4))+2.04e-8*((t(4)+t(3))/2)^3*(t(3)-t(4))-1.005*6.83e5*(8.82*(t(4)-10)/t(4))^0.5*(10-12);<BR>

[ 本帖最后由 lxq 于 2006-12-26 00:03 编辑 ]

happy 发表于 2006-3-6 09:04

回复：（houshang）非线性方程组的解法

试一下fsolve

grta 发表于 2006-4-4 12:04

fsolve对初值的敏感性太强了吧

happy 发表于 2006-4-4 21:50

回复：（houshang）非线性方程组的解法

那你就自己找一种算法写程序吧，比如延拓法，呵呵<br>对初值要求很低，不过程序很复杂<br><br>fortran有现成的延拓法程序
[此贴子已经被作者于2006-4-4 21:51:07编辑过]

huqingwan 发表于 2006-12-25 15:41

我也想问同样的问题，如果用fsolve 或者 多维的ode45 都要求有好的初值， 对于延拓方法，如何构造延拓函数，有什么要求，
高手给介绍一些相关书籍，谢谢！

dingd 发表于 2006-12-25 21:21

1stOpt不需初值，下面是一组结果：

t1: -174.411218902856
t2: 483.892193989488
t3: -712.057566657294
t4: 10.0000000776654

dingd 发表于 2006-12-25 21:39

还有一组：

t1= -331.73893388449
t2= -398.067008191837
t3= -433.653484507166
t4= 10.0000001665828

huqingwan 发表于 2006-12-28 11:12

1stOpt不需初值，

请问这句话是什么意思啊？解释一下好么？
:handshake

dingd 发表于 2006-12-28 11:27

就是不需象Matlab那样要猜或赋初值，直接算就行。

wenzicao 发表于 2008-1-25 12:58

回复 9楼 的帖子

解方程目标函数不是零对吗?

无水1324 发表于 2008-1-25 15:51

回复 8楼 的帖子

1stOpt也是一个数学软件，它求解你的问题不需要初值
页: [1]
查看完整版本: 非线性方程组的解法