Нарисуйте оба сигнала и отметьте, где каждый из них отличен от нуля.
Отразите h и сдвигайте её: для каждого момента выхода найдите диапазон, где оба сигнала перекрываются.
Перемножайте и складывайте (дискретный случай) или перемножайте и интегрируйте (непрерывный случай) только на этом перекрытии.
Проверьте результат: выход начинается в сумме моментов начала, заканчивается в сумме моментов окончания, а его площадь или сумма равна произведению площадей или сумм входов.
2Тренировочные задачи
1.Свёртка и линейные стационарные системыСредняя
Импульс x(t) = 2 при 0 ≤ t < 1 s (в остальное время нуль) подан на систему с импульсной характеристикой h(t) = 1 при 0 ≤ t < 3 s (в остальное время нуль). Найдите выход y(3,5) при t = 3,5 s.
x(t) = 2, 0 ≤ t < 1 sh(t) = 1, 0 ≤ t < 3 st = 3,5 s
y(3.5) =
2.Свёртка и линейные стационарные системыСложная
Система подчиняется уравнению y[n] = 0,4·y[n − 1] + 2·x[n] и начинает с состояния покоя. На входе единичная ступенька x[n] = u[n]. Найдите y[2].
y[n] = 0,4y[n − 1] + 2x[n]x[n] = u[n]
y[2] =
3.Преобразования ЛапласаСредняя
Найдите f(t) для F(s) = (6s + 9)/(s² + 8s + 16) и вычислите f(0,25).
F(s) = (6s + 9)/(s² + 8s + 16)t = 0,25 s
f(0.25) =
4.Преобразования ЛапласаСложная
Система описывается дифференциальным уравнением y″ + 7y′ + 10y = 4u′ + 20u, где u(t) — вход, y(t) — выход. Найдите её передаточную функцию G(s) = Y(s)/U(s) (при нулевых начальных условиях) и вычислите импульсную характеристику g(t) при t = 0,75 s. Найдите g(0,75).
y″ + 7y′ + 10y = 4u′ + 20ut = 0,75 s
g(0.75) =
5.Переходная характеристикаСредняя
Система имеет передаточную функцию G(s) = 35,5/(s² + 6s + 71) и получает на вход единичную ступеньку. Полюсы равны s = σ ± jωd. Найдите вещественную часть σ.
G(s) = 35,5/(s² + 6s + 71)
σ =
6.Переходная характеристикаСложная
У контура с единичной отрицательной обратной связью G(s) = K/(s(s + 5)). Найдите коэффициент K, при котором переходная характеристика имеет перерегулирование 15 %.
G(s) = K/(s(s + 5))%OS = 15 %
K =
7.Диаграммы Боде и запасыСредняя
Найдите фазу ∠G(jω) функции G(s) = 27(s + 20)/(s(s + 6)) при ω = 15 rad/s в градусах.
G(s) = 27(s + 20)/(s(s + 6))ω = 15 rad/s
∠G(jω) = °
8.Диаграммы Боде и запасыСложная
Контур с единичной обратной связью имеет передаточную функцию разомкнутой системы L(s) = K/(s(s + 1)(s + 9)). Найдите коэффициент усиления K, который даёт запас по усилению 20 dB.
L(s) = K/(s(s + 1)(s + 9))GM = 20 dB
K =
3Растворы
1.Свёртка и линейные стационарные системы
Ответ: y(3.5) = 1
Непрерывная свёртка: y(t) = ∫ x(τ)·h(t − τ) dτ. Отразите h, сдвиньте в положение t и проинтегрируйте произведение перекрывающихся частей. Для причинных сигналов интеграл берётся от 0 до t. Полезные факты: площадь y равна площади x, умноженной на площадь h; если x длится от t₁ до t₂, а h — от t₃ до t₄, то y длится от t₁ + t₃ до t₂ + t₄. Для h(t) = e−atu(t) и ступеньки высотой A получаем y(t) = (A/a)(1 − e−at).
Перекрытие отражённого импульса
h(t − τ) — это импульс h, отражённый и сдвинутый так, что он заканчивается в точке τ = t.
τ ∈ [0,5, 1] → длина 0,5 s
Интеграл
y(3,5) = 2·1·0,5 = 1
Ответ
y(3,5) = 1
y(t) имеет форму трапеции: она растёт от 0, держится на уровне 2 и снова падает до 0 при t = 4 s.
2.Свёртка и линейные стационарные системы
Ответ: y[2] = 3,12
Импульсная характеристика и ряд Фурье: импульсная характеристика h(t) полностью описывает линейную стационарную систему: для y′ + a·y = b·x она равна h(t) = b·e−atu(t), а для y[n] = a·y[n − 1] + b·x[n] — h[n] = b·aⁿu[n]. Синусоида cos ωt проходит с коэффициентом передачи |H(jω)| и фазой ∠H(jω), где H — преобразование Фурье функции h. Периодический сигнал — это сумма гармоник на частотах k/T: меандр с амплитудой ±A содержит только нечётные гармоники с амплитудами 4A/(kπ), треугольная волна с пиком A — 8A/(π²k²), а последовательность импульсов высотой A и коэффициентом заполнения D имеет среднее значение A·D и гармоники (2A/(kπ))·|sin(kπD)|.
Итерации
y[0] = 2
y[1] = 2,8
y[2] = 3,12
Результат
y[2] = 2·(1 − (0,4)3)/(1 − 0,4) = 3,12
Ответ
y[2] = 3,12
3.Преобразования Лапласа
Ответ: f(0.25) = 0,8277
Преобразование Лапласа: F(s) = ∫₀∞ f(t)·e−st dt превращает производные в умножение на s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Стандартные пары: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), а сдвиг s → s + σ умножает функцию времени на e−σt. Чтобы найти оригинал рациональной F(s), разложите её на простые дроби: вычет простого полюса p равен A = [(s − p)·F(s)] при s = p. Теорема о конечном значении: f(∞) = lims→0 s·F(s), верна только когда все полюсы s·F(s) лежат в левой полуплоскости. Теорема о начальном значении: f(0⁺) = lims→∞ s·F(s).
Разложение знаменателя на множители
s² + 8s + 16 = (s + 4)²
Простые дроби
6s + 9 = 6(s + 4) + −15
F(s) = 6/(s + 4) − 15/(s + 4)²
Обратное преобразование
f(t) = 6·e−4t − 15·t·e−4t
f(0,25) = (6 − 15·0,25)·e−4·0,25 = 0,8277
Ответ
f(0,25) = 0,8277
4.Преобразования Лапласа
Ответ: g(0.75) = 0,8925
Преобразование Лапласа: F(s) = ∫₀∞ f(t)·e−st dt превращает производные в умножение на s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Стандартные пары: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), а сдвиг s → s + σ умножает функцию времени на e−σt. Чтобы найти оригинал рациональной F(s), разложите её на простые дроби: вычет простого полюса p равен A = [(s − p)·F(s)] при s = p. Теорема о конечном значении: f(∞) = lims→0 s·F(s), верна только когда все полюсы s·F(s) лежат в левой полуплоскости. Теорема о начальном значении: f(0⁺) = lims→∞ s·F(s).
Передаточная функция
G(s) = Y(s)/U(s) = (4s + 20)/(s² + 7s + 10)
Простые дроби
s² + 7s + 10 = (s + 2)(s + 5)
A = 4, B = 0
Импульсная характеристика
g(t) = 4·e−2t
g(0,75) = 0,8925
Ответ
g(0,75) = 0,8925
5.Переходная характеристика
Ответ: σ = −3
Переходная характеристика: система первого порядка G(s) = K/(τs + 1) отвечает на ступеньку величиной U сигналом y(t) = K·U·(1 − e−t/τ): за одну постоянную времени τ она достигает 63,2 % конечного значения K·U, а за 4τ — 98 %. Система второго порядка G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) имеет полюсы s = −ζωₙ ± jωd, где ωd = ωₙ√(1 − ζ²). При 0 < ζ < 1 переходная характеристика перерегулируется: перерегулирование %OS = 100·e−ζπ/√(1−ζ²), время пика tp = π/ωd, а время установления по критерию 2 % ts ≈ 4/(ζωₙ). Сравните знаменатель с s² + 2ζωₙs + ωₙ², чтобы найти ωₙ и ζ. Для системы масса–пружина–демпфер ωₙ = √(k/m) и ζ = c/(2√(km)); для последовательной RLC-цепи ωₙ = 1/√(LC) и ζ = (R/2)·√(C/L).
Сравните со стандартной формой
Знаменатель должен совпадать с s² + 2ζωₙs + ωₙ².
ωₙ² = 71 → ωₙ = 8,426 rad/s
2ζωₙ = 6 → ζ = 6/(2·8,426) = 0,356
Результат
σ = −ζωₙ = −6/2 = −3
Ответ
σ = −3
6.Переходная характеристика
Ответ: K = 23,39
Переходная характеристика: система первого порядка G(s) = K/(τs + 1) отвечает на ступеньку величиной U сигналом y(t) = K·U·(1 − e−t/τ): за одну постоянную времени τ она достигает 63,2 % конечного значения K·U, а за 4τ — 98 %. Система второго порядка G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) имеет полюсы s = −ζωₙ ± jωd, где ωd = ωₙ√(1 − ζ²). При 0 < ζ < 1 переходная характеристика перерегулируется: перерегулирование %OS = 100·e−ζπ/√(1−ζ²), время пика tp = π/ωd, а время установления по критерию 2 % ts ≈ 4/(ζωₙ). Сравните знаменатель с s² + 2ζωₙs + ωₙ², чтобы найти ωₙ и ζ. Для системы масса–пружина–демпфер ωₙ = √(k/m) и ζ = c/(2√(km)); для последовательной RLC-цепи ωₙ = 1/√(LC) и ζ = (R/2)·√(C/L).
Передаточная функция замкнутой системы
T(s) = K/(s² + 5s + K)
Коэффициент демпфирования по перерегулированию
ζ = −ln(0,15)/√(π² + ln²(0,15)) = 0,5169
Коэффициент усиления
ωₙ = 5/(2ζ) = 4,836 rad/s
K = ωₙ² = 23,39
Ответ
K = 23,39
7.Диаграммы Боде и запасы
Ответ: ∠G(jω) = −121,3 °
Частотная характеристика: подставьте s = jω в передаточную функцию. Модуль в децибелах равен 20·log₁₀|H(jω)|, а фаза — сумме углов множителей числителя минус сумма углов множителей знаменателя. У фильтра первого порядка RC или RL частота среза (−3 dB) равна ωc = 1/(RC) или R/L, и fc = ωc/(2π); фильтр нижних частот даёт |H| = 1/√(1 + (ω/ωc)²) и φ = −atan(ω/ωc). На асимптотической (кусочно-линейной) диаграмме Боде каждый полюс на своей частоте среза загибает модуль вниз на 20 dB/decade, каждый нуль загибает его вверх на 20 dB/decade, а интегратор 1/s даёт −20 dB/decade на всех частотах. Частота среза усиления ωc: |L(jω)| = 1 (0 dB); запас по фазе PM = 180° + ∠L(jωc). Частота перехода фазы ω180: ∠L(jω) = −180°; запас по усилению GM = −20·log₁₀|L(jω180)| dB. Замкнутая система устойчива, когда оба запаса положительны.
Частотная характеристика: подставьте s = jω в передаточную функцию. Модуль в децибелах равен 20·log₁₀|H(jω)|, а фаза — сумме углов множителей числителя минус сумма углов множителей знаменателя. У фильтра первого порядка RC или RL частота среза (−3 dB) равна ωc = 1/(RC) или R/L, и fc = ωc/(2π); фильтр нижних частот даёт |H| = 1/√(1 + (ω/ωc)²) и φ = −atan(ω/ωc). На асимптотической (кусочно-линейной) диаграмме Боде каждый полюс на своей частоте среза загибает модуль вниз на 20 dB/decade, каждый нуль загибает его вверх на 20 dB/decade, а интегратор 1/s даёт −20 dB/decade на всех частотах. Частота среза усиления ωc: |L(jω)| = 1 (0 dB); запас по фазе PM = 180° + ∠L(jωc). Частота перехода фазы ω180: ∠L(jω) = −180°; запас по усилению GM = −20·log₁₀|L(jω180)| dB. Замкнутая система устойчива, когда оба запаса положительны.