6 svar
132 visningar
naytte behöver inte mer hjälp
naytte 8271 – Moderator
Postad: 25 aug 23:11 Redigerad: 25 aug 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 759
Postad: 26 aug 00:48 Redigerad: 26 aug 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 8271 – Moderator
Postad: 26 aug 01:01 Redigerad: 26 aug 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 759
Postad: 26 aug 07:11 Redigerad: 26 aug 07:15
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. det är samma cc och mm för hela talföljden oavsett kk. (Notera dock att mm inte är givet entydigt eftersom man alltid kan välja ett större värde på mm.)

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 med |ak|c·mk|a_k| \le c\cdot m^k, 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 \dfrac{m}{|z|}. Då får du att Q<1Q<1 ifall |z|>m|z| > m.

naytte 8271 – Moderator
Postad: 26 aug 15:21

Ah okej, då är jag med!

Men vad är det examinatorn gör för att inverstransformera? Jag hänger itne riktigt medp å vad som händer i facit...

LuMa07 759
Postad: 27 aug 11:40 Redigerad: 27 aug 12:41
naytte skrev:

Men vad är det examinatorn gör för att inverstransformera? Jag hänger itne riktigt medp å vad som händer i facit...

Hon höll nog på med att vända och vrida på funktionen tills hon lyckades med att få ett enklare uttryck som motsvarade en geometrisk summa.

Det är följande steg hon tagit, fast i omvänd ordning:

k=0z-k= \displaystyle \sum_{k=0}^\infty z^{-k}= {}[geom. summa]=11-z-1=zz-1{} = \dfrac{1}{1-z^{-1}} = \dfrac{z}{z-1} .

Man kan sedan notera att

(zz-1)'=-1(z-1)2\displaystyle (\frac{z}{z-1})^\prime = \frac{-1}{(z-1)^2}

Om den geometriska serien deriveras termvis, så får man alltså att

k=1-kz-k-1=(zz-1)'=-1(z-1)2 \displaystyle \sum_{k=1}^\infty -k z^{-k-1} = (\frac{z}{z-1})^\prime = \frac{-1}{(z-1)^2} .

När man nu multiplicerar båda leden med -2z-2z, så får man att

k=12kz-k=2z(z-1)2 \displaystyle \sum_{k=1}^\infty 2k z^{-k} = \frac{2z}{(z-1)^2} .

 

Det är en sådan lösning som man återskapar när man redan vet vad som är svaret.

naytte 8271 – Moderator
Postad: 27 aug 15:14

Ja okej.

Jag tror att strategin på tentamen då helt enkelt är att be en bön om att man kan manipulera serien till något man känner till. Om det är så som examinatorn gör ska väl det duga.

Tack för hjälpen!

Svara
Close