实验一离散时间信号与系统时域分析
离散时间信号分析
实验一
学院：电气工程学院专业：测控技术与仪器班级：测仪101
实验二
学院：电气工程学院专业：测控技术与仪器班级：测仪101
实验三
学院：电气工程学院专业：测控技术与仪器班级：测仪101
实验四
学院：电气工程学院专业：测控技术与仪器班级：测仪101
实验五
学院：电气工程学院专业：测控技术与仪器班级：测仪101
%双线性变换法设计ButterWorth数字滤波器[n,Wn]=buttord(0.2,0.3,1,25,’s’);
[b,a]=butter(n,Wn,’s’);
freqs(b,a)
[bz,az]=bilinear(b,a,1);
通过本次实验,我基本掌握了双线性变换法及脉冲相应不变法设计
实验六
学院：电气工程学院专业：测控技术与仪器班级：测仪101。
信号与系统分析第六章 离散时间信号与系统的时域分析
应用上述性质, 可以将任意离散信号f(k)表示为单位序
列的延时加权和,
f ( k ) f ( 1 ) ( k 1 ) f ( 0 ) ( k ) f ( 1 ) ( k 1 )
f (n)(k n) n
同样, 根据单位序列δ(k)的特点,
(6.5)
f(k)(k) f(0)
第六章 离散时间信号与系统的时域分析
2. 单位阶跃序列ε(k) 单位阶跃序列ε(k)
(k) 10
k 0 k 0
(6.8)
ε(k)的波形如图6.3所示。 单位阶跃序列ε(k)类似于
连续时间系统的单位阶跃信号ε(t), 但应注意, ε(t)在t=0点
处发生跳变, 在此处不定义或定义为 定义为1。
, 而1 ε(k)在k=0处 2
实际处理时, 常把信号存放在处理器的存储单元 中, 随时取用, 也可以先记录数据后分析或短时间内存 入, 数据在较长时间内完成处理过程。 考虑到上述因 素, 离散时间信号f(kTs)可以不必以kTs为变量, 而可以 直接用f(k)表示离散信号, k为信号出现的序号。 用f(k) 表示离散信号不仅简便而且具有更为普遍的意义, 即 离散变量k可以不限于代表时间。 通常, 离散时间信 号也称为序列, 可以把它看成是一组序列值的集合。
可以看出, 任意信号与单位序列δ(k)相乘得到的仍然是 一个δ(k)序列, 只不过序列的幅度不再为1而是被f(0)加 权,δ(k)的这个性质称之为“加权性”, 或“取样性”。 推广后可以得到, 对于任意延时的单位序列δ(k-n),
f(k)δ(k-n)=f(n)δ(k-n) (6.4)
第六章 离散时间信号与系统的时域分析
k 0 k 0
(6.1)
第六章 离散时间信号与系统的时域分析
离散时间系统的时域特性分析实验报告
信号、系统与信号处理实验报告实验一、离散时间系统的时域特性分析姓名:学号:班级：专业：一．实验目的线性时不变（LTI）离散时间系统在时域中可以通过常系数线性差分方程来描述，冲激响应列可以刻画时域特性。
本次实验通过使用MATLAB函数研究离散时间系统的时域特性，以加深对离散时间系统的差分方程、冲激响应和系统的线性和时不变性的理解。
二．基本原理一个离散时间系统是将输入序列变换成输出序列的一种运算。
离散时间系统中最重要、最常用的是“线性时不变系统”。
1.线性系统满足叠加原理的系统称为线性系统，即若某一输入是由N个信号的加权和组成的，则输出就是系统对这几个信号中每一个输入的响应的加权和。
即那么当且仅当系统同时满足和时，系统是线性的。
在证明一个系统是线性系统时，必须证明此系统同时满足可加性和比例性，而且信号以及任何比例系数都可以是复数。
2.时不变系统系统的运算关系在整个运算过程中不随时间（也即序列的先后）而变化，这种系统称为时不变系统（或称移不变系统）。
若输入的输出为，则将输入序列移动任意位后，其输出序列除了跟着位移外，数值应该保持不变，即则满足以上关系的系统称为时不变系统。
3.常系数线性差分方程线性时不变离散系统的输入、输出关系可用以下常系数线性差分方程描述：当输入为单位冲激序列时，输出即为系统的单位冲激响应。
当时，是有限长度的，称系统为有限长单位冲激响应（FIR）系统；反之，则称系统为无限长单位冲激响应（IIR）系统。
三．实验内容及实验结果1.实验内容考虑如下差分方程描述的两个离散时间系统：系统1：系统2：输入：（1）编程求上述两个系统的输出，并画出系统的输入与输出波形。
（2）编程求上述两个系统的冲激响应序列，并画出波形。
（3）若系统的初始状态为零，判断系统2是否为时不变的？是否为线性的？2.实验结果（1）编程求上述两个系统的输出和冲激响应序列，并画出系统的输入、输出与冲激响应波形。
clf;n=0:300;x=cos((20*pi*n)/256)+cos((200*pi*n)/256);num1=[0.5 0.27 0.77];den1=[1];num2=[0.45 0.5 0.45];den2=[1 -0.53 0.46];y1=filter(num1,den1,x);y2=filter(num2,den2,x);subplot(3,1,1);stem(n,x);xlabel('时间信号');ylabel('信号幅度');title('输入信号');subplot(3,1,2);stem(y1);xlabel('时间信号n');ylabel('信号幅度');title('输出信号');subplot(3,1,3);stem(y2);xlabel('时间序号n ');ylabel('信号幅度');title('冲激响应序列');（2）N=40;num1=[0.5 0.27 0.77];den1=[1];num2=[0.45 0.5 0.45];den2=[1 -0.53 0.46];y1=impz(num1,den1,N);y2=impz(num2,den2,N);subplot(2,1,1);stem(y1);xlabel('时间信号n ');ylabel('信号幅度');title('³冲激响应');subplot(2,1,2);stem(y2);xlabel('时间信号n ');ylabel('信号幅度');title('³冲激响应');1.应用叠加原理验证系统2是否为线性系统：clear allclcn = 0 : 1 : 299;x1 = cos(20 * pi * n / 256);x2 = cos(200 * pi * n / 256);x = x1 + x2;num = [0.45 0.5 0.45];den = [1 -0.53 0.46];y1 = filter(num, den, x1);y2 = filter(num, den, x2);y= filter(num, den, x);yt = y1 + y2;figuresubplot(2, 1, 1);stem(n, y, 'g');xlabel('时间信号n');ylabel('信号幅度');axis([0 100 -2 2]);grid;subplot(2, 1, 2);stem(n, yt, 'r');xlabel('时间信号n');ylabel('信号幅度');axis([0 100 -2 2]);grid;2.应用时延差值来判断系统2是否为时不变系统。
离散时间信号和系统的时域分析
离散时间信号和系统的时域分析河南工业大学实验报告课程名称：数字信号处理开课实验室：6316实验报告撰写要求：认真总结实验，规范撰写实验报告。
实验报告内容应包括实验目的、实验要求、实验过程、实验总结，其中实验过程应附必要的截图，给出详细说明，对本实验中自行完成的较复杂网络拓扑的配置实现，应用表格给出各设备的主要参数配置（见下表），最后，对实验中遇到的问题和解决进行描述和剖析，总结收获。
并完成思考题。
实验一：离散时间信号和系统的时域分析一、实验目的：掌握用Matlab分析离散时间信号和系统的时域特性的方法。
二、实验环境：1. 运行Windows 2000 / 2007 / XP操作系统的PC一台；2. Matlab仿真环境；三、实验内容与要求：用Matlab在时域中产生一些基本的离散时间信号，并对这些信号进行一些基本运算,用Matlab仿真一些简单的离散时间系统，研究它们的时域特性。
四、实验步骤：Q1.23 产生长度为50、频率为0.08、振幅为2.5、相移为90度的正弦序列并显示它n = 0:50;f = 0.08;phase =90;A = 2.5;arg = 2*pi*f*n - phase;x = A*cos(arg);clf; % Clear old graphstem(n,x); % Plot the generated sequenceaxis([0 40 -3 3]);grid;title('Sinusoidal Sequence');xlabel('Time index n');ylabel('Amplitude');axis;Q 2.20 修改程序2.1，clf;N = 45;num = [0.9 -0.45 0.35 0.002];den = [1 0.71 -0.46 -0.62];y = impz(num,den,N);% Plot the impulse responsestem(y);xlabel('Time index n'); ylabel('Amplitude');title('Impulse Response'); grid;五、实验结果Q1.23运行结果Q2.1运行结果六、实验心得。
6.离散时间信号与系统的时域分析
0, n 1 1 z ( n) x ( n) y ( n) , n 1 2 1 n 1 ( 2 )( n 1)( 2 ) , n 0
6 线性时不变离散系统的时域分析
5. 累加 设某一序列为x(n),则x(n)的累加序列 y(n)定义为
y ( n)
k
x(k ) x(n) * u(n)
n
根据上述性质可以推得以下结论：
f (n n1 ) * (n n2 ) f (n n1 n2 )
6 线性时不变离散系统的时域分析
例 已知 x1 (n) (n) 3 (n 1) 2 (n 2) x2 (n) u(n) u(n 3) 试求信号 x (n) ，它满足 x(n) x1 (n) x2 (n) 解：可利用上面讲述的性质求解。
1 1/ 2 1/4 -2 -1 0 1 1/8 ... 2
n
x(-n) 1 1/2 1/8 1/4 ... -2 -1 0
1
2
n
6 线性时不变离散系统的时域分析
3.序列的加减 两序列的加、减是指同序号(n)的序列值逐项对 应相加得一新序列。
6 线性时不变离散系统的时域分析
例：
x(n) 1 1/2 1/4 -2 -1 0 y(n) 2 1 1/4 1/2 1 2 n …
6 线性时不变离散系统的时域分析
2.单位阶跃序列
u（n）
1, u ( n) 0,
n0 n0
u(n)
...
-1 0 1 2 3 n
(n) u (n) u (n) u (n 1)
m 0
u (n) (n m) (n) (n 1) (n 2)
实验一 时域离散信号、系统及系统响应
四、 思考题
• 1 在分析理想采样序列特性的实验中， 采样频率不同时， 相应 在分析理想采样序列特性的实验中， 采样频率不同时， 理想采样序列的傅里叶变换频谱的数字频率度量是否都相同? 理想采样序列的傅里叶变换频谱的数字频率度量是否都相同 它 们所对应的模拟频率是否相同? 为什么? 们所对应的模拟频率是否相同 为什么 • 2 在卷积定理验证的实验中， 如果选用不同的频域采样点数 值， 在卷积定理验证的实验中， 如果选用不同的频域采样点数M值 例如， 例如， 选M=10和M=20， 分别做序列的傅里叶变换， 求得 和 ， 分别做序列的傅里叶变换，
• 3 调通并运行实验程序， 完成下述实验内容： 调通并运行实验程序， 完成下述实验内容： 分析采样序列的特性。 ① 分析采样序列的特性。 a. 取采样频率 s=1 kHz, 即T=1 ms。 取采样频率f 。 b. 改变采样频率 fs=300 Hz， 观察 改变采样频率, 的变化， ， 观察|X(ejω)|的变化， 并 的变化 做记录(打印曲线 打印曲线)； 进一步降低采样频率， 做记录 打印曲线 ； 进一步降低采样频率， fs=200 Hz， ， 观察频谱混叠是否明显存在， 说明原因， 并记录(打印 打印) 观察频谱混叠是否明显存在， 说明原因， 并记录 打印 这时的|X(ejω)|曲线。 曲线。 这时的 曲线 • ② 时域离散信号、 系统和系统响应分析。 时域离散信号、 系统和系统响应分析。 a. 观察信号 b(n)和系统 b(n)的时域和频域特性； 利用 观察信号x 和系统h 的时域和频域特性； 和系统 的时域和频域特性 线性卷积求信号x 通过系统h 的响应y(n)， 比较 线性卷积求信号 b(n)通过系统 b(n)的响应 通过系统 的响应 ， 所求响应y(n)和hb(n)的时域及频域特性， 注意它们之 的时域及频域特性， 所求响应 和 的时域及频域特性 间有无差别， 绘图说明， 并用所学理论解释所得结果。 间有无差别， 绘图说明， 并用所学理论解释所得结果。 b. 观察系统 a(n)对信号 c(n)的响应特性。 观察系统h 对信号x 的响应特性。 对信号 的响应特性 ③ 卷积定理的验证
离散时间信号与系统的时域分析实验报告
离散时间信号与系统的时域分析实验报告报告⼆：⼀、设计题⽬1．绘制信号)()(1k k f δ=和)2()(2-=k k f δ的波形2.绘制直流信号)()(1k k f ε=和)2(2-=k f ε的波形3绘制信号)()(6k G k f =的波形⼆实验⽬的1．掌握⽤MATLAB 绘制离散时间信号（序列）波形图的基本原理。
2．掌握⽤MATLAB 绘制典型的离散时间信号（序列）。
3．通过对离散信号波形的绘制与观察，加深理解离散信号的基本特性。
三、设计原理离散时间信号(也称为离放序列)是指在时间上的取值是离散的，只在⼀些离放的瞬间才有定义的,⽽在其他时间没有定义，简称离放信号(也称为离散序列) 序列的离散时间间隔是等间隔(均匀)的，取时间间隔为T.以f(kT)表⽰该离散序列,k 为整数(k=0,±1.±2，...)。
为了简便，取T=1.则f(kT)简记为f(k)， k 表⽰各函数值在序列中出现的序号。
序列f(k)的数学表达式可以写成闭合形式,也可逐⼀列出f(k)的值。
通常,把对应某序号K0的序列值称为序列的第K0个样点的“样点值”。
四、设计的过程及仿真1clear all; close all; clc;k1=-4;k2=4;k=k1:k2;n1=0;n2=2;f1=[(k-n1)==0];f2=[(k-n2)==0];subplot(1,2,1)stem(k,f1,'fill','-k','linewidth',2);xlabel('k');ylabel('f_1(k)');title('δ(k)')axis([k1,k2,-0.1,1.1]);subplot(1,2,2)stem(k,f2,'filled','-k','linewidth',2);ylabel('f_2(k)');title('δ(k-2)')axis([k1,k2,-0.1,1.1]);程序运⾏后，仿真绘制的结果如图所⽰:2c lear all; close all; clc;k1=-2;k2=8;k=k1:k2;n1=0;n2=2; %阶跃序列开始出现的位置f1=[(k-n1)>=0]; f2=[(k-n2)>=0];subplot(1,2,1)stem(k,f1,'fill','-k','linewidth',2);xlabel('k');ylabel('f_1(k)');title('ε(k)')axis([k1,k2+0.2,-0.1,1.1])subplot(1,2,2)stem(k,f2,'filled','-k','linewidth',2);xlabel('k');ylabel('f_2(k)');title('ε(k-2)')axis([k1,k2+0.2,-0.1,1.1]);程序运⾏后，仿真绘制的结果如图所⽰:3clear all; close all; clc;k1=-2;k2=7;k=k1:k2; %建⽴时间序列n1=0;n2=6; f1=[(k-n1)>=0];f2=[(k-n2)>=0];f=f1-f2;stem(k,f,'fill','-k','linewidth',2);xlabel('k');ylabel('f(k)');title('G_6(k)')axis([k1,k2,-0.1,1.1]);程序运⾏后，仿真绘制的结果如图所⽰:五、设计的结论及收获实现了⽤matlab绘制离散时间信号, 通过对离散信号波形的绘制与观察，加深理解离散信号的基本特性。
数字信号处理实验离散时间 LTI 系统的时域分析与 Z 域分析
实验一离散时间LTI系统的时域分析与Z域分析一、实验目的1、掌握用MATLAB求解离散时间系统的零状态响应、单位脉冲响应和单位阶跃响应；2、掌握离散时间系统系统函数零极点的计算方法和零极点图的绘制方法，并能根据零极点图分析系统的稳定性。
二、实验原理1、离散时间系统的时域分析（1）离散时间系统的零状态响应离散时间LTI系统可用线性常系数差分方程来描述，即MATLAB中函数filter可对式（1-1）的差分方程在指定时间范围内的输入序列所产生的响应进行求解。
函数filter的语句格式为：y=filter(b,a,x)其中，x为输入的离散序列；y为输出的离散序列；y的长度与x的长度一样；b与a分别为差分方程右端与左端的系数向量。
（2）离散时间系统的单位脉冲响应系统的单位脉冲响应定义为系统在 (n)激励下系统的零状态响应，用h(n)表示。
MATLAB求解单位脉冲响有两种方法：一种是利用函数filter；另一种是利用函数impz。
impz函数的常用语句格式为impz(b,a,n)，其中b和a的定义见filter，n表示脉冲响应输出的序列个数。
（3）离散时间系统的单位阶跃响应系统的单位阶跃响应定义为系统在ε(n)激励下系统的零状态响应。
MATLAB求解单位脉冲响应有两种方法：一种是利用函数filter，另一种是利用函数stepz。
stepz函数的常用语句格式为stepz(b,a,N)其中，b和a的定义见filter，N表示脉冲响应输出的序列个数。
2、离散时间系统的Z域分析(1)系统函数的零极点分析离散时间系统的系统函数定义为系统零状态响应的z变换与激励的z变换之比，即如果系统函数H(z)的有理函数表示式为那么，在MATLAB中系统函数的零极点就可通过函数roots得到，也可借助函数tf2zp得到。
roots的语法格式为：Z=roots(b)%计算零点b=[b1b2…bmbm+1]P=roots(a)%计算极点a=[a1a2…anan+1]tf2zp的语句格式为[Z,P,K]=tf2zp(b,a)其中，b与a分别表示H(z)的分子与分母多项式的系数向量。
信号与系统实验报告
信号与系统实验报告一、信号的时域基本运算1．连续时间信号的时域基本运算两实验之一实验分析：输出信号值就等于两输入信号相加（乘）。
由于b=2，故平移量为2时，实际是右移1，符合平移性质。
两实验之二心得体会：时域中的基本运算具有连续性，当输入信号为连续时，输出信号也为连续。
平移，伸缩变化都会导致输出结果相对应的平移伸缩。
2.离散时间信号的时域基本运算两实验之一实验分析：输出信号的值是对应输入信号在每个n值所对应的运算值，当进行拉伸变化后，n值数量不会变，但范围会拉伸所输入的拉伸系数。
两实验之二心得体会：离散时间信号可以看做对连续时间信号的采样，而得到的输出信号值，也可以看成是连续信号所得之后的采样值。
二、连续信号卷积与系统的时域分析1.连续信号卷积积分两实验之一实验分析：当两相互卷积函数为冲激函数时，所卷积得到的也是一个冲激函数，且该函数的冲激t值为函数x，函数y冲激t值之和。
两实验之二心得体会：连续卷积函数每个t值所对应的卷积和可以看成其中一个在k值取得的函数与另外一个函数相乘得到的一个分量函数，并一直移动k值直至最后，最后累和出来的最终函数便是所得到的卷积函数。
3.RC电路时域积分两实验之一实验分析：全响应结果正好等于零状态响应与零输入响应之和。
两实验之二心得体会：具体学习了零状态，零输入，全响应过程的状态及变化，与之前所学的电路知识联系在一起了。
三、离散信号卷积与系统的时域分析1.离散信号卷积求和两实验之一实验分析：输出结果的n值是输入结果的k号与另一个n-k的累和两实验之二心得体会：直观地观察到卷积和的产生，可以看成连续卷积的采样形式，从这个方面去想，更能深入地理解卷积以及采样的知识。
2.离散差分方程求解两实验之一实验分析：其零状态响应序列为0 0 4 5 7.5，零输入响应序列为2 4 5 5.5 5.75，全状态响应序列为2 4 9 10.5 13.25，即全状态=零输入+零状态。
两实验之二心得体会：求差分方程时，可以根据全状态响应是由零输入输入以及零状态相加所得，分开来求，同时也加深了自己对差分方程的求解问题的理解。
信号与系统实验报告实验一 信号与系统的时域分析
实验一信号与系统的时域分析一、实验目的1、熟悉和掌握常用的用于信号与系统时域仿真分析的MA TLAB函数；2、掌握连续时间和离散时间信号的MA TLAB产生，掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程；3、牢固掌握系统的单位冲激响应的概念，掌握LTI系统的卷积表达式及其物理意义，掌握卷积的计算方法、卷积的基本性质；4、掌握利用MA TLAB计算卷积的编程方法，并利用所编写的MA TLAB程序验证卷积的常用基本性质；掌握MATLAB描述LTI系统的常用方法及有关函数，并学会利用MA TLAB求解LTI 系统响应，绘制相应曲线。
基本要求：掌握用MA TLAB描述连续时间信号和离散时间信号的方法，能够编写MATLAB程序，实现各种信号的时域变换和运算，并且以图形的方式再现各种信号的波形。
掌握线性时不变连续系统的时域数学模型用MATLAB描述的方法，掌握卷积运算、线性常系数微分方程的求解编程。
二、实验原理信号（Signal）一般都是随某一个或某几个独立变量的变化而变化的，例如，温度、压力、声音，还有股票市场的日收盘指数等，这些信号都是随时间的变化而变化的，还有一些信号，例如在研究地球结构时，地下某处的密度就是随着海拔高度的变化而变化的。
一幅图片中的每一个象素点的位置取决于两个坐标轴，即横轴和纵轴，因此，图像信号具有两个或两个以上的独立变量。
在《信号与系统》课程中，我们只关注这种只有一个独立变量（Independent variable）的信号，并且把这个独立变量统称为时间变量（Time variable），不管这个独立变量是否是时间变量。
在自然界中，大多数信号的时间变量都是连续变化的，因此这种信号被称为连续时间信号（Continuous-Time Signals）或模拟信号（Analog Signals），例如前面提到的温度、压力和声音信号就是连续时间信号的例子。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验一离散时间信号与系统时域分析实验目的1学习MATLAB语言编程和调试技巧2学会简单的矩阵输入和图形表示法3掌握简单的绘图命令
一实验目的
1学习MATLAB语言编程和调试技巧
2学会简单的矩阵输入和图形表示法
3掌握简单的绘图命令
二、实验原理
本实验主要为了熟悉MATLAB环境，重点掌握简单的矩阵（信号）输入和绘图命令，特别是绘图命令tem()和plot()。
实验内容中涉及到信号的无失真采样、离散卷积运算和差分方程求解这三个主要的问题。
其基本原理分别如下：
对一个模拟信号某(t)进行采样离散化某(n)，为了不失真地从采样信号某(n)中恢复原始信号某(t)，采样时必须满足采样定理，即采样频率必须大于等于模拟信号中最高频率分量的2倍。
一个离散时间系统，输入信号为某(n)，输出信号为y(n)，运算关系用T[﹒]表示，则输入与输出的关系可表示为y(n)＝T[某(n)]。
（1）线性时不变(LTI)系统的输入输出关系可通过h(n)表示：
y(n)=某(n)某h(n)=
式中某表示卷积运算。
（2）LTI系统的实现
可物理实现的线性时不变系统是稳定的、因果的。
这种系统的单位脉冲响应是因果的（单边）且绝对可和的，即：
h(n)0,n0;nh(n)0在MATLAB语言中采用conv实现卷积运算，即：
Y=conv(某,h),它默认从n=0开始。
常系数差分方程可以描述一个LTI系统，通过它可以获得系统的结构，也可以求信号的瞬态解。
利用MATLAB 自带的filter（），可以代替手工迭代运算求解系统的
差分方程，求解的过程类似于对输入信号进行滤波处理。
三、实验内容
1、试画出如下序列的波形
（1）某(n)3(n3)(n2)2(n1)4(n1)2(n2)3(n3)（2）某(n)0.5R10(n)解：用MATLAB描述波形
1(1)某=[3120-42-3];%矩阵输入某
n=-3:1:3;%输入自变量n,以间隔为1从-3到3变化n
实验目的1学习MATLAB语言编程和调试技巧2学会简单的矩阵输入和图形表示法3掌握简单的绘图命令
tem(n,某);%tem()函数绘制火柴杆图，注意n,某元素个数必须相等某label('n');%横坐标显示n
ylabal('某(n)');%纵坐标显示某(n)
grid;%绘制网格
1(2)n=0:9;
某=0.5.^n;
tem(n,某);
某label('n');
ylabel('某(n)');
gri
实验目的1学习MATLAB语言编程和调试技巧2学会简单的矩阵输入和图形表示法3掌握简单的绘图命令
2、用MATLAB计算序列{-201–13}和序列{120-1}的离散卷积，即计算
某(n)2(n)(n2)(n3)3(n4)与h(n)(n)2(n1)(n3)解：用MATLAB描述波形。
程序：
某=[-201-13];
h=[120-1];
c=conv(某,h);%conv()函数计算序列卷积
M=length(c)-1;%计算序列c的有效长度
n=0:1:M;
tem(n,c);
某label('n');
ylabel('幅度');
实验目的1学习MATLAB语言编程和调试技巧2学会简单的矩阵输入和图形表示法3掌握简单的绘图命令
10t),进行理想采样。
3、已知模拟信号为某(n)3co(2t)2co(6t)co(（1）由MATLAB仿真程序画出某(t)的波形。
（2）为了不失真地从采样信号中恢复原始信号，求奈奎斯特采样频率。
（3）分别以采样周期T=0.2和0.1对模拟信号进行采样，求采样输出序列解：（1）某(t)波形仿真程序如下：
t=0:1/256:2;
某=3某co(2某pi某t)-2某co(6某pi某t)+co(10某pi某t);%MATLAB表示某
plot(t,某);
grid;
实验目的1学习MATLAB语言编程和调试技巧2学会简单的矩阵输入和图形表示法3掌握简单的绘图命令
（2）信号某(t)的最高频率分量fma某=5Hz，由此奈奎斯特采样频率=10Hz=2fma某。
（3）由MATLAB仿真程序画出波形及其采样点，下面可取T=0.2或0.1,可分别获得所要求的信号。
T=0.2;
t=0:T:2;
某=3某co(2某pi某t)-2某co(6某pi某t)+co(10某pi某t);
tem(t,某);
grid;
实验目的1学习MATLAB语言编程和调试技巧2学会简单的矩阵输入和图形表示法3掌握简单的绘图命令
4、用MATLAB计算差分方程
当输入序列为某(n)(n)时的输出结果y(n),0n40
解MATLAB程序如下
N=41;
a=[0.8-0.440.360.22];%按顺序输入某(n)的系数矩阵
b=[10.7-0.45-0.6];%按顺序输入y(n)的系数矩阵
某=[1zero(1,N-1)];%仿真输入信号δ(n)
n=0:1:N-1;
y=filter(a,b,某);%利用filter函数求解差分方程，结果保存在矩阵y中
tem(n,y);
某label('n');
ylabel('幅度');
实验目的1学习MATLAB语言编程和调试技巧2学会简单的矩阵输入和图形表示法3掌握简单的绘图命令。