09.05.2014, 05:18 | #1 (permalink) |
Новичок
Регистрация: 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, 05:18 | |
Helpmaster
Member
Регистрация: 08.03.2016
Сообщений: 0
|
По содержанию ваш топик аналогичен данным темам Метод Рунге-Кутта 2-го порядка, Паскаль Метод Рунге Кутта Решить дифференциальное уравнение 4-го порядка в PASCAL Метод Ньютона-Рафсона на Паскаль Методом Рунге-Кутта найти решение дифференциальных уравнений |
09.05.2014, 08:28 | #2 (permalink) |
Специалист
Регистрация: 27.08.2008
Адрес: Санкт-Петербург
Сообщений: 27,807
Сказал(а) спасибо: 340
Поблагодарили 583 раз(а) в 208 сообщениях
Репутация: 113184
|
Пожалуйста. У нас разбиралась абсолютно аналогичная задача, только функция другая и метод интегрирования надо поправить: там - Симпсон, у Вас - трапеции. Метод Рунге Кутта
От поста #5 и дальше. |
09.05.2014, 22:36 | #4 (permalink) | |
Специалист
Регистрация: 27.08.2008
Адрес: Санкт-Петербург
Сообщений: 27,807
Сказал(а) спасибо: 340
Поблагодарили 583 раз(а) в 208 сообщениях
Репутация: 113184
|
Цитата:
Код:
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. |
|
10.05.2014, 00:13 | #5 (permalink) |
Новичок
Регистрация: 09.05.2014
Сообщений: 3
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
|
Прибольшущее тебе спасибо, да пусть в твоем доме интернет будет без трабл, и потерь пакетов, Говарю тебе как работник одной из компании провайдера интернета
|
Ads | |
Member
Регистрация: 31.10.2006
Сообщений: 40200
Записей в дневнике: 0
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 55070
|
10.05.2014, 08:28 | #6 (permalink) |
Специалист
Регистрация: 27.08.2008
Адрес: Санкт-Петербург
Сообщений: 27,807
Сказал(а) спасибо: 340
Поблагодарили 583 раз(а) в 208 сообщениях
Репутация: 113184
|
Спасибо. И да, поправил одну неточность в части интегрирования: сложение квадратов первой и последней точек (было сложение самих токов). Приношу извинения - поторопился.
|
Ads | |
Member
Регистрация: 31.10.2006
Сообщений: 40200
Записей в дневнике: 0
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 55070
|
|
|