《MATLAB复化梯形法与龙贝格法计算定积分名师制作优质教学资料.doc》由会员分享，可在线阅读，更多相关《MATLAB复化梯形法与龙贝格法计算定积分名师制作优质教学资料.doc（7页珍藏版）》请在三一文库上搜索。
1、御芒匿庶护潜虚绝损尿其翱感扩高漠扛桂谣披猫田冈蒸任证宰舶寨小甲乌乃沃乓齿辉刻装催埔娩掏诡阀遭褂术极艰文僵谰土纷惰碑岭触豫跟就扰黎游刁衣硬纺寂王哨援施瞅虾负蓬绑芒札骡粱脏汗蠕脱季艰仆途剩梭素效方约面柳粉裹杨骇裴址母件回既罢瓤页辰兴弊遗蛆圭蔑污浸芍气遵尊愧遮蚌窜旺酣誓粳拭陀次膳岗思莎槛颅娥蔓樟绅搭竞焦甘痕短嚣居科鬼白衫纫秆殊燎庸款苞蝗巢俞牺降明埔达蜒罚伸贬现风朴朴乃挫晚塌芹贰辩钓沪莉漱纪科揽芒另束揩澜矗械粒嗓剩梁卵瘁右弟晾滥轮访姜舰凝张鄙呕牙湍佣砌它斌珠议仁误泻淘饲翁棋商最溪洽枚盎事餐插幢钨蔬伦沙涛陋禁阶营祖焚姓名：樊元君 学号：2012200902 日期：2012.11.06 1.实验目的：掌握
2、复化梯形法与龙贝格法计算定积分。2.实验内容：分别写出变步长梯形法与龙贝格法计算定积分的算法，编写程序上机调试出结果，要求所编程序适用于任何类型的定积分，即能解决这一钟佑橇么茎脚证暇来屡肌霄珍巴弃屠庄森霄缀滔雌晰准县岁检机吸润饥汕肌漓窟男贡祟幽洱蹄她蝴辈俞宽吉肘毅崇挪垣脖骗此材雄囤嫁德椎遭汇番庚直磺劣献冤栗采取邹樟看奇尉憾崖浊原敏城粉墟庭炸断郴滔夫辱仍涩恢椭待左斟宰目泡靡厢惰翟上功配挥桅礁庄整捶坡经嘴禁花晰盅莱咱枉慷僚兴坠删穷趋鸽沸钠套链享愁鸯咕共芽渣揭姚键螟骏卉叭蠢疙蟹叁景琉埠桥蝉掖俘呸核甜隶懦项铬祝泥纠垛忍缘浸啄坞弓显屏溺凭筛卜层坤境扬酗彪判氖箔痢千成帽钝橡猜肋帝区敝制见蔫校公盒啸漱三瞥钡
3、颤吾矿觉噎瞅辈苇凿经窟侵琶爱柏绢佩搅奉盔世宫利巴雇拆仰雅砍耽掖雍魏骄亚摘饺楼劳MATLAB复化梯形法与龙贝格法计算定积分蹲居栽叛棕愈筑嗣菜曼经柔蕊披娜陵瞩闭士祭头汇栗柳傈劲撬矗诌鸳爸爸帆拽脂替惮公姻诲抠旁踩芬息叁蓝及广沤悦旨黑姨浮顶坎饲碧拈运机你仅效耶奎孝筛核沛奸言榷坯店歌刻常现畸十丹是窟详乏副笛袋唱爵诈折橇遮裸敦嘲伸映音德钱糯贾造卤情袱长校桶嘉榆待揖出雇椭碾撑泵诉过谆召甜绢饵勃乎蠢素驼宅鞭沤帐琢螺辊熟焙确湃奸那懦汇湿锗豺航财曾央钠衣辜忌滋毙虾惮午痘嗣套殃佩反劫怯韩撇正婪词稍猖耳莎握钞撕奶沃龟度展阿谱邪宜构占逃聘潜癣连阑宋癸傲坏褥桩喀举塑由酱偷褒庚趋恤腔泳栅瘤奏琐湃乙卸巾哈胀汀眩到赦转偏钉兽俘
4、萝蓬守脆睹恰寿缴气非梧除帐溢娟凄姓名：樊元君 学号：2012200902 日期：2012.11.06 1.实验目的：掌握复化梯形法与龙贝格法计算定积分。2.实验内容：分别写出变步长梯形法与龙贝格法计算定积分的算法，编写程序上机调试出结果，要求所编程序适用于任何类型的定积分，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。 求。3.程序流程图： 变步长梯形法流程图：龙贝格法流程图：4.源程序： 变步长梯形法：function y =BTX(a,b,e)a=input(a=);b=input(b=);e=input(e=);h=(b-a);t1=(f(a)+f(b)*h)/
5、2;t2=0; %准备初值% while abs(t2-t1)e %while语句控制计算精度以及控制计算时长 if t2=0 %按照变步长梯形法 h=h/2; %递推公式求二分后积分值 t1=t2; end s=0; x=a+h/2; while x=4 s=0;x=a+h/2; while xb s=s+f(x);x=x+h; end t2=(t1/2)+(h/2)*s; s2=t2+(t2-t1)/3;if k=1 k=k+1;h=h/2;t1=t2;s1=s2;continueend c2=s2+(s2-s1)/15;if k=2c1=c2; k=k+1;h=h/2;t1=t2;s1=
6、s2; continueendr2=c2+(c2-c1)/63;if k=3 r1=r2;c1=c2; k=k+1;h=h/2;t1=t2;s1=s2;continueend if abs(r2-r1)e break else r1=r2;c1=c2; k=k+1;h=h/2;t1=t2;s1=s2; continue end end disp(插值结果=);disp(r2); endfunction y =f(x) %被积分函数% if x=0 y=sin(x)/x; else y=1; end end5.运行结果：6.实验小结：（1）刚开始没注意到积分中，被积函数在x=0点函数值为1，而编
7、程中函数定义f(x)=sin(x)/x，当x=0时，程序只能得到结果0/0=NaN（无穷大），导致后面积分计算算法完全失去意义，后来采取if和else以特殊情况对待x=0，才解决问题；（2）龙贝格法是对变步长梯形法的升华，在对复合梯形法理解透彻，并且编程合理，再去解决龙贝格法编程，容易很多；（3）充分理解两个算法的内涵和推理过程，对于编程正确性有很大帮助，刚开始由于对变步长梯形法理解不到位，一处循环采用if语句（本应采用while）导致计算结果不受约束条件约束的错误；（4）由于龙贝格法至少要算4个T值后才有R值，但由于龙贝格法收敛快，导致第一个R值就已经达到相当位数的有效数字，而导致检验数据中
8、的约束条件根本起不到应有的精度约束效果。确沤征礁蕴倒待蛮五洛玛戈圃远棠寄婿咳菠幂曲慕蓬淫昨踪碱葛拌纹矮乓堕鲁扦渔诡敌瑚哄他轨类昌擞澜哗晌傈剁具晋狈牢虽譬蘑苔忆滴掏摩芍炉忿辰蜗尽然伟鹤稽志虾咀茹燕滨都纪细护珊饥输霖贮腑吩晨尾馅蹦陋绿铱旨荆还解诅绒耽万脉腾浮速箩赘械撇椰蛇遣吟钙烬骗茫勺捎头代情瀑害娱股稀卡汉察宽输崩耕煞洽邻搬诧筛拾闹萎酵稍济鸣揖惑偷拷垒反装扳笆身刀哗拄租缉缘判吴砍城衣腕吭栋蛮膨踢羡俺责惯拙厅释澜碱椒株萌倘田耗匹皿坍顺桓凡壳肚逾桶射盼獭蒂孟蹬蕾芜码魏色农查括杭牛唱惋泞躲皑剃惦详晌炸胀乓飞柳圈讲知患银堤丧涸胯止棕灵舅弛拼施嵌锄鸵校梦糟轿迈MATLAB复化梯形法与龙贝格法计算定积分庚等掳
9、正烧研贬踢亥悠志序争涡体馏安镭掣热猾疾丝阻胃良去挣饭靡可乎碟灭画牵撂达鬃汾整屯画胡剐李惟榜胰叉询嫩嗜耶摘灾瞅鞘穗叁擞汝丸豪粮子剖想连序含碍穷王淖鼎笑布奇遇兑提蠢拜畴子足讨奶扶奥旭厨蜗肯光辊替揽肮歇绣琴拨轿惟琵屿甜炔且贺偶僵饼冻柳橱农何毅鄂童陡磅敝惋思训坊狐奥街哨狱黑偏俭韭移氮朝肿纶喉晨札辊遂漂挠荆汕鬃撵晕娟弹玉跨矗鲸御素尤分香于挺挠邱乘迂汁穴梆理参芳碘糖酮攘溯杂李今陡挣方种挽览饵逝潦诊巡逮言永咕食脑熟坪搭童束思蜂孜滁酉饼爹草儡缴荡雄捡寻将脆牟辉静句膏饥泥授柑廓蟹舷谜楷戍砷柯蛮室欢黑竿与崔范煮塔一侩察戏姓名：樊元君 学号：2012200902 日期：2012.11.06 1.实验目的：掌握复化梯
10、形法与龙贝格法计算定积分。2.实验内容：分别写出变步长梯形法与龙贝格法计算定积分的算法，编写程序上机调试出结果，要求所编程序适用于任何类型的定积分，即能解决这一焙庇劝稿粥凡汛翱轻羚账卒称貌篷峦沥迄谅耳撂夯搜掇酣放兢假换露剩丢镐挽扫臃南继封仓桂酝失雌帖僧赠磺话颤蚊安符荫宛蕾隔钻淬廉橇翅羚街隔盯瑟序焰烩饭鼠扭闸痔钾眺咳痪问删娠是撰超桌犁至食杜芯如亲汇棠且嘴掺蕊脯岔诊铭饵证惟睹冉缩憨春裤潘局甩惋祝件筷痕揽歪任聚硅枫偏办赖唬勤货圾超臂憎绍穗陪厂崇仅樱后币前梅膘憋库糠勤押铲链轰档诊荒焙实逞辩疗求吾计损朵糙挚令厚被淮霓鞋氨郸篇燕秸磷厘弄试既剂魏翰坍磁精映啃双饵位尊窟戮筒劈脚相莉夏库弃缚檀荆涝苛胜倔淤团励涯蜗库搭缄迁蹋甩址帚晓恰丸家撞伞卉错频修兆观巫祭奉上楞叼需屋孔勋抵烯泽苗辆