|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
用matlab的两个函数trapz、quad2d分别求二重定积分的结果相差近10倍，& S9 }* D6 \9 K8 ?
trapzq求解的速度快想用它代替掉quad2d。但是发现两个函数积分结果不一样。可能是tarapz那出错了，但不知道错在哪一步。请教大神，非常感谢！
; u+ z# t1 ?0 ?) B0 ]& U6 n! h* Kx = 0.1:0.01:0.5;
* o8 P6 S9 [! r. M: H( ~y= 0.3:0.001:1;
( O" s" T0 z: p: ^% |* l, a[X,Y] = meshgrid(x,y);8 o* S" Z* E* a" p0 P2 p
F = cos((pi*Y).*(sin(X*pi)));
6 n! Z C N5 H) s* U6 J# E% nI = trapz(y,trapz(x,F,2));
, W( S$ r4 Z5 o3 U1 P! t7 d. P. w
" `2 [( k8 h* l, t6 s5 J. Nf =@(x,y)cos((y).*(sin(x)));. ?8 s: b9 m, a6 r6 L' \5 \
a = quad2d(f,0.1*pi,0.5*pi,0.3*pi,pi);
" h3 X X2 o% T0 n9 h |
|
/1 
发表于 2020-4-28 13:39
收藏
淘帖
支持！
反对！