basic | Van der pol |
페이지 정보
작성자 cemtool 작성일14-03-31 13:55 조회3,499회 댓글0건본문
Van der pol
Van der Pol 2차 비선형 미분방정식의 해를 구해보자.
2차미분이기 때문에 ode23 혹은 ode45 함수를 사용할 수 있도록 해야한다. 아래의 식과 같이 1차 미분 형태로 변환한다.
1차 미분의 형태로 변환되었으면, CEM 모드와 M 모드에서 아래와 같은 사용자정의 함수로 구현한다. |
function xdot <> t,x xdot = zeros(2,1); xdot(1) = x(1)*(1-x(2)^2)-x(2); xdot(2) = x(1); return; // demo ////////////////////////////// t0 = 0; tf = 20; x0 = [0 0.25]; [t,x] = ode45("vdpol",t0,tf,x0); plot(t,x(:,1)); hold on; plot(t,x(:,2)); xtitle("t [sec]"); ytitle("x(t)"); title("Van der pol");
댓글목록
등록된 댓글이 없습니다.