Ode45 и решение системы ДУ, ошибки

@Punktir · Регистрация: 01.12.2009

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

Здравствуйте!
Имеется система ДУ, которую нужно промоделировать в Матлабе, посмотреть график, но постоянно что-то да не так, застрял на моменте описания ДУ-ошибки и хоть убей. Поправьте, пожалуйста, если компетентны?
Прилагаю картинку - система ДУ в привычном виде. Прочел соответствующие темы на форуме, но эффект 0

Данные таковы:
Временные границы: 0-1000
Tist(0)=20;
Tg(0)=20;

Вот так описываю функцию:

Matlab M
1
2
3
4
5
6
function [ dT ] = difur( t, Tist, Tg )
 
dT(1) = (0-0.5*(Tist-Tg))/0.5;
dT(2) = (0.01*(Tist-Tg)-0.005*(Tg-20))/0.5;
 
end

А вот так пытаюсь увидеть изображение:

Matlab M
1
2
[T,Y] = ode45(@difur, [0,1000], [20 20]);
plot (T,Y);

@Punktir · 03.07.2015, 02:31 **[ТС]**

Чуть-чуть прорвался, но теперь встретил иную проблему.
Код теперь такой:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
function dt = enot( t,Tist,Tg )
Pel  = 0.1;
Gist = 0.5;
Cist = 0.5;
Gvn  = 0.005;
Tvn  = 20;
Cg   = 0.5;
 
dt   = zeros(2,1);
dt(1)= (Pel-Gist*(Tist(t)-Tg(t)))/Cist;
dt(2)= (Gist*(Tist(t)-Tg(t))-Gvn*(Tg(t)-Tvn))/Cg;
 
end

Вывожу так:

Matlab M
1
2
3
4
5
clear all, clc
T = [1 1000]; 
y0 = [20 20]; 
[t,Y] = ode45('enot',T,y0); 
 plot(t,Y);

Но, при это вылазиет ошибка:

Matlab M
1
2
3
4
Attempted to access Tg(1); index out of bounds because numel(Tg)=0.
 
Error in enot (line 10)
dt(1)= (Pel-Gist*(Tist(t)-Tg(t)))/Cist;

На форумах почитал, что дело в размерности вектора, но как бы я её не менял-все по-прежнему...Кто-нибудь в курсе что такое?

@SSC · 03.07.2015, 07:42

Вот код программы

Matlab M
1
2
3
4
%global k m1 x0 Lst Sst w mu0
 
[T,Y] = ode45(@diff_ur_sys, [0 1000], [25 28]);
plot (T,Y);

и функции

Matlab M
1
2
3
4
5
6
7
8
function dy = diff_ur_sys(t,y)
%global k m1 x0 Lst Sst w mu0
% y(1) - Tist
% y(2) - Tg
dy = zeros(2,1);    % a column vector
dy(1)=(0-0.5*(y(1)-y(2)))/0.5;
dy(2)=(0.01*(y(1)-y(2))-0.005*(y(2)-20))/0.5;
end

Постоянные оформи через глобальные переменные (в заготовке закоментированы)
И начальные условия 20 20 ведут к нулевым производным поэтому в решении прямые линии, при других значениях более интересно

Новые блоги и статьи Все статьи Все блоги /
Автозаполнение реквизита при выборе элемента справочника Maks 27.03.2026 Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. При выборе "Спецтехники" (Тип Справочник. Спецтехника), заполняется. . .	Сумматор с применением элементов трёх состояний. Hrethgir 26.03.2026 Тут. https:/ / fips. ru/ EGD/ ab3c85c8-836d-4866-871b-c2f0c5d77fbc Первый документ красиво выглядит, но без схемы. Это конечно не даёт никаких плюсов автору, но тем не менее. . . всё может быть. . .	Автозаполнение реквизитов при создании документа Maks 26.03.2026 Программный код из решения ниже размещается в модуле объекта документа, в процедуре "ПриСозданииНаСервере". Алгоритм проверки заполнения реализован для исключения перезаписи значения реквизита,. . .	Команды формы и диалоговое окно Maks 26.03.2026 1. Команда формы "ЗаполнитьЗапчасти". Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. В качестве источника данных. . .
Кому нужен AOT? DevAlt 26.03.2026 Решил сделать простой ланчер Написал заготовку: dotnet new console --aot -o UrlHandler var items = args. Split(":"); var tag = items; var id = items; var executable = args;. . .	Отправка уведомления на почту при создании или изменении элементов справочника Maks 24.03.2026 Программная отправка письма электронной почты на примере типового справочника "Склады" в конфигурации БП3. Перед реализацией необходимо выполнить настройку системной учетной записи электронной. . .	модель ЗдравоСохранения 5. Меньше увольнений- больше дохода! anaschu 24.03.2026 Теперь система здравосохранения уменьшает количество увольнений. 9TO2GP2bpX4 a42b81fb172ffc12ca589c7898261ccb/ https:/ / rutube. ru/ video/ a42b81fb172ffc12ca589c7898261ccb/ Слева синяя линия -. . .	Midnight Chicago Blues kumehtar 24.03.2026 Такой Midnight Chicago Blues, знаешь?. . Когда вечерние улицы становятся ночными, а ты не можешь уснуть. Ты идёшь в любимый старый бар, и бармен наливает тебе виски. Ты смотришь на пролетающие. . .