[Matlab] practice_7

% simulation execises : 미방을 푼다
% initial value problem(IVP) = ODE + IC problem
% matlab -> state equation / simulink -> integrator method
% state equation : ODE 의 특별한 형태
% 조건 : 1차 미방, 하나의 근을 가져야한다.
% N 차 -> state eq -> N 개의 1차 미방
% pendulum problem

pend1(0.1, [1 2])
% [ tout, xout ] = ode45(@pend1, [0, 30], [1 0]);
whos tout xout
% 0-30초, 1rad-0rad 로 변할 때의 output
plot(tout, xout(:, 1)) % 각도
plot(tout, xout(:, 2)) % 각속도
plot(tout, ones(size(tout)), '.')

% phase plane
plot(xout(:, 1), xout(:, 2))
grid

% IC change
[ tout, xout ] = ode45(@pend1, [0, 30], [0 10]);
plot(xout(:, 1), xout(:, 2))
grid