Matlab在常微分方程求解中的应用教程.pptVIP

17世纪：初等解法 18世纪：初等解法和无穷级数方法 19世纪：解的存在性、奇点理论、定性理论、稳定性理论 Van der Pol微分方程 1928年荷兰的范德波耳(Van der Pol)为描述LC回路的电子管振荡器建立了著名的vanderPol方程.它在自激振荡理论中有着重要的意义,一直作为数学物理方程中的一个基本方程.这是一个具有可变非线性阻尼的微分方程,代表了一类极为典型的非线性问题.和其他非线性微分方程在数学上无法精确求解一样,人们一直在努力寻找求解这类方程近似解析解的方法,并乐于用Van der Pol方程来检验求解方法的有效性. odeGetData01.m [t, y] = ode45(vdp1, [0 25], [3 3]); plot(t, y(:,1), t, y(:,2), :); xlabel(Time t); ylabel(Solution y(t) and y(t)); legend(y(t), y(t)); 也可以画出 及 在相位平面（Phase Plane）的运动情况 odePhastPlot01.m [t, y] = ode45(vdp1, [0 25], [3 3]); plot(y(:,1), y(:,2), -o); xlabel(y(t)); ylabel(y(t)); 当 值越来越大时，Van der Pol方程就渐变为一个Stiff系统，此时若要求解该类问题，就必须用对应于Stiff系统的指令 将 值改成 1000， ODE 文件改写成（vdp2.m）： type vdp2.m function dy = vdp2(t, y) mu = 1000; dy = [y(2); mu*(1-y(1)^2)*y(2)-y(1)]; 选用专门用来解决 Stiff 系统问题的 ODE 指令，例如 ode15s，來求解此系统并作图显示 ode15s01.m [t, y]= ode15s(vdp2, [0 3000], [2 1]); subplot(2,1,1); plot(t, y(:,1), -o); xlabel(Time t); ylabel(y(t)); subplot(2,1,2); plot(t, y(:,2), -o); xlabel(Time t); ylabel(y(t)); 由上图可知， 的变化相当剧烈（超过 ），这也就是 Stiff 系统的特点 绘制二维平面相位图： ode15s02.m 若要产生在某些特定时间点的状态变量值，则调用 ODE 指令的格式可改成： [t,y] = solver(odeFile, [t0,t1, …,tn], y0) 其中[t0, t1, …, tn]即是特定时间点所形成的向量 [t, y]= ode15s(vdp2, [0 3000], [2 1]); subplot(1,1,1); plot(y(:, 1), y(:, 2), -o); xlabel(y(t)); ylabel(y(t)) ODE 指令的选项 ODE 指令可以接受第四个输入变量，代表积分过程用到的各种选项（Options），此种ODE指令的格式为： [t,y] = solver(odeFile, [t0,tn], y0, options) 其中 options 是由 odeset 指令來控制的结构变量 结构变量即包含了积分过程用到的各种选项 odeset的一般格式如下： options = odeset(name1, value1, name2, value2, …) name1的值为value1，name2的值为value2… 也可以只改变一个现有的options中某栏位的值，格式如下： newOptions = odeset(options, name, value); 若要获取某栏位的值，可用 odeget，格式如下： value = odeget(otpions, name); 当使用 odeset 指令時，若无输入任何变量名，则 odeset 会显示所有变量名及其值，并以大括号代表预设值 odeset AbsTol: [ positive scalar or vector {1e-6} ] RelTol: [ positive scalar {1e-3} ] NormControl: [ on | {off} ] NonNegative: [ vector