目录
一、公式
频域变量:z 时域变量:n
二、代码
以sin(nπ)为例
# 正弦序列
def sinZ(x,z):
'''
:param x: 时间序列
:param z: 频率序列
:return:
'''
sinZValues = []
for zz in z:
sinZValue = 0
for xx in x:
sinZValue += np.sin(np.pi*xx)*pow(zz,float(-xx))
sinZValues.append(sinZValue)
return sinZValues
if __name__ == '__main__':
from matplotlib import pyplot as plt
x = np.arange(0,20,1)
z = np.arange(0,20,0.1)
x = np.arange(0,20,1)
sinZV = sinZ(x,z)
plt.subplot(223)
plt.plot(z,sinZV)
三、结果