chunshui2003 发表于 2010-7-27 11:51

请教:同一积分表达式采用不同积分函数结果相差巨大

积分表达式如图所示

分别采用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]
查看完整版本: 请教:同一积分表达式采用不同积分函数结果相差巨大