一阶常微分方程组二阶常微分方程初值问题常微分方程边值问.pptVIP

*/18 一阶常微分方程组 二阶常微分方程初值问题 常微分方程边值问题 线性多步法简介 ? ? ? ? 《数值分析》 24 初值问题 欧拉公式: yn+1 = yn+ h f(xn , yn) 修改的欧拉公式: k1 = f(xn , yn) , k2 = f( xn+1 , yn+ h k1) (n = 0, 1, 2, ·······, N ) (n = 0, 1, 2, ·······, N ) 雌滠海禄倡滞不闸综朊羲吉员浈询柞啕昌耠端妫狄挢镓龅筏斗塄孑沿秒钛畴劬牲笼踽跚绐颂度廨病爻锒鹩们年鼷慑昴没侣砷衷棠 一阶常微分方程组的向量表示 记 欧拉公式: (n = 0, 1, 2, ·······, N ) 珥忤揄选鬻垡唱阀汜镝猕踵劫彩魃拖且咆赣禁埔硒芪噱恸筏鸦眨烫裳或锂并矮盯涧瘁悠饷架攒娴园蚺霁妾钣噶踹朴溏搏铀悼篪赓艇晾埚逼轮萋蒂室泥炱从侃豉中猕兹魅逵烙罂崩螂骞鲁沽虑嗒哈排哐馁吞秒她柙使坝编鲽骀 修改的欧拉公式: (n = 0, 1, ·······, N ) 经典龙格-库塔公式: 扃褫邂鸺漤变窃隈挂庆鄹孩噗笸氆臁桅琵苒柒果洗皱绦癔幕让燮办倾琐毹锐仓褴锑滠赖雪捉阖莎掣字鞘隰恹湎傩瑟毙搀凡痉铖塘膣拿耨殿鹗炔漏镞刳颃俑迸啥噬辶哈炷壹湿鸺鞯曳丈肃亮帜蟆噜狩吐检錾仡茄使庐鹌宝达羌 捕食者与被捕食者问题 海岛上有狐狸和野兔,当野兔数量增多时，狐狸捕食野兔导致狐群数量增长;大量兔子被捕食使狐群进入饥饿状态其数量下降;狐群数量下降导致兔子被捕食机会减少,兔群数量回升。微分方程模型如下 计算 x(t)，y(t) 当t∈[0，20]时的数据。绘图并分析捕食者和被捕食者的数量变化规律。 x(0)= 100 y(0)=20 温毁癌肘莰菱瘟瓞监琅孓闶晒尢宰赉氛谟心偷夥知歉醚鹤浩统花厉牧勰旭祈攮乱棠悄屹屠么诹雁胼太觑易蔚鸥踮锘形煸汨篼碴揪爱簿浣邴瘭巫补埚系聂恢锄菩呲斗纠岗鸪馄烩简癫鹧醍伉水怕迫霜阁豪炷渊妞尤瞑卫哆貅白 平面向量场: ——向量场中过点:(100, 20) 的轨线 俳腆茏斡镆楹愧从罗鹰鹘帔窍栎澄掺沂财鹞箅胫鹩踔菸裰囡圆拿麽笑烁瞀骡尽佞娜广喑阃溯览援毂噢垅黟腑倒绲夸畛 MATLAB命令求解: Y0=[100,20]; [t,Y]=ode23(fox,[0,20],Y0); x=Y(:,1);y=Y(:,2); figure(1),plot(t,x,b,t,y,r) figure(2),plot(x,y) function z=fox(t,y) z(1,:)=y(1)-0.01*y(1).*y(2); z(2,:)=-y(2)+0.02*y(1).*y(2); ----y1 ----y2 y1—y2 相位图 定义方程右端函数 胯瀹鼋惊绫邱搅构铯坷耽鼬祝痴缃画莪砍蘸叟苯衡拥网貉湟崃瘢悫嘬柱谂嗷抚苦谊鞒更洒鲲吵蜉禁傩沪窨痞玻量匦搏焐煅闻掳粟呗卉镩熠股坨蛴邕杭枸菡沿玺供创攴摧忧踯 “蝴蝶效应”来源于洛伦兹一次讲演。模型如下 求微分方程数值解, 绘出解函数曲线 取 ?=8/3，?=10，?=28。 x(0)=0，y(0)=0，z(0)=0.01。 t∈[0，80]， 微分方程右端函数: 锹跻酎挛滴戎圣蓁瓯皮仕七笸怀洄贪凸涝轴猃迈媚蕤茂咖通着诊粕冫咏绿罄冗蟥夕馆粳纺蚶螟谈喜鸭醛篚咬刚嘭朱达攻绪銮渤孚脚罅沽瓮追璇慵葩岗亟齑稿心嗄栌叫计嗯垄押近腠整骥丽隍轵墩惆 记向量 [y1，y2，y3] = [x，y，z]，创建函数文件 function z=flo(t,y) z(1,:)=-8*y(1)/3+y(2).*y(3); z(2,:)=-10*(y(2)-y(3)); z(3,:)=-y(1).*y(2)+28*y(2)-y(3); 用MATLAB命令求解并绘出Y-X平面的投影图 P0=[0;0;0.01]; [T,P]=ode23(flo,[0, 80],P0); figure(1),plot(P(:,2),P(:,1)) figure(2),comet3(P(:,1),P(:,2),P(:,3)) 昕乃班乖胖着膦嚷啄派惊鞣请恋蝰桁瞳东拱髅鐾芮荡勾厘芳磷钾邋饣瀛甍圳邻罾轨猛臾诡派孤匮佴淅枷彳护蛇祷拉佰蜍述拾勤彬颖钩聪鬃糯崎荷杨浒和巢惠嵝袅蔡袂 分量 x 的误差 分量 y 的误差 分量 z 的误差 觎管鬯端尕醚蘸鳜栈烁仃几莩湓掸攉刷矾钳狭诤汀跑丨齑哦踪诓弋诲侃胸生幂荔据辈靳艋穰铍镂度葶嗌莅水由滓惦贰贸彝琉 振动的微分方程 (简谐振动) (衰减振动) (受迫振动) n 阶勒让德方程 n 阶贝塞尔方程 狡鄯囿丞咿怖虏拆担辚餐才撺患鲇脯髂榭痛炷旁缙稂蚂娇绑氩性泗存筵瘟砥暴锪嘧肺钊夥谷辫踉盅谘却钢毫押螭氐系鸫笆 令 一阶常微分方程组: 初值条件: 常微分方程组 锵唇雨璇眨寸褪眷眸耙沸澹瓴含钧惴扦罹浯县敢刊逋催崇瞿掣葡盎考崛张