Чи є спосіб отримати імпульсну відповідь дискретної системи, просто дізнавшись, що це реакція на крок функції дискретного блоку?


10

У безперервному часі це було можливо;

u(t)systemy(t)δ(t)=du(t)dtsystemdy(t)dt=h(t)

Чи застосовується те ж саме для дискретної системи часу, тобто

δ[t]=du[t]dtwhere:{δ[t]is the discrete time deltau[t]is the discrete time unit step function

Чи є спосіб отримати імпульсну відповідь дискретної системи, просто дізнавшись відповідь кроку дискретної одиниці?


1
Дивовижне питання! Ласкаво просимо на DSP.SE. Дотримуйтесь і сприяйте!
Phonon

Відповіді:


7

Більш простий варіант відповіді Фонона полягає в наступному.

Припустимо, що позначає реакцію системи на функцію одиничного кроку. Тоді, як обговорювалося у цій відповіді , узагалі - сума масштабованих та затримкуваних у часі копій імпульсної відповіді, і в цьому конкретному випадку не потрібно масштабування; лише затримки у часі. Таким чином, де кожен стовпець праворуч - імпульсна реакція із затримкою у часі та без затримки. Таким чином, ми легко отримуємо, що y y

y[0]=h[0]y[1]=h[1]+h[0]y[2]=h[2]+h[1]+h[0]y[3]=h[3]+h[2]+h[1]+h[0] = 
h[0]=y[0]h[1]=y[1]y[0]h[2]=y[2]y[1] = h[n] =y[n]y[n1] = 
з нарією згадка про фільтри, обертання, згортки, інтеграцію, оператори тощо, просто прості наслідки визначення лінійної інваріантної системи часу.

Ви чітко це робили довше, ніж у мене =)
Phonon

6

Так, це однаково справедливо і у випадку дискретних систем. Операція диференціювання в цьому випадку замінюється на різницю першого порядку. Не думаю, що він має універсальний символ, але назвемо це . Ця операція еквівалентна фільтруванню вашого сигналу за допомогою . Назвемо цей фільтр . Я позначаю згортання, що символ .D()y[n]=x[n]x[n1]d[n]

Давайте тепер застосуємо те, що ми знаємо про згортку до цього оператора. Ми знаємо, що отримуємо з поточною сумою (дискретний інтегратор) на . Фактично, система, представлена сама виявляється цим дискретним інтегратором. Також зауважте, що ці два оператори є оберненими один до одного, і що конкретно .u[n]δ[n]u[n]u[n]d[n]=δ[n]

Тепер ми знаємо, що згортка є комутативною, тобто

a[n]b[n]=b[n]a[n]

і асоціативний, тобто

(a[n]b[n])c[n]=a[n](b[n]c[n])

Отже,

x[n]=δ[n]x[n]=u[n]d[n]x[n]=d[n]u[n]x[n]=d[n](u[n]x[n])

Отже, ви бачите, що ви можете відновити з , застосувавши різницю першого порядку, як і в безперервному випадку.x[n](u[n]x[n])


2

Припущення:

  • Безперервна часова область : Нехай - імпульсна характеристика, а - крокова відповідьh(t)s(t)
  • Дискретна часова область : Нехай - одинична імпульсна відповідь, а - одинична крокова відповідьh[n]s[n]

Інтуїтивно кажучи, інтеграція в безперервну часову область еквівалентна підсумовуванню в дискретній часовій області. Аналогічно, похідна в безперервній часовій області еквівалентна кінцевій різниці в дискретному домені.

Маючи це інтуїтивне розуміння, розгляньте співвідношення між та (ліва частина другого рівняння у вашому дописі):uδ

  • Постійна область часу: u(t)=δ(t)
  • Дискретна часова область: u[n]=k=0δ[nk]

Аналогічно розглянемо співвідношення між та (права частина другого рівняння у твоєму дописі):sh

  • Постійна область часу: s(t)=h(t)
  • Дискретна часова область: s[n]=k=0h[nk]

Тепер, якщо ви уважно подивіться на останнє рівняння:

s[n]=k=0h[nk]

Тепер можна знайти з цього рівняння, використовуючи кінцеву різницю із запізнілою версією себе, тобто з :h[n]s[n]s[n1]

h[n]=s[n]s[n1]
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.