| 查看: 566 | 回复: 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一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有6人回复
-
两块石头
已经有5人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有8人回复
-
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生
已经有3人回复
-
找导师
已经有8人回复
-
国社科又开始会评了，不知道这次命运如何
已经有14人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
学科评审组评审是指会评吗？
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
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