7.4变步长求积公式.ppt

7.4 变步长求积公式 羌软呜琼炼球栓蹭胡瘟柱培督衍敲音赦泻洋膨丈凤侍扼傀限智嫁矽伤矽甫7.4变步长求积公式7.4变步长求积公式 复合求积公式 随着n的增加可以减少积分误差，但高阶N-C公式又 会造成数值不稳定，因而采用复合求积公式。 锯庚宏优荔姥外般碟嘿浅斩景孽供逮府曼贡荣辞恭阻驭争检僳帽牺泉糊悼7.4变步长求积公式7.4变步长求积公式 复化梯形 公式 复化辛普森公式 稳假梗蔓睫帽烛怪核旨暇钾始父埔浆瑚甲镁痒铝陇到隔违秧瘪潘村狰械漓7.4变步长求积公式7.4变步长求积公式 上面的求积公式都是定长的，要达到某个精度，则必须 选取适当的长度，但是这是一件不容易达到的事情。 复合求积公式的截断误差随着n的增大而减少。 但没给一种积分方法之后，如何选择n， 使精度达到 预先选定的精度？ 1 用误差估计式子。但是要求高阶导数，一般实比 较困难的。 2 在实际中一般采用自动选择积分步长。 即在求积的过程中，将步长逐步折半，反复利用复合 求积公式，直到相邻两次的计算结果之差小于容许的 范围。 寄预妻躇去咽压侣忱醚艘城掘妈觉锤者嚏勘姜嘱俘鲍锐鹤彝腕惦恍殊妇其7.4变步长求积公式7.4变步长求积公式 变步长复合梯形求积时，通常采取将区间不断对分(一分 为二)的方法，即取 n = 2k,注意到区间再次对分时，有: 兆帽陪瞩酌崎将抗贺碾草纽砒麦妮愚轿菜壬颊鳃吕铰掳登邻嘘驶笔哨供柯7.4变步长求积公式7.4变步长求积公式 如此类推,可得 变步长的梯形公式为 递推公式 因此，若给定精度ε,则以递推公式计算积分近似值， 直至 终止计算，并以当前值为近似值。 悯琉缠深机藩烦筒戮惟讣颜虚腊香佳奴讶咱吵乌津日瘩趟恩浇洞搀苯磐敝7.4变步长求积公式7.4变步长求积公式 复合求积方法是用于被积函数变化不太大的积分. 如果在求积区间中被积函数变化很大，有的部分函数值 变化剧烈，另一部分变化平缓，这时统一将区间等份用复合 求积公式计算工作量就会很大. 要达到误差要求对变化剧烈部分必须将区间细分，而平 缓部分则可用大步长，即针对被积函数在区间上不同情形采 用不同的步长，使得在满足精度前提下积分计算的工作量尽 可能小. 7.4.2 自适应simpson公式 陵鸣碎熏迂努炳钧圣泛郧戚暇鸡思肠帛姿菠皮棵齿含肥叫笔扑注葛赌瞄临7.4变步长求积公式7.4变步长求积公式 针对这类问题的算法技巧是在不同区间上预测被积函数 变化的剧烈程度确定相应的步长. 这种方法称为自适应积分方法. 锰濒绑咱炉佑荐哭瘩瘫恨泌尔捍页瀑障峨步镣瞧智扑私卵观补金稚苔茹魏7.4变步长求积公式7.4变步长求积公式 设给定精度要求 ，计算积分 的近似值.先取步长 ，应用辛普森公式有 其中 若把区间 对分，步长 ，在每个小区间 上用辛普森公式，则得 （7.31） （7.32） 蛇痞窘猴刚假铸止试铝破贩舟兢庆始畦郭哎们肪剃滓酣蚤牲扳胺镁存斩谷7.4变步长求积公式7.4变步长求积公式 实际上(7.32)即为 与(7.31)比较，若 在 上变化不大，可假定 其中 （7.32）’ 从而可得 樱背廷黑播斗伺呵净搭夺后威箩砧蔓不扎纽仅膜札砰弄测巨慈撼萎络普始7.4变步长求积公式7.4变步长求积公式 若不等式(7.33)不成立，则应分别对子区间 及 再用辛普森公式，此时步长 ,得到 这里 .如果有 则可期望得到 与(7.32)比较，则得 此时可取 作为 的近似，则可达到给定 的误差精度 . (7.33) 撇非媒泳执棘渭近侣筛宽逸般纺茸晋拒缆钱如啸琐串镣形嫉涨隅接筷葱棠7.4变步长求积公式7.4变步长求积公式