| 查看: 547 | 回复: 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的表达式。 |
回复此楼» 猜你喜欢
听说今天filecode变了
已经有34人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
基金中了
已经有6人回复
-
filecode
已经有14人回复
-
固定端突然变了，今天
已经有7人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有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
50