求高手用matlab解决这个二重积分，跪谢啦！

genixingfu

呼唤高手！在附件里

rocwoods

把相应参数值给出来，用dblquad函数,不要用符号积分。自己看看dblquad的函数帮助，用法写得很清楚。

[ 本帖最后由 rocwoods 于 2009-2-12 10:15 编辑 ]

genixingfu

在附件中，已经改了，麻烦了！

genixingfu

高手们好，我也在编了一下，但是不完善，应该是那几个参数都是已知的，就是说最后的程序是这样的，在键盘上随便输入那几个参数，比如，u1=2呀就是一组值，最后出现的结果应该是在0~1之间的数。我的程序如下，仅供抛砖引玉！
clear;
syms a1 a2 b1 b2 r fai R
a=1/4*(1/a2^2-1/a1^2);
b=1/4*(1/a2^2+1/a1^2);
d=1/a1^2*exp(-1/2*((b1/a1)^2+(b2/a2)^2));
jfsi=r*exp(-b*r^2+a*r^2*cos(2*fai)+r*(b1/a1^2*cos(fai)+b2/a2^2*sin(fai)));%积分式
jg=int(jfsi,fai,0,2*pi);%第一重积分，
P=d/(2*pi)*int(jg,r,0,R);%第二重积分，


结果：

P =

1/2/a1^2*exp(-1/2*b1^2/a1^2-1/2*b2^2/a2^2)/pi*int(int(r*exp((-1/4/a2^2-1/4/a1^2)*r^2+(1/4/a2^2-1/4/a1^2)*r^2*cos(2*fai)+r*(b1/a1^2*cos(fai)+b2/a2^2*sin(fai))),fai = 0 .. 2*pi),r = 0 .. R)