声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 5035|回复: 8

请教:关于MATLAB中带参数的数值积分

[复制链接]
发表于 2005-10-29 18:20 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
MATLAB中有数值和解析的积分,现在有一个找不到解析解的被积函数,其中包含多个变量。将其中的一个变量设为积分变量,其余的设为参数。现在只能退而求其次,要求对参数求解析解并且对积分变量求数值解,即:积出的函数是保留完整参数的。<BR>请教诸位前辈,这样的问题该如何处理?<BR>不胜感激!

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2005-10-30 09:53 | 显示全部楼层

回复:(tongjiayan)请教:关于MATLAB中带参数的数值...

<DIV class=quote><B>以下是引用<I>tongjiayan</I>在2005-10-29 18:20:33的发言:</B><BR>MATLAB中有数值和解析的积分,现在有一个找不到解析解的被积函数,其中包含多个变量。将其中的一个变量设为积分变量,其余的设为参数。现在只能退而求其次,要求对参数求解析解并且对积分变量求数值解,即:积出的函数是保留完整参数的。<BR>请教诸位前辈,这样的问题该如何处理?<BR>不胜感激!</DIV>
<P>找不到解析解的被积函数,这是什么意思?</P>
 楼主| 发表于 2005-10-30 10:49 | 显示全部楼层

回复:(suffer)回复:(tongjiayan)请教:关于MATL...

您好!<BR>我的意思是由于被积函数的复杂性,MATLAB找不到解析解。
发表于 2005-10-31 22:49 | 显示全部楼层

回复:(tongjiayan)请教:关于MATLAB中带参数的数值...

我用过这样的方法,不知道能否解决你的问题.<BR>比如需要求f=a*exp(x^2)+b*x*exp(x^2)/c关于x的在[0 1]上的积分.<BR>我这样做的----------&gt;&gt;syms a b c x;<BR>&gt;&gt;f=a*exp(x^2)+b*x*exp(x^2)/c;<BR>&gt;&gt;q=int((f),x,0,1)<BR>q=-1/2*(i*a*pi^(1/2)*erf(i)*c-b*exp(1)+b)/c<BR>&gt;&gt;Q=vpa(q)<BR>Q=-.500000*(-2.925303*a*c-1.718281*b)/c
发表于 2005-11-1 09:18 | 显示全部楼层

回复:(tongjiayan)回复:(suffer)回复:(tongji...

<DIV class=quote><B>以下是引用<I>tongjiayan</I>在2005-10-30 10:49:14的发言:</B><BR>您好!<BR>我的意思是由于被积函数的复杂性,MATLAB找不到解析解。</DIV>
<P>找出不可及的自变量,将其作为常量处理<BR><BR>剩下的就是常规符号积分的操作了</P>
 楼主| 发表于 2005-11-4 19:00 | 显示全部楼层

回复:(brianlei)回复:(tongjiayan)请教:关于MA...

<DIV class=quote><B>以下是引用<I>brianlei</I>在2005-10-31 22:49:46的发言:</B><br>我用过这样的方法,不知道能否解决你的问题.<br>比如需要求f=a*exp(x^2)+b*x*exp(x^2)/c关于x的在[0 1]上的积分.<br>我这样做的----------&gt;&gt;syms a b c x;<br>&gt;&gt;f=a*exp(x^2)+b*x*exp(x^2)/c;<br>&gt;&gt;q=int((f),x,0,1)<br>q=-1/2*(i*a*pi^(1/2)*erf(i)*c-b*exp(1)+b)/c<br>&gt;&gt;Q=vpa(q)<br>Q=-.500000*(-2.925303*a*c-1.718281*b)/c</DIV>您好!<br>非常感谢您。<br>
[此贴子已经被作者于2005-11-10 12:46:43编辑过]

发表于 2009-5-25 15:36 | 显示全部楼层
置顶,这方法特别好!谢谢
发表于 2011-4-1 15:42 | 显示全部楼层
楼主我遇到同样的问题 请问你后来是怎么解决的呢?
发表于 2012-5-2 08:54 | 显示全部楼层
好好好!!!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-5-19 06:40 , Processed in 0.053635 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表