Дифференциальное уравнение численными методами

@lalumieredelune · Регистрация: 01.11.2015

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

Пробовала решить по первым условиям, но сомневаюсь

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
>> solinit = bvpinit(linspace(0,2),[0 0]);
>> sol = bvp4c(@myfun,@bc2,solinit);
>> plot(sol.x,sol.y(1,:)); 
hold on
 
function f = myfun( t, y )
f=[y(2);sin(t)-t*y(1)];
end
 
function res=bc2(ya,yb) 
res = [ya(2)
     yb(2)2];
end

Зосима · 17.11.2015, 16:45

а какие у тебя граничные условия?

@lalumieredelune · 17.11.2015, 20:49 **[ТС]**

Зосима, спасибо, что откликнулся!
нужно решить при условиях y(0)=y(2)=0 при $https://www.cyberforum.ru/cgi-bin/latex.cgi?t\in$ [0,2] и при условиях y(0)=1, y`(0)=0 и сравнить результаты

Зосима · 18.11.2015, 12:35

Ты почти угадала с граничными условиями

гляди, ya относится к начальному моменту времени (t=0), yb относится к конечному моменту времени (t=2), причем каждая переменная ya,yb в данной задаче содержит два элемента: первый описывает функцию, второй - первую производную

То есть:

ya(1) = y(0),
ya(2) = y'(0),
yb(1) = y(2),
yb(2) = y'(2)

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

)

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
function temp1
 
solinit = bvpinit(linspace(0,2),[0 0]);
 
% расчет с первыми граничными условиями:
sol = bvp4c(@myfun,@bc1,solinit);
plot(sol.x,sol.y(1,:),'b');
hold on
% расчет со вторыми граничными условиями:
sol = bvp4c(@myfun,@bc2,solinit);
plot(sol.x,sol.y(1,:),'r');
 
end
 
function f = myfun( t, y )
f = [y(2); sin(t)-t*y(1)];
end
 
function res=bc1(ya,yb)
res = [ya(1) % y(0)=0
    yb(1)]; % y(2)=0
end
 
function res=bc2(ya,yb)
res = [ya(1)-1 % y(0)=1
    ya(2)]; % y'(0)=0
end

Дифференциальное уравнение численными методами

@lalumieredelune · 18.11.2015, 22:38 **[ТС]**

Зосима, спасибо за ответ и подробное объяснение

Новые блоги и статьи Все статьи Все блоги /
Сумматор с применением элементов трёх состояний. 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, знаешь?. . Когда вечерние улицы становятся ночными, а ты не можешь уснуть. Ты идёшь в любимый старый бар, и бармен наливает тебе виски. Ты смотришь на пролетающие. . .	SDL3 для Desktop (MinGW): Вывод текста со шрифтом TTF с помощью библиотеки SDL3_ttf на Си и C++ 8Observer8 24.03.2026 Содержание блога Финальные проекты на Си и на C++: finish-text-sdl3-c. zip finish-text-sdl3-cpp. zip

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

@lalumieredelune 0 / 0 / 0 Регистрация: 01.11.2015 Сообщений: 4
	17.11.2015, 20:49 [ТС]
	Зосима, спасибо, что откликнулся! нужно решить при условиях y(0)=y(2)=0 при $https://www.cyberforum.ru/cgi-bin/latex.cgi?t\in$ [0,2] и при условиях y(0)=1, y`(0)=0 и сравнить результаты 0

@lalumieredelune 0 / 0 / 0 Регистрация: 01.11.2015 Сообщений: 4
	18.11.2015, 22:38 [ТС]
	Зосима, спасибо за ответ и подробное объяснение 0