Краевая задача. (bvp4c, bvpinit.)

@qviksy · Регистрация: 10.11.2014

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

Чем дальше в MatLab, тем больше непонятных ошибок.
программа не работает и не понятно - почему.

$https://www.cyberforum.ru/cgi-bin/latex.cgi?y''-{exp}^{x}y'=x;$
y(1.5) = 0.2;
y(2.5) = 2.5;

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
function example
solinit = bvpinit(linspace(1.5,2.5,10),[0.2 2.5]);
% ВЫЗЫВАЕМ СОЛВЕР BVP4C:
sol = bvp4c(@ode2, @bc2, solinit);
% CТРУКТУРА sol:  
% sol.x – ВЕКТОР, В КОТОРМ СОДЕРЖАТЬСЯ КООРДИНАТЫ УЗЛОВ СЕТКИ
% sol.y – МАТРИЦА, СОСТОЯЩАЯ ИЗ ДВУХ СТРОК,
% sol.y(1, :) – СООТВЕТСТВУЕТ ЗНАЧЕНИЯМ ФУНКЦИИ y1 
% sol.y(2, :) - СООТВЕТСТВУЕТ ЗНАЧЕНИЯМ ФУНКЦИИ y2
plot(sol.x, sol.y)
grid on
legend('y(t)','x(t)')
xlabel('t')
end
 
function dydx = ode2(x, y)
dydx = [y(2)
    exp(x).*y(1)+x];
%y1 = y
%y2 = y'
end
 
function res = bc2(ya, yb)
res = [ya(1.5)-0.2
       yb(2.5)-2.5];
end

Добавлено через 16 минут
Я нашла ошибку! Программа работает, но почему должны быть единицы, я не понимаю.

Matlab M
1
2
3
function res = bc2(ya, yb)
res = [ya(1)-0.2
       yb(1)-2.5];

Зосима · 08.12.2014, 09:13

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

почему должны быть единицы, я не понимаю.

гляди, в функции граничных условий ya, yb - это два массива (а не функции! поэтому запись типа ya(-3.24) сразу отпадает). длинна каждого из них равна порядку уравнения (или количеству переменных системы).
В данном случае это 2: ya(1) и yb(1) - описывают поведение y(x), а ya(2) и yb(2) - описывают поведение dy/dx.
Массив ya описывает поведение функций в начальный момент x=a, а массив yb описывает поведение функций в конечный момент x=b.

Ясненько?

@qviksy · 11.12.2014, 19:18 **[ТС]**

Ох, спасибо огромное. Я вроде поняла смысл того, что вы хотели донести до меня.

Новые блоги и статьи Все статьи Все блоги /
Символьное дифференцирование igorrr37 13.02.2026 / * Логарифм записывается как: (x-2)log(x^2+2) - означает логарифм (x^2+2) по основанию (x-2). Унарный минус обозначается как ! */ #include <iostream> #include <stack> #include <cctype>. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»
SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .	Установка Qt-версии Lazarus IDE в Debian Trixie Xfce volvo 10.02.2026 В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .

@qviksy 0 / 0 / 0 Регистрация: 10.11.2014 Сообщений: 7
	11.12.2014, 19:18 [ТС]
	Ох, спасибо огромное. Я вроде поняла смысл того, что вы хотели донести до меня. 0

Краевая задача. (bvp4c, bvpinit.)

Решение