|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用matlab的两个函数trapz、quad2d分别求二重定积分的结果相差近10倍,
9 t( @+ N' \, N% w# jtrapzq求解的速度快想用它代替掉quad2d。但是发现两个函数积分结果不一样。可能是tarapz那出错了,但不知道错在哪一步。请教大神,非常感谢!
: O a [( V/ O, [9 C+ P$ |- hx = 0.1:0.01:0.5;4 c/ B) ~+ ?; M+ r' ~/ d1 d2 |
y= 0.3:0.001:1;1 ?* E1 Z% ^: C0 f; y' o6 a; {
[X,Y] = meshgrid(x,y);
6 O& p) _+ ?5 G7 XF = cos((pi*Y).*(sin(X*pi)));, d0 p! @% {; N/ @
I = trapz(y,trapz(x,F,2));! T& s$ Q7 w0 w2 k3 ~3 J: g
7 X A& l( j! g. D$ Qf =@(x,y)cos((y).*(sin(x)));+ S- j/ `* T) O9 f, v. z/ _ ?0 }: ~$ R
a = quad2d(f,0.1*pi,0.5*pi,0.3*pi,pi);) ]$ N* U/ L. {# ~9 l4 ]
|
|