Skip to main content

Difference equation with given input sequence

Short answer

Solve recursively: x(0)=0, x(1)=0.3679, x(2)=0.8463. For k>=3, x(k)=1.4284 - 5.4375(0.3679)^k, settling at 1.4284.

Consider the difference equation: x(k+2)-1.3679x(k+1)+0.3679x(k)=0.3679u(k+1)+0.2642u(k), where x(k) is the output and x(k)=0 for k≤ 0, where u(k) is the input and is given by u(k)=0, k<0; u(0)=1; u(1)=0.2142; u(2)=0.2142; u(k)=0, k=3,4,5,6,… Determine the output x(k).

The output is found by solving the equation recursively, then confirmed with a closed form from the z-transform.

Step 1: Initial values

x(k)=0x(k) = 0 for k≤0k \le 0, so x(0)=0x(0) = 0. Put k=−1k = -1 in the equation, with x(−1)=0x(-1) = 0, u(−1)=0u(-1) = 0, u(0)=1u(0) = 1:

x(1)−1.3679x(0)+0.3679x(−1)=0.3679u(0)+0.2642u(−1)  ⇒  x(1)=0.3679x(1) - 1.3679x(0) + 0.3679x(-1) = 0.3679u(0) + 0.2642u(-1) \;\Rightarrow\; x(1) = 0.3679

Step 2: Pulse transfer function

With zero initial conditions (the x(1)x(1) and u(0)u(0) terms cancel), taking z-transforms gives:

G(z)=X(z)U(z)=0.3679z+0.2642z2−1.3679z+0.3679=0.3679z+0.2642(z−1)(z−0.3679)G(z) = \frac{X(z)}{U(z)} = \frac{0.3679z + 0.2642}{z^2 - 1.3679z + 0.3679} = \frac{0.3679z + 0.2642}{(z-1)(z-0.3679)}

The input is a finite sequence:

U(z)=1+0.2142z−1+0.2142z−2=z2+0.2142z+0.2142z2U(z) = 1 + 0.2142z^{-1} + 0.2142z^{-2} = \frac{z^2 + 0.2142z + 0.2142}{z^2}

Step 3: Recursive solution

x(k+2)=1.3679x(k+1)−0.3679x(k)+0.3679u(k+1)+0.2642u(k)x(k+2) = 1.3679x(k+1) - 0.3679x(k) + 0.3679u(k+1) + 0.2642u(k)
kkCalculation of x(k+2)x(k+2)x(k+2)x(k+2)
01.3679(0.3679)−0+0.3679(0.2142)+0.2642(1)1.3679(0.3679) - 0 + 0.3679(0.2142) + 0.2642(1)0.8463
11.3679(0.8463)−0.3679(0.3679)+0.3679(0.2142)+0.2642(0.2142)1.3679(0.8463) - 0.3679(0.3679) + 0.3679(0.2142) + 0.2642(0.2142)1.1576
21.3679(1.1576)−0.3679(0.8463)+0+0.2642(0.2142)1.3679(1.1576) - 0.3679(0.8463) + 0 + 0.2642(0.2142)1.3288
31.3679(1.3288)−0.3679(1.1576)1.3679(1.3288) - 0.3679(1.1576)1.3918
41.3679(1.3918)−0.3679(1.3288)1.3679(1.3918) - 0.3679(1.3288)1.4149
51.3679(1.4149)−0.3679(1.3918)1.3679(1.4149) - 0.3679(1.3918)1.4234

Step 4: Closed form for k≥3k \ge 3

X(z)=G(z)U(z)=(0.3679z+0.2642)(z2+0.2142z+0.2142)z2(z−1)(z−0.3679)X(z) = G(z)U(z) = \dfrac{(0.3679z+0.2642)(z^2+0.2142z+0.2142)}{z^2(z-1)(z-0.3679)}. For k≥3k \ge 3 only the poles at z=1z = 1 and z=0.3679z = 0.3679 contribute to the residues of X(z)zk−1X(z)z^{k-1}:

Resz=1=(0.6321)(1.4284)1 (0.6321)=1.4284Resz=0.3679=(0.3995)(0.4284)(0.3679)2(−0.6321)(0.3679)k−1=−5.4375(0.3679)k\begin{aligned} \text{Res}_{z=1} &= \frac{(0.6321)(1.4284)}{1\,(0.6321)} = 1.4284 \\ \text{Res}_{z=0.3679} &= \frac{(0.3995)(0.4284)}{(0.3679)^2(-0.6321)}(0.3679)^{k-1} = -5.4375(0.3679)^k \end{aligned} x(k)=1.4284−5.4375(0.3679)k,k≥3x(k) = 1.4284 - 5.4375(0.3679)^k, \quad k \ge 3

Check: k=3k = 3 gives 1.4284−5.4375(0.04979)=1.15761.4284 - 5.4375(0.04979) = 1.1576, which matches.

Result

kk0123456∞\infty
x(k)x(k)00.36790.84631.15761.32881.39181.41491.4284

The final value from the theorem is lim⁡z→1(1−z−1)X(z)=0.6321(1+0.2142+0.2142)0.6321=1.4284\lim_{z\to1}(1-z^{-1})X(z) = \dfrac{0.6321(1 + 0.2142 + 0.2142)}{0.6321} = 1.4284.

Answer: x(0)=0x(0) = 0, x(1)=0.3679x(1) = 0.3679, x(2)=0.8463x(2) = 0.8463, and x(k)=1.4284−5.4375(0.3679)kx(k) = 1.4284 - 5.4375(0.3679)^k for k≥3k \ge 3. The output settles at 1.4284.

(With u(2)=−0.2142u(2) = -0.2142, the textbook deadbeat case, the output is exactly 1 for all k≥3k \ge 3.)

✓ Answer checked by Claude Opus 5.5 (high effort)

All Digital Control System chapter 2 questions →