| 查看: 488 | 回复: 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的表达式。 |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有8人回复
-
26申博
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有22人回复
-
最失望的一年
已经有4人回复
-
国自然申请面上模板最新2026版出了吗？
已经有19人回复
-
请教限项目规定
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
推荐一本书
已经有13人回复
-
疑惑？
已经有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