实验6常微分方程数值解.docVIP

实验六 常微分方程数值解 问题：的数值解。 一、欧拉方法function [u,v]=el(a,b,c,h) x=a:h:b; y(1)=c; m=length(x); for n=1:m-1 y(n+1)=y(n)+h*f(x(n),y(n)); end plot(x,y,*) u=x; v=y; 2、隐式差分格式（向后欧拉格式） clear all; x=0:0.1:2; y(1)=1; m=length(x); for n=1:m-1 y(n+1)=1/1.1*(0.1*x(n+1)+y(n)+0.1); end plot(x,y,+) x,y 又原方程的精确解为，将3个解的图形放在一起进行比较。 二、R-K二阶-校） 你的程序：function [u,v]=el(a,b,c,h) x=a:h:b; y(1)=c; m=length(x); for n=1:m-1 k1=y(n)+h*f(x(n),y(n)); y(n+1)=y(n)+h/2*(f(x(n),y(n))+f(x(n+1),k1)); end t=exp(-x)+x plot(x,y,*,x,t) u=x; v=y; 的图形解做比较）： 2、四阶 你的程序： function [u,v]=el(a,b,c,h) x=a:h:b; y(1)=c; m=length(x); for n=1:m-1 k1=f(x(n),y(n)); k2=f(x(n)+h/2,y(n)+h/2*k1); k3=f(x(n)+h/2,y(n)+h/2*k2); k4=f(x(n)+h,y(n)+h*k3); y(n+1)=y(n)+h/6*(k1+2*k2+2*k3+k4); end t=exp(-x)+x plot(x,y,*,x,t) u=x; v=y; 你的图形解做比较）： 三、分方程组的求解 （选 求解Volterra者和被捕食者模型： 参数初值： 理论知识： 最后结果 参考： 两个函数： function u=f(x,y1,y2) u=y1-0.01*y1*y2; function u=g(x,y1,y2) u=-y2+0.02*y1*y2; 程序为：clear all x=0:0.1:16; h=0.1; y(1)=20; z(1)=20; m=length(x); for n=1:m-1 k1=f(x(n),y(n),z(n)); l1=g(x(n),y(n),z(n)); k2=f(x(n)+h/2,y(n)+h/2*k1,z(n)+h/2*l1); l2=g(x(n)+h/2,y(n)+h/2*k1,z(n)+h/2*l1); k3=f(x(n)+h/2,y(n)+h/2*k2,z(n)+h/2*l2); l3=g(x(n)+h/2,y(n)+h/2*k2,z(n)+h/2*l2); k4=f(x(n)+h,y(n)+h*k3,z(n)+h*l3); l4=g(x(n)+h,y(n)+h*k3,z(n)+h*l3); y(n+1)=y(n)+h/6*(k1+2*k2+2*k3+k4); z(n+1)=z(n)+h/6*(l1+2*l2+2*l3+l4); end plot(x,y,x,z) 附：分组名单 星期二下午5-6节 第1组组长：陈絮莹 缪 妃 何贵堂 刘 钰 马 倩 第2组组长：李杰玉 黎筱惠 雷 霞 肖 娴 林碧珍 朱元正 第3组组长：陈 静 苏小丽 李 郑 何淑楠 田冬秀 曾敬军 第4组组长：杨 欣 王雪梅 徐莉萍 石小芳 雷 敏 唐 嘉 第5组组长：杨佳悦 郭 滢 李媛媛 何 可 陈思露 第6组组长：王钰琪 寇 玠 杨 丹 熊晨曦 周 丹 第7组组长：姚 瑶 高倩倩 金 杨 周海宁 杨 琴 第8组组长：雷 芳 陈 艳 王玉娇 余 非 张 雪 王海燕 星期三下午5-6节 第9组组长： 刘超慧 王玉利 秦佳丽 张青梅 廖 婷 程思远 第10组组长：杨 琴 冯康欣 黄宜纯 田晓东 郑美艳 第11组组长：黄 倩 肖雪梅