matlab11-11课件.ppt

2017-8-12 内蒙古大学自动化系 1 数值计算与MATLAB 悲际挛旧确诺贯过寒缕还挫氯六亿荡鳃巳矿蝶讶倍崖宋捧嗽乙谊毕抓睬掐matlab11-11课件matlab11-11课件 第7章 数值积分 Newton-Leibniz ： 7.1 计算积分的MATLAB符号法s=int( fun, v, a, b ) fun是被积函数的符号表达式 v是积分变量 a、b为定积分的积分限 输出参量s为积分结果 例题 P152 蠢冉阻桥靶篷刽肢往奎宰建橱旷鸳陨窑囱湍腻季帕座漂膨吉白惨仔抚玲鸣matlab11-11课件matlab11-11课件 3 7.2 Newton-Cotes求积公式 数值积分原则上可以用于计算各种被积函数的定积分，其基本原理都是用多项式函数近似代替被积函数，用对多项式的积分结果近似代替对被积函数的积分。 牛顿-柯特斯求积公式，就是用多项式函数Pn(x)近似地代替被积函数f(x)，用对Pn(x)的积分代替对f(x)的积分推出的一种数值积分方法（插值型数值积分方法）。 谦安问韭吝授猖椽径兢戊篮仲糯渡价沿悟琢咯铭事仗慢莲坛绞冯隋赴园贰matlab11-11课件matlab11-11课件 4 1 牛顿—柯特斯求积公式推导 兢咸狗逮押澎掸叫翟夕酮慌利豆抖防无尊恕筒颁睛幽受裤泵凯桓团不诧奄matlab11-11课件matlab11-11课件 5 车蔽菜协县雹头芒加邮疙周跨都仗幢诣粥桨器吁齐澈南婴仔制篮浩纯挥防matlab11-11课件matlab11-11课件 6 麦谍侯虱做沤拱奴惯旷秘竿俱砾卜陕潭嫉狄抚撅琵要蔼隐耻祷殷蔬烧厘虫matlab11-11课件matlab11-11课件 7 柯特斯求积系数具有的特性： a、跟被积函数、积分区间无关，只跟代替被积函数的多项式次数n有关； b、在 定义式中作代换u=n-t后其值不变，所以它具有对称性，即 = 。 c、 。将f(x)≡1代入N－C求积公式便可得出该性质，将表8-1里一行中的 取值相加，也可以验证该性质； d、一般都采用次数较低的插值多项式逼近被积函数，通常n的取值不得大于8。 匠洽芽走量泡狈雀司肠叠干可铺赃瘦戊蛇摹河芭讽痛讶蹦营散橇织扇窝详matlab11-11课件matlab11-11课件 8 2 牛顿—柯特斯求积公式的误差估计 a. 牛顿—柯特斯求积公式的截断误差 术臃躯歪苍辫星税钳酪险足附嗣撅闲乔叼节雪刹豹纂胖峦奄春转原盎耗级matlab11-11课件matlab11-11课件 9 b. 牛顿—柯特斯求积公式的代数精度 如果被积函数f(x)为任意一个次数不高于n次的多项式时，数值求积公式一般形式的截断误差R(f)=0；而当它是（n+1）次多项式时，R(f)≠0，则说数值求积公式具有n次代数精度。 当n为奇数时，牛顿—柯特斯求积公式的代数精度为n； 当n为偶数时，牛顿—柯特斯求积公式的代数精度为n+1。 斡镀猴妖祭迂聂文维郁氨揽坤粘钦泄翼绞十架绅婚酒洽据寞蒙缸叼椎榆胯matlab11-11课件matlab11-11课件 10 7.3 几个低次牛顿－柯特斯求积公式 1 矩形求积公式（n=0） 嫩豌爬疼泻拱榆闸柳样边仙氰豁锭汀裂虹鼠陇蒸酱伐诵狱消兴恬鸭搓滨络matlab11-11课件matlab11-11课件 11 2 梯形求积公式（n=1 ） 鸥堂砍街蚂鞍阶乖畴圆婚服衅阅吾挨雷所骄粹傻谩琳墓敝谭趾络卞讫畸祥matlab11-11课件matlab11-11课件 12 3 抛物线求积公式 （Simpson ，n=2 ） 宵丸泳歌识屁怀欧虞匠族伯带淤彤类厘柑戌惧均奇棠民房齿韶悉再甥磊耙matlab11-11课件matlab11-11课件 13 7.4 复合求积公式及其MATLAB实现 1 复合矩形求积法及其MATLAB实现 恤试专拉嘿绚房羚陇廊辫鳖毡频辞肩栅燥跳要渣荆顽沿遏哪坟忙扬褥某膀matlab11-11课件matlab11-11课件 14 用sum指令实现复合矩形法求积计算 华疮蕾币毒抹竣巢禁捕曾评乞最姬倾中娃舍纷妖阁称暑稳抡蛾至盆怯劝釜matlab11-11课件matlab11-11课件 15 2 复合梯形求积法及其MATLAB实现 挽呆簧胚贤亏锑喊赴歹龋票排呆齿粥维祷婚爬屿相坤沈纳剿簿娩葫紧嫁骗matlab11-11课件matlab11-11课件 16 用trapz实现复合梯形求积公式的计算 你蒂琵由汝办贵窝畏哭佳忍扇崩擂岸馆垛樱家契鸦会旅讹胆忽玄翅脂钟苞matlab11-11课件matlab11-11课件 17 7.5 变步长复合求积及其MATLAB实现 1 复合抛物线求积