序列的移位和周期延拓课程设计
课程设计拓扑序列求解

课程设计拓扑序列求解一、教学目标本课程的教学目标是使学生掌握拓扑序列求解的基本概念和方法,能够运用所学知识解决实际问题。
具体目标如下:1.理解拓扑序列的定义及其性质。
2.掌握拓扑排序的算法及其应用。
3.了解拓扑序列求解在图论和其他学科中的应用。
4.能够运用拓扑排序算法对给定的有向无环图进行拓扑序列的求解。
5.能够运用所学知识解决实际问题,如任务调度、项目规划等。
情感态度价值观目标:1.培养学生的逻辑思维能力和问题解决能力。
2.激发学生对图论和算法研究的兴趣。
3.培养学生团队合作和自主学习的精神。
二、教学内容本课程的教学内容主要包括拓扑序列的定义和性质、拓扑排序算法及其应用。
具体安排如下:1.拓扑序列的定义和性质:介绍拓扑序列的概念,解释其性质和特点,举例说明其应用。
2.拓扑排序算法:讲解拓扑排序的基本思想和算法步骤,分析算法的正确性和时间复杂度。
3.拓扑序列求解的应用:通过实例介绍拓扑序列求解在图论和其他学科中的应用,如任务调度、项目规划等。
三、教学方法为了激发学生的学习兴趣和主动性,本课程将采用多种教学方法相结合的方式进行教学。
具体方法如下:1.讲授法:通过讲解拓扑序列的定义、性质和拓扑排序算法,使学生掌握基本概念和方法。
2.案例分析法:通过分析实际问题案例,使学生了解拓扑序列求解的应用和意义。
3.实验法:安排实验课程,让学生动手实践,加深对拓扑排序算法的理解和掌握。
四、教学资源为了支持教学内容和教学方法的实施,丰富学生的学习体验,将选择和准备以下教学资源:1.教材:选择一本与拓扑序列求解相关的教材,作为学生学习的主要参考资料。
2.参考书:推荐一些与拓扑序列求解相关的参考书籍,供学生深入研究。
3.多媒体资料:制作课件和教学视频,以图文并茂的形式展示拓扑序列求解的概念和方法。
4.实验设备:准备相应的实验设备,如计算机、网络设备等,以便进行实验教学。
五、教学评估本课程的评估方式将包括平时表现、作业和考试等,以全面客观地评估学生的学习成果。
mt—课程设计

1、 课程设计目的:综合运用本课程的理论知识进行语音信号的频谱分析,通过理论推导得出相应结论,并利用MATLAB 作为编程工具进行计算机实现,从而复习巩固课堂所学的理论知识,提高对所学知识的综合应用能力,并从实践上初步实现对数字信号的处理。
2、课程设计题目及要求(略)3、原理设计与分析在本次课程设计中,不管是连续信号还是离散信号,我们都采用FFT 对其进行分析。
若是连续信号先对其进行采样,由采样定理知道,若对信号进行理想采样,那么在频域其频谱将会发生周期延拓,那么得到采样信号的频谱后,经过低通滤波器,进行滤波(滤波频率 Ωm 应满足Ωc <Ωm <Ωs -Ωc ),这样就可以得到原信号的频谱,也即可以将原信号恢复。
这也即是,为什么可以对连续信号分析时用离散信号的分析方法对其进行分析。
当对连续信号进行采样后得到离散信号x(n ),我们作傅里叶变换得到X(e j ω)对其频谱进行分析。
然后,在频域我们经过对X(e j ω)进行采样,就得到频域的离散频谱X(k )。
对离散信号作DFT 时,我们可以通过利用MATLAB 里面的库函数fft(x(n ),N),这样就可以省去我们大部分的时间去做算法设计,在观察栅栏效应时对于同一个离散信号x(n ),我们作不同的点数的DFT 时,我们只需要对函数fft(x(n ),N)中的N 进行赋不同的之即可,省去了我们对信号的后部分补零的麻烦,方便快捷。
最后,我们还需对信号作卷积,其中包括线性卷积和循环卷积。
针对线性卷积h(n)=x(n )*y(m ),我们直接使用MATLAB 里边的库函数conv(x(n ),y(m ))来做,很方便。
然而对于循环卷积就没有现成的函数可以使用,这就需要我们自己做算法分析了,下面给出两种做法:(1)、算法设计)(1n x 与)(2n x 的N 点循环卷积定义为:)(]))(()([)(112n R m n x m x n y N N m N ∑-=-=循环卷积的实现步骤:)()(21m x m x Nm x ))((2Nm x ))((2-)())((2m R m x N N -)())((2m R m n x N N -当1,,2,1,0-=N n 时,分别将)(1m x 与)())((2m R m n x N N -相乘,并在0=m 到1-N 区间内求和,便得到)(1n x 与)(2n x 的循环卷积)(n y 。
《离散傅里叶变换-第三章》

n0 0 = kn 8 7
3
3
2π − j kn 8
3 − j kπ 8
(2) 设变换区间N=16, 则
X(k) = ∑ x(n)W
n= 0
3π k −j 16
π
N= 0 = n0 0
2 = ∑ e, k = 0,1, ⋅ ⋅ ⋅, 7 π N =0 sin( k ) 8
2. 时域循环移位定理 设x(n)是长度为N的有限长序列,y(n)为x(n)的循环移位,即: y(n)=x((n+m))NRN(n) 则: Y(k)=DFT[y(n)]=W-kmNX(k) 其中:X(k)=DFT[x(n)], 0≤k≤N-1
kn 证明: Y ( k ) = DFT [ y (n )] = x (( n + m )) N RN (n )WN ∑ N− 令n+m=n′,则有1 n =0 N −1
~
~ ∞
x (n ) =
m =−∞
∑
x ( n + mN )
(3.1.5)
(3.1.6) ••
~
x (n ) ••
0
••
N-1
•
n
x (n ) = x ( n ) ⋅ RN (n )
~
~
••
••
~(n ) x
•• •
0
••
•
••
•• •
~
••
N-1
•
n
一般定义周期序列 x(n) 中从n=0到N-1的第一个周期为 x(n)的主 n) x(n) (3.1.7) x( 值区间,而主值区间上的序列称为x(n) 的主值序列。(3.1.7) x(n)
数字信号除处理课程设计-DFT在信号频谱分析中的应用[1]
![数字信号除处理课程设计-DFT在信号频谱分析中的应用[1]](https://img.taocdn.com/s3/m/4aa75deaf8c75fbfc77db2ba.png)
湘南学院课程设计课程名称数字信号处理系别:计算机科学系专业班级:通信一班学号: 04 06 02 03 10 36 姓名:胡金霞、肖雅青、许芬、李真真、薛明、蒋小松题目: DFT在信号频谱分析中的应用完成日期: 2010年 12月 23日指导老师:樊洪斌2010年 12月 23日目录1、设计题目 (3)2、设计目的 (3)3、设计原理 (3)4、实现方法 (3)5、设计内容及结果 (6)6、改进建议 (12)7、思考题及解答 (15)8、设计体会 (15)9、参考文献 (16)Ⅰ.设计题目DFT 在信号频谱分析中的应用Ⅱ.设计目的掌握离散傅里叶变换的有关性质,利用Matlab 实现DFT 变换。
了解DFT 应用,用DFT 对序列进行频谱分析,了解DFT 算法存在的问题及改进方法。
学习并掌握FFT 的应用。
Ⅲ.设计原理所谓信号的频谱分析就是计算信号的傅里叶变换。
连续信号与系统的傅里叶分析显然不便于直接用计算机进行计算,使其应用受到限制,而DFT 是一种时域和频域均离散化的变换,适合数值运算,成为分析离散信号和系统的有力工具。
工程实际中,经常遇到的连续信号Xa(t),其频谱函数Xa(jW)也是连续函数。
数字计算机难于处理,因而我们采用DFT 来对连续时间信号的傅里叶变换进行逼近,进而分析连续时间信号的频谱。
Ⅳ.实现方法离散傅里叶变换是有限长序列的傅里叶变换,它相当于把信号的傅里叶变换进行等频率间隔采样,并且有限长序列的离散傅里叶变换和周期序列的离散傅里叶级数本质是一样的。
快速傅里叶变换(FFT )并不是一种新的变换,它是离散傅里叶变换的一种快速算法,并且主要是基于这样的思路而发展起来的:(1)把长度为N 的序列的DFT 逐次分解成长度较短的序列的DFT 来计算。
(2)利用WN(nk)的周期性和对称性,在DFT 运算中适当的分类,以提高运算速度。
(对称性nkNnk N W W N-=+2,12-=NN W ;周期性nkN nk N nrN N k rN n N W W W W ---==)(,r 为任意整数,1=nrNNW )离散傅里叶变换的推导:离散傅里叶级数定义为nk j N k p p ek x Nn x N21)(1)(π∑-==(1-1)将上式两端乘以nm j Neπ2-并对n 在0~N-1求和可得⎥⎦⎤⎢⎣⎡==∑∑∑∑∑-=---=-=-=---=-10)(110101)(1N2N2N2)()(1)(N n m k n j N N k p N n N k m k n j pN n nm j pe k X ek XNen xπππ 因为{m k 1mk 0)(N )(10)(N 2N2N2-1-1N 11=≠---=-==∑m k j m k j N n m k n je eeNπππ所以∑∑-=-=--=110)()()(N2N k p N n nm j pm k k X en xδπ 这样∑-=-=10N2)()(N n nm j p p en x m X π用k 代替m得∑-=-=10N2)()(N n nk j pP en xk X π(1-2)令N2πj NeW -=则(1-2)成为DFS []∑-===1)()()(N n nkNpp p W n xk X n x (1-3) (1-1)成为IDFS []∑-=-==10)(1)()(N n nkN p p p W k X N n x k X (1-4) 式(1-3)、(1-4)式构成周期序列傅里叶级数变换关系。
求拓扑排序序列课程设计

求拓扑排序序列课程设计一、教学目标本节课的教学目标是让学生掌握拓扑排序序列的概念、性质和求解方法,能够运用拓扑排序解决实际问题。
具体分为以下三个部分:1.知识目标:学生能够理解拓扑排序的定义,了解拓扑排序的性质和应用场景,掌握拓扑排序的求解方法。
2.技能目标:学生能够运用拓扑排序解决实际问题,如任务调度、项目规划等,提高问题解决的效率。
3.情感态度价值观目标:培养学生对计算机科学和图论的兴趣,培养学生的逻辑思维和创新能力。
二、教学内容本节课的教学内容主要包括以下三个方面:1.拓扑排序的定义和性质:介绍拓扑排序的定义,解释拓扑排序的性质,如无环性、唯一性等。
2.拓扑排序的求解方法:讲解拓扑排序的求解方法,如DFS、BFS等,并通过实例进行演示。
3.拓扑排序的应用:介绍拓扑排序在实际问题中的应用,如任务调度、项目规划等,并通过实例进行讲解。
三、教学方法为了提高学生的学习兴趣和主动性,本节课将采用以下教学方法:1.讲授法:讲解拓扑排序的定义、性质和求解方法,为学生提供系统的知识结构。
2.案例分析法:通过分析实际问题,让学生了解拓扑排序的应用,提高学生的解决问题的能力。
3.实验法:让学生动手实践,求解实际问题,培养学生的实际操作能力。
四、教学资源为了支持教学内容和教学方法的实施,丰富学生的学习体验,我们将准备以下教学资源:1.教材:为学生提供系统的知识结构,方便学生课后复习。
2.多媒体资料:通过动画、图片等形式,直观地展示拓扑排序的概念和性质,提高学生的学习兴趣。
3.实验设备:为学生提供实际操作的机会,培养学生的实际操作能力。
五、教学评估本节课的评估方式包括以下几个方面:1.平时表现:通过观察学生在课堂上的参与程度、提问回答等情况,评估学生的学习态度和理解程度。
2.作业:布置相关的拓扑排序题目,评估学生对拓扑排序概念和求解方法的理解和应用能力。
3.考试:设计考试题目,全面考察学生对拓扑排序的定义、性质、求解方法和应用的掌握程度。
时间序列分析(课程设计)教学大纲

《时间序列分析(课程设计)》教学大纲一、课程设计基本信息课程设计环节代码:110951课程设计环节名称:时间序列分析课程设计英文名称:TimeSeriesAna1ysis课程设计周数:一周学分:1适用对象:统计学专业本科学生先修课程与环节:统计学原理、抽样调查、多元统计分析等二、时间序列分析课程设计目的和任务课程设计是时间序列分析的同步课程,是上机实习课。
课程中,学生通过上机,学习揣摩这些方法模型的思想和构造,评价的优劣性。
另外,通过上机,加强学生的动手能力,为今后从事和软件开发打下良好的基础。
主要是通过对一些典型、通用时间序列分析方法的练习,加深学生对方法中的基本概念和基本理论的认识,使得学生掌握时间序列分析的基本方法和技巧介绍时间序列的基本知识、常用的建模和预测方法,在内容上强调平稳序列的特性和非平稳时间序列的处理方法。
通过学习本课程,使学生理解并掌握时间序列分析中的基本原理和方法,并具备建立合适的时间序列模型对时间序列数据进行拟合的基本技巧;训练学生介绍用软件对时间序列进行建模、预测;培养学生自行处理常规时间序列分析问题的能力和综合运用知识分析、解决问题的能力,达到理论与实践的统一。
三、时间序列分析课程设计方式先分组,每组确定一个题目,设计报告中要出现时间序列分析方法,对具体的问题通过网络及其他渠道收集相关数据,再对数据信息进行科学的分类、分析,然后在收集的数据基础之上,建立实证模型,并对模型进行检验和修正,完成整个课程设计任务。
四、时间序列分析指导方法与要求时间序列分析课程设计主要由学生上机自主完成,教师答疑的方式进行。
设计中要求综合运用所学的时间序列理论方面的知识,根据统计任务,对调查得来的(或者收集得到的二手资料)原始资料进行科学的分类、综合与加工分析,并利用常用的EVieWS和SAS等统计软件对数据进行描述性统计分析,画出相关统计图表,构建合理的统计模型进行模型分析,再对结果加以检验并修正。
信号系统MATLAB实验报告

信号与系统实验报告桂林理工大学信息科学与工程学院 电子信息工程实验二 信号及其表示【实验目的】了解各种常用信号的表达方式掌握部分绘图函数【实验内容】一、绘出连续时间信号x(t)=t e 707.0 sin 32t 关于t 的曲线,t 的范围为 0~30s ,并以递增。
MATLAB 源程序为:t=0::30; %对时间变量赋值x=exp*t).*sin(2/3.*t); %计算变量所对应得函数值 plot(t,x);grid; %绘制函数曲线ylabel('x(t)');xlabel('Time(sec)')二、产生周期为的方波。
MATLAB源程序为:Fs=100000;t=0:1/Fs:1;x1=square(2*pi*50*t,20);x2=square(2*pi*50*t,80);subplot(2,1,1),plot(t,x1),axis([0,,,]); subplot(2,1,2),plot(t,x2),axis([0,,,]);三、产生sinc(x)函数波形。
MATLAB源程序为:x=linspace(-4,4);y=sinc(x);plot(x,y)四、绘制离散时间信号的棒状图。
其中x(-1)=-1,x(0)=1,x(1)=2,x(2)=1,x(3)=0,x(4)=-1,其他时间x(n)=0。
MATLAB源程序为:n=-3:5; %定位时间变量x=[0,0,-1,1,2,1,-1,0,0];stem(n,x);grid; %绘制棒状图line([-3,5],[0,0]); %画X轴线xlabel('n');ylabel('x[n]')五、单位脉冲序列δ(n-0n )={00...1...0n n n n =≠直接实现:x=zeros(1,N);x(1,n0)=1;函数实现:利用单位脉冲序列)(0n n -δ的生成函数impseq,即 function[x,n]=impseq(n0,ns,nf)n=[ns:nf];x=[(n-n0)==0];plot(n,x);stem(n,x);输入参数:impseq(0,0,9)——连续图形012345678900.10.20.30.40.50.60.70.80.91输入参数:impseq(0,0,9)——离散图形六、单位阶跃序列ε(n-0n )={00...1...0n n n n ≥<直接实现:n=[ns:nf];x=[(n-n0)>=0];函数实现:利用单位阶跃序列)(0n n -ε的生成函数stepseq ,即 Function[x,n]=stepseq(n0,ns,nf)n=[ns:nf];x=[(n-n0)>=0];plot(n,x);七、实指数序列=,∀)(x n∈,Ranna直接实现:n=[ns:nf]:x=a.^n;函数实现:利用实指数序列n a(的生成函数rexpseq,即n)x=Function[x,n]=rexpseq(a,ns,nf)n=[ns:nf];x=a,^n:八、复指数序列n e n x n j ∀=+,)()(ωδ直接实现:n=[ns:nf];x=exp((sigema+jw)*n);函数实现:利用复指数序列n j e n x )()(ωδ+=的生成函数cexpseq,即 Function[x,n]=cexpseq(sigema,w,ns,nf)n=[ns:nf];x=exp((sigema+j*w)*n);0123456789-3000-2000-1000100020003000400050006000九、正(余)弦序列n wn n x ∀+=),cos()(θ直接实现:n=[ns:nf];x=cos(w*n+sita);函数实现:利用正(余)弦序列x(n)=cos(wn+θ)的生成函数cosswq,即Function[x,n]=cosseq(w,ns,nf,sita)n=[ns:nf];x=cos(w*n+sita);输入参数:cosseq,0,9,30)——连续信号0123456789-0.2-0.15-0.1-0.0500.050.10.150.2输入参数:cosseq,0,9,30)——离散信号0123456789实验三信号的运算【实验目的】了解信号处理的基本操作。
dsp.1-2复习教案

4.逆Z变换的域内环绕原点的一条逆时针闭合单围线.
可由留数定理去求该围线积分:
为c内的第k个极点, 为c外的第m个极点,Res[ ]表示极点处的留数。
当Zr为一阶极点时的留数:
(2)部分分式法
(3)幂级数展开法
5.Z变换的性质与定理
6.系统函数
(1)定义:单位采样响应的Z变换
(2)系统的因果稳定性的判断:
稳定条件:系统函数的收敛域包含单位圆。
因果稳定的条件:系统函数的所有极点都在单位圆内。
7.系统函数的零极点与系统幅频特性函数
一般情况下系统函数可表示成有理分式:
模:
相角:
小结
理解序列的概念及几种典型序列,掌握序列的运算,掌握线性卷积过程,会判断序列的周期性
4.用单位抽样序列表示任意序列
任意序列可表示成单位抽样序列的位移加权和.
二、线性移不变系统
1.系统的定义
一个离散时间系统 定义为:将输入序列 转换成输出序列 的一种运算或者变换即为 。
2.线性系统:满足叠加原理的系统。
3.移不变系统:系统对输入信号的处理不随时间变化。
如T[x(n)]=y(n),则T[x(n-m)]=y(n-m)。
三、时域离散线性移不变系统的描述
(1)时域描述:
①单位采样响应: 表示系统对 的响应输出,它可以将输出输入信号联系起来( ),所以可以用 表示系统对输入信号的处理功能.
②差分方程:
(2)频域描述:
①系统函数:
②频率响应函数:
四、时域离散线性移不变系统的输入/输出关系
(1)时域: 或差分方程
(2)频域:
一、时间离散信号——序列
1.序列的定义:在数字信号处理中,一般用x(n)表示时域表示离散信号(序列)。x(n)可看作是对模拟信号 的采样,即 。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
摘要复杂的序列通常可由基本序列通过运算和组合构成的,序列的移位和周期延拓是序列的基本运算。
序列的移位是序列的每一个样本都向右或向左移动k个单位,形成另一个序列。
周期延拓是把一个周期序列x(n),拓展为有K个周期的新序列。
MATLAB是“矩阵实验室”(MATrix LABoratoy)的缩写,是一种科学计算软件,主要适用于矩阵运算及控制和信息处理领域的分析设计。
本课题利用MATLAB的元素集运算和矩阵的运算实现了对序列移位和周期延拓运算的软件实现。
关键词:MATLAB;序列;移位;周期延拓目录1 课题描述 (1)2 设计原理 (1)3 设计过程 (2)3.1软件介绍 (2)3.2设计内容 (3)3.3设计步骤 (3)4程序运行结果及分析 (4)总结 (6)参考文献 (7)1 课题描述时域离散信号用x(n)表示,时间变量n(表示采样位置)只能取整数。
因此,x(n)是一个离散序列,以后简称序列。
序列适合计算机存储与处理。
序列的基本运算包括相加、相乘、移位、周期延拓等。
MATLAB是MATrix LABoratory的缩写,早期主要用于现代控制中复杂的矩阵、向量的各种运算。
MATLAB以矩阵作为基本编程单元,它提供了各种矩阵的运算与操作,并有较强的绘图功能。
本课题是利用MATLAB元素集运算和矩阵的运算实现了对序列移位和周期延拓运算的软件实现。
开发工具: MATLAB2设计原理设计原理如下:移位:在这个序列运算中,x[n]的每一个样本都移动(即延迟)k个采样周期,设移位后的序列为y(n)。
当k >0时每一个样本向右移动,称为x(n)的延时序列;当k<0时,每一个样本向左移动,称为x(n)的超前序列。
y(n)=x(n-k) (2.1)在MATLAB中,如果原始的序列用x和nx表示,移位后的序列用y和yn 表示,移位运算并不影响向量x的值,因此y=x。
移位体现为位置向量的改变。
ny的每个元素都比nx加了一个k,即ny=nx+k。
y和ny就是移位后的向量的表述,说明y取k拍前的x值。
向左移位可令k取负号,意味着y取k拍后的xz-进行标注,它被称为迟延算子,表示把输入序列右移一值。
在系统框图中用1位;用z进行标注,它是左移运算是右移算子的逆运算。
实际上迟延算子取的是序列过去的值,具有物理可实现性;而左移算子是提前算子,它要知道序列未来z-算子。
的值,物理上无法实现。
所以数字信号处理中通常都用1满足:周期延拓:如果对于所有的n,序列x[n](2.2)x[n]=x[n+kN]称为周期序列。
满足上述关系的最小N称为基本周期。
用MATLAB 则x[n]把一个周期序列{x(n),0<n<N-1},拓展为有K个周期的序列xtide(n)一般有两种方法。
简单复制法:设x是一个已赋值的行向量,其长度为N=length(x),把它复制K次,得到的xtide长度为K*N。
因此有,xtide=[X,X,……X];nxtide=0:K:N-1。
这个方法在K太大时容易数错,不是很好。
余函数mod法:函数n1=(n mod N)。
这个算式把大于等于N的n值,减去N的整倍数,使余数n1在0与N-1之间。
对于小于等于0的n值,则加以N 的整倍数。
也使n1在0与N-1之间。
把这一运算用到位置向量上,就可以方便的实现有限序列的周期延拓。
先设置位置向量,要复制K个x,则新向量的长度应为K*N。
其中N为x的长度,也就是所取的周期。
假如起始位置为0,则可用下列语句:nxtide=0:K*N-1; %设置延拓序列的位置向量xtide=x(mod(nxtide,N)+1) %确定位置向量各点对应的x值3设计过程3.1软件介绍MATLAB是一种科学计算软件,适用于工程应用个领域的分析设计与复杂计算,它使用方便,输入简捷,运算高效且内容丰富,很容易由用户自行扩展。
因此,当前已成为美国和其他发达国家大学教学和科学研究中最常用且必不可少的工具。
MATLAB以矩阵作为基本编程单元,它提供了各种矩阵的运算与操作,并有较强的绘图功能。
MATLAB集科学计算、图像处理、声音处理于一身,是一个高度的集成系统,有良好的用户界面,并有良好的帮助功能。
MATLAB不仅流行于控制界,在机械工程、生物工程、语音处理、图像处理、信号分析、计算机技术等各行各业中都有极广泛的应用。
MATLAB 语言的特点1.编程效率高 2.用户使用方便 3.扩充能力强 4.语句简单,内涵丰富 5.高效方便的矩阵和数组运算 6.方便的绘图功能数字信号处理的主要对象是数字信号,且是采用运算的方法达到处理目的,因此,其实现方法不同于模拟信号的实现方法,基本上可以分为两种实现方法,即软件实现方法和硬件实现的方法。
而硬件实现指的是选用合适的DSP 芯片,配有适合芯片语言及任务要求的软件,实现某种信号处理功能的一种方法。
数字信号处理的特点1.灵活性2.高精度和高稳定性3.便于大规模集成4.对数字信号可以存储、运算;系统可以获得高性能指标3.2设计内容已知()()n R n x n 88.0=,利用MATLAB 生成()n x ,()m n x -,()()()n R n x N 8和()()()n R m n x N 8-。
3.3设计步骤建模:取序列的观察区间为24。
利用MATLAB 的矩阵乘法和冒号运算可使周期延拓程序简短明了。
假如取三个延拓周期,则 x=[1 2 3 4]; y=x`*ones(1,3); =1 1 1 2 2 2 3 3 3 4 4 4y1=(y(:))’=[1 2 3 4 1 2 3 4 1 2 3 4]另一种更好的方法是采用MATLAB 求余函数mod ,y =x(mod(n ,M)+1)可实现对x(n)以M 为周期的周期延拓,其中求余后加1是因为MATLAB 向量下标只能从1开始,这样使程序更为简洁。
设计程序如下:N=24;M=8; m=3;n=0:N-1;x1=(0.8).^n;x2=[(n>=0)&(n<M)];%产生x (n ) xn=x1.*x2;xm=zeros(1,N);%设定Mxm 的长度 for k=m+1:m+M xm(k)=xn(k-m); endxc=xn(mod(n,8)+1);%产生x (n )的周期延拓xcm=xn(mod(n-m,8)+1);%产生x (n )移位后的周期延拓 subplot(2,2,1);stem(n,xn,'.'); xlabel('n');ylabel('xn'); axis([0,length(n),0,1]); subplot(2,2,2);stem(n,xm,'.'); xlabel('n');ylabel('xm'); axis([0,length(n),0,1]); subplot(2,2,3);stem(n,xc,'.'); xlabel('n');ylabel('xc'); axis([0,length(n),0,1]); subplot(2,2,4);stem(n,xcm,'.'); xlabel('n');ylabel('xcm'); axis([0,length(n),0,1]);4程序运行结果及分析运行结果如下图所示,显示出()n x ,()m n x -,()()()n R n x N 8和()()()n R m n x N 8-的图形。
5101520nx (n )05101520nx (n -m )5101520n周期延拓5101520n移位后周期延拓图4.1 序列()n x ,()m n x -,()()()n R n x N 8,()()()n R m n x N 8-的图形分析:图中x(n)是原序列;x(n-m)是原序列向右移位3得到的。
x(n)下面的图是对原序列进行3个周期的延拓;x(n-m )下面的图是对x(n-m)这一移位后的序列进行3个周期的延拓得到的。
本题用的是求余函数法实现序列的周期延拓的。
总结通过该课程设计,我基本掌握了MATLAB这一软件,还有MATLAB编译程序的基本过程,和构造工具及其相关的技术。
课本上的知识是有限的,并且是抽象的,枯燥的,难以理解的。
在本次课程设计,我有很大的收获,这不仅仅是理论知识上的完善,而且在自学能力和动手能力上有了新的飞跃!设计中,我感觉到自己知识极度匮乏,网查阅资料,翻阅各类相关书籍,完全靠自己设计,这让我的思维逻辑更加清晰。
在操作MATLAB中我基本握了一些基本程序的设计,为以后进一步的深入学习做好了铺垫。
在我设计好之后,老师对我进行指导,使得我的课程设计进一步完善,更加完美。
通过这次课程设计我还对数字信号处理中的关于序列的移位和周期延拓有了更深刻的认识,通过自己编程生成的图形更直接的看到了序列的移位和周期延拓的产生,对今后在数字信号处理方面和通信方面的学习有很大的帮助。
在这次课程设计过程中,我发现了自己综合应用能力的欠缺,也认识到计算机软件的强大功能。
以后,我会更加重视用软件编程,应用计算机来对处理信号。
参考文献[1] 丁玉美,高西全.《数字信号处理》(第三版).西安电子科技大学出版社,2000.[2] 陈怀堔,吴大正,高西全.《MATLAB及在电子信息课程中的应用》,电子工业出版社,2006.[3] 求是科技.《MATLAB7.0从入门到精通》.人民邮电出版社,2006.[4]高西全,丁玉美.《数字信号处理(第三版)》学习指导,西安科技大学出版社,2001.。