菜鸟求助：双重积分报错

WindLee

ay=0.77,az=0.1,uc=0.49,r=0.25,w=5;
bj=@(ky,kz)4*ay*az/[(ay^2+uc^2*ky.^2/w^2)*(az^2+uc^2*kz.^2/(w-1)^2)]*[2*besselj(1,(ky.^2+kz.^2)^0.5*r)/((ky.^2+kz.^2)^0.5*r)]^2

s=dblquad(bj,1,2,1,2);
报错：
??? Error using ==> mrdivide
Matrix dimensions must agree.
Error in ==> quad at 62
y = f(x, varargin{:});
Error in ==> dblquad>innerintegral at 88
    Q(i) = feval(quadf, intfcn, xmin, xmax, tol, trace, y(i), varargin{:});
Error in ==> quad at 62
y = f(x, varargin{:});
Error in ==> dblquad at 64
Q = feval(quadf, @innerintegral, ymin, ymax, tol, trace, intfcn, ...
刚刚用Matlab，不知道什么原因，希望高手指导！感激不尽！

如果上面的w是参数，我想对ky，kz分别在负无穷大正无穷积分，怎么搞？w一般怎么处理？

[ 本帖最后由 WindLee 于 2009-8-21 15:48 编辑 ]

WindLee

郁闷！自己解决了。又有一个问题。如果上面的w是参数，我想对ky，kz分别在负无穷大正无穷积分，最后得到关于w的函数。怎么搞？w一般怎么处理？

[ 本帖最后由 WindLee 于 2009-8-21 15:49 编辑 ]