输入函数必须返回 'double' 或 'single' 值。找到 'sym'。

syms theta0 z0;

miu=4*pi*10^(-7);

a=0.015;

h=0.03;

x=0:0.005:0.03;

z=0.002;

Z=z+0.015-z0;

M=8.68*10^5;

I=M;

theta=atan(x./(z+0.002));

r=(x.^2+(z+0.015)^2).^0.5;

dBx=miu.*M.*a.*cos(theta0).*Z./(4.*pi.*(x.^2+a^2-2*a.*r.*cos(theta-theta0)+Z.^2).^(3/2));

f=@(theta0,z0)dBx;

Bx=integral2(f,0,2*pi,0,h);

错误使用 integral2Calc>integral2t/tensor (line 231)

输入函数必须返回 'double' 或 'single' 值。找到 'sym'。

出错 integral2Calc>integral2t (line 55)

[Qsub,esub] = tensor(thetaL,thetaR,phiB,phiT);

出错 integral2Calc (line 9)

[q,errbnd] = integral2t(fun,xmin,xmax,ymin,ymax,optionstruct);

出错 integral2 (line 106)

Q = integral2Calc(fun,xmin,xmax,yminfun,ymaxfun,opstruct);