blob: 10045067b0a41c227f997a74a13c6ac381fc1407 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
import numpy as np
import matplotlib.pyplot as plt
from scipy import signal
# numH = []
# denH =
# sysH = signal.TransferFunction(numH, denH)
# impulse = signal.impulse(sysH)
num = 1
den = [1, 2, 2]
sys = signal.TransferFunction(num, den)
t = [0.0001 * t for t in range(60000)]
y = signal.step(sys, t)
plt.plot(t, y)
plt.xlabel('t (sec)')
plt.ylabel('y(t)')
plt.title('MAE 171A - Homework 3.3\n Step Response')
plt.show()
|
