西安邮电数字信号处理实验报告
数字信号处理实验报告
实验一 信号、系统及系统响应一、实验目的1、熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对时域采样定理的理解。
2、熟悉离散信号和系统的时域特性。
3、熟悉线性卷积的计算编程方法：利用卷积的方法，观察、分析系统响应的时域特性。
4、掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离散信号、系统及其系统响应进行频域分析。
二、 实验原理1.理想采样序列：对信号x a (t)=A e −αt sin(Ω0t )u(t)进行理想采样，可以得到一个理想的采样信号序列x a (t)=A e −αt sin(Ω0nT ),0≤n ≤50，其中A 为幅度因子，α是衰减因子，Ω0是频率，T 是采样周期。
2.对一个连续时间信号x a (t)进行理想采样可以表示为该信号与一个周期冲激脉冲的乘积，即x ̂a (t)= x a (t)M(t),其中x ̂a (t)是连续信号x a (t)的理想采样；M(t)是周期冲激M(t)=∑δ+∞−∞(t-nT)=1T ∑e jm Ωs t +∞−∞,其中T 为采样周期，Ωs =2π/T 是采样角频率。
信号理想采样的傅里叶变换为X ̂a (j Ω)=1T ∑X a +∞−∞[j(Ω−k Ωs )],由此式可知：信号理想采样后的频谱是原信号频谱的周期延拓，其延拓周期为Ωs =2π/T 。
根据时域采样定理，如果原信号是带限信号，且采样频率高于原信号最高频率分量的2倍，则采样以后不会发生频率混叠现象。
三、简明步骤产生理想采样信号序列x a (n),使A=444.128，α=50√2π，Ω0=50√2π。
（1） 首先选用采样频率为1000HZ ，T=1/1000，观察所得理想采样信号的幅频特性，在折叠频率以内和给定的理想幅频特性无明显差异，并做记录；（2） 改变采样频率为300HZ ，T=1/300，观察所得到的频谱特性曲线的变化，并做记录；（3） 进一步减小采样频率为200HZ ，T=1/200，观察频谱混淆现象是否明显存在，说明原因，并记录这时候的幅频特性曲线。
数字信号处理实验报告
数字信号处理实验报告黎美琪通信一、实验名称：（快速傅里叶变换）的探究二、实验目的.学习理解的基本实现原理(注：算法主要有基时间抽取法和基频域抽取法，此实验讨论的是基频率抽取算法，课本上主要讲解的是基时间抽取算法).编写代码实现基频率抽取算法三、实验条件机四、实验过程(一)基础知识储备.基频率抽取( )算法基本原理：输入[]前后分解，输出[]奇偶分解。
设序列的点数为^，为整数（公式中的、定义不一样，打印后统一改正）将输入的[]按照的顺序分成前后两段：对输出的[]进行奇偶分解()、()和()之间可以用下图所示的蝶形运算符表示：的一次分解流图：的二次分解流图：最后完整的分解流图（^一共分解了三次）：的运算过程规律。
）^点的共进行级运算，每级由个蝶形运算组成。
同一级中，每个蝶形的两个输入数据只对计算本蝶形有用，而且每个蝶形的输入、输出数据结点又同在一条水平线上，也就是说计算完一个蝶形后，所得输出数据可立即存入原输入数据所占用的存储单元。
这样，经过级运算后，原来存放输入序列数据的个存储单元中便依次存放（）的个值。
（注：这种利用同一存储单元存储蝶形计算输入、输出数据的方法称为原位计算。
原位计算可节省大量内存，从而使设备成本降低。
））旋转因子的变化规律 ：以点的为例，第一级蝶形，，，，；第二级蝶形，；第三级的蝶形，。
依次类推，对于级蝶形，旋转因子的指数为∙^(−)，，，，，……，^()这样就可以算出每一级的旋转因子。
）蝶形运算两节点之间的“距离” ：第一级蝶形每个蝶形运算量节点的“距离”为，第二级每个蝶形运算另节点的“距离”为，第三级蝶形每个蝶形运算量节点的“距离”为。
依次类推：对于等于的次方的，可以得到第级蝶形每个蝶形运算量节点的“距离”为的次方。
.旋转因子 的性质1) 周期性 2) 对称性mk N N mk N W W -=+2 )可约性为整数/,//n N W W n mk n N mk N =．频率抽取（）基算法和时间抽取（）基算法比较：两种算法是等价的，其相同之处：（）与两种算法均为原位运算。
数字信号处理综合实验报告
综合实验1. 实验目的能综合利用信号处理的理论和Matlab 工具实现对信号进行分析和处理(1)熟练对信号进行时域和频域分析；(2)熟练进行滤波器设计和实现；(3)掌握对信号的滤波处理和分析。
2．实验原理设计并实现滤波器对信号进行分析和处理是信号处理课程学习的主要内容。
通过对信号进行频谱分析，能发现信号的频率特性，以及组成信号的频率分量。
对信号进行滤波处理，能改善信号的质量，或者为数据处理(如传输，分类等)提供预处理，等。
本次实验是对特定信号进行分析并进行滤波处理，需要综合应用之前的实验内容，主要有以下几个方面。
(1)离散时间信号与系统的时域分析Matlab 为离散时间信号与系统的分析提供了丰富且功能强大的计算函数和绘图分析函数，便于离散时间信号和系统的时域表示和分析。
(2)信号的频域分析信号处理课程主要学习了离散信号和系统的频域分析方法与实现，以及滤波器的设计与实现。
离散信号与系统的频域分析包括DTFT DFT Z变换等，FFT则是DFT的快速实现。
用Matlab分析信号的频谱可以用freqz函数或者FFT函数。
(3)滤波器设计滤波器的设计首先要确定滤波器的类型，即低通、高通、带通还是带阻。
滤波器的边缘频率可以通过对信号的频谱分析得到，滤波器的幅度指标主要有阻带最小衰减As 和通带最大衰减Ap。
一般来说，As越大，对截止通过的频率分量的衰减越大；Ap越小，对需要保留的频率分量的衰减越小。
因此，As 越大，Ap 越小，滤波器的性能越好，但随之而来，滤波器的阶数越大，实现的代价(包括计算时间和空间)越大。
由此，滤波器的设计需要对滤波器性能和实现代价进行均衡考虑。
另外根据冲激响应的长度可以分为IIR 和FIR 两种类型。
两种类型的滤波器各有特点。
用FIR 滤波器可以设计出具有严格线性相位的滤波器，但在满足同样指标的条件下，FIR 滤波器的阶数高于IIR 滤波器。
Matlab 为各种类型的滤波器的设计提供了丰富的函数，可以借助这些函数方便地设计出符合要求地滤波器。
数字信号处理实验报告
数字信号处理实验报告引言数字信号处理（Digital Signal Processing，DSP）是一门研究数字信号的获取、分析、处理和控制的学科。
在现代科技发展中，数字信号处理在通信、图像处理、音频处理等领域起着重要的作用。
本次实验旨在通过实际操作，深入了解数字信号处理的基本原理和实践技巧。
实验一：离散时间信号的生成与显示在实验开始之前，我们首先需要了解信号的生成与显示方法。
通过数字信号处理器（Digital Signal Processor，DSP）可以轻松生成和显示各种类型的离散时间信号。
实验设置如下：1. 设置采样频率为8kHz。
2. 生成一个正弦信号：频率为1kHz，振幅为1。
3. 生成一个方波信号：频率为1kHz，振幅为1。
4. 将生成的信号通过DAC（Digital-to-Analog Converter）输出到示波器上进行显示。
实验结果如下图所示：（插入示波器显示的正弦信号和方波信号的图片）实验分析：通过示波器的显示结果可以看出，正弦信号在时域上呈现周期性的波形，而方波信号则具有稳定的上下跳变。
这体现了正弦信号和方波信号在时域上的不同特征。
实验二：信号的采样和重构在数字信号处理中，信号的采样是将连续时间信号转化为离散时间信号的过程，信号的重构则是将离散时间信号还原为连续时间信号的过程。
在实际应用中，信号的采样和重构对信号处理的准确性至关重要。
实验设置如下：1. 生成一个正弦信号：频率为1kHz，振幅为1。
2. 设置采样频率为8kHz。
3. 对正弦信号进行采样，得到离散时间信号。
4. 对离散时间信号进行重构，得到连续时间信号。
5. 将重构的信号通过DAC输出到示波器上进行显示。
实验结果如下图所示：（插入示波器显示的连续时间信号和重构信号的图片）实验分析：通过示波器的显示结果可以看出，重构的信号与原信号非常接近，并且能够还原出原信号的形状和特征。
这说明信号的采样和重构方法对于信号处理的准确性有着重要影响。
数字信号处理实验报告_完整版
实验1 利用DFT 分析信号频谱一、实验目的1.加深对DFT 原理的理解。
2.应用DFT 分析信号的频谱。
3.深刻理解利用DFT 分析信号频谱的原理，分析实现过程中出现的现象及解决方法。
二、实验设备与环境 计算机、MATLAB 软件环境 三、实验基础理论1.DFT 与DTFT 的关系有限长序列 的离散时间傅里叶变换 在频率区间 的N 个等间隔分布的点 上的N 个取样值可以由下式表示：212/0()|()()01N jkn j Nk N k X e x n eX k k N πωωπ--====≤≤-∑由上式可知，序列 的N 点DFT ,实际上就是 序列的DTFT 在N 个等间隔频率点 上样本 。
2.利用DFT 求DTFT方法1：由恢复出的方法如下：由图2.1所示流程可知：101()()()N j j nkn j nN n n k X e x n eX k W e N ωωω∞∞----=-∞=-∞=⎡⎤==⎢⎥⎣⎦∑∑∑ 由上式可以得到：IDFTDTFT( )12()()()Nj k kX e X k Nωπφω==-∑ 其中为内插函数12sin(/2)()sin(/2)N j N x eN ωωφω--= 方法2：实际在MATLAB 计算中，上述插值运算不见得是最好的办法。
由于DFT 是DTFT 的取样值，其相邻两个频率样本点的间距为2π/N ，所以如果我们增加数据的长度N ，使得到的DFT 谱线就更加精细，其包络就越接近DTFT 的结果，这样就可以利用DFT 计算DTFT 。
如果没有更多的数据，可以通过补零来增加数据长度。
3.利用DFT 分析连续信号的频谱采用计算机分析连续时间信号的频谱，第一步就是把连续信号离散化，这里需要进行两个操作：一是采样，二是截断。
对于连续时间非周期信号,按采样间隔T 进行采样，阶段长度M ，那么：1()()()M j tj nT a a a n X j x t edt T x nT e ∞--Ω-Ω=-∞Ω==∑⎰对进行N 点频域采样，得到2120()|()()M jkn Na a M kn NTX j T x nT eTX k ππ--Ω==Ω==∑因此，可以将利用DFT 分析连续非周期信号频谱的步骤归纳如下： （1）确定时域采样间隔T ，得到离散序列（2）确定截取长度M ，得到M 点离散序列,这里为窗函数。
西电电院数字信号处理上机实验报告六
实验六、FIR数字滤波器设计及其网络结构班级: 学号: 姓名: 成绩:1实验目得(1)熟悉线性相位FIR数字滤波器得时域特点、频域特点与零极点分布;(2)掌握线性相位FIR数字滤波器得窗函数设计法与频率采样设计法;(3)了解IIR数字滤波器与FIR数字滤波器得优缺点及其适用场合。
2 实验内容(1)设计计算机程序,根据滤波器得主要技术指标设计线性相位FIR数字低通、高通、带通与带阻滤波器;(2)绘制滤波器得幅频特性与相频特性曲线,验证滤波器得设计结果就是否达到设计指标要求;(3)画出线性相位FIR数字滤波器得网络结构信号流图。
3实验步骤(1)设计相应得四种滤波器得MATLAB程序;(2)画出幅频相频特性曲线;(3)画出信号流图。
4 程序设计%% FIR低通f=[0、2,0、35];m=[1,0];Rp=1;Rs=40;dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);dat2=10^(-Rs/20);rip=[dat1,dat2];[M,f0,m0,w]=remezord(f,m,rip);M=M+2;hn=remez(M,f0,m0,w);w=0:0、001:pixn=[0:length(hn)-1];H=hn*exp(-j*xn'*w);figuresubplot(2,1,1)plot(w/pi,20*log10(abs(H)));gridon;xlabel('\omega/\pi'),ylabel('|H(e^j^w)|/dB')subplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel('\omega/\pi'),ylabel('\phi(\omega)/\pi') %% FIR高通f=[0、7,0、9];m=[0,1];Rp=1;Rs=60;dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);dat2=10^(-Rs/20);rip=[dat2,dat1];[M,f0,m0,w]=remezord(f,m,rip);hn=remez(M,f0,m0,w);w=0:0、001:pixn=[0:length(hn)-1];H=hn*exp(-j*xn'*w);figuresubplot(2,1,1)plot(w/pi,20*log10(abs(H)));gridon;xlabel('\omega/\pi'),ylabel('|H(e^j^w)|/dB')subplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel('\omega/\pi'),ylabel('\phi(\omega)/\pi') %% FIR带通f=[0、2,0、35,0、65,0、8];m=[0,1,0];Rp=1;Rs=60;dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);dat2=10^(-Rs/20);rip=[dat2,dat1,dat2];[M,f0,m0,w]=remezord(f,m,rip);M=M+3hn=remez(M,f0,m0,w);w=0:0、001:pixn=[0:length(hn)-1];H=hn*exp(-j*xn'*w);figuresubplot(2,1,1)plot(w/pi,20*log10(abs(H)));gridon;xlabel('\omega/\pi'),ylabel('|H(e^j^w)|/dB')subplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel('\omega/\pi'),ylabel('\phi(\omega)/\pi') %% FIR带阻f=[0、2,0、35,0、65,0、8];m=[1,0,1];Rp=1;Rs=60;dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);dat2=10^(-Rs/20);rip=[dat1,dat2,dat1];[M,f0,m0,w]=remezord(f,m,rip);hn=remez(M,f0,m0,w);w=0:0、001:pixn=[0:length(hn)-1];H=hn*exp(-j*xn'*w);figuresubplot(2,1,1)plot(w/pi,20*log10(abs(H)));gridon;xlabel('\omega/\pi'),ylabel('|H(e^j^w)|/dB')subplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel('\omega/\pi'),ylabel('\phi(\omega)/\pi') 5实验结果及分析(1)FIR低通滤波器自动得到得M值不满足要求,故我们将M加上2 在w=0、2π时,H=-0、5dB;w=0、35π时,H=-41dB。
数字信号处理(西电上机实验)
数字信号处理实验报告实验一：信号、系统及系统响应一、实验目的：(1) 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。
(2) 熟悉时域离散系统的时域特性。
(3) 利用卷积方法观察分析系统的时域特性。
(4) 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。
二、实验原理与方法：(1) 时域采样。
(2) LTI系统的输入输出关系。
三、实验内容、步骤(1) 认真复习采样理论、离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容，阅读本实验原理与方法。
(2) 编制实验用主程序及相应子程序。
①信号产生子程序，用于产生实验中要用到的下列信号序列：a. xa(t)=A*e^-at *sin(Ω0t)u(t)A=444.128;a=50*sqrt(2)*pi;b. 单位脉冲序列：xb(n)=δ(n)c. 矩形序列：xc(n)=RN(n), N=10②系统单位脉冲响应序列产生子程序。
本实验要用到两种FIR系统。
a. ha(n)=R10(n);b. hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3)③有限长序列线性卷积子程序用于完成两个给定长度的序列的卷积。
可以直接调用MATLAB语言中的卷积函数conv。
conv用于两个有限长度序列的卷积，它假定两个序列都从n=0 开始。
调用格式如下：y=conv (x, h)四、实验内容调通并运行实验程序，完成下述实验内容：①分析采样序列的特性。
a. 取采样频率fs=1 kHz, 即T=1 ms。
b. 改变采样频率，fs=300 Hz，观察|X(ejω)|的变化，并做记录(打印曲线)；进一步降低采样频率，fs=200 Hz，观察频谱混叠是否明显存在，说明原因，并记录(打印)这时的|X(ejω)|曲线。
②时域离散信号、系统和系统响应分析。
a. 观察信号xb(n)和系统hb(n)的时域和频域特性；利用线性卷积求信号xb(n)通过系统hb(n)的响应y(n)，比较所求响应y(n)和hb(n)的时域及频域特性，注意它们之间有无差别，绘图说明，并用所学理论解释所得结果。
数字信号处理实验报告
一、实验目的1. 理解数字信号处理的基本概念和原理。
2. 掌握离散时间信号的基本运算和变换方法。
3. 熟悉数字滤波器的设计和实现。
4. 培养实验操作能力和数据分析能力。
二、实验原理数字信号处理（Digital Signal Processing，DSP）是利用计算机对信号进行采样、量化、处理和分析的一种技术。
本实验主要涉及以下内容：1. 离散时间信号：离散时间信号是指时间上离散的信号，通常用序列表示。
2. 离散时间系统的时域分析：分析离散时间系统的时域特性，如稳定性、因果性、线性等。
3. 离散时间信号的变换：包括离散时间傅里叶变换（DTFT）、离散傅里叶变换（DFT）和快速傅里叶变换（FFT）等。
4. 数字滤波器：设计、实现和分析数字滤波器，如低通、高通、带通、带阻滤波器等。
三、实验内容1. 离散时间信号的时域运算（1）实验目的：掌握离散时间信号的时域运算方法。
（2）实验步骤：a. 使用MATLAB生成两个离散时间信号；b. 进行时域运算，如加、减、乘、除等；c. 绘制运算结果的时域波形图。
2. 离散时间信号的变换（1）实验目的：掌握离散时间信号的变换方法。
（2）实验步骤：a. 使用MATLAB生成一个离散时间信号；b. 进行DTFT、DFT和FFT变换；c. 绘制变换结果的频域波形图。
3. 数字滤波器的设计和实现（1）实验目的：掌握数字滤波器的设计和实现方法。
（2）实验步骤：a. 设计一个低通滤波器，如巴特沃斯滤波器、切比雪夫滤波器等；b. 使用MATLAB实现滤波器；c. 使用MATLAB对滤波器进行时域和频域分析。
4. 数字滤波器的应用（1）实验目的：掌握数字滤波器的应用。
（2）实验步骤：a. 采集一段语音信号；b. 使用数字滤波器对语音信号进行降噪处理；c. 比较降噪前后的语音信号，分析滤波器的效果。
四、实验结果与分析1. 离散时间信号的时域运算实验结果显示，通过MATLAB可以方便地进行离散时间信号的时域运算，并绘制出运算结果的时域波形图。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
西安邮电大学通信与信息工程学院实验报告（2016/2017学年第2学期）实验名称：数字信号处理实验学生姓名：***专业：通信工程学号：********指导教师：**2017年06月01日西安邮电大学通院《数字信号处理实验》独立设课实验过程考核表西安邮电大学通院《数字信号处理实验》独立设课实验成绩鉴定表1 熟悉MATLAB1.1实验目的1、了解MATLAB程序设计语言的基本特点，熟悉MATLAB软件的运行环境。
2、掌握各种信号的建模方式。
3、掌握各种信号的托兴表示方法。
4、掌握变量等有关概念，具备初步的将一般数学模型转化为对应的计机模型并进行处理的能力。
1.2 实验原理1、图形保持hold on/off命令控制是保持原有图形还是刷新原有图形，不带参数的hold命令在两种状态之间进行切换。
2、绘图修饰命令title(图形名称)xlabel(x轴说明)ylabel(y轴说明)text(x,y,图形说明)legend1.3 实验内容以及步骤单位采样序列代码：单位阶跃序列代码：矩形序列代码：三角波代码：方波代码：锯齿波代码：Sinc函数代码：1.4 实验结果分析学习使用MATLAB7.0软件；学习信号的图形表示方法，掌握各种各种信号的的建模方式；实现单位采样序列、单位阶跃序列、举行序列、三角波、方波、锯齿波、Sinc函数。
1.5 实验总结以及体会构成经过本次试验，我对于用MATLAB来作出一些数字信号中的一些基本波形的图像有了深刻的认识，可以熟练而又准确地以MATLAB为工具做出诸如单位阶跃序列，矩形序列等多种序列的图形，不仅如此，在作图中，我对于MATLAB的一些作图的基本指令，以及需要作出的序列有了深刻的认识与见解；本节课做出这些图形并不是重点，重点是了解MATLAB程序设计语言的基本特点，掌握这一个重要的工具，并且可以将数学模型熟练地转化为计算机模型，并进行处理，这才是本次实验真正的作用与意义，也是作为当代大学生的我们所必须具备的基本能力。
1.6 指导老师评语2 信号的表示2.1 实验目的1、了解MATLAB 程序设计语言的基本特点，熟悉MATLAB 软件运行环境。
2、掌握各种信号的建模方式。
3、掌握各种信号的图形表示方法。
4、掌握变量等有关概念，具备初步的将一般数学模型转化为对应的计算机模型并进行处理的能力2.2 实验原理常用的 MATLAB 绘图语句有figure、plot、subplot、stem 等，图形修饰语具有title、axis、text 等。
2.2.1 figure 语句figure 有两种用法。
当只有一句figure 命令时，程序会创建一个新的图形窗口，并返回一个整数型的窗口编号。
当采用figure(n)时，表示将第n 个图形窗口作为当前的图形窗口，将其显示在所有窗口的最前面。
如果该图形窗口不存在，则新建一个窗口，并赋以编号n。
2.2.2 plot 语句线形绘图函数。
用法为plot(x,y,’s’)。
参数x 为横轴变量，y 为纵轴变量，s 用以控制图形的基本特征如颜色、粗细等，通常可以省略，常用方法如表1-1 所示。
表 1-1 plot 命令的参数及其含义2.2.3 subplot 语句subplot(m,n,i)是分割显示图形窗口命令，它把一个图形窗口分为m 行n 列共m×n 个小窗口，并指定第i 个小窗口为当前窗口。
2.2.4 二维统计分析图在MATLAB 中，二维统计分析图形很多，常见的有条形图、阶梯图、杆图和填充图等，所采用的函数分别是：bar(x,y,选项)stairs(x,y,选项)stem(x,y,选项)fill(x1,y1,选项1,x2,y2,选项2,…)2.2.5 图形保持hold on/off 命令控制是保持原有图形还是刷新原有图形，不带参数的hold 命令在两种状态之间进行切换。
2.2.6 绘图修饰命令title(图形名称)xlabel(x 轴说明)ylabel(y 轴说明)text(x,y,图形说明)legend(图例1,图例2,…)2.2.7 MATLAB 常用信号生成函数（1） ZEROS• 功能：产生全零阵列• 调用格式：X=ZEROS (N) %产生N 行N 列的全零矩阵• X=ZEROS(M,N) %产生M 行N 列的全零矩阵（2） ONES• 功能：产生全 1 阵列• 调用格式：X=ONES(N) %产生N 行N 列的全1 矩阵• X=ONES(M,N) %产生M 行N 列的全1 矩阵（3） SINC• 功能：辛格函数• 调用格式：Y= SINC(X) %（4） RECTPULS• 功能：产生矩形脉冲信号• 调用格式：Y= RECTPULS (T) %产生高度为1、宽度为1、关于T＝0 对称的矩形脉冲• Y= RECTPULS (T，W) %产生高度为1、宽度为W、关于T＝0 对称的矩形脉冲（5） RAND• 功能：产生伪随机序列• 调用格式：Y= RAND (1,N) %产生[0,1]上均匀分布的随机序列• Y= RANDN (1,N) %产生均值为0，方差为1 的白噪声序列（6） SAWTOOTH• 功能：产生周期锯齿波或三角波• 调用格式：Y= SAWTOOTH (T) %产生幅值为＋1，－1，以2 为周期的方波• Y=SAWTOOTH(T,WIDTH)% 产生幅值为＋1，－1，以WIDTH *2 为周期的方波（7） SQUARE• 功能：产生方波• 调用格式：Y= SQUARE (T) %产生幅值为＋1，－1，以2 为周期的锯齿波• Y= SQUARE(T,DUTY) % 产生幅值为＋1，－1，以占空比为DUTY 的方波• 例： t = 0:.0001:.0625;• y = SQUARE(2*pi*30*t，80); plot(t,y)％产生一个占空比为80%的方波（8） FLIPLR• 功能：序列左右翻转• 调用格式：Y = FLIPLR(X)• % X = 1 2 3 翻转后 3 2 14 5 6 6 5 4（9） CUMSUM 、SUM• 功能：计算序列累加• 调用格式：Y = CUMSUM(X) % 向量X 元素累加，记录每一次的累加结果,而SUM只记录最后的结果。
2.3：实验内容以及步骤1、分别以条形图、阶梯图、杆图和填充图形式绘制曲线y=2sin(x)。
程序如下：2、采用图形保持，在同一坐标内绘制曲线y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)。
程序如下：3、在0≤x≤2π区间内，绘制曲线y1=2e-0.5x 和y2=cos(4πx)，并给图形添加图形标注。
程序如下：4、实现单位采样序列δ(n)、单位阶跃序列u(n)、矩形序列( ) N R n ，并用图形显示。
写出程序及输出图形5、实现三角波、方波、锯齿波、Sinc 函数，并用图形显示。
写出程序及输出图形2.4：实验结果分析1、分别以条形图、阶梯图、杆图和填充图形式绘制曲线y=2sin(x)。
2、采用图形保持，在同一坐标内绘制曲线y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)。
3、在0≤x≤2π区间内，绘制曲线y1=2e-0.5x 和y2=cos(4πx)，并给图形添加图标注。
4、实现单位采样序列δ (n)、单位阶跃序列u(n)、矩形序列 ( ) N R n ，并用图形显示，写出程序及输出图形。
5、实现三角波、方波、锯齿波、Sinc 函数，并用图形显示。
写出程序及输出图形2.5：实验总结以及体会构成本次实验主要是熟悉了如何使用MATLAB 表示信号，通过MATLAB 的函数来实现是主要方法，其次也可以自行编写函数，表示自己需要的特定信号，来适应不同的应用场合。
在信号表示中学会画图是非常重要的一件事，几乎所有的信号都可以也需要用图像表示出来，这样对分析和实验都是有着好处的，熟练掌握画图函数，这些便都可以迎刃而解，画图函数有很多种，其参数也是形形色色，记住每个参数的位置意义是非常重要的。
因为是对数字信号的处理，因此免不了会用到很多的理论知识，这时就是理论和时间结合的时候，MATLAB 是一个可以仿真信号的工具，可以模拟实际，需要用理论知识来设计解决问题所需要的模型，然后用MATLAB 将其仿真出来，看看实际结果是否是我们预期的，这个结果也将是我们在现实中实施该方案所能看到的结果。
2.6：指导老师评语3 FFT 频谱分析及应用3.1 实验目的1、通过实验加深对FFT 的理解；2、熟悉应用FFT 对典型信号进行频谱分析的方法。
3.1 实验原理在各种信号序列中，有限长序列占重要地位。
对有限长序列可以利用离散傅立叶变换(DFT)进行分析。
DFT 不但可以很好的反映序列的频谱特性，而且易于用快速算法(FFT)在计算机上进行分析。
有限长序列的DFT 是其z 变换在单位圆上的等距离采样，或者说是序列傅立叶的等距离采样，因此可以用于序列的谱分析。
FFT 是DFT 的一种快速算法，它是对变换式进行一次次分解，使其成为若干小数据点的组合，从而减少运算量。
在MATLAB 信号处理工具箱中的函数fft（x,n）,可以用来实现序列的N 点快速傅立叶变换。
经函数fft求得的序列一般是复序列，通常要求出其幅值和相位。
MATLAB中提供了求复数的幅值和相位的函数：abs、angle，这些函数一般和fft同时使用。
3.3.3 实验内容以及步骤被噪声污染得信号，比较难看出所包含得频率分量，如一个由50Hz 和120Hz正弦信号构成的信号，受零均值随机噪声的干扰，数据采样率为1000Hz，试用FFT 函数来分析其信号频率成分，要求：○1 画出时域波形；○2 分析信号功率谱密度。
3.3.4 实验结果分析3.3.5实验总结以及体会构成本次实验只要是针对有限长序列，我们可以通过MATLAB对其进行离散傅里叶变化，是非常方便的，因为工具内部已经将函数写好，我们只需要将变化的信号交给函数就可以，可以说大大地减小了手工计算量，并且能过直观得将图画出来。
其中需要注意的问题也有很多，比如因为操作失误，将函数名称写错，软件机会报错，这个问题还是比较容易查找的，其次是在表示信号的时候，将一些数值写错，到会最后的变换不是我们需要的，画出的图像有明显的问题，这些都需要细心才能避免，当然还有一些比较头疼的错误，比如函数不会使用，所给的参数错误，或者不符合标准。
3.3.6指导老师评语4 信号的运算-卷积4.1 实验目的1、掌握信号的线性卷积运算。
2、掌握信号的循环卷积运算。
3、掌握信号循环卷积计算线性卷积的条件。
4.2 实验原理数字信号卷积定理4.3 实验内容以及步骤1、假设卷积下面信号X(n)=9.0n 0≤n<130 elseh(n)= 1 0≤n<120 else选定循环卷积的长度为N = 21。