Matlab程序代码

%%
clear;clc;
deltax = 0.01;
x = 0;
n = 101;
x = zeros(1,n);
y = zeros(1,n);
x(1) = 0;
y(1) = 0;

for i = 2:n
    x(i) = x(i-1) + deltax;
    y(i) = y(i-1) + 2 * x(i-1) * deltax;
end
plot(x,y, '-o')

%%
clear;clc;
deltat = 0.01;
n = ceil(2*pi/deltat) + 1;
t = zeros(1,n);
x1 = zeros(1,n);
x2 = zeros(1,n);
t(1) = 0;
x1(1) = 0;
x2(1) = 1;

for i = 2:n
    t(i) = t(i-1) + deltat;
    x1(i) = x1(i-1) + x2(i-1) * deltat;
    x2(i) = x2(i-1) - x1(i-1) * deltat;
end
plot(t, x1, '-', t, x2,'-o');