| 查看: 1041 | 回复: 17 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】请高手帮忙解下方程，命令过程请留下，谢谢！
|
|||
|
2*exp(-x)*(x^3+18*x^2+88*x+96)/(x*(x^4+20*x^3+120*x^2+240*x+120))=exp(-y)*(y^2+18*y^2+88*y+96)/(y*(y^4+20*y^3+120*y^2+240*y+120)) （1） 447.33*log(x)=464.56*log(y) （2） 怎么求解啊，哪种方法最简便，占用计算机资源最少，谢谢了！ [ Last edited by spring180 on 2009-8-2 at 10:26 ] |
回复此楼» 猜你喜欢
梦想
已经有4人回复
-
中青基了要发朋友圈吗？
已经有6人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上意见出来了
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
找导师
已经有7人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
spring180(金币+1):谢谢参与
sunxiao(金币+2,VIP+0):谢谢参与，欢迎常来 8-3 00:47
spring180(金币+5,VIP+0):谢谢！ 8-3 09:13
hitzhang(金币+2,VIP+0): 8-3 10:03
spring180(金币+1,VIP+0):不过我感觉你非线性拟合和fsolve一起是不是重复了？有什么意义吗？ 8-3 11:13
spring180(金币+1):谢谢参与
sunxiao(金币+2,VIP+0):谢谢参与，欢迎常来 8-3 00:47
spring180(金币+5,VIP+0):谢谢！ 8-3 09:13
hitzhang(金币+2,VIP+0): 8-3 10:03
spring180(金币+1,VIP+0):不过我感觉你非线性拟合和fsolve一起是不是重复了？有什么意义吗？ 8-3 11:13
|
function aatest clear all clc format long x0=30*rand(2,1); options=optimset('TolX',1e-16,'TolFun',1e-12,'LargeScale','on','MaxFunEvals',100000,'MaxIter',100000); [X,FVAL,EXITFLAG,OUTPUT]=fsolve(@myfun,x0,options) [X,RESNORM,RESIDUAL,EXITFLAG] = lsqnonlin(@myfun,x0,[],[],options) function F = myfun(x) f1=2*exp(-x(1))*(x(1)^3+18*x(1)^2+88*x(1)+96)/(x(1)*(x(1)^4+20*x(1)^3+120*x(1)^2+240*x(1)+120))-... exp(-x(2))*(x(2)^2+18*x(2)^2+88*x(2)+96)/(x(2)*(x(2)^4+20*x(2)^3+120*x(2)^2+240*x(2)+120)); f2=447.33*log(x(1))*100-464.56*log(x(2))*100; F=[f1;f2]; 查看EXITFLAG =1时的结果，但是即便是如此也与初值有关。 |
13楼2009-08-02 19:00:19
2楼2009-08-02 09:31:16
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
3楼2009-08-02 09:48:56
4楼2009-08-02 09:57:56
投票:
30