| 查看: 501 | 回复: 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的表达式。 |
回复此楼» 猜你喜欢
“人文社科而论，许多学术研究还没有达到民国时期的水平”
已经有4人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有11人回复
-
什么是人一生最重要的？
已经有4人回复
-
版面费该交吗
已经有3人回复
-
今年春晚有几个节目很不错，点赞！
已经有12人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有12人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有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