- 1、本文档共9页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字信号处理 课程设计
指导书
南通大学 电子信息学院
信息工程系
2013年11月
一、课程设计要求
1.熟练掌握Matlab语言的编程方法;
2.熟悉用于一维数字信号处理的Matlab主要函数的应用;
3.记录实验结果(包括波形和数据),撰写课程设计报告。
二、课程设计内容:
序列的产生
编写程序产生以下序列: x1(n)=(0.8)nu(n)、x3(n)=δ(n-4), x3(n)= R4(n),并画出波形;
参考程序:
n= -50:50; %生成时间向量n
x1=(n=0); %产生阶跃序列
%阶跃序列也可以用赋值的方法产生:x1=[zeros(1,50),ones(1,51)]
x21=(0.8).^n; x2=x21.*x1;
subplot(3,1,1);stem(n,x1,.)
subplot(3,1,2);stem(n,x2,.)
注意点乘和乘的区别。
求序列的傅里叶变换(DTFT)
(1) 编写序列傅里叶变换的函数
设x(n)是有限长的因果序列,编写求x(n)傅立叶变换的函数:function [X]=dtft(x,n,w),其中X代表,x代表x(n),w代表ω,n代表离散时间。
参考函数:
function [X]=dtft(x, n, w)
X=x*exp(-j* n *w); %通过行向量和列向量的矩阵乘,实现求和
将以上定义的函数存成名为dtft.m的文件,以备(2)中调用。
(2)利用傅里叶变换求序列的频谱
编写程序,在程序中调用(1)编写的dtft函数,求序列x(n)=u(n+2)-u(n-2)的幅度频谱和相位频谱。
说明:求幅度频谱和相位频谱,即求dtft复数值的幅值和辐角,可分别调用matlab函数库中函数abs和angle。
参考程序:
n= -100:100;
x=(n=-2)(n2); %或者x=[zeros(1,98),1,1,1,1,zeros(1,99)];
w=-2*pi : 0.01 : 2*pi;
X=dtft(x,n,w);
plot(w,abs(X));%求幅度频谱
以及的频谱,与x(n)频谱比较,验证傅里叶变换的频移性质和时移性质。
3.运用离散傅里叶变换(DFT)分析有限长序列的离散频谱
(1)对矩形窗序列x1(n)=R8(n)进行频谱分析,截取长度N=64,记录幅度频谱和相位频谱;
提示:运用matlab函数库中的fft函数计算序列的DFT,示例:y=fft(x,N)表示求x序列N点的DFT。
注意:由于需要绘出离散谱线,所以建议运用stem函数画图。
参考程序:
N=64;
n=0:N-1;
x=(n=0)(n8);
X=fft(x,N);
stem(n,abs(X),.) %注意了解该句中.的含义
2)选取截取长度N为周期序列周期的整数倍,对以下周期序列进行频谱分析:,请画出x2(n)的幅度频谱特性,记下最大峰点的高度和位置;改变N值(仍为周期的整数倍),观测峰点位置及高度的变化,并与理论结果相比较。
(3)将截取长度N改为周期的非整数倍,再次分析x2(n)幅频特性,与(2)的结果作比较,说明有何区别。
(4)观测截取长度N=250时、在分别对x2(n)加矩形窗和加海明窗两种情况下x2(n)的幅度频谱,记录幅度频谱的大致形状,并分析两种情况下频谱泄漏(单根谱线宽度加大)的程度。
提示:调用matlab函数hamming(如x=hamming(N)的形式),可以生成元素个数为N的海明窗列向量x。
4.抽样定理的验证
(1)生成连续信号,时间t取值区间为-5ms到5ms,画出该连续信号的波形;
提示:当时间点之间的间隔取很小数值时,所生成的离散信号可以用来近似表示连续信号。
(2)选择取样频率分别为20KHz、5KHz、1KHz对xa(t)进行取样,得到离散信号x1(n)、x2(n)、 x3(n),画出各离散信号的波形以及对应的幅度频谱, 离散时间的取值区间为-500到500,角频率的取值区间为到;
提示:生成离散信号的方法是:,其中取样周期Ts=1/Fs,Fs为取样频率。
(3)比较(2)中不同取样频率下离散信号的频谱,你得出什么结论?
参考程序:
dt=0.0005;
t=-0.005:dt:0.005;
xa=exp(-1000*abs(t));
figure(1)?;
plot(t,xa);
n=-500:500;
w=-3*pi:0.001:3*pi;
Ts=1/20000; %取样的频率为20KHz
x1=exp(-1000*abs(n*Ts));
X1=dtft(x1,n,w);
figure(2);
subplot(3,1,1);%
文档评论(0)