找回密码
 注册
关于网站域名变更的通知
查看: 443|回复: 1
打印 上一主题 下一主题

不知道错在哪一步怎么办?

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-4-28 13:39 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

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
  • TA的每日心情
    开心
    2019-11-19 15:19
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2020-4-28 15:10 | 只看该作者
    可以随便更改 积分区间?简单放缩可以?4 A% x$ N7 s  j" d- W! S) I
    假设 f(x)=x,[0,2],可以简单化为  f(x)=2*x,[0,1]????  Y7 A1 S( z% n7 M' Y) }
    # M+ F3 s9 ]) b+ o( v/ R: l. j
    a=quad2d(@(X,Y)cos((pi*Y).*sin(X*pi)),0.1,0.5,0.3,1)' a5 y( c5 F0 ~* S; Z9 e! k1 q7 [; w
    6 {. V7 y5 ^- b6 O5 u
    a = 0.0102
    " t' ?& `; Y5 K! O1 ?5 J8 P与I接近。
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-11-24 05:37 , Processed in 0.140625 second(s), 23 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表