找回密码
 注册
Simdroid-非首页
查看: 60|回复: 0

请帮忙看下这个复杂的积分

[复制链接]
发表于 2010-2-22 15:19:37 | 显示全部楼层 |阅读模式 来自 北京海淀
本帖最后由 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

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-10-7 01:34 , Processed in 0.035604 second(s), 14 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表