2015/4/14 回答同學之疑問¶
問題：y[n]=(x[n])^n 比較delay前後之差別
In [89]:
%matplotlib inline
import numpy as np
import matplotlib.pyplot as plt
n=np.arange(6)
x=np.array([0,1,2,3,4,0])
plt.figure(1)
plt.stem(n,x)
plt.xlim((-2,7))
plt.ylim((0,10))
plt.title('x[n]')
Out[89]:
<matplotlib.text.Text at 0xb504e90>
In [90]:
n=np.arange(6)
x_n_minus_n0=np.array([0,0,1,2,3,4,])
plt.figure(1)
plt.stem(n,x_n_minus_n0)
plt.xlim((-2,7))
plt.ylim((0,10))
plt.title('x[n-n0]')
Out[90]:
<matplotlib.text.Text at 0xa737a10>
再送入系統
In [91]:
w=x_n_minus_n0**n
#plt.stem(n,y)
print n
print w
[0 1 2 3 4 5] [ 1 0 1 8 81 1024]
In [94]:
y=x**n
y_n_minus_n0=np.insert(y,0,0)
print 'n'
print n
print 'x'
print x
print 'x_n_minus_n0'
print x_n_minus_n0
print 'w'
print w
print 'y'
print y
print 'y_n_minus_n0'
print y_n_minus_n0
n [0 1 2 3 4 5] x [0 1 2 3 4 0] x_n_minus_n0 [0 0 1 2 3 4] w [ 1 0 1 8 81 1024] y [ 1 1 4 27 256 0] y_n_minus_n0 [ 0 1 1 4 27 256 0]
結論： w[n] ==> 第3個時間點為8，是x[2]的3次方 y[n-1] ==> 第3個時間點為4，是x[2]的2次方 或許能回答你的問題
In [ ]: