非线性方程组的解法

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

请教有那位高手能帮忙解下下面这个方程: 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)); 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)); 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)); 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); 

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

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

试一下fsolve

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

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

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

那你就自己找一种算法写程序吧，比如延拓法，呵呵 对初值要求很低，不过程序很复杂[em01] 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

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

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

1stOpt也是一个数学软件，它求解你的问题不需要初值

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 非线性方程组的解法

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

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

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

回复 9楼的帖子

回复 8楼的帖子

[编程技巧] 非线性方程组的解法

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

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

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

回复 9楼 的帖子

回复 8楼 的帖子

回复 9楼的帖子

回复 8楼的帖子