tongjiayan 发表于 2005-10-29 18:20

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

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

suffer 发表于 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>

tongjiayan 发表于 2005-10-30 10:49

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

您好!<BR>我的意思是由于被积函数的复杂性,MATLAB找不到解析解。

brianlei 发表于 2005-10-31 22:49

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

我用过这样的方法,不知道能否解决你的问题.<BR>比如需要求f=a*exp(x^2)+b*x*exp(x^2)/c关于x的在上的积分.<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

suffer 发表于 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>

tongjiayan 发表于 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的在上的积分.<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编辑过]

lzxin1234 发表于 2009-5-25 15:36

置顶,这方法特别好!谢谢

weicheng519 发表于 2011-4-1 15:42

楼主我遇到同样的问题 请问你后来是怎么解决的呢?

qq1522491397 发表于 2012-5-2 08:54

好好好!!!
页: [1]
查看完整版本: 请教:关于MATLAB中带参数的数值积分