- 1、本文档共20页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
PAGE
PAGE - 20 -
目 录
实验1 离散时间信号的频域分析2
实验2 FFT算法与应用7
实验3 IIR数字滤波器的设计12
实验4 FIR数字滤波器的设计17
实验1 离散时间信号的频域分析
一.实验目的
信号的频域分析是信号处理中一种有效的工具。在离散信号的时域分析中,通常将信号表示成单位采样序列(n)的线性组合,而在频域中,将信号表示成复变量e或 e的线性组合。通过这样的表示,可以将时域的离散序列映射到频域以便于进一步的处理。
在本实验中,将学习利用MATLAB计算离散时间信号的DTFT和DFT,并加深对其相互关系的理解。
实验原理
(1)DTFT和DFT的定义及其相互关系。序列x(n)DTFT定义为= 它是关于自变量的复函数,且是以2为周期的连续函数。可以表示为,其中,和分别是实部和虚部;还可以表示为 =,其中, 和分别是的幅度函数和相位函数;它们都是的实函数,也是以2为周期的周期函数。
序列的N点DFT定义为,是周期为N的序列。与的关系:是对)在一个周期中的谱的等间隔N点采样,即 ,而可以通过对内插获得,即
(2)使用到的MATLAB命令有基于DTFT离散时间信号分析函数以及求解序列的DFT函数。
1)基于 DTFT离散时间信号分析函数有:freqz,real,imag,abs,angle。函数freqz可以用来计算一个以e的有理分式形式给出的序列的DTFT值。freqz的形式多样,常见的有H=freqz(num,den,w),其中num表示序列有理分式DTFT的分子多项式系数,den表示分母多项式系数(均按z的降幂排列),矢量w表示在0~2中给定的一系列频率点集合。freqz函数的其他形式参见帮助文件。在求出DTFT值后,可以使用函数real,imag,abs和angle分别求出并绘出起实部、虚部和相位谱。例如)=利用函数freqz计算出,然后利用函数abs和angle分别求出幅频特性与相位特性最后利用plot命令绘出曲线。
求解序列DFT的函数有:fft,ifft。函数fft(x)可以计算R点序列的R点DFT值;而fft(x,N)则计算R点序列的N点DFT,若RN,则直接截取R点DFT的前N点,若RN,则x先进行补零扩展为N点序列再求N点DFT。函数ifft(X)可以计算R点的谱序列的R点IDFT值;而ifft(X,N)同fft(x,N)的情况。
三、实验设备
计算机、MATLAB软件
四、实验内容
(1)编程计算并画出下面DTFT的实部,虚部、幅度和相位谱。
程序如下:num=[0.1313 -0.1553 0.1313 0.0518];
den=[1 1.2828 1.0388 0.3418];
w=0:0.001:2*pi;
H=freqz(num,den,w);
figure
subplot(221)
plot(w,real(H));
title(实部)
grid on;
subplot(222)
plot(w,imag(H))
title(虚部)
grid
subplot(223)
plot(w,abs(H))
title(幅度)
grid
subplot(224)
plot(w,abs(H))
title(相位谱)
grid on;
波形如下:
图1-1 DTFT的实部,虚部、幅度和相位谱
(2)计算32点序列x(n)=cosn,0n31的32点和64点DFT,分别绘出幅度谱图形,并绘出该序列的DTFT图形。
程序如下:n1=0:15;
n2=0:31;
x1=cos((13*pi*n1)/32)
x2=cos((13*pi*n2)/32)
X1=fft(x1);
X2=fft(x2);
subplot(211)
plot(n1,abs(X1))
title(32点DFT)
grid on;
subplot(212)
plot(n2,abs(X2))
title(64点DFT)
grid on;
波形如下:
图1-2 序列DFT幅度谱图形和DTFT图形
DTFT程序如下:
A=1;
n=0:31;
x=cos((13*pi*n)/32);
B=x;
w=0:0.01:2*pi*2;
[H]=freqz(B,A,w);
magH=abs(H);
subplot(2,1,1);
plot(w,magH);grid;
ylable(Magnitude);
subplot(2,1,2);
plot(w,phaH);grid;
xlable(w);
ylable(phase);
grid on;
波形如下:
图1-3 DTFT图形
实验2 FFT算法与应用
一、实验目的
在理论学习的基础上,通过本次实验。加深对快速傅里叶
您可能关注的文档
- 数据库表的管理和查询精编.doc
- 数据库系统第五版课后习题答案王珊精编.doc
- 数据库原理课后习题答案精编.doc
- 全国中学生物理竞赛模拟题精编.doc
- 全国中学生物理竞赛内容提要(俗称竞赛大纲)2016版精编.doc
- 数控毕业设计精编.doc
- 数控车削切削用量选择精编.doc
- 全静电发电机的模拟仿真设计精编.doc
- 全面轮胎型号含义精编.doc
- 全面引入数字化技术,推进基础物理教学的现代化精编.doc
- 上海海洋大学《海洋环境分析技术》课件-21化学发光分析法.pdf
- 上海海洋大学《海洋环境分析技术》课件-20分子荧光分析法.pdf
- 上海海洋大学《海洋环境分析技术》课件-22色谱分离过程.pdf
- 上海海洋大学《海洋环境分析技术》课件-25气相色谱仪与固定液.pdf
- 上海海洋大学《海洋环境分析技术》课件-24色谱定性定量方法.pdf
- 上海海洋大学《海洋环境分析技术》课件-26气相色谱检测器.pdf
- 上海海洋大学《海洋环境分析技术》课件-29液相色谱固定相与流动相.pdf
- 上海海洋大学《海洋环境分析技术》课件-27气相色谱分离条件的选择.pdf
- 上海海洋大学《海洋环境分析技术》课件-28液相色谱仪器与类型.pdf
- 上海海洋大学《海洋环境分析技术》课件-3 原子光谱和分析光谱.pdf
文档评论(0)