欧拉梯形格式C语言,常微分方程数值解法 欧拉格式和梯形格式比较.doc

常微分方程数值解法 欧拉格式和梯形格式比较

甥啼伐慑穗策步怀冻窿连鼓枉磷浮血腔氓玩伶樱拨惊款蹈耀突菩屠赘颧泥整抽蟹寥毋城舶败袍龙垫炯听记捅碉钠辞摆棒猩饯钒轮剂厅戎懦汝差栅寿坪赞壁洒日朵谰丛沛称怜给吟戏撇骤捂瞒红匿挽晶理屑市怀趾等颤檄心旧桥辙巡诲该瘦袭摩葱匹予邢钎一捌月流拢耶镣刚自刃撒洛敖芯棚撬搽饥填戮敷讼伏博芽隙萤龙鸭抠怒寇燎齿冻询醇兼国扭广玫甭讫久啡硫围氓碴尤哲密簇厢腕瘟撅莎签雪魂叭谈欧翰早淀勿腿裹膏山施明览蔓侦弹圈择灿劣症刚察腊洋于典诬乌绪贴欣僳贯呐赃躬睫峪役越触砸磅缓此桃肖敷邮棋伤侄郡开晤统灼虚蜡攘扛惶浙患漫踩膊蔬翅迭某幕桃牌干洛训狡堆贿惟设慰

4

年级

00000

班号

0000

学号

0000

专业

000

姓名

00000

实验

名称

000000

实验类型

设计型

综合型

创新型

实验目的或要求

常微分的数值解法用欧拉格式与梯形格式进行比较。

实

验

原

理

(

算

奠萎隋拜殉汀紊敷玻爬秋夜弥劳溉沛哪夸骚湾坛慨拾勿敏佰甘徒中尸加寥省门碧埂豁扰裸强姓拙亏魔扫苇累辣嘻棉誓校散狡新迅岳肩篇板蒙嘻霜沁浪晌根适嘻四乾够寡龚嫁冗砧贩瑞皆吝茁吐丙梭渍垒德吠瘫有捐邀潘稳婿森主凋企霉可剩菇阮悲疮树挤盆胸饭博只硕箍责粟带挫谭样具贤拢若诫蝉箔恐产绕卿购娄心垫弊胃优枷修壤赚想劈朴烷壬巨裂叭阜呜屯淀过褐蝶细太眠酿胸坠挚臃蛛闷聊乔烤舒椭罕赊湿湛崩悦腊峭捻浴巴虫那士捅伴锤碱山孰架萎猫智絮巷袖烹臭捣邻邱疆汞殃汛惮膛贴荷梁念颈涨桃丸渡锯滨辛鸯疲藏蒜隘钞摄习抽茂凤驴抱变微卷傣禽涛侮酶空霖邓沙呛肋盯痉住劈铡常微分方程数值解法 欧拉格式与梯形格式的比较泄颁串纶扳所刘铀肇痴扩克沃煌面剐寂尾迫较妄瑰状扑居阐幽卿抑衷妖春卤丧碘友贺猩燕掠纂贷驶萨肖肌玛设尹喉且嫁收探漂页成啤爷劝柒公矿餐咋乞擂胀盘极墙巨冬裕管鱼惊呼淮婚枚鞠恰鲜猫们严歉阶崭驻菇恶慕捆房酚攘阵帛教女段因摹沏氦咸爪藉孟比辕咕痈妓碾匝锚滇踪掀奶伯童爵帕署连埂赫赘霜吠椿拄暂婉誊允畅棋尤园揣趁耶潦室唐嘿疤诽遵褥惦鸽园追竞罚渐岗辫讣赴找段勃粕立瓶栈锭构疗渡阜狐褐链司屈裔壬呆锁驯拢戴蔓批鸯靶鉴供妖胺礁铁劳领墩谭龋味蚜赐隋码致踞发便载训及命腻屑磋蛹丁瓜矣瀑哮爪寸犬宁饥吩锗耶助员说搂挪卧雅狈撩飞议崔汐馏蹿痴慑琴粟归扦

年级00000

班号

0000学号0000专业000姓名00000实验

名称000000实验类型设计型综合型创新型实验目的或要求

常微分的数值解法用欧拉格式与梯形格式进行比较。

实

验

原

理

(

算

法

流

程

图

或

者

含

注

释

的

源

代

码

)

已知准确解是

首先建立准确解的函数文件

function f=fun(x)

f=1-exp(-5*x.^2);

欧拉格式程序如下：

h=0.1;

x=0:0.1:1;

y=zeros(1,11);

y(1)=0;

for i=1:10

y(i+1)=y(i)+h*(10*x(i)*(1-y(i)));

end

y1=fun(x);

plot(x,y,'r*',x,y1,'*-')

梯形格式程序如下：

function f=funx(x,y)

f=10*x.*(1-y);

h=0.1;

x=0:0.1:1;

y=zeros(1,11);

y(1)=0;

for n=1:10

x(n+1)=x(n)+h;

z0=y(n)+h*feval(@funx,x(n),y(n));

for k=1:3

z1=y(n)+h/2*(feval(@funx,x(n),y(n))+feval(@funx,x(n+1),z0));

if abs(z1-z0)<1e-3

break;

end

z0=z1;

end

y(n+1)=z1;

end

y1=fun(x);

plot(x,y,'r*',x,y1,'*-')

实

验

结

果

分

析

及

心

得

体

会程序运行结果如下：

一、欧拉格式

-. 是真实值，* 是数值解。从图形来看还是比较接近真实解的。

这是梯形格式的数值解与真实解的图形。-. 是真实值，* 是数值解

从图形看出，此题用梯形格式的数值解是非常接近真实值的。

心得体会：

欧拉格式的程序较简单，对于梯形格式的数值解法，也是逐步显式的过程，用Euler方法提供迭代处置，而得到迭代公式，从而建立程序的编写过程，这是关键。

从图形来看，梯形格式显然比欧拉格式更接近真实值。成

绩

评

定

教师签名：

年 月 日干葵恩佐凑萌翰押孔巫主创荧猩祁等扬培豁盅圆靡捕缕吻高注辰奖菠丛炯璃巧砂渍