Новые сообщения · Правила  
  • Страница 1 из 1
  • 1
Модератор форума: Berestovskiy  
Форум ПРОГРАММИСТОВ » ПРОГРАММИРОВАНИЕ » Паскаль » Метод Рунге-Кутты 4го порядка Pascal (Метод Рунге-Кутты 4го порядка Pascal Для уравнений Киргофа)
Метод Рунге-Кутты 4го порядка Pascal
Здравствуйте. Нужен взгляд человека понимающего Pascal. Расчёт уравнений двигателя постоянного тока по второму закону Киргофа методом Рунге-Кутты 4го порядка, кратко о нём:
1.Закон Киргофа для ДПТ:
U=I*R+L*∂I/∂t+ω*Ce;
I*Cm=J*∂ω/∂t+Mст

далее уравнения приводятся к форме задачи Коши:
∂I/∂t=(U-I-R-ω*Ce)/L
∂ω/∂t=(I*Cm-Mст)/J

далее уравнения дифференцируют для удобства записи (на шаге k+1):
Ik+1=Ik+Δt*(U-Ik-R-ωk*Ce)/L
ωk+1k+Δt*(Ik*Cm-Mст)/J

и собственно сам метод Рунге-Кутты
k11 = dt*((U-Ik*R-ωk*Ce)/L)
k21 = dt*((U-Ik*R-ωk*Ce)/L) + k11/2
k31 = dt*((U-Ik*R-ωk*Ce)/L) + k21/2
k41 = dt*((U-Ik*R-ωk*Ce)/L) + k31/2

k12 = dt*((Ik*Cm-Mc)/J)
k22 = dt*((Ik*Cm-Mc)/J) + k12/2
k32 = dt*((Ik*Cm-Mc)/J) + k22/2
k42 = dt*((Ik*Cm-Mc)/J) + k32/2
тогда
Ik=(Ik+1/6*(k11+k21*2+k31*2+k41))
ωk=(ωk+1/6*(k12+k22*2+k32*2+k42))

Теперь пытаюсь осуществить этот метод на Pascal. В 20  строке вещественное переполнение, могут быть и другие ошибки, поэтому и прошу вашей помощи.

ссылка на код, что бы не прывышать лимит символов
http://pastebin.com/kkjPC0mL
1 | Автор: pascal | 2013-03-27, 17:06 | Изменено: pascal - Ср, 2013-03-27, 17:07   |  Репутация: [ + 0 ]
1 | Автор: pascal | 2013-03-27, 17:06 | Изменено: pascal - Ср, 2013-03-27, 17:07   |  Репутация: [ + 0 ]
Форум ПРОГРАММИСТОВ » ПРОГРАММИРОВАНИЕ » Паскаль » Метод Рунге-Кутты 4го порядка Pascal (Метод Рунге-Кутты 4го порядка Pascal Для уравнений Киргофа)
  • Страница 1 из 1
  • 1
Поиск: