1343阅读
0回复

[求助]关于syms定义的问题（二重积分）

购货方

UID ：57635

注册：2010-04-20
登录：2014-05-25
发帖：28
等级：仿真新人

0楼发表于: 2013-12-30 15:58:01

此悬赏帖已过期
最佳答案：3 rf币，热心助人剩余点数: 1 rf币。

问题描述：二重积分算四个量： Yx1x1，Yx1y1，Yy1x1，Yy1y1。如果只计算Yx1x1，程序则没错，但是如果要同时算四个量则报错。由于是新手，不知道如何解决，还请各位知道的告知指导下。万分感激。报错如下（后面附有涉及到的程序）：

Error using symengine>@()0.0
Too many input arguments.

Error in @(t,v)fx1y1_2(t.*sin(v),t.*cos(v)).*t

Error in integral2Calc>integral2t/tensor (line 229)
        Z = FUN(X,Y);  NFE = NFE + 1;

Error in integral2Calc>integral2t (line 56)
[Qsub,esub] = tensor(thetaL,thetaR,phiB,phiT);

Error in integral2Calc (line 10)
    [q,errbnd] = integral2t(fun,xmin,xmax,ymin,ymax,optionstruct);

Error in integral2 (line 107)
    Q = integral2Calc(fun,xmin,xmax,yminfun,ymaxfun,opstruct);

Error in M_MOM_SE_OB_test (line 209)
            fx1y1_22 =   integral2(fx1y1_2_polar,0,k0,0,2*pi,'AbsTol',1e-3) ;
-------------------------------------以下是涉及到的部分程序----------------------------------------------------------
syms kx ky
            f_p = L(1)*W(1)/4/1j*(exp(1j*p*pi/2)*sin(p*pi/2-kx*L(1)/2)/(p*pi/2-kx*L(1)/2)-exp(-1j*p*pi/2)*sin(p*pi/2+kx*L(1)/2)/(p*pi/2+kx*L(1)/2))*(exp(1j*q*pi/2)*sin(q*pi/2-ky*W(1)/2)/(q*pi/2-ky*W(1)/2)+exp(-1j*q*pi/2)*sin(q*pi/2+ky*W(1)/2)/(q*pi/2+ky*W(1)/2));
            f_p1 = L(1)*W(1)/4/1j*(exp(1j*p1*pi/2)*sin(p1*pi/2+kx*L(1)/2)/(p1*pi/2+kx*L(1)/2)-exp(-1j*p1*pi/2)*sin(p1*pi/2-kx*L(1)/2)/(p1*pi/2-kx*L(1)/2))*(exp(1j*q1*pi/2)*sin(q1*pi/2+ky*W(1)/2)/(q1*pi/2+ky*W(1)/2)+exp(-1j*q1*pi/2)*sin(q1*pi/2-ky*W(1)/2)/(q1*pi/2-ky*W(1)/2));

            f_f = L(1)*W(1)/4/1j*(exp(1j*p*pi/2)*sin(p*pi/2-kx*L(1)/2)/(p*pi/2-kx*L(1)/2)+exp(-1j*p*pi/2)*sin(p*pi/2+kx*L(1)/2)/(p*pi/2+kx*L(1)/2))*(exp(1j*q*pi/2)*sin(q*pi/2-ky*W(1)/2)/(q*pi/2-ky*W(1)/2)-exp(-1j*q*pi/2)*sin(q*pi/2+ky*W(1)/2)/(q*pi/2+ky*W(1)/2));
            f_f1 = L(1)*W(1)/4/1j*(exp(1j*p1*pi/2)*sin(p1*pi/2+kx*L(1)/2)/(p1*pi/2+kx*L(1)/2)+exp(-1j*p1*pi/2)*sin(p1*pi/2-kx*L(1)/2)/(p1*pi/2-kx*L(1)/2))*(exp(1j*q1*pi/2)*sin(q1*pi/2+ky*W(1)/2)/(q1*pi/2+ky*W(1)/2)-exp(-1j*q1*pi/2)*sin(q1*pi/2-ky*W(1)/2)/(q1*pi/2-ky*W(1)/2));
  %-------------------------------------Yx1x1------------------------------------------------- %
            fx1x1_1 =   f_p*f_p1*(k0^2-kx^2)/sqrt(k0^2-kx^2-ky^2);
            fx1x1_2 =   matlabFunction(fx1x1_1) ;
            fx1x1_2_polar =  @(t,v) fx1x1_2(t.*sin(v),t.*cos(v)).*t ;
            fx1x1_22 =   integral2(fx1x1_2_polar,0,k0,0,2*pi,'AbsTol',1e-3) ;

            fx1x1_3 =   f_p*f_p1*(k0^2-kx^2)/sqrt(kx^2+ky^2-k0^2)/(-1j);
            fx1x1_4 =   matlabFunction(fx1x1_3) ;
            fx1x1_4_polar = @(t,v) fx1x1_4(sin(v)./t,cos(v)./t)./t.^3 ;
            fx1x1_44 =   integral2(fx1x1_4_polar,0,1/k0,0,2*pi,'AbsTol',1e-3);

            Yx1x1_22 =  Yx1x1_22 -omega*epsilon0/(4*pi^2*k0^2)*(fx1x1_22+fx1x1_44) ;%the another part of (65)

   %--------------------------------------Yx1y1------------------------------------------------ %
            fx1y1_1 =   f_p1*f_f*(-kx*ky)/sqrt(k0^2-kx^2-ky^2);
            fx1y1_2 =   matlabFunction(fx1y1_1) ;
            fx1y1_2_polar =  @(t,v) fx1y1_2(t.*sin(v),t.*cos(v)).*t ;
            fx1y1_22 =   integral2(fx1y1_2_polar,0,k0,0,2*pi,'AbsTol',1e-3) ;%（line 209）%

            fx1y1_3 =   f_p1*f_f*(-kx*ky)/sqrt(kx^2+ky^2-k0^2)/(-1j);
            fx1y1_4 =   matlabFunction(fx1y1_3) ;
            fx1y1_4_polar = @(t,v) fx1y1_4(sin(v)./t,cos(v)./t)./t.^3 ;
            fx1y1_44 =   integral2(fx1y1_4_polar,0,1/k0,0,2*pi,'AbsTol',1e-3);

            Yx1y1_22 =   Yx1y1_22 + omega*epsilon0/(4*pi^2*k0^2)*(fx1y1_22+fx1y1_44) ;%the another part of (65)

   %--------------------------------------Yy1x1------------------------------------------------ %
            fy1x1_1 =   f_p*f_f1*(-kx*ky)/sqrt(k0^2-kx^2-ky^2);
            fy1x1_2 =   matlabFunction(fy1x1_1) ;
            fy1x1_2_polar =  @(t,v) fy1x1_2(t.*sin(v),t.*cos(v)).*t ;
            fy1x1_22 =   integral2(fy1x1_2_polar,0,k0,0,2*pi,'AbsTol',1e-3) ;

            fy1x1_3 =   f_p*f_f1*(-kx*ky)/sqrt(kx^2+ky^2-k0^2)/(-1j);
            fy1x1_4 =   matlabFunction(fy1x1_3) ;
            fy1x1_4_polar = @(t,v) fy1x1_4(sin(v)./t,cos(v)./t)./t.^3 ;
            fy1x1_44 =   integral2(fy1x1_4_polar,0,1/k0,0,2*pi,'AbsTol',1e-3);

            Yy1x1_22 =   Yy1x1_22 -omega*epsilon0/(4*pi^2*k0^2)*(fy1x1_22+fy1x1_44) ;%the another part of (65)

   %-----------------------------------------Yy1y1--------------------------------------------- %
            fy1y1_1 =   f_f*f_f1*(k0^2-ky^2)/sqrt(k0^2-kx^2-ky^2);
            fy1y1_2 =   matlabFunction(fy1y1_1) ;
            fy1y1_2_polar =  @(t,v) fy1y1_2(t.*sin(v),t.*cos(v)).*t ;
            fy1y1_22 =   integral2(fy1y1_2_polar,0,k0,0,2*pi,'AbsTol& ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

发帖回复