wjlvdp.m

(有源代码)数值分析作业,本文主要包括两个部分,第一部分是常微分方程(ODE)的三个实验题,第二部分是有关的拓展讨论,包括高阶常微分的求解和边值问题的求解(BVP).文中的算法和算例都是基于Matla

%function wjlvdp
% van der pol equation 

%
clear
[t1,Y1] = ode45('ODEfun5',[-8:0.08:16],[-2,-2]);
[t2,Y2] = ode45('ODEfun5',[-10:0.1:20],[-0.7,-0.7]);
[t3,Y3] = ode45('ODEfun5',[1:0.1:31],[2,0]);
%[2,0] [2;0] 似乎效果是一样的
%注意这里的点的构造数是一样多的,要不下面一句就会错,如0.1变为0.2
numsol5 = [t1,Y1,t2,Y2,t3,Y3];
figure
plot(t1,Y1(:,1))
hold on
plot(t1,Y1(:,2),'--r')
figure

subplot(3,1,1),plot(t1,Y1,'-g');
subplot(3,1,2),plot(t2,Y2,'-b');
subplot(3,1,3),plot(t3,Y3,'-r');
figure;
plot(Y1(:,1),Y1(:,2),'-g');
hold on;
plot(Y2(:,1),Y2(:,2),'-b');
plot(Y3(:,1),Y3(:,2),'-r');