Решение системы ДУ

@maest · Регистрация: 11.11.2010

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

Добрый день, решаю систему ДУ:

Matlab M
1
2
3
4
5
function dydt = RLC(t,y,E,R,R1,R2,Rm,L,L1,L2,Lm,C)
dydt =[ ( -y(1)*R1+Lm*( ( (y(1)-y(2))*Rm-y(2)*R2-y(4) )/(L2+Lm) )-(y(1)-y(2))*Rm+E )/( (L1+Lm)*( 1+Lm^2/(L2+Lm) ) );
        ( (y(1)-y(2))*Rm+Lm*( (-(y(1)-y(2))*Rm-y(1)*R1+E )/(L2+Lm) )-y(2)*R2-y(4) )/( (L2+Lm)*( 1+Lm^2/(L1+Lm) ) );
        y(4)/L-y(3)*R/L;
        y(2)/C-y(3)/C ];

Matlab M
1
2
3
4
tspan = [ 0 1 ];
y0 = [ 0; 0; 0; 0 ];
[ t, y ] = ode45(@RLC, tspan, y0,[],E,R,R1,R2,Rm,L,L1,L2,Lm,C);
plot(t,y(:,4));grid;

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

Matlab M
1
2
syms T;
E = 220*sin(w*T);

@SSC · 07.12.2016, 14:38

Вы смешиваете численные расчеты с символьными, а это делать нельзя.
попробуйте

Matlab M
1
2
3
4
5
6
function dydt = RLC(t,y,E,R,R1,R2,Rm,L,L1,L2,Lm,C)
Em = E*sin(2*pi*50*t);
dydt =[ ( -y(1)*R1+Lm*( ( (y(1)-y(2))*Rm-y(2)*R2-y(4) )/(L2+Lm) )-(y(1)-y(2))*Rm+Em )/( (L1+Lm)*( 1+Lm^2/(L2+Lm) ) );
        ( (y(1)-y(2))*Rm+Lm*( (-(y(1)-y(2))*Rm-y(1)*R1+Em )/(L2+Lm) )-y(2)*R2-y(4) )/( (L2+Lm)*( 1+Lm^2/(L1+Lm) ) );
        y(4)/L-y(3)*R/L;
        y(2)/C-y(3)/C ];

Новые блоги и статьи Все статьи Все блоги /
10 пpимет, которые всегда сбываются Maks 31.03.2026 1. Чтобы, наконец, пришла маршрутка, надо закурить. Если сигарета последняя, маршрутка придет еще до второй затяжки даже вопреки расписанию. 2. Нaдоели зима и снег? Не надо переезжать. Достаточно. . .	Перемещение выделенных строк ТЧ из одного документа в другой Maks 31.03.2026 Реализация из решения ниже выполнена на примере нетипового документа "ВыдачаОборудованияНаСпецтехнику" с единственной табличной частью "ОборудованиеИКомплектующие" разработанного в конфигурации КА2. . . .	Functional First Web Framework Suave DevAlt 30.03.2026 Sauve. IO Апнулись до NET10. Из зависимостей один пакет, работает одинаково хорошо как в режиме проекта так и в интерактивном режиме. из сложностей - чисто функциональный подход. Решил. . .	Автоматическое создание документа при проведении другого документа Maks 29.03.2026 Реализация из решения ниже выполнена на нетиповых документах, разработанных в конфигурации КА2. Есть нетиповой документ "ЗаявкаНаРемонтСпецтехники" и нетиповой документ "ПланированиеСпецтехники". В. . .
Настройка движения справочника по регистру сведений Maks 29.03.2026 Решение ниже реализовано на примере нетипового справочника "ТарифыМобильнойСвязи" разработанного в конфигурации КА2, с целью учета корпоративной мобильной связи в коммерческом предприятии. . . .	Автозаполнение реквизита при выборе элемента справочника Maks 27.03.2026 Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. При выборе "Спецтехники" (Тип Справочник. Спецтехника), заполняется. . .	Сумматор с применением элементов трёх состояний. Hrethgir 26.03.2026 Тут. https:/ / fips. ru/ EGD/ ab3c85c8-836d-4866-871b-c2f0c5d77fbc Первый документ красиво выглядит, но без схемы. Это конечно не даёт никаких плюсов автору, но тем не менее. . . всё может быть. . .	Автозаполнение реквизитов при создании документа Maks 26.03.2026 Программный код из решения ниже размещается в модуле объекта документа, в процедуре "ПриСозданииНаСервере". Алгоритм проверки заполнения реализован для исключения перезаписи значения реквизита,. . .