- 1、本文档共18页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
第三章离散时间系统的频域分析
例3.1计算离散时间傅立叶变换
%程序P3_1
%离散时间傅里叶变换的求取
clf;
%计算离散时间傅立叶变换的频率样本
w=-4*pi:8*pi/511:4*pi;
num=[21];den=[1-0.6];
h=freqz(num,den,w);
%PlottheDTFT
subplot(2,1,1)
plot(w/pi,real(h));grid
title(title(的实部)
xlabel(xlabel();
ylabel(振幅);
subplot(2,1,2)
plot(w/pi,imag(h));grid
title(title(的虚部)
xlabel(xlabel();
ylabel(振幅);
pause
subplot(2,1,1)
plot(w/pi,abs(h));grid
title(title(幅度谱)
xlabel(xlabel();
ylabel(振幅);
subplot(2,1,2)
plot(w/pi,angle(h));grid
title(相位谱相位谱)
xlabel(xlabel();
ylabel(以弧度为单位的相位);
Q3.1在程序P3.1中,计算离散时间傅里叶变换的原始系列是什么?MATLAB
命令pause的作用是什么?
原始序列:H(jw)=(2+e^(jw))/(1-0.6e^(-jw))
pause的作用:程序暂停,直至用户按任意一个按键.
Q3.2运行程序P3.1,求离散时间傅里叶变换的实部,虚部以及幅度和相位谱。
离散时间福利叶变换是w的周期函数吗?若是,周期是多少?描述这四个图形
表示的对称性。
离散时间福利叶变换是w的周期函数,周期为2π。
例3.2离散时间傅立叶变换的性质
%p3.2
%离散时间傅里叶变换的时移性质
clf;
w=-pi:2*pi/255:pi;
w0=0.4*pi;
D=10;
num=[123456789];
h1=freqz(num,1,w);
h2=freqz([zeros(1,D)num],1,w);
subplot(2,2,1);
plot(w/pi,abs(h1));grid
title(原序列的幅度谱);
subplot(2,2,2);
plot(w/pi,abs(h2));grid
title(时移后序列的幅度谱);
subplot(2,2,3)
plot(w/pi,angle(h1));grid
title(原序列的相位谱)
subplot(2,2,4)
plot(w/pi,angle(h2));grid
title(时移后序列的相位谱)
Q3.7运行修改后的程序并讨论你的结果。
根据结果可知当D增大时,序列的幅度谱不变;相位谱变密,周期缩短。
Q3.8选取不同的时移值重做习题Q3.7。
D=50
D=3
Q3.9选取两个改变了长度的序列以及两个不同的时移值,重做习题Q3.7。
D=30;num=[123456];
D=100;num=[12345678910];
%程序P3_3
%离散时间傅里叶变换的频移性质
clf;
w=-pi:2*pi/255:pi;wo=0.4*pi;
num1=[1357911131517];
L=length(num1);
h1=freqz(num1,1,w);
n=0:L-1;
num2=exp(wo*i*n).*num1;
h2=freqz(num2,1,w);
subplot(2,2,1)
plot(w/pi,abs(h1));grid
title(原序列的幅度谱)
subplot(2,2,2)
plot(w/pi,abs(h2));grid
title(频移后序列的幅度谱)
subplot(2,2,3)
plot(w/pi,angle(
文档评论(0)