deval

Оцените структуру решения для дифференциального уравнения

Описание

пример

y = deval(sol,x) и y = deval(x,sol) оцените решение sol из проблемы дифференциального уравнения в точках, содержавшихся в x.

пример

y = deval(___,idx) возвращает только компоненты решения с индексами, перечисленными в векторном idx. Можно использовать любую из ранее перечисленных комбинаций входных аргументов.

пример

[y,yp] = deval(___) также возвращает yp, который является первой производной числового решения, произведенного решателем.

Примеры

свернуть все

Этот пример решает уравнение DDE y' = ddex1de(t,y,Z) использование dde23, затем строит решение.

Решите систему с помощью dde23.

sol = dde23(@ddex1de, [1 0.2], @ddex1hist, [0 5]);

Оцените решение в 100 точках в интервале [0 5].

x = linspace(0,5);
y = deval(sol,x);

Постройте решение.

plot(x,y)

Этот пример решает систему y' = vdp1(t,y) использование ode45, затем строит первый компонент решения.

Решите систему с помощью ode45.

sol = ode45(@vdp1, [0 20], [2 0]);

Оцените первый компонент решения в 100 точках в интервале [0 20].

x = linspace(0,20,100);
y = deval(sol,x,1);

Постройте решение.

plot(x,y)

Решите простой ОДУ y' = t^2 с начальным условием y0 = 0 в интервале [0,3] использование ode23.

sol = ode23(@(t,y) t^2, [0 3], 0);

Оцените решение в семи точках. Структура решения sol содержит функцию интерполяции что deval использование, чтобы произвести непрерывное решение в этих точках. Задайте второй выходной аргумент с deval чтобы также возвратить производную интерполяции функционируют в заданных точках.

x = linspace(0,3,7);
[y,yp] = deval(sol,x)
y = 1×7

         0    0.0417    0.3333    1.1250    2.6667    5.2083    9.0000

yp = 1×7

         0    0.2500    1.0000    2.2500    4.0000    6.2500    9.0000

Входные параметры

свернуть все

Структура решения в виде структуры, возвращенной одним из этих решателей для дифференциальных уравнений.

  • Решатель задачи с начальными значениями — ode45, ode23, ode113, ode15s, ode23s, ode23t, ode23tb, ode15i

  • Решатель дифференциальных уравнений с запаздывающим аргументом — dde23, ddesd, или ddensd

  • Решатель для краевой задачи — bvp4c или bvp5c

Пример: sol = ode45(@myode,tspan,y0)

Типы данных: struct

Оценка указывает в виде вектора. x задает точки, в которых вы хотите значение решения. Элементы x должен содержаться в исходном интервале интегрирования, [sol.x(1) sol.x(end)]. Для каждого индекса i, решение y(:,i) соответствует x(i).

Пример: 0:0.1:1

Пример: [2 3 5 8]

Типы данных: single | double

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

Пример: y = deval(sol,x,[2 3]) возвращает только вторые и третьи компоненты решения.

Типы данных: single | double

Выходные аргументы

свернуть все

Интерполированное решение, возвращенное как вектор или матрица. Количество строк в y равно количеству возвращаемых компонентов решения.

Для многоточечных краевых задач, решение, полученное bvp4c или bvp5c может быть прерывистым в интерфейсах. Для интерфейсной точки xc, deval функция возвращает среднее значение пределов слева и права на xc. Чтобы получить предельные значения, установите значение x быть немного больше или меньшим, чем xc.

Производная непрерывного решения производится sol, возвращенный как вектор или матрица. yp одного размера с y и указывает на наклон функции интерполяции, используемой sol в каждом решении указывают в y.

Смотрите также

| | | | | | | | | | | |

Представлено до R2006a
Для просмотра документации необходимо авторизоваться на сайте