Решить систему двух нелинейных ОДУ - Matlab - Страница 2

Зосима · 21.02.2015, 16:33

Студворк — интернет-сервис помощи студентам

значит они небезрамерные!

*не запускал расчеты со значениями констант из твоего предыдущего сообщения?

@Dimarik21 · 21.02.2015, 17:09 **[ТС]**

Зосима, вопрос на миллион, как вывести график первой функции дебита? что надо изменить в

Matlab M
1
2
3
[L U] = ode45(@func1,[0 100],[0 20e6]);
n = 130; %кол-во отсчетов
plot(L(1:n),U(1:n,2))

сейчас вроде кое что получается

Добавлено через 10 минут
Понял как. Буду разбираться с людьми, которые эти формулы и придумали теперь. Может оин подскажут, что где да как с размерностями тут.

Зосима · 21.02.2015, 17:17

гляди, если у тебя давление в МПа, то в нач. условиях (внутри ode45) должно быть [0, 20], без "е6"

@Dimarik21 · 21.02.2015, 17:23 **[ТС]**

Да, я так и сделал. Когда рассчитываю распределение давления на 100 метров, оно слишком быстро падает. Видимо проблема в коэффициентах, надо с ними играть. Спрошу у человека, кто эти уравнения выводил, какие размерности он брал, потом буду продолжать разбираться.

@Dimarik21 · 24.02.2015, 21:49 **[ТС]**

Альфа \frac{кг*МПа}{м^7} = 0,06
Q \frac{м^3}{сек}
P МПа = 20,7
D метр = 0,2
Ag \frac{МПа^2*сут}{тыс.м^3} 0.06
Bg \frac{МПа^2*сут}{(тыс.м^3)^2}

Никак не могу разобраться, почему вылетают комплексные корни, которых быть не должно(

Зосима · 25.02.2015, 15:16

Неа, единицы все-равно не сходятся

@Dimarik21 · 25.02.2015, 16:16 **[ТС]**

Зосима, я чуть чуть ошибся. Все оба равнения сходятся.
Альфа \frac{кг*МПа}{м^7} = 40
Q \frac{м^3}{сек}
P МПа = 20,7
D метр = 0,2
Ag \frac{МПа^2*сут}{тыс.м^3} 0.06
Bg \frac{МПа^2*сут^2}{(тыс.м^3)^2} 0,0003

Добавлено через 43 секунды
при таких данных получаются комплексные числа, поччему???

Новые блоги и статьи Все статьи Все блоги /
Запустил конкурс "тем и промптов для текстовых квестов созданных почти чисто ИИ" Adler 06.10.2026 Всем привет! За последние три-четыре дня я создал более 16 текстовых квестовых игр используя преимущественно по одному запросу к ИИ на игру. Мне так понравилось смотреть все ветки/ сцены во всех. . .	ИИ не может найти нужный язык в списке Supersumestria 05.10.2026 Я ему даю вот такое изображение и прошу найти и подчеркнуть немецкий язык. Возвращает он вот это: https:/ / i. **********/ vqBWLe2. png Нужную строчку в 3й колонке просто выдумал. . Это. . .	Новая последняя моя музыка в SUNO zorxor 05.10.2026 Здравствуйте, дорогие мои друзья! С большой радостью я хотел бы представить вам свою новую последнею музыку, которую сгенерировала мне по моей просьбе нейросеть SUNO. С уважением, zorxor. Это. . .	Программный домашний кинотеатр russiannick 27.09.2026 Сподобился на программный домашний кинотеатр. В качестве ЯВУ по традиции выбрал js. В помощники взял Яндекс-Алису. Было создано три зала на разные интересы. исторические и ретро сериал Хичкок. . .
Беседа с ИИ о программистах, недопускающих к созданию и правке кода генеративные ИИ и причины этого zorxor 21.09.2026 Раньше я радовался или получал некоторые эмоции, пусть небольшие, но всё же, от самого процесса написания кода, рекомпиляции и запуска, видя постепенное развитие программы и прочее. А теперь лень. . .	Мобильное приложение ColorStep pavlinmavlin 17.09.2026 Реализовал приложение Красный, Зеленый, Синий в Unity3d + c#. Название изменил на ColorStep. Приложение прошло модерацию и теперь доступно для скачивания. Делал его сам, шаг за шагом — и вот,. . .	Запрет дублирования строк в табличной части Maks 13.09.2026 Реализация из решения ниже выполнена на нетиповом справочнике "Нормы ТО" с табличной часть "Виды ТО", разработанного в КА2, со следующими реквизитами: - ВидТО (СправочникСсылка. ВидыТО); - ВидГСМ. . .	Скрипты Tampermonkey для CyberForum, ChatGPT, Claude и пр. Jin X 06.09.2026 Скрипты Tampermonkey для CyberForum, ChatGPT, Claude и пр. Работая с форумом и нейросетями в браузере часто хочется что-то подкорректировать или добавить какого-то функционала. Ниже прикреплён. . .

Зосима 5246 / 3574 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	21.02.2015, 16:33
	значит они небезрамерные! *не запускал расчеты со значениями констант из твоего предыдущего сообщения? 0

Зосима 5246 / 3574 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	21.02.2015, 17:17
	гляди, если у тебя давление в МПа, то в нач. условиях (внутри ode45) должно быть [0, 20], без "е6" 0

@Dimarik21 0 / 0 / 0 Регистрация: 18.02.2015 Сообщений: 18
	21.02.2015, 17:23 [ТС]
	Да, я так и сделал. Когда рассчитываю распределение давления на 100 метров, оно слишком быстро падает. Видимо проблема в коэффициентах, надо с ними играть. Спрошу у человека, кто эти уравнения выводил, какие размерности он брал, потом буду продолжать разбираться. 0

@Dimarik21 0 / 0 / 0 Регистрация: 18.02.2015 Сообщений: 18
	24.02.2015, 21:49 [ТС]
	Альфа \frac{кгМПа}{м^7} = 0,06 Q \frac{м^3}{сек} P МПа = 20,7 D метр = 0,2 Ag \frac{МПа^2сут}{тыс.м^3} 0.06 Bg \frac{МПа^2*сут}{(тыс.м^3)^2} Никак не могу разобраться, почему вылетают комплексные корни, которых быть не должно( 0

Зосима 5246 / 3574 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	25.02.2015, 15:16
	Неа, единицы все-равно не сходятся 0

@Dimarik21 0 / 0 / 0 Регистрация: 18.02.2015 Сообщений: 18
	25.02.2015, 16:16 [ТС]
	Зосима, я чуть чуть ошибся. Все оба равнения сходятся. Альфа \frac{кгМПа}{м^7} = 40 Q \frac{м^3}{сек} P МПа = 20,7 D метр = 0,2 Ag \frac{МПа^2сут}{тыс.м^3} 0.06 Bg \frac{МПа^2сут^2}{(тыс.м^3)^2} 0,0003 Добавлено через 43 секунды* при таких данных получаются комплексные числа, поччему??? 0