请帮忙看下这个复杂的积分
本帖最后由 pper 于 2010-2-22 15:22 编辑图中的表达式是我要求的表达式，具体的程序如下：
syms x y
xx1=sym('xx1')
yy1=sym('yy1')
D=110262.614;l1=0.2;l2=0.3;kk=20000;...
vv=60;t=0.1;pa=3.14;wo=12.5;pz=75;cc=1000;mm=840;
www1=pz*exp(wo*i*t)/4/pa^2
www2=sin(x*l1)*sin(y*l2)/x/y/l1/l2
www3=exp((x*xx1+y*yy1)*i)
www31=exp(-(x*x1+y*y1)*i)
www32=exp(-(x*x1-y*y1)*i)
www33=exp((x*x1+y*y1)*i)
www34=exp((x*x1-y*y1)*i)
www4=www3*(www31+www32+www33+www34)
www5=D*(x^2+y^2)^2-mm*vv^2*y^2+kk-cc*(vv*y*i)
www6=www2/www5*www4
www6=char(www6)
www6=strrep(www6,'*','.*');
www6=strrep(www6,'/','./');
www6=strrep(www6,'\','.\');
www6=strrep(www6,'^','.^')
www7=dblquad(www6,-1000,1000,-1000,1000)
www8=www1.*www7
xx1=-100:0.1:100;yy1=-100:0.1:100
plot3(xx1,yy1,www8)
出现错误
??? Error using ==> inline.subsref
Not enough inputs to inline function.
Error in ==> dblquad>innerintegral at 78
fcl = intfcn(xmin, y(1), varargin{:}); %evaluate only to get the class below
Error in ==> quad at 63
y = f(x, varargin{:});
Error in ==> dblquad at 58
Q = quadf(@innerintegral, ymin, ymax, tol, trace, intfcn, ...
感觉也没有错在哪儿啊，请各位帮我看一下，
红色的字母是常数，蓝色的是变量
程序里面xx1=x'；yy1=y' ；x=kexi；y=yita
页:
[1]