3 svar
13 visningar
naytte 8247 – Moderator
Postad: Igår 23:11 Redigerad: Igår 23:13

Komplex analys - använda Z-transform för att lösa differensekvation

Hej!

Jag håller på att lära mig om Z-transformen för att lösa differensekvationer. Nedan följer en exempeluppgift som jag sitter med:

Här tänker jag att man kan se hela VL och HL som varsin följd med k0k\ge 0, så om vi tillämpar Z-transformen på båda led får vi:

Ƶ{ak+1-ak}k=0=Ƶ({2}k=0)\displaystyle \mathcal{Z}\left( \{ a_{k+1}-a_k \}_{k=0}^{\infty}\right)=\mathcal{Z}(\{ 2 \}_{k=0}^\infty)

Kalla A(z):=Ƶ({ak}k=0)A(z):=\mathcal{Z}( \{a_k \}_{k=0}^{\infty}). I så fall får vi efter förenkling att

zAz-za0-Az=2k0z-k\displaystyle zA\left(z\right)-za_0-A\left(z\right)=2\sum_{k\ge 0}z^{-k}

Nu kommer min första fråga. Jag känner igen att HL är en geometrisk summa som endast konvergerar då |z|>1|z|>1. Vet vi någonting om zz som låter oss göra den omskrivningen? Jag vet inte mycket mer om Z-transformen än definitionen.

Låt säga att vi sedan tillåter oss själva att skriva om serien som en sluten formel och stoppar in a0=0a_0=0. Då får vi:

Az=2z(z-1)2\displaystyle A\left(z\right)=\frac{2z}{(z-1)^2}

Nu skulle det ju passa bra att inverstransformera. Men jag vet inte riktigt hur det fungerar. Jag hittar följande formel online

akk=0=Ƶ-1Az=12πiC2z(z-1)2zk-1dz\displaystyle \left\{ a_k \right\}_{k=0}^\infty=\mathcal{Z}^{-1}\left(A\left(z\right)\right)=\frac{1}{2\pi i}\oint_{C}\frac{2z}{(z-1)^2}z^{k-1}dz

Där CC är någon speciell kurva. Jag är ganska lost på detta. Dokumentet som det står att vi ska läsa i på Canvas ligger inte på Canvas längre så... ja det är lite problematiskt det här. Jag skulle behöva hjälp med denna uppgift, men mina konkreta frågor är också:

  • Vad ställer vi för krav på zz i transformen? Kan det vara vilket komplext tal som helst?
  • Vad ställer vi för krav på kurvan CC då vi söker inverstransformen?
  • Kan man inverstransformera utan att integrera?
LuMa07 Online 755
Postad: Idag 00:48 Redigerad: Idag 00:53

I kursen i Transformteori respektive Signaler och system (och dylika kurser) skulle jag nog säga att man aldrig inverstransformerar m.h.a. integralen. Just nu kommer jag på tre grundläggande sätt för att invertera z-transformen:

  • Tabeller med räknelagar för z-transformen och standardtransformerna
  • Lång division ifall A(z) är ett rationellt uttryck
  • Vanliga maclaurinutvecklingar

Funktionen A(z)A(z) som du fått finns bland standardtransformerna och den motsvarar talföljden ak=2ka_k=2k, kk\in\mathbb{N}.

Du ställer dock frågan inom kursen i komplex analys och därmed är jag osäker om ni kanske vill räkna ut integralen.


När talföljden är av exponentiell ordning, d.v.s. det finns c,m>0c,m>0 sådana att akc·mka_k \le c\cdot m^k för varje kk \in \mathbb{N}, så kommer z-transformen konvergera utanför en cirkel av radien mm. När man arbetar med zz-transformen, så förutsätter man att |z|>R|z| > R för något R0R \ge 0.

Om du börjar med talföljden, så bestäms RR i enlighet med talföljdens exponentiella ordning.

Om du börjar med en funktion i zz, så bestäms RR så att den stämmer överens med den singularitet som ligger längst bort från origo.


Om vi tittar på ditt A(z)A(z) så kan den inverteras m.h.a. binomialutvecklingen (d.v.s. maclaurinserien för w(1+w)αw\mapsto (1+w)^\alpha):

A(z)=2z(z-1)2=2z·1(1-z-1)2=2z-1(1-z-1)-2=2z-1·k=0-2k(-z)-k\displaystyle {A(z)} = \dfrac{2z}{(z-1)^2} = \dfrac{2}{z} \cdot \dfrac{1}{(1-z^{-1})^2} = 2\,z^{-1} {(1-z^{-1})}^{-2} = 2z^{-1} \cdot \sum_{k=0}^\infty \binom{-2}{k} {(-z)}^{-k}

=k=02·-2k(-1)kz-k-1=k=12·-2k-1(-1)k-1z-k\displaystyle = \sum_{k=0}^\infty 2\cdot \binom{-2}{k} {(-1)}^k z^{-k-1} = \sum_{k=1}^\infty 2 \cdot\binom{-2}{k-1} {(-1)}^{k-1} z^{-k}

Därmed är a0=0a_0 = 0, och sedan ak=2·-2k-1(-1)k-1\displaystyle a_k = 2\cdot \binom{-2}{k-1} {(-1)}^{k-1}k1k\ge 1. Det återstår att förenkla binomialkoefficienten.


Tillägg: 26 aug 2026 06:57

Jag tappade bort absolutbeloppet när jag skrev om exponentiell ordning. Villkoret borde lyda |ak|c·mk|a_k| \le c\cdot m^k

naytte 8247 – Moderator
Postad: Idag 01:01 Redigerad: Idag 01:02

Hej, tack för svar!

Vi ska nog inte faktiskt behöva sitta med några integraler. Nedan följer facit till uppgiften, men jag förstår inte riktigt vad som händer. Men jag tror allmänt att det är tänkt att man ska använda listiga "trick" för att inverstransformera igen. Vi får inga formelblad så jag tror inte att man ska behöva kunna "tabelltransformer".

Dessutom: vad gäller att RR bestäms av talföljdens exponentiella ordning, hur går den att se i vårt fall? Och vad gäller c,mc,m, är dessa fixerade eller beror de på kk\in \mathbb{N}?

LuMa07 Online 755
Postad: Idag 07:11 Redigerad: Idag 07:13
naytte skrev:

...

Dessutom: vad gäller att RR bestäms av talföljdens exponentiella ordning, hur går den att se i vårt fall? Och vad gäller c,mc,m, är dessa fixerade eller beror de på kk\in \mathbb{N}?

Konstanterna cc och mm gäller likformigt, d.v.s. hela talföljden oavsett kk.

Eftersom uttrycket A(z)=2z(z-1)2A{(z)} = \dfrac{2z}{(z-1)^2} har en pol i punkten z=1z=1 vars avstånd till origo är R:=1R := 1 och det finns ingen pol med längre avstånd till origo, så kommer den tillhörande talföljden uppfylla kravet |ak|c·mk|a_k| \le c\cdot m^k för varje m>Rm > R(=1) med ett lämpligt val på c=c(m)c=c(m). I allmänhet kan det hända att m=Rm=R(=1) också skulle funka, dock så är inte fallet här. Jag är ganska säker att m=Rm=R inte funkar ifall polen har högre multiplicitet. För tillfället är jag dock osäker hur det är med enkelpoler.

Om du har en talföljd, så bilda "potensserien" enligt definitionen av zz-transformen och undersök dess konvergens m.h.a. rotkriteriet:

Q=limsupk|akz-k|km/|z|Q = \limsup_{k\to\infty} \sqrt[k]{|a_k\,z^{-k}|} \le m/|z|. Då får du att Q<1Q<1 ifall |z|>m|z| > m.

Svara
Close