原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
数字信号处理(程佩青)课后习题解答(1)--第1页
未知驱动探索，专注成就专业
数字信号处理(程佩青)课后习题解答(1)
题目1
题目描述：
给定一个长度为N的实数序列x(n)，编写一段Python代
码，实现求该序列的傅里叶变换，并将结果绘制成幅度谱和相
位谱。
解答：
importnumpyasnp
importmatplotlib.pyplotasplt
defDFT(x):
N=len(x)
X=np.zeros(N,dtype=complex)
forkinrange(N):
forninrange(N):
X[k]+=x[n]*np.exp(-2j*np.pi*k
*n/N)
returnX
输入实数序列
#x(n)
x=[1,2,3,4,5]
#计算傅里叶变换
1
数字信号处理(程佩青)课后习题解答(1)--第1页
数字信号处理(程佩青)课后习题解答(1)--第2页
未知驱动探索，专注成就专业
X=DFT(x)
#计算幅度谱和相位谱
X_mag=np.abs(X)
X_phase=np.angle(X)
#绘制幅度谱
plt.figure()
plt.stem(range(len(X_mag)),X_mag)
plt.xlabel(Frequency)
plt.ylabel(Magnitude)
plt.title(MagnitudeSpectrum)
#绘制相位谱
plt.figure()
plt.stem(range(len(X_phase)),X_phase)
plt.xlabel(Frequency)
plt.ylabel(Phase)
plt.title(PhaseSpectrum)
plt.show()
题目2
题目描述：
根据以下信号的采样频率和采样点数，计算并绘制它的离
散傅里叶变换（DFT）的幅度谱和相位谱。
采样频率：8kHz
2
数字信号处理(程佩青)课后习题解答(1)--第2页
数字信号处理(程佩青)课后习题解答(1)--第3页
未知驱动探索，专注成就专业
采样点数：1024
信号：sin(2πf0n)+0.5sin(2πf1n)
其中，f0=500Hz，f1=1000Hz。
解答：
importnumpyasnp
importmatplotlib.pyplotasplt
defDFT(x):
N=len(x)
X=np.zeros(N,dty
文档评论（0）