exa20.m

来自「市面上电力系统的范例 直得收藏学习的ㄧ本书」· M 代码 · 共 16 行

t0= 0;   tfinal =5;               % time interval
x0 = [1, 0];                    % initial conditions
% [t,x] = ode23('pendulum', t0, tfinal, x0);   % use with MATLAB 4
tspan = [t0, tfinal];                          % use with MATLAB 5
[t,x] = ode23('pendulum', tspan, x0);          % use with MATLAB 5
subplot(2,1,1),plot(t,x)
title('Time response of pendulum on rigid rod')
xlabel('t,  sec')
text(3.3,3.1,'Velocity')
text(3.3,1.25,'Angle, rad')
th= x(:,1);   w = x(:,2);
subplot(2,1,2),plot(th, w)
title('Phase plane plot of pendulum')
xlabel('Position,  rad')
subplot(111)