| 查看: 561 | 回复: 1 | ||
[求助]
在matlab中用牛顿迭代法
|
|
方程如下：289.95*x^2+60*x-60*x*(exp((x-20)/32.3081)+1)^(-1)-1938.486*log(exp((x-20)/32.3081)+1)+1043.8243-J*y*10^6=0 将此方程用牛顿迭代法求出x关于y的表达式，其中J当做未知常量，这个方程比较复杂，求完之后对再求积分int（x，y，0，7.7*10^(-6))，最后的结果应该是一个关于J的表达式。 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
zjz311
木虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 2209.1
- 红花: 8
- 帖子: 139
- 在线: 80.5小时
- 虫号: 1555164
- 注册: 2011-12-28
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
xzhdty: 金币+2, 谢谢参与 2012-09-21 20:59:11
感谢参与，应助指数 +1
xzhdty: 金币+2, 谢谢参与 2012-09-21 20:59:11
|
为什么用牛顿迭代？ y=(289.95*x.^2+60*x-60*x.*(exp((x-20)/32.3081)+1).^(-1)-1938.486*log(exp((x-20)/32.3081)+1)+1043.8243)*10^(-6)./J; matlab代码： y=@(J) @(x) (289.95*x.^2+60*x-60*x.*(exp((x-20)/32.3081)+1).^(-1)-1938.486*log(exp((x-20)/32.3081)+1)+1043.8243)*10^(-6)./J; integ=@(J) quadl(y(J),0,7.7*10^(-6)); J=rand(1); integ(J); |
2楼2012-09-16 12:35:58
投票:
10