Система дифференциальных уравнений в форме Коши

@М_Скофилд · Регистрация: 14.12.2013

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

Здравствуйте! Вынужден обратиться к вам, так как сам не смог справиться с задачей

необходимо решить систему ОДУ в нормальной форме Коши
начальные значения не даны по заданию, поэтому их выбираю произвольно

вот код функции

Matlab M
1
2
3
4
5
6
7
8
9
10
function dy=pr8(t,y)
dy=zeros(8,1);
dy(1)=-1000.*y(1)+1000.*y(3)-1000.*y(5)-1000;
dy(2)=-333,33.*y(2)-333,33.*y(3)+333,33.*y(4)+333,33.*y(5);
dy(3)=111,11.*y(1)-200.*y(2)-311,11.*y(3)+200.*y(4)+311,11.*y(5)+111,11;
dy(4)=142,86.*y(2)+142,86.*y(3)-157,14.*y(4)-157,14.*y(5)+14,29.*y(6)+14,29.*y(7);
dy(5)=-61,73.*y(1)+111,11.*y(2)+172,84.*y(3)-728,4.*y(4)-200,3.*y(5)+11,11.*y(6)+27,45.*y(7)-61,73;
dy(6)=-111,11.*y(4)-111,11.*y(5)-144,78.*y(6)-144,78.*y(7)+33,67.*y(8);
dy(7)=-142,86.*y(4)+352,94.*y(5)-186,15.*y(6)-1110,52.*y(7)+43,29.*y(8);
dy(8)=60,61.*y(6)+60,61.*y(7)+60,61.*y(8)];

а вот и ее вызов

Matlab M
1
2
3
4
5
>> clear all;
>> t=[0 5];
>> y0=[0, 0, 0, 0, 0, 0, 0, 0];
>> 
>> [T, Y]=ode45(@pr8, t, y0);

на это система ругается следующим

Matlab M
1
2
3
4
5
6
7
8
??? Error: File: C:\MATLAB6p5\work\pr8.m Line: 12 Column: 42
Missing operator, comma, or semicolon.
 
Error in ==> C:\MATLAB6p5\toolbox\matlab\funfun\private\odearguments.m
On line 104  ==> f0 = feval(ode,t0,y0,args{:});
 
Error in ==> C:\MATLAB6p5\toolbox\matlab\funfun\ode45.m
On line 155  ==> [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, args, ...

установлен MatLab 6.5, возможно по этой причине ошибка, а может и нет. В любом случае посмотрите пожалуйста код и помогите решить. Заранее благодарю!!!

@nullpointer · 14.12.2013, 00:45

М_Скофилд, вы пропустили какой-то оператор в вашей функции, но так как код не полностью представлени точнее сказать не могу

@М_Скофилд · 14.12.2013, 09:56 **[ТС]**

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

М_Скофилд, вы пропустили какой-то оператор в вашей функции, но так как код не полностью представлени точнее сказать не могу

Мне необходимо решить систему уравнений, которая приведена выше. И нужно решить используя функцию ode23 или ode45. Я уже сталкивался с несколькими видами синтаксисов этих функций и уже прилично запутался

@nullpointer · 14.12.2013, 11:23

М_Скофилд, я понимаю что вам нужно решить ту систему, но вы покажите полностью код функции pr8, потому что там в 12 строке ошибка, а я ее не вижу. Вам Matlab написал же "C:\MATLAB6p5\work\pr8.m Line: 12 Column: 42" - это значит что в файле pr8.m в 12 строке на 42 символе отсутствует оператор, запятая или точка с запятой.

@М_Скофилд · 14.12.2013, 13:00 **[ТС]**

в том и дело, что у меня после 10 строки ничего там нет) ни 11-й строчки, ни тем более 12-й

@Hofstadter · 14.12.2013, 13:06

Сообщение от М_Скофилд

вот код функции

жесть:
1. убери точки в строках 3-10.
2. замени запятые строках 3-10 на точки.
3. удали квадратную скобку в конце строки 10.
4. систему нужно править т.к. решение не сходится.

P.S. сорри ТС. М_Скофилд, темы свои делать нужно.

@М_Скофилд · 14.12.2013, 13:16 **[ТС]**

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

жесть:
1. убери точки в строках 3-10.
2. замени запятые строках 3-10 на точки.
3. удали квадратную скобку в конце строки 10.
4. систему нужно править т.к. решение не сходится.

P.S. сорри ТС. М_Скофилд, темы свои делать нужно.

спасибо за замечания и советы. Насчет того, что новую тему создавать не успел прочитать FAQ)
вечером опробую, надеюсь получится

@М_Скофилд · 14.12.2013, 21:41 **[ТС]**

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

производная пропорциональна времени, значит функция будет расти в квадрате

вот она и растет.
выбери диапазон поменьше.

я конечно прошу прощения, но не мог ли бы вы показать свой код, так как графики функций у нас с вами расходятся) все исправления, о которых вы говорили внес и получилось следующее

@Hofstadter · 15.12.2013, 01:44

все то же самое.

Зосима · 16.12.2013, 15:09

М_Скофилд,

Сообщение от М_Скофилд

не мог ли бы вы показать свой код

Matlab M
1
2
3
4
5
6
7
8
9
10
function dy=pr8(t,y)
dy=zeros(8,1);
dy(1) = -1000.*y(1) + 1000.*y(3) - 1000.*y(5) - 1000;
dy(2) = -333.33.*y(2) - 333.33.*y(3) +333.33.*y(4) + 333.33.*y(5);
dy(3) = 111.11.*y(1) - 200.*y(2) - 311.11.*y(3) + 200.*y(4) + 311.11.*y(5) + 111.11;
dy(4) = 142.86.*y(2) + 142.86.*y(3) - 157.14.*y(4) - 157.14.*y(5) + 14.29.*y(6) + 14.29.*y(7);
dy(5) = -61.73.*y(1) + 111.11.*y(2) + 172.84.*y(3) - 728.4.*y(4) - 200.3.*y(5) + 11.11.*y(6) + 27.45.*y(7) - 61.73;
dy(6) = -111.11.*y(4) - 111.11.*y(5) - 144.78.*y(6) - 144.78.*y(7) + 33.67.*y(8);
dy(7) = -142.86.*y(4) + 352.94.*y(5) - 186.15.*y(6) - 1110.52.*y(7) + 43.29.*y(8);
dy(8) = 60.61.*y(6) + 60.61.*y(7) + 60.61.*y(8);

*основная ошибка в том, что разделитель дробной части не запятая, а точка!

Проверь внимательно систему! Может где-то знак перепутал

@М_Скофилд · 16.12.2013, 15:13 **[ТС]**

Спасибо за отзывчивость, я кажется стану вашим форумчанином теперь))

Новые блоги и статьи Все статьи Все блоги /
Памятка для бота и "визитка" для читателей "Semantic Universe Layer (Слой семантической вселенной)" Hrethgir 19.04.2026 Сгенерировано для краткого описания по случаю сборки и компиляции скелета серверного приложения. И пусть после этого скажут, что статьи сгенерированные AI - туфта и не интересно. И это не реклама -. . .	Запрет удаления строк ТЧ документа при определенном условии Maks 19.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "Аккумуляторы", разработанного в конфигурации КА2. У данного документа есть ТЧ, в которой в зависимости от прав доступа. . .	Модель заражения группы наркоманов alhaos 17.04.2026 Условия задачи сформулированы тут Суть: - Группа наркоманов из 10 человек. - Только один инфицирован ВИЧ. - Колются одной иглой. - Колются раз в день. - Колются последовательно через. . .	Мысли в слух. Про "навсегда". kumehtar 16.04.2026 Подумалось тут, что наверное очень глупо использовать во всяких своих установках понятие "навсегда". Это очень сильное понятие, и я только начинаю понимать край его смысла, не смотря на то что давно. . .
My Business CRM MaGz GoLd 16.04.2026 Всем привет, недавно возникла потребность создать CRM, для личных нужд. Собственно программа предоставляет из себя базу данных клиентов, в которой можно фиксировать звонки, стадии сделки, а также. . .	Знаешь почему 90% людей редко бывают счастливыми? kumehtar 14.04.2026 Потому что они ждут. Ждут выходных, ждут отпуска, ждут удачного момента. . . а удачный момент так и не приходит.	Фиксация колонок в отчете СКД Maks 14.04.2026 Фиксация колонок в СКД отчета типа Таблица. Задача: зафиксировать три левых колонки в отчете. Процедура ПриКомпоновкеРезультата(ДокументРезультат, ДанныеРасшифровки, СтандартнаяОбработка) / / . . .	Настройки VS Code Loafer 13.04.2026 { "cmake. configureOnOpen": false, "diffEditor. ignoreTrimWhitespace": true, "editor. guides. bracketPairs": "active", "extensions. ignoreRecommendations": true, . . .

@nullpointer 48 / 48 / 14 Регистрация: 30.03.2009 Сообщений: 523
	14.12.2013, 00:45
	М_Скофилд, вы пропустили какой-то оператор в вашей функции, но так как код не полностью представлени точнее сказать не могу 0

@nullpointer 48 / 48 / 14 Регистрация: 30.03.2009 Сообщений: 523
	14.12.2013, 11:23
	М_Скофилд, я понимаю что вам нужно решить ту систему, но вы покажите полностью код функции pr8, потому что там в 12 строке ошибка, а я ее не вижу. Вам Matlab написал же "C:\MATLAB6p5\work\pr8.m Line: 12 Column: 42" - это значит что в файле pr8.m в 12 строке на 42 символе отсутствует оператор, запятая или точка с запятой. 0

@М_Скофилд 0 / 0 / 0 Регистрация: 14.12.2013 Сообщений: 6
	14.12.2013, 13:00 [ТС]
	в том и дело, что у меня после 10 строки ничего там нет) ни 11-й строчки, ни тем более 12-й 0

@Hofstadter 42 / 42 / 0 Регистрация: 27.10.2012 Сообщений: 228
	15.12.2013, 01:44
	все то же самое. 0

@М_Скофилд 0 / 0 / 0 Регистрация: 14.12.2013 Сообщений: 6
	16.12.2013, 15:13 [ТС]
	Спасибо за отзывчивость, я кажется стану вашим форумчанином теперь)) 0

Опции темы