作者wahaluna (绿狒)
看板MATLAB
标题[问题] 有关多重积分
时间Mon Nov 26 20:46:00 2012
大家好,我是matlab新手,目前遇到一些问题,希望大家可以帮助我
我要算的式子是三层积分,要求的未知数c在最外层的积分下界,
最里面是两个积分相乘,
g(w,y1)= S f(v1,y1) dv1 * S f(v2,y1,w) dv2 范围都是0~inf,
之後把g对y1积分 wp(w)=S g(w,y1) dy1 范围-inf~inf
最外层是 S wp(w) dw 范围c~inf,且全部积分完要等於alpha
以下是我的程式码
g=int(exp(((-9*n1*(v1*y1-p1)^2)/2)+(((-n1+1)/2)*v1^2))*v1^(n1-1),v1,0,inf)*
int(exp(((-9*n2*(v2*(w+y1)-p2)^2)/2)+(((-n2+1)/2)*v2^2))*v2^(n2-1),v2,0,inf)
wp=A*int(g,y1,-inf,inf)
eq=int(wp,w,c,inf)-alpha
其中n1=n2为已知,p1=p2也是已知,A为一个常数,alpha为已知
我不知道哪边使用错误
g跑出来的结果跟我想的不一样,当中包含了v1和v2
类似这样
limit(-1/19683*(-327115800000*y1^7*pi^(1/2)*erf(30*y1/(20*y1^2+2)^(1/.....
wp跑出来也是一样很复杂
2489620740387015/4398046511104*int(limit(-1/19683*(-327115800000*y1^7*pi^(1/..
然後就坏掉跑不出来了...
我的问题是,对某个变数积分後,该变数不是应该要消失吗?
为什麽还是会一直出现很多变数
ex.g(w,y1)还是有v1,v2,wp(w)还是有y1...
还有像这样变数比较多的情况下我积分该怎麽写?
要求的c值,是不是不能用fzero来解?
我知道我的问题很多很复杂
不过还是希望大家可以给我一点意见跟想法,谢谢~
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.108.142
1F:→ vaakaa:同学你是在解什麽色散的问题吗?@@ 11/27 01:07
2F:→ vaakaa:另外 解根只靠fzero是很容易失败的 11/27 01:11
3F:→ wahaluna:不是欸...那对积分的部分有什麽建议吗 11/27 08:51
4F:→ vaakaa:一定只能数值算?积分可不可以进一步简化? 11/27 17:45