Метод Рунге-Кутта 2-го порядка, Паскаль

Xstyle74 · 09.05.2014, 05:18

Помогите с заданием, плиииззз, вообще ниче понять не могу,
Задание
Ток в электрической цепи описывается дифференциальным уравнением (см. стр. 24), которое необходимо решить, найдя y(0.1), y(0.2) ... y(1) с помощью метода Рунге-Кутта второго порядка с коррекцией по средней производной (в дифференциальном уравнении k = 6 По найденным значениям y определить количество тепла, выделяющееся на единичном сопротивлении за единицу времени, по формуле с помощью метода трапеций.
Уравнение со стр.24

Vladimir_S · 09.05.2014, 08:28

Цитата:

Сообщение от Xstyle74

Помогите с заданием, плиииззз, вообще ниче понять не могу

Пожалуйста. У нас разбиралась абсолютно аналогичная задача, только функция другая и метод интегрирования надо поправить: там - Симпсон, у Вас - трапеции. Метод Рунге Кутта
От поста #5 и дальше.

Xstyle74 · 09.05.2014, 21:29

Если бы я разбирался еще...если не сложно то помогите плиз,ну паскаль вообще не дается мне....

Vladimir_S · 09.05.2014, 22:36

Цитата:

Сообщение от Xstyle74

Если бы я разбирался еще...если не сложно то помогите плиз,ну паскаль вообще не дается мне....

Очень это грустно. На самом деле переделка той программы заняла полторы-две минуты. Ладно, вот вам для тупого бездумного скатывания и сдавания.

Код:

Const
 h=0.1;
 N=10;

Var
 I:Array[0..N] of Real;
 j:integer;
 t,Z,Sum:Real;

Function F(x,y:real):real;
begin
 F:=(6.0-y*y)*Cos(x)+y*2;
end;

Begin
 I[0]:=0.3;
 t:=0;
 writeln('t= ',t:3:1,'   I= ',I[0]:6:3);
 for j:=1 to N do
  begin
   t:=t+h;
   Z:=I[j-1]+h*F(t-h,I[j-1]);
   I[j]:=I[j-1]+h*(F(t-h,I[j-1])+F(t,Z))/2;
   writeln('t= ',t:3:1,'   I= ',I[j]:6:3);
  end;
 writeln;
 Sum:=(Sqr(I[0])+Sqr(I[N]))/2*h;
 for j:=1 to N-1 do
  Sum:=Sum+h*Sqr(I[j]);
 Writeln('Heat: ', Sum:0:5);
 Readln
End.

Xstyle74 · 10.05.2014, 00:13

Прибольшущее тебе спасибо, да пусть в твоем доме интернет будет без трабл, и потерь пакетов, Говарю тебе как работник одной из компании провайдера интернета

Vladimir_S · 10.05.2014, 08:28

Цитата:

Сообщение от Xstyle74

Прибольшущее тебе спасибо, да пусть в твоем доме интернет будет без трабл, и потерь пакетов, Говарю тебе как работник одной из компании провайдера интернета

Спасибо. И да, поправил одну неточность в части интегрирования: сложение квадратов первой и последней точек (было сложение самих токов). Приношу извинения - поторопился.

09.05.2014, 05:18	#1 (permalink)
Xstyle74 Новичок Регистрация: 09.05.2014 Сообщений: 3 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 10	Метод Рунге-Кутта 2-го порядка, Паскаль Помогите с заданием, плиииззз, вообще ниче понять не могу, Задание Ток в электрической цепи описывается дифференциальным уравнением (см. стр. 24), которое необходимо решить, найдя y(0.1), y(0.2) ... y(1) с помощью метода Рунге-Кутта второго порядка с коррекцией по средней производной (в дифференциальном уравнении k = 6 По найденным значениям y определить количество тепла, выделяющееся на единичном сопротивлении за единицу времени, по формуле с помощью метода трапеций. Уравнение со стр.24 Миниатюры

09.05.2014, 21:29	#3 (permalink)
Xstyle74 Новичок Регистрация: 09.05.2014 Сообщений: 3 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 10	Если бы я разбирался еще...если не сложно то помогите плиз,ну паскаль вообще не дается мне....

10.05.2014, 00:13	#5 (permalink)
Xstyle74 Новичок Регистрация: 09.05.2014 Сообщений: 3 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 10	Прибольшущее тебе спасибо, да пусть в твоем доме интернет будет без трабл, и потерь пакетов, Говарю тебе как работник одной из компании провайдера интернета

Ads
Яндекс Member Регистрация: 31.10.2006 Сообщений: 40200 Записей в дневнике: 0 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 55070