马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
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 编辑 ] |