Решить систему диф.уравнений первого порядка в аналитическом виде и вывести графики найденых функций

@Wlk · Регистрация: 07.12.2011

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

Собсно сабж. Нужно решить систему диф.уравнений первого порядка в аналитическом виде (в символьном) и вывести графики найденых функций.

Matlab M
1
2
3
4
5
6
7
clear all
L=1e-2;
M=1;
T=50000;
t=0:T/100:T;
[P1, P2, P3]= dsolve('DP1=(-L)*P1+M*P2', 'DP2=L*P1-M*P2-L*P2', 'DP3=L*P2','P1(0)=1','P2(0)=0','P3(0)=0','t');
plot(t,P1+P2);

В итоге, символьно решается система диф. уравнений (правильно или нет, не знаю, так как не могу проверить, из-за того, что не могу увидеть график). Значения L и M в выражения для P1,P2,P3 не подставляются. Мне нужны будут как выражения в символьном виде так и уже с подставлеными значениями L и M. Пробовал играться с eval(), не дало эффекта (наверно, не правильно использую, пробовал применять для параметров L и М (в виде eval('L=1e-2')), так же для Р1 и Р2 (в виде eval(Р1) или eval('Р1') или Р1=eval('Р1') или Р01=eval('Р1') или eval(Р01='P1') или eval('Р01=Р1'); не дало эффекта)). Далее при попытке вывода графика коммандой plot(t,P1); получаю следующюю ошибку.

Error using plot
Conversion to double from sym is not possible.

Error in KCE (line 7)
plot(t,P1)

Пробовал играться с ezplot(P1) - в некоторых случаях строило графики, но не подходящие для работы, так как нигде не указано какую переменную по абсциссе ставить. А в других случаях просто выдавало ошибку.

Error using ezplot>ezparam (line 393)
Cannot plot parametrized surfaces. Try ezsurf.

Error in ezplot (line 164)
hp = ezparam(cax, f{1}, f{2}, vars, labels, args{2:end});

Error in sym/ezplot (line 50)
h = ezplot(char(f),char(y));

Error in KCE (line 7)
ezplot(t,P1+P2);

Кажется все cвязано с переводом символьных даных в чисельные. Кто может помочь?

@Wlk · 26.01.2012, 16:34 **[ТС]**

Получилось подставить чисельные данные в символьное выражение, кроме независимой переменной используя

Matlab M
1
2
3
[P1, P2, P3]= dsolve('DP1=(-L)*P1+M*P2', 'DP2=L*P1-M*P2-L*P2', 'DP3=L*P2','P1(0)=1','P2(0)=0','P3(0)=0','t');
P=[P1 P2 P3];
Pt=subs(P,[L M],[1e-2 1]);

Остался вопрос как построить график. Ezplot() строит и показывает только характер кривой, мне же нужны значения фунцкии. Сейчас програма выглядит так:

Matlab M
1
2
3
4
5
6
7
8
clear all
syms L M
T=50000;
t=0:T/100:T;
[P1, P2, P3]= dsolve('DP1=(-L)*P1+M*P2', 'DP2=L*P1-M*P2-L*P2', 'DP3=L*P2','P1(0)=1','P2(0)=0','P3(0)=0','t');
P=[P1 P2 P3];
Pt=subs(P,[L M],[1e-2 1]);
ezplot(Pt(1));

По-прежнему не знаю как построить график подставив в символьные выражения для решений диф. уравнений, которые хранятся в переменных Pt(1),Pt(2),Pt(3), значения t.

Добавлено через 20 часов 56 минут
Добился некоего эффекта, но строит только первые дветочки графика, в чем пробелма может быть? Матлаб 2011б

@Wlk · 10.02.2012, 18:46 **[ТС]**

Решено. Проблема была в матлабе. В 2010б все строит нормально.

Новые блоги и статьи Все статьи Все блоги /
Инструменты COM: Сохранение данный из VARIANT в файл и загрузка из файла в VARIANT bedvit 28.01.2026 Сохранение базовых типов COM и массивов (одномерных или двухмерных) любой вложенности (деревья) в файл, с возможностью выбора алгоритмов сжатия и шифрования. Часть библиотеки BedvitCOM Использованы. . .	Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 28.01.2026 Содержание блога SDL3 имеет собственные средства для загрузки и отображения PNG-файлов с альфа-каналом и базовой работы с ними. В этой инструкции используется функция SDL_LoadPNG(), которая. . .	Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL3_image 8Observer8 27.01.2026 Содержание блога SDL3_image - это библиотека для загрузки и работы с изображениями. Эта пошаговая инструкция покажет, как загрузить и вывести на экран смартфона картинку с альфа-каналом, то есть с. . .	Влияние грибов на сукцессию anaschu 26.01.2026 Бифуркационные изменения массы гриба происходят тогда, когда мы уменьшаем массу компоста в 10 раз, а скорость прироста биомассы уменьшаем в три раза. Скорость прироста биомассы может уменьшаться за. . .
Воспроизведение звукового файла с помощью SDL3_mixer при касании экрана Android 8Observer8 26.01.2026 Содержание блога SDL3_mixer - это библиотека я для воспроизведения аудио. В отличие от инструкции по добавлению текста код по проигрыванию звука уже содержится в шаблоне примера. Нужно только. . .	Установка Android SDK, NDK, JDK, CMake и т.д. 8Observer8 25.01.2026 Содержание блога Перейдите по ссылке: https:/ / developer. android. com/ studio и в самом низу страницы кликните по архиву "commandlinetools-win-xxxxxx_latest. zip" Извлеките архив и вы увидите. . .	Вывод текста со шрифтом TTF на Android с помощью библиотеки SDL3_ttf 8Observer8 25.01.2026 Содержание блога Если у вас не установлены Android SDK, NDK, JDK, и т. д. то сделайте это по следующей инструкции: Установка Android SDK, NDK, JDK, CMake и т. д. Сборка примера Скачайте. . .	Использование SDL3-callbacks вместо функции main() на Android, Desktop и WebAssembly 8Observer8 24.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .

@Wlk 90 / 90 / 48 Регистрация: 07.12.2011 Сообщений: 215
	10.02.2012, 18:46 [ТС]
	Решено. Проблема была в матлабе. В 2010б все строит нормально. 0