
- •Практическая работа в matlab
- •1.1. Запуск matlab
- •1.2. Простейшие арифметические действия
- •Задание 1
- •1.3. Формирование одномерного и двумерного массивов
- •Начальное_значение : Шаг : Конечное_значение
- •1.4. Построение и форматирование графиков функций
- •1.4.1. Двумерные графики функций
- •Задание 2
- •1.4.2.Трехмерные графики функций
- •1.5. Символьные вычисления
- •1.5.1. Классический метод решения дифференциальных уравнений
- •Задание 4
- •1.5.2. Создание символьных переменных, выражений, матриц
- •1.5.8. Вычисление сумм рядов
- •1.5.9. Символьное дифференцирование
- •1.5.10. Символьное интегрирование
- •Задание 5
Задание 2
Согласно выданному варианту и, используя соотношения y3, y4, …, y14, построить и отформатировать графики функций принадлежности, приведенных ниже.
при a = 2, d1 = 4, d2 = 10, d3 = 16, b = 20.
при a = 2, d1 = 4, d2 = 10, d3 = 16, b = 20.
п
ри
a
= 2, d1
= 4, d2
= 10, d3
= 16, b' =
20, d4
= 22, b"
= 26.
при a = 2, d1 = 4, d2 = 10, d3 = 16, d4 = 20, b = 24.
при a = 2, d1 = 4, d2 = 10, d3 = 16, d4 = 20, b = 24.
при a = 2, d1 = 4, d2 = 10, d3 = 14, d4 = 18, b' = 22, d5 = 24, d6 = 28, b" = 30.
при a = 2, d1 = 4, d2 = 10, d3 = 16, b = 18.
при a = 2, d1 = 4, d2 = 10, d3 = 16, b = 18.
при a = 2, d1 = 4, d2 = 10, d3 = 16, b' = 18, d4 =22, b" = 24.
при a = 1, d1 = 2, d2 = 4, d3 = 6, d4 = 8, d5 = 10, d6 = 12, b = 14.
при a = 1, d1 = 2, d2 = 4, d3 = 6, d4 = 8, d5 = 10, d6 = 12, b = 14.
при a = 1, d1 = 2, d2 = 4, d3 = 6, d4 = 8, d5 = 10, d6 = 12, b' = 14, d7 = 14, d8 = 16,
b" = 18.
при d1 = 4.5, d2 = 5, d3 = 5.5, Δd = 0.5, а = 3.
п
ри
d1
= 4.5, d2
= 5, d3
= 5.5, Δd
= 0.5, b
= 8
при d1 = 4.5, d2 = 5, d3 = 5.5 , d4 = 5.3, d5 = 5.8, , d6 = 6.3, Δd = 0.5.
при a = 2, d1 = 4, d2 = 10, d3 = 16, b' = 20, d4 =22, b" = 26.
1.4.2.Трехмерные графики функций
1.5. Символьные вычисления
Для выполнения символьных вычислений предназначен Symbolic Math Toolbox. Ниже приводятся примеры выполнения символьных вычислений.
1.5.1. Классический метод решения дифференциальных уравнений
Классический метод решения будем применять к неоднородному линейному дифференциальному уравнению
,
приведенному к виду
путем деления на a2 левой и правой части исходного дифференциального уравнения. Тогда a = a0/a2; b = a0/a2; с = b0/a2. Рассмотрим пример решения приведенного неоднородного дифференциального уравнения, у которого корни характеристического уравнения являются отрицательными и действительными.
Пример 1.5.1. Найти решение неоднородного линейного дифференциального уравнения
(1.2)
при
нулевых начальных условиях y(0)
= 0,
(0)
= 0 и единичном ступенчатом входном
сигнале u(t)
= 1(t).
Решение. Решение неоднородного дифференциального уравнения всегда равно сумме y(t) = yо(t) + yч(t), где yо(t) – общее решение соответствующего однородного дифференциального уравнения, полученного приравниванием нулю правой части (1.2)
,
(1.3)
а yч(t) – частное решение, определяемое видом правой части (1.2). Решение однородного дифференциального уравнения (1.3) начинается с задания соответствующего характеристического уравнения
p2 + 8p + 12 = 0
и определения его корней
p1,2
= –4 ±
,
p1 = – 4 + 2 = –2; p2 = – 4 – 2 = – 6.
Общее решение однородного дифференциального уравнения запишется так
,
где С1, С2 – постоянные интегрирования.
Перейдем к определению частного решения yч. Правая часть уравнения (1.2) – константа. Частное решение ищется в форме yч = B, где B – константа. Для нахождения B подставим yч = B в уравнение (1.2)
12B = 48, B = 4.
Таким образом, полное решение неоднородного уравнения (1.2) имеет вид
.
(1.4)
Исходя из нулевых начальных условий y(0) = 0, (0) = 0, запишем (1.4) в виде двух уравнений
y(0) = С1 + С2 + 4 = 0,
(0) = –2С1 – 6С2 = 0
и решим их относительно С1 и С2
С1 = – 6, С2 = 2.
Подставляя найденные значения С1 и С2 в (1.4), получим искомое решение
.
(1.5)
Найдите в MATLAB и сравните с (1.5) решение неоднородного дифференциального уравнения, используя оператор dsolve. Наберите программу решения дифференциального уравнения (1.2) с нулевыми начальными условиями.
>> y=dsolve('D2y+8*Dy+12*y=48','y(0)=0','Dy(0)=0')
y = 2*exp(-6*t)-6*exp(-2*t)+4
Последнее выражение и будет решением линейного неоднородного дифференциального уравнения (1.2).
Постройте график переходного процесса при u(t) = 1(t) по программе
>> t = 0:0.1:10;
>> y = 2*exp(-6*t)-6*exp(-2*t)+4;
>> plot(t,y,'k'), grid on