VPF::[C]Метод Рунге-Кутты 4го порядка

Lenk

Дата 18.6.2008, 00:27 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 2
Регистрация: 18.6.2008

Репутация: нет
Всего: нет

нужно разработать программу численного интегрирования  уравнения :
d^2x/dt^2=-dU/dx,
где U(x)=e^(-x)  +x-1.
при   t=0    x=0, v=dx/dt=v0

не понимаю, как найти зависимость x(t) и поэтому не получается.

Код


#include <stdlib.h>
#include <stdio.h>
#include <math.h>
main()
{double v1,v2,v3,x1,x2,h,t,k11,k12,k13,k14,k21,k22,k23,k24,a;
int i;
FILE *inp,*out, *out1;
inp=fopen("input.txt","r");
out=fopen("output.txt","w");
out1=fopen("output1.txt","w");
fscanf(inp,"%lf %lf",&x1,&v1);
h=0.01; 
for(i = 0;i <= 5000;i++)
{
  k11=exp(-x1)-1;
  k12=exp(-x1-h/2*k11)-1;
  k13=exp(-x1-h/2*k12)-1;
  k14=exp(-x1-h*k13)-1;
  v2=v1+h/6*(k11+2*k12+2*k13+k14);
  k21=v2;
  k22=v2+h*k21/2;
  k23=v2+h*k22/2;
  k24=v2+h*k23;
  x2=x1+h/6*(k21+2*k22+2*k23+k24);
  x1=x2;
  v1=v2;
  fprintf(out,"%lf %lf\n",x1,v1);
  fprintf(out1,"%lf %lf\n",t,x1);
}

fclose(inp);
fclose(out);
}

помогите, пожалуйста!

GIK

Дата 18.6.2008, 15:16 (ссылка)

(нет голосов)

Загрузка ...

Добрый человек

Профиль
Группа: Участник
Сообщений: 985
Регистрация: 3.6.2005
Где: я только не небыв ал

Репутация: 4
Всего: 14

Чет я смысла немогу понять smile

а когда то соображал в интегрировании...
Может тебе цикл надо создать, в котором изменяеться t? или я туплю smile

А вообще интегрирование несет в себе прикладной характер, это интересно и занимательно smile

--------------------

Математика=>пиво=> програмирование, три вещи последовательны и совместимы !!!
Программирование - это не деятельнось! Программирование - это состояние души!
Бог - самый крутой программист.

kBepTu

Дата 18.6.2008, 19:30 (ссылка)

(нет голосов)

Загрузка ...

Шустрый

Профиль
Группа: Участник
Сообщений: 68
Регистрация: 6.12.2007

Репутация: 2
Всего: 2

Цитата

не понимаю, как найти зависимость x(t) и поэтому не получается

Задачу надо переформулировать следующим образом: обознач dx/dt через y, численно решаешь уравнение dy/dt = -dU/dx. В цикле находишь y (методом Рунге Кута - численно решаешь дифур). Потом просто у тебя есть производная - численно ее интыгрируешь sum(y*delta(t)).
Алгоритм почти готов smile

пользуйся

Это сообщение отредактировал(а) kBepTu - 18.6.2008, 19:32

Lenk

Дата 18.6.2008, 20:44 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 2
Регистрация: 18.6.2008

Репутация: нет
Всего: нет

спасибо большое, щас буду пробовать!

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »