|
Lorenz 방정식의 상태-공간 표현 (4th Runge Kutta Methode)
방정식 초기값은 t = 0
dx/dt = -10x + 10y x = 5, y = 5, z = 5
dy/dt = 28x - y - xz
dz/dt = -2.6666667z + xy
Table 작성
t x y z k11 k12 k13
0.0 5.000000 5.000000 5.000000 0.000000 110.000000 11.666665
0.1 9.781470 17.079464 10.439
|
|
|
f) = %lf u2(%1.1f) = %lf error = %lf\n", t, w2[i], t, u2(t), u2(t)-w2[i]);
}
}
double f1(double t, double w1, double w2){
return 3*w1 + 2*w2 - (2*t*t+1)*exp(2*t);
}
double f2(double t, double w1, double w2){
return 4*w1 + w2 + (t*t + 2*t - 4)*exp(2*t);
}
double u1(double t){
return exp(5*t)/3 - exp(
|
|
|
매틀랩 M_FILE
Matlab Curvefitting Gragh
*Reference
Matlab의 입문과 활용, 김용수저, 높이깊이 -Runge-Kutta C언어 소스
-결과 출력값
-매틀랩으로 그래프 출력
-선형회귀분석 C언어 소스
-결과 출력값
-매틀랩으로 그래프 출력
|
|
|
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");
// 도
|
|
메뉴펼치기
