|
|
Регистрация Восстановить пароль |
Регистрация | Задать вопрос |
Заплачу за решение |
Новые сообщения |
Сообщения за день |
Расширенный поиск |
Правила |
Всё прочитано |
|
Опции темы | Поиск в этой теме |
14.04.2014, 10:31 | #1 |
Пользователь
Регистрация: 24.09.2010
Сообщений: 22
|
В ответе inf и nan (решений уравнений методом Эйлера).
Здравствуйте!
При реализации задачи на С++, в ответе получила inf (как я понимаю, бесконечность) и nan. Возможно, влияет то, что нужно для начала инициализировать все переменные (включая Ucbf, Uc1f...)? Или необходимо выбирать более мелкий шаг (dt)? Решение легко нашлось в системе Wolfram Mathematica, по ним строились графики. Помогите, пожалуйста, найти проблему! Думаю, что метод Эйлера реализован правильно (решение проверялось на закомментированной системе //Ucf = Uci + Ili/C*(tf-ti); //Ilf = Ili + (E-Uci)/L*(tf-ti) Исходные данные: E2 = 5; Cb = 2 10^-12; C1 = 10^-6; C2 = 10^-6; i = 10^-12; R = 10^6; Rb = 20; R1 = 1000; R2 = 1000; R3 = 1000; L1 = 10^-5; E1[t] := 10 Sin[20000 Pi t]; Система: UCb'[t] == (1/Cb) ((1/Rb) (E1[t] - E2 - UC1[t] - UCb[t]) - i (Exp[UCb[t]/0.026 - 1]) - (1/R) UCb[t]), UC1'[t] == (1/C1) (IL1[t] + (1/Rb) (E1[t] - E2 - UC1[t] - UCb[t]) - (1/R1) UC1[t]), UC2'[t] == (1/C2) (IL1[t] - (1/R3) UC2[t]), IL1'[t] == (-1/L1) (UC1[t] + UC2[t] - E1[t] + R2 IL1[t]), UCb[0] == 0, UC1[0] == 0, UC2[0] == 0, IL1[0] == 0. Код:
t ucb uc1 uc2 il E1 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000010 219463.130731 0.438926 0.000000 5.877853 5.877853 0.000020 -inf -109670.316484 58.778525 -5862.903032 9.510565 0.000030 nan inf -58570.839575 5966661.176993 9.510565 0.000040 nan nan 59608626.638747 -inf 5.877853 0.000050 nan nan -inf nan -0.000000 0.000060 nan nan nan nan -5.877853 0.000070 nan nan nan nan -9.510565 0.000080 nan nan nan nan -9.510565 0.000090 nan nan nan nan -5.877853 0.000100 nan nan nan nan 0.000000 И т.д. Последний раз редактировалось KEIego; 14.04.2014 в 11:40. Причина: Это редактирование дополнит суть проблемы. |
Похожие темы | ||||
Тема | Автор | Раздел | Ответов | Последнее сообщение |
проблема с системой дифференциальных уравнений методом Эйлера | Anubys | Помощь студентам | 0 | 24.04.2011 13:53 |
РЕШЕНИЕ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ методом Эйлера | ruslan 91 | Паскаль, Turbo Pascal, PascalABC.NET | 1 | 10.01.2011 22:12 |
решение систем диф.уравнений методом Эйлера??? в VisualBasic | goog2 | Помощь студентам | 0 | 14.12.2010 21:50 |
Задача Методом Эйлера и методом Рунге-Кутта. Прошу помочь. Очень срочно. | BeNeDiKT | Паскаль, Turbo Pascal, PascalABC.NET | 0 | 12.05.2009 13:14 |