|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用matlab的两个函数trapz、quad2d分别求二重定积分的结果相差近10倍,
' Y+ r1 }% [4 j9 I- J5 u I5 |trapzq求解的速度快想用它代替掉quad2d。但是发现两个函数积分结果不一样。可能是tarapz那出错了,但不知道错在哪一步。请教大神,非常感谢!
0 G+ O/ w* V& m. fx = 0.1:0.01:0.5;% P2 |% [! l" k0 g T% |& p
y= 0.3:0.001:1;
7 l G/ b0 ^- F& ?9 Q' n[X,Y] = meshgrid(x,y);
. d+ p( d6 ~( Q, X; q0 hF = cos((pi*Y).*(sin(X*pi)));
\$ b" j! H8 o$ A) x+ t7 i( L0 y: C2 `I = trapz(y,trapz(x,F,2));
6 {* ]) b9 [8 P8 t2 M1 h2 s2 S: Q' o8 r) i! r$ Y! O
f =@(x,y)cos((y).*(sin(x)));
) u2 j- g6 v( E* R5 \a = quad2d(f,0.1*pi,0.5*pi,0.3*pi,pi);
+ Y5 G# v8 n% S' }9 q5 }- Q/ [8 T |
|