|
Runge-Kutta 적분기
void rk4(double* y, double x,
void (*f)(double* dy, double* y, double x),
double h, int n);
// 운동방정식에 대한 연립 1차 미분방정식
void derv4(double* dz, double* z,double t);
#define DIM 4
void main()
{
int n;
double *y, h, t;
FILE *fp;
fp=fopen("problem8_2.dat", "w");
// 도
|
|
|
스프링상수와 감쇠계수를 2배 한 경우(실험자 임의)
◆ 종합 고찰 및 느낀점
◆ 심화 학습
1) 오일러법
2) 개선된 오일러법
3) Runge Kutta법
◆ 이상적인 현가장치
1) 완전한 노면점착(perfect road holding) - 딱딱한 스프링 시스템
2)
|
|
|
x2dot];
2)
function xdot=young(t,x)
global c n
c = 10;
n = 5;
x1dot = 3*(x(1)+x(2)-(x(1)^3/3)-1.3);
x2dot = (c*x(1) - n*x(2)^2);
xdot=[x1dot;x2dot]; ◎ 안정도 판별
◎ 선형 비선형 모든 미분방정식 풀수 있는 방법:Runge-kutta법
◎ 2계 미분 방정식(R-L-C시스템)
|
|
|
x2dot];
2)
function xdot=young(t,x)
global c n
c = 10;
n = 5;
x1dot = 3*(x(1)+x(2)-(x(1)^3/3)-1.3);
x2dot = (c*x(1) - n*x(2)^2);
xdot=[x1dot;x2dot]; ◎ 안정도 판별
◎ 선형 비선형 모든 미분방정식 풀수 있는 방법:Runge-kutta법
◎ 2계 미분 방정식(R-L-C시스템)
|
|
|
Springer Verlag, 2006. ISBN 0-387-26102-8.
http://en.wikipedia.org/wiki/Runge%E2%80%93Kutta_methods
http://en.wikipedia.org/wiki/Runge-Kutta_methods 1.상미분 방정식
1) 1계 상미방
2) Euler Method
① 전방 Euler
② 수정된 Euler
③ 후방 Euler
④ Euler 법의 정확도
⑤ Euler 법
|
|
메뉴펼치기
