楼主: ivylee_777
4552 8

[问答] 如何用matlab计算下面积分的上限 [推广有奖]

  • 23关注
  • 4粉丝

硕士生

99%

还不是VIP/贵宾

-

威望
0
论坛币
1631 个
通用积分
1.0600
学术水平
0 点
热心指数
0 点
信用等级
5 点
经验
601 点
帖子
168
精华
0
在线时间
239 小时
注册时间
2005-7-6
最后登录
2017-3-17

楼主
ivylee_777 发表于 2013-8-21 02:03:33 |AI写论文
5论坛币
大家好,我想计算 这个积分的上限a. 请问用matlab应该怎么求?十分感谢!

最佳答案

TimeT 查看完整内容

我想到的笨方法是利用TAILOR展开来求积分。 1. 令x+0.5=y,这样exp(x^2+x)=exp(y^2-0.25)=exp(y^2)/exp(0.25)。 所以,只要求得出exp(y^2)的积分,你的问题就能解决。 2. 利用TAILOR展开 exp(y^2)=1+y^2+y^4/2!+y^6/3!+y^8/4!+... 3. 对这个展开后的多项式求积分=y + 1/3*y^3+1/5*y^5/2!+ 1/7*y^7/3!+1/9*y^9/4!+... 只要展开项足够多,就能求出相当精确值。 4. 再用上限(a,注意变换y=x+0.5, 所以应用a+0.5代入)代入就能把你 ...
关键词:MATLAB atlab matla Mat Lab matlab 如何

本帖被以下文库推荐

沙发
TimeT 发表于 2013-8-21 02:03:34
ivylee_777 发表于 2013-8-27 05:44
不一定matlab,能解出来就行,能教我用excel怎么解吗?谢谢!
我想到的笨方法是利用TAILOR展开来求积分。
1. 令x+0.5=y,这样exp(x^2+x)=exp(y^2-0.25)=exp(y^2)/exp(0.25)。 所以,只要求得出exp(y^2)的积分,你的问题就能解决。
2. 利用TAILOR展开 exp(y^2)=1+y^2+y^4/2!+y^6/3!+y^8/4!+...
3. 对这个展开后的多项式求积分=y + 1/3*y^3+1/5*y^5/2!+ 1/7*y^7/3!+1/9*y^9/4!+...
只要展开项足够多,就能求出相当精确值。
4. 再用上限(a,注意变换y=x+0.5, 所以应用a+0.5代入)代入就能把你的方程变为:含a的多项式=b的形式(即一元高次方程)。
5. 把这个方程输入到EXCEL,EXCEL中有SOLVER(或GOAL SEEK)功能来解一元方程,迅速由b就能解出a来。
具体的就请你自己实现吧。
已有 1 人评分学术水平 热心指数 收起 理由
ivylee_777 + 1 + 1 分析的有道理

总评分: 学术水平 + 1  热心指数 + 1   查看全部评分

藤椅
ivylee_777 发表于 2013-8-21 02:05:08
不知道为什么公式显示不出加号。积分的公式应该是 exp(x+x^2)。谢谢!

板凳
ivylee_777 发表于 2013-8-21 23:58:27
还没有解决,希望能够在这里得到解答,十分感谢!

报纸
TimeT 发表于 2013-8-24 18:26:46
一定要用MATLAB?我倒有个用EXCEL求的很笨很繁的方法(当B不是很大的情况下)。

地板
ivylee_777 发表于 2013-8-27 05:44:12
不一定matlab,能解出来就行,能教我用excel怎么解吗?谢谢!

7
ivylee_777 发表于 2013-8-28 01:02:43
十分感谢你的回答,这是个很活的解法啊,一点都不笨。还有个问题,第二步中,你的Taylor展开我看得出是围绕0展开的,你是先求exp(y)的展开式,然后将y^2带到y的位置,是吧?因为我直接求exp(y^2)的展开,没有得到你这个结果。谢谢你!

8
TimeT 发表于 2013-8-29 22:55:35
ivylee_777 发表于 2013-8-28 01:02
十分感谢你的回答,这是个很活的解法啊,一点都不笨。还有个问题,第二步中,你的Taylor展开我看得出是围绕 ...
我是先求exp(y)的展开式,然后将y^2带到y的位置,因为似乎这样方便些。:-)

9
ivylee_777 发表于 2013-8-30 04:28:11
TimeT 发表于 2013-8-29 22:55
我是先求exp(y)的展开式,然后将y^2带到y的位置,因为似乎这样方便些。:-)
十分感谢!我按照这个方法解解试试。

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-26 21:30