请教:同一积分表达式采用不同积分函数结果相差巨大
积分表达式如图所示分别采用quadgk quadl quad 来进行计算,具体参数如下:
%%%% 电磁刚度参数数值 %%%%%
delta_0 = 18*10^-3;
Rg = 6.24;
L1 = 2.1; %转子有效长度m
k_u = 1.102; %饱和度
mu_0 = 4*pi*10^(-7); %空气导磁系数
f = 50;
theta = 30.64/180*pi;
phi = acos(0.875); %功率因数角
p = 40; %磁极对数
Fsm = 19210;
Fjm = 24214;
Lambda_0 = mu_0./(k_u*delta_0); %发电机均匀气隙磁导
sigma = k_u*delta_0; %饱和度与均匀气隙大小的乘积
omega_f = 2*pi*f./p; %发电机同步转速
积分表达式如下:
coeff_1 = Rg.*L1.*Lambda_0./(2*sigma.^2);
B = @(alpha,t) Fsm.*cos(omega_f.*t-p.*alpha);
C = @(alpha,t) Fjm.*cos(omega_f.*t-p.*alpha+theta+phi+pi/2);
C11 = @(alpha) 1+cos(2.*alpha);
K11_exp2 = @(alpha,t) coeff_1.*C11(alpha).*(B(alpha,t)+C(alpha,t)).^2;
K11_quadgk = quadgk(@(alpha) K11_exp2(alpha,0),0,2*pi);
K11_quadl = quadl(@(alpha) K11_exp2(alpha,0),0,2*pi);
K11_quad = quad(@(alpha) K11_exp2(alpha,0),0,2*pi);
计算结果分别为
5.0701e^8 3.1042e^8 3.8713e^8
看了一些关于积分函数的介绍,其中quadgk几乎是一个万能的而且精度较高。
三个函数得到三个不同的结果,请问根据这个积分表达式哪一个积分函数得出的结果更加可信呢?
页:
[1]