matlab解常和微分方程的特别方法.pptVIP

运行命令文件runf3.m [t,x]=ode45(function3,[0,20],[30;20]); plot(t,x); xlabel(time t0=0,tt=20); ylabel(x values x1(0)=30,x2(0)=20); 例5 创建函数function4，存在function4.m中 function f = function4(t,x) global U f = [x(2);U*(1-x(1)^2)*x(2)-x(1)]; 再运行命令文件runf4.m： global U U = 7; Y0=[1;0]; [t,x] = ode45(function4,0,40,Y0); x1=x(:,1);x2=x(:,2); plot(t,x1,t,x2) 数学科学与技术学院 胡金燕 lionfr@ 前沿： 科学技术和工程中许多问题是用微分方程的形 式建立数学模型，因此微分方程的求解有很实 际的意义。 一、常微分方程（组）的符号解 二、常微分方程（组）数值解 一、常微分方程（组）的符号解 函数 dsolve 格式： r = dsolve(eq1,eq2,…’,cond1,cond2,…,v) 说明 (1)对给定的常微分方程（组）eq1,eq2,… 中指定的符号自变量v，与给定的边界条件 和初始条件cond1,cond2