QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4732|回复: 1
打印 上一主题 下一主题

[问题求助] 关于数值积分的程序

[复制链接]
字体大小: 正常 放大
753085848 实名认证       

8

主题

3

听众

448

积分

升级  49.33%

  • TA的每日心情
    开心
    2011-10-7 20:45
  • 签到天数: 3 天

    [LV.2]偶尔看看I

    自我介绍
    目前正在数模培训中
    跳转到指定楼层
    1#
    发表于 2010-11-22 21:08 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    function [I,n]=TX(f,a,b,eps)
    ( V1 [9 y6 M  G  Y0 d  dh=b-a;5 `5 G4 w4 K' _& t# X
    n=1;: `2 a4 L: @  K7 b. D$ v
    eps=1e-3;. w' a. ]9 }0 `3 Y4 q9 ]3 e
    T1=0;
    % ]/ J+ r: `& ^! _8 r" xT2=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));6 @( Q1 {8 _/ G2 p! n9 f% S& Q
    while abs(T2-T1)>eps
    2 _) M2 _7 Y$ C  K5 P% W; M    h=h/2;/ z$ t4 |; k, I3 k2 I& n
        T1=T2;+ |5 D7 ~5 ~" Y  K8 I9 a8 T, S
        x=a+h;% q' W3 i5 L. N  u! H2 y
        s=0;" k, p$ m1 Y; Y3 T# X) o; d/ K
        while x<b* v- N# R& A+ s0 \) c( Q
            s=s+subs(sym(f),findsym(sym(f)),x);$ i' j/ x& a+ \9 q8 M: l% z
            x=x+h;   
    6 s+ @5 s  u0 P) R' C    end& _& D' G  I8 b: v) Q
        T2=1/2*(T1+h*s);
    6 {  |- n4 e: X    n=n+1;/ u/ @- `, L' G2 ~) W, m5 I$ @
    end7 Z- {& q8 f5 f! O& v
    I=T2+1/3*(T2-T1);
    ' j  `# X! g8 I! V' ^2 h 这是一个关于变步长梯形积分的程序,但是怎么调试就是得不出正确的结果,请各位帮忙看看。帮忙修改一下,顺便指出我的错误,感激不尽!
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    753085848 实名认证       

    8

    主题

    3

    听众

    448

    积分

    升级  49.33%

  • TA的每日心情
    开心
    2011-10-7 20:45
  • 签到天数: 3 天

    [LV.2]偶尔看看I

    自我介绍
    目前正在数模培训中
    自己修改过的,是正确的!4 e" X" ~4 y* M, D0 Z
    function [I,n]=FHTX(f,a,b,ep)
    ' O$ }' s: k" a3 o0 R) Va=a+eps;
    ! `7 q2 C# L! X/ r5 Wh=b-a;
    8 R8 ^2 ~8 C, o$ M3 w# cn=1;" N3 z% }6 a2 C4 f& n) g
    ep=1e-6;
    7 e2 B! Z  A- `3 l7 W) A. ST1=1/2*h*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
    ' S% Y+ d% o6 c8 \3 etol=1;
    # |% f4 I! Q* }# G; `6 J) xwhile tol>ep! W7 s/ w( E; {7 N3 q
        u=h/2;5 O2 _/ w# `; g# V, k
        x=a+u;% ]* H3 ]/ W" y
        H=0;
    : Q; x% W4 J! p+ |2 k    while x<b
    9 R6 H0 a8 u/ J+ R        H=H+subs(sym(f),findsym(sym(f)),x);
    7 ^% k# y# {" S0 y  Z/ c9 S' K) }        x=x+h;    . ?" v, j7 [- k+ F
        end
    8 a/ I/ G" x: Y/ S; u: X7 Y    T2=1/2*(T1+h*H);: ~' Y/ D0 I6 w" e& J1 f
        tol=abs(T2-T1);
    ' Q5 W% ?# z% V+ w, j    h=u;9 X7 G! c% b) P
        T1=T2;6 i. ?# q3 l. D, v; k
        n=n+1;# }/ `- }3 _* B5 k
    end
    ' {+ G6 V8 D: ]1 E0 i  w' I% KI=T2+1/3*(T2-T1);
    # L# V3 \  u. E* E. s4 I
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-5-29 22:51 , Processed in 0.439556 second(s), 56 queries .

    回顶部