蒋朝晖数字信号处理实验
蒋朝晖数字信号处理实验.
《数字信号处理》实验指导书编写蒋朝辉中南大学信息科学与工程学院2014年4月目录实验一 MATLAB仿真软件的基本操作命令和使用方法 (2)实验二常见离散信号的MATLAB产生和图形显示 (8)实验三离散时间系统的时域分析 (12)实验四离散时间信号的DTFT (16)实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验内容1、帮助命令使用 help 命令，查找 sqrt （开方）函数的使用方法；2、MATLAB 命令窗口（1）在MATLAB 命令窗口直接输入命令行计算31)5.0sin(21+=πy 的值；（2）求多项式 p(x) = x3 + 2x+ 4的根；3、矩阵运算（1）矩阵的乘法已知A=[1 2;3 4]，B=[5 5;7 8]，求A^2*B（2）矩阵的行列式已知A=[1 2 3;4 5 6;7 8 9]，求A（3）矩阵的转置及共轭转置已知A=[1 2 3;4 5 6;7 8 9]，求A'已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B'（4）特征值、特征向量、特征多项式已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式；（5）使用冒号选出指定元素已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素；4、Matlab 基本编程方法（1）编写命令文件：计算1+2+…+n<2000 时的最大n 值；（2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。
5、MATLAB基本绘图命令（1）绘制余弦曲线 y=cos(t)，t∈[0，2π]（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π]（3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求：（a）线形为点划线、颜色为红色、数据点标记为加号；（b）坐标轴控制：显示范围、刻度线、比例、网络线（c）标注控制：坐标轴名称、标题、相应文本；>> clear;t=0:pi/10:4*pi;y=10*sin(t);plot(t,y);plot(t,y,'-+r');grid>> xlabel('X'),ylabel('Y');>> title('Plot:y=10*sin(t)');>> text(14,10,'完整图形');实验二常见离散信号的MATLAB产生和图形显示实验内容与步骤１. 写出延迟了np个单位的单位脉冲函数impseq，单位阶跃函数stepseq, n=ns:nf function [x,n]=impseq[np,ns,nf];function [x,n]=stepseq[np,ns,nf];2. 产生一个单位样本序列x1(n)，起点为ns= -10, 终点为nf=20, 在n0=0时有一单位脉冲并显示它。
数字信号 实验四
一、 实验目的和要求：（１）进一步掌握线性卷积的计算机编程方法，利用卷积的方法观察系统响应的时域特性。
（２）掌握循环卷积的计算机编程方法，并比较与线性卷积的差别，验证二者之间的关系。
利用循环卷积的方法观察、分析系统响应的时域特性。
二、 实验内容与原理：1.实验原理：（１）线性卷积：线性时不变系统（Linear Time-Invariant System, or L. T. I 系统）输入、输出间的关系为：当系统输入序列为)(n x ，系统的单位脉冲响应为)(n h ，输出序列为)(n y ，则系统输出为：∑∞-∞=-=*=m m n h m x n h n x n y )()()()()(；上式称为线性卷积。
（２）循环卷积设两个有限长序列)(1n x 和)(2n x ，长度分别为1N 和2N ，)()(11k X n x D FTN −−−→←点 )()(22k X n x D F T N −−−→←点。
如果)()()(21k X k X k X ⋅=，则∑---==121)())(()()]([)(N m N N n R m n x m x k X IDFT n x上式称为)(1n x 和)(2n x 的循环卷积。
（３）两个有限长序列的线性卷积序列)(1n x 和)(2n x ，长度分别为L 点和P 点，)(3n x 为这两个序列的线性卷积，则)(3n x 为∑∞-∞=-=*=m m n xm x n x n x n x )()()()()(21213且线性卷积)(3n x 的非零值长度为1-+P L 点。
（４）循环卷积与线性卷积的关系序列)(1n x 为L 点长，序列)(2n x 为P 点长，若序列)(1n x 和)(2n x 进行N 点的循环卷积)(n x c ，其结果是否等于该两序列的线性卷积)(n x l ，完全取决于循环卷积的长度。
由教材相关推导，得∑∞-∞=+=q Nlc n RqN n x n x )()()(，也就是说，循环卷积是线性卷积的周期延拓序列再取主值区间。
数字信号处理实验(设计性实验修改)
《数字信号处理》实验指导书（实验报告）Digital Signal Processing Laboratory湛柏明编蒋伟荣审班级:姓名:湖北汽车工业学院电子信息科学系二〇〇六年十二月修订前言《信号与系统》、《数字信号处理》是电子信息类专业的两门主要技术基础课程，是电子信息类专业本科生的必修课程，也是电子信息类专业硕士研究生入学必考课程。
该课程的任务在于研究信号与系统理论的基本概念和基本分析方法，使学生初步认识如何建立信号与系统的数学模型，如何经适当的数学分析求解，并对所得结果给予物理解释，赋予物理意义。
该课程的基本理论和方法大量用于计算机信息处理的各个领域，特别是通信、数字语音处理、数字图像处理、数字信号分析等领域应用更为广泛。
通过实验，配合《信号与系统》和《数字信号处理》课程的教学、加强学生对信号与系统理论的感性认识、提高学生的综合能力具有重要的意义。
长期以来，《信号与系统》和《数字信号处理》课程一直采用黑板式的单一教学方式，学生仅依靠做习题来巩固和理解教学内容，对课程中大量的应用性较强的内容不能实际动手设计、调试、分析，严重影响和制约了教学效果。
由于黑板式教学，课程中大量的信号分析结果缺乏可视化的直观表现，学生自己设计系统也不能直观地得到系统特性的可视化测试结果，学生将大量的时间和精力用于繁杂的手工数学运算，而未真正理解所得结果在信号处理中的实际意义。
近年来，计算机多媒体教序手段的运用逐步普及，大量优秀的科学计算和系统仿真软件不断涌现，为我们实现计算机辅助教学和学生上机实验提供了很好的平台。
通过对这些软件的分析和对比，我们选择MATLAB语言作为辅助教学工具，借助MATLAB强大的计算能力和图形表现能力，将《信号与系统》和《数字信号处理》中的概念、方法和相应的结果，以图形的形式直观地展现给学生，大大的方便学生迅速掌握和理解教学内容。
然而，我们意识到，按照之前的《信号与系统》和《数字信号处理》课程的各8个实验学时进行实验，实验效果比较不尽如人意，由于实验学时数太少，没有给学生更的时间先去了解MATLAB语言，以至于使实验课流于形式，由于实验学时太少，也导致我们无法安排更为细致的具有综合型和设计型的实验项目。
数字信号处理实验报告(二)
数字信号处理第二次实验报告学院：信息工程学院班级：2012级电子信息工程*班姓名：学号：20125507**指导老师：实验四：IIR数字滤波器设计及软件实现一、实验目的1、熟悉双线性变换设计IIR滤波器的原理与方法2、掌握IIR滤波器的MATLAB实现方法二、实验原理简述IIR数字滤波器间接法基本设计过程：1、将给定的数字滤波器的指标转换成过渡模拟滤波器的指标；2、设计过渡模拟滤波器；3、将过渡模拟滤波器系统函数转换成数字滤波器的系统函数三、程序与图形1、%-----------------信号产生函数mstg---------------function st=mstg %功能函数的写法%产生信号序列向量st,并显示st的时域波形和频谱%st=mstg 返回三路调幅信号相加形成的混合信号，长度N=1600N=1600 %N为信号st的长度。
Fs=10000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz，Tp为采样时间t=0:T:(N-1)*T;k=0:N-1;f=k/Tp;fc1=Fs/10; %第1路调幅信号的载波频率fc1=1000Hz,fm1=fc1/10; %第1路调幅信号的调制信号频率fm1=100Hzfc2=Fs/20; %第2路调幅信号的载波频率fc2=500Hzfm2=fc2/10; %第2路调幅信号的调制信号频率fm2=50Hzfc3=Fs/40; %第3路调幅信号的载波频率fc3=250Hz,fm3=fc3/10; %第3路调幅信号的调制信号频率fm3=25Hzxt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); %产生第1路调幅信号xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); %产生第2路调幅信号xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %产生第3路调幅信号st=xt1+xt2+xt3; %三路调幅信号相加fxt=fft(st,N); %计算信号st的频谱%-------绘制st的时域波形和幅频特性曲线-----subplot(2,1,1)plot(t,st);grid;xlabel('t/s');ylabel('s(t)');axis([0,Tp/8,min(st),max(st)]);title('(a) s(t)的波形')subplot(2,1,2)stem(f,abs(fxt)/max(abs(fxt)),'.');grid;title('(b) s(t)的频谱') axis([0,Fs/5,0,1.2]);xlabel('f/Hz');ylabel('幅度')-10123t/ss (t )(b) s(t)的频谱f/Hz幅度2、%-------实验4-2--------- clear all;close allFs=10000;T=1/Fs; %采样频率%调用信号产生函数mstg 产生由三路抑制载波调幅信号相加构成的复合信号st st=mstg;fp=280;fs=450; %下面wp,ws,为fp,fs 的归一化值范围为0-1wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF 指标（低通滤波器的通、阻带边界频）[N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord 计算椭圆DF 阶数N 和通带截止频率wp[B,A]=ellip(N,rp,rs,wp); %调用ellip 计算椭圆带通DF 系统函数系数向量B 和A[h,w]= freqz(B,A);y1t=filter(B,A,st); %滤波器软件实现 figure(2);subplot(2,1,1); plot(w,20*log10(abs(h))); axis([0,1,-80,0]) subplot(2,1,2);t=0:T:(length(y1t)-1)*T; plot(t,y1t);%axis([0,1,-80,0])-10123t/ss (t )(b) s(t)的频谱f/Hz幅度-80-60-40-20000.020.040.060.080.10.120.140.16-1-0.500.511.53、%-------实验4-3---------fpl=440;fpu=560;fsl=275;fsu=900;wp=[2*fpl/Fs,2*fpu/Fs];ws=[2*fsl/Fs,2*fsu/Fs];rp=0.1;rs=60;[N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord 计算椭圆DF 阶数N 和通带截止频率wp[B,A]=ellip(N,rp,rs,wp); %调用ellip 计算椭圆带通DF 系统函数系数向量B 和A[h,w]= freqz(B,A); y2t=filter(B,A,st);figure(3);subplot(2,1,1);plot(w,20*log10(abs(h))); axis([0,1,-80,0]) subplot(2,1,2);t=0:T:(length(y2t)-1)*T; plot(t,y2t);00.20.40.60.81-80-60-40-20000.020.040.060.080.10.120.140.16-2-10124、%-------实验4-4--------- fp=900;fs=550;wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF 指标（低通滤波器的通、阻带边界频）[N,wp]=ellipord(wp,ws,rp,rs);%调用ellipord 算椭圆DF 阶数N 通带截止频率 [B,A]=ellip(N,rp,rs,wp,'high'); %调用ellip 计算椭圆带通DF 系统函数系数向量B 和A[h,w]= freqz(B,A); y3t=filter(B,A,st);figure(4);subplot(2,1,1); plot(w,20*log10(abs(h))); axis([0,1,-80,0]) subplot(2,1,2);t=0:T:(length(y3t)-1)*T; plot(t,y3t);-80-60-40-20000.020.040.060.080.10.120.140.16-2-1012四、实验结果分析由图可见，三个分离滤波器指标参数选取正确，损耗函数曲线达到所给指标。
《数字信号处理》上机实验指导书
《数字信号处理》上机实验指导书一、引言“数字信号处理”是一门理论和实验密切结合的课程，为了深入地掌握课程内容，应当在学习理论的同时，做习题和上机实验。
上机实验不仅可以帮助学生深入地理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。
所以，根据本课程的重点要求编写了四个实验。
第一章、二章是全书的基础内容，抽样定理、时域离散系统的时域和频域分析以及系统对输入信号的响应是重要的基本内容。
由于第一、二章大部分内容已经在前期《信号与系统》课程中学习完，所以可通过实验一帮助学生温习以上重要内容，加深学生对“数字信号处理是通过对输入信号的一种运算达到处理目的”这一重要概念的理解。
这样便可以使学生从《信号与系统》课程顺利的过渡到本课程的学习上来。
第三章、四章DFT、FFT是数字信号处理的重要数学工具，它有广泛的使用内容。
限于实验课时，仅采用实验二“用FFT对信号进行谱分析”这一实验。
通过该实验加深理解DFT的基本概念、基本性质。
FFT是它的快速算法，必须学会使用。
所以，学习完第三、四章后，可安排进行实验二。
数字滤波器的基本理论和设计方法是数字信号处理技术的重要内容。
学习这一部分时，应重点掌握IIR和FIR两种不同的数字滤波器的基本设计方法。
IIR滤波器的单位冲激响应是无限长的，设计方法是先设计模拟滤波器，然后再通过S~Z平面转换，求出相应的数字滤波器的系统函数。
这里的平面转换有两种方法，即冲激响应不变法和双线性变换法，后者没有频率混叠的缺点，且转换简单，是一种普遍应用的方法。
学习完第六章以后可以进行实验三。
FIR滤波器的单位冲激响应是有限长的，设计滤波器的目的即是求出符合要求的单位冲激响应。
窗函数法是一种基本的，也是一种重要的设计方法。
学习完第七章后可以进行实验四。
以上所提到的四个实验，可根据实验课时的多少恰当安排。
例如：实验一可根据学生在学习《信号与系统》课程后，掌握的程度来确定是否做此实验。
若时间紧，可以在实验三、四之中任做一个实验。
(完整版)数字信号处理实验三
3.41；3.42 由教材可知： ，即序列的偶部分的傅立叶变换是序列的傅立叶变换的实部。
5、实验步骤
1、进行本实验，首先必须熟悉matlab的运用，所以第一步是学会使用matlab。
2、学习相关基础知识，根据《数字信号处理》课程的学习理解实验内容和目的。
plot(w/pi,angle(h1));grid
xlabel('\omega/\pi');ylabel('以弧度为单位的相位');
title('原序列的相位谱')
subplot(2,2,4)
plot(w/pi,angle(h2));grid
xlabel('\omega/\pi');ylabel('以弧度为单位的相位');
grid;
title('相位谱arg[H(e^{j\omega})]');
xlabel('\omega/\pi');
ylabel('以弧度为单位的相位');
3.4
clf;
w=-4*pi:8*pi/511:4*pi;
num1=[1 3 5 7 9 11 13 15 17];
h=freqz(num,1,w);
Q3.32 通过加入合适的注释语句和程序语句，修改程序P3.8，对程序生成的图形中的两个轴加标记。时移量是多少？
Q3.33 运行修改后的程序并验证离散傅里叶变换的圆周时移性质。
Q3.36 运行程序P3.9并验证离散傅里叶变换的圆周卷积性质。
Q3.38 运行程序P3.10并验证线性卷积可通过圆周卷积得到。
数字信号处理上机实验答案(第三版,第十章)[自己整理完善的]
![数字信号处理上机实验答案(第三版,第十章)[自己整理完善的]](https://img.taocdn.com/s3/m/ecd1b05f33687e21af45a9b3.png)
第十章 上机实验数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。
上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。
本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。
实验一 系统响应及系统稳定性。
实验二 时域采样与频域采样。
实验三 用FFT 对信号作频谱分析。
实验四 IIR 数字滤波器设计及软件实现。
实验五 FIR 数字滤波器设计与软件实现实验六 应用实验——数字信号处理在双音多频拨号系统中的应用任课教师根据教学进度，安排学生上机进行实验。
建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR 数字滤波器设计及软件实现在。
学习完第六章进行；实验五在学习完第七章后进行。
实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。
10.1 实验一: 系统响应及系统稳定性1.实验目的（1）掌握 求系统响应的方法。
（2）掌握时域离散系统的时域特性。
（3）分析、观察及检验系统的稳定性。
2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。
已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。
在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB 语言的工具箱函数filter 函数。
也可以用MATLAB 语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。
系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。
重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。
系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。
或者系统的单位脉冲响应满足绝对可和的条件。
数字信号处理实验2
数字信号处理实验2 ——离散系统频率响应和零极点分布姓名：李倩学号：班级：通信四班指导教师：周争一．实验原理离散时间系统的常系数线性差分方程：求一个系统的频率响应:H(e^jw)是以2pi为周期的连续周期复函数，将其表示成模和相位的形式：H(e^jw)=|H(e^jw)|*e^(jarg[H(e^jw)])其中|H(e^jw)|叫做振幅响应（幅度响应），频率响应的相位arg[H(e^jw)]叫做系统的相位响应。
将常系数线性差分方程的等式两边求FT,可以得到系统的频率响应与输入输出的频域关系式：H(e^jw)=Y(e^jw)/X(e^jw)将上式中的e^jw用z代替，即可得系统的系统函数：H(z)=Y(z)/X(z)H(z)=∑h(n)*z^(-n)（n的取值从负无穷到正无穷）将上式的分子、分母分别作因式分解，可得到LTI系统的零极点增益表达式为：H(z)=g∏(1-zr*z^(-1))/∏(1-pk*z^(-1))其中g为系统的增益因子，pk(k=1,2,3,…,N)为系统的极点，zr(r=1,2,3,…,M)为系统的零点。
通过系统的零极点增益表达式，可以判断一个系统的稳定性，对于一个因果的离散时间系统，若所有的极点都在单位圆内，则系统是稳定的。
二．实验内容一个LTI离散时间系统的输入输出差分方程为y(n)-1.6y(n-1)+1.28y(n-2)=0.5x(n)+0.1x(n-1)(1)编程求此系统的单位冲激响应序列，并画出其波形。
(2)若输入序列x(n)=＆(n)+2＆(n-1)+3＆(n-2)+4＆(n-3)+5＆(n-4),编程求此系统输出序列y(n)，并画出其波形。
(3)编程得到系统频响的幅度响应和相位响应并画图。
(4)编程得到系统的零极点分布图，分析系统的因果性和稳定性。
三．程序与运行结果(1)编程求此系统的单位冲激响应序列，并画出其波形。
程序：clear;N=100;b=[0.5 0.1];a=[1 -1.6 1.28];h1=impz(b,a,N); %计算系统的冲激响应序列的前N个取样点x1=[1 zeros(1,N-1)]; %生成单位冲激序列h2=filter(b,a,x1); %计算系统在输入单位冲激序列时的输出subplot(2,1,1);stem(h1);xlabel('时间序号n');ylabel('单位冲激响应序列值');title('单位冲激响应序列h1(n)');subplot(2,1,2);stem(h2);xlabel('时间序号n');ylabel('单位冲激响应序列值');title('单位冲激响应序列h2(n)');运行结果：0102030405060708090100-10125时间序号n单位冲激响应序列值单位冲激响应序列h1(n)0102030405060708090100-10125时间序号n单位冲激响应序列值单位冲激响应序列h2(n)结果说明：可以用impz 函数直接求出系统的单位冲激响应序列，也可输入单位冲激序列，用filter 函数求出系统的单位冲激响应序列，两者求得的结果相同。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《数字信号处理》实验指导书编写蒋朝辉中南大学信息科学与工程学院2014年4月目录实验一 MATLAB仿真软件的基本操作命令和使用方法 (2)实验二常见离散信号的MATLAB产生和图形显示 (8)实验三离散时间系统的时域分析 (12)实验四离散时间信号的DTFT (16)实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验内容1、帮助命令使用 help 命令，查找 sqrt （开方）函数的使用方法；2、MATLAB 命令窗口（1）在MATLAB 命令窗口直接输入命令行计算31)5.0sin(21+=πy 的值；（2）求多项式 p(x) = x3 + 2x+ 4的根；3、矩阵运算（1）矩阵的乘法已知A=[1 2;3 4]，B=[5 5;7 8]，求A^2*B（2）矩阵的行列式已知A=[1 2 3;4 5 6;7 8 9]，求A（3）矩阵的转置及共轭转置已知A=[1 2 3;4 5 6;7 8 9]，求A'已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B'（4）特征值、特征向量、特征多项式已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式；（5）使用冒号选出指定元素已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素；4、Matlab 基本编程方法（1）编写命令文件：计算1+2+…+n<2000 时的最大n 值；（2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。
5、MATLAB基本绘图命令（1）绘制余弦曲线 y=cos(t)，t∈[0，2π]（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π]（3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求：（a）线形为点划线、颜色为红色、数据点标记为加号；（b）坐标轴控制：显示范围、刻度线、比例、网络线（c）标注控制：坐标轴名称、标题、相应文本；>> clear;t=0:pi/10:4*pi;y=10*sin(t);plot(t,y);plot(t,y,'-+r');grid>> xlabel('X'),ylabel('Y');>> title('Plot:y=10*sin(t)');>> text(14,10,'完整图形');实验二常见离散信号的MATLAB产生和图形显示实验内容与步骤１. 写出延迟了np个单位的单位脉冲函数impseq，单位阶跃函数stepseq, n=ns:nf function [x,n]=impseq[np,ns,nf];function [x,n]=stepseq[np,ns,nf];2. 产生一个单位样本序列x1(n)，起点为ns= -10, 终点为nf=20, 在n0=0时有一单位脉冲并显示它。
修改程序，以产生带有延时11个样本的延迟单位样本序列x2(n)= x1(n-11)，并显示它。
>> clear;>> ns=-10;nf=20;n0=0;>> [x1,n1]=impseq(n0,ns,nf);>> subplot(1,2,1),stem(n1,x1);title('n0=0时的单位脉冲')>> np=11;>> [x2,n2]=impseq(np,ns,nf);>> subplot(1,2,2),stem(n2,x2);title('延迟11个样本后')3．产生一个序列X(n)= n(u(n)-u(n-8)), 0<=n<=20,并显示。
>> clear>> n=[0:20];>> x=n.*(stepseq(0,0,20)-stepseq(8,0,20));>> stem(n,x);4.编写序列相加，相乘，以及序列翻转、移位的函数文件 function [y,ny] = seqadd(x1,n1,x2,n2)；function [y,ny] = seqmult(x1,n1,x2,n2)；function [y,ny] = seqfold(x,nx)；function [y,ny] = seqshift(x,nx,k)；5．已知序列x=[0,1,2,3,4,3,2,1,0]，n= -5:3, 产生一个序列y(n) =2*x(n+3)+x(-n);并显示它。
>> x=[0,1,2,3,4,3,2,1,0];>> n=[-5:3];>> y=2*seqshift(x,n,3)+seqfold(x,n);stem(x,y)>> stem(n,y)6．复杂信号的产生：复杂的信号可以通过在简单信号上执行基本的运算来产生 试产生一个振幅调制信号 ，并显示出来。
)1.02cos())01.02cos(4.01()2cos())2cos(1()(n n n f n f m n y H L ⨯⋅⨯⋅+=⋅⋅+=ππππn=0:100>> n=[0:100];>> y=(1+0.4*cos(2*pi*0.01*n)).*cos(2*pi*0.1*n); >> stem(n,y)实验三 离散时间系统的时域分析实验内容与步骤１. 假定一因果系统为y(n)-0.4y(n-1)+0.75y(n-2)=2.2403x(n)+2.4908x(n-1)+2.2403x(n-2) 用MA TLAB 程序仿真该系统，输入三个不同的输入序列：)1.02c o s()(1n n x ⋅=π，)4.02cos()(2n n x ⋅=π，)(3)(221n x n x x -= 计算并并显示相应的输出)(1n y ， )(2n y 和)(n y 。
>> n=0:40; a=2; b=-3;x1=cos(2*pi*0.1*n); x2=cos(2*pi*0.4*n); x=a*x1+b*x2;num=[2.2403 2.4908 2.2403]; den=[1 -0.4 0.75];y1=filter(num,den,x1); %计算出y1(n) y2=filter(num,den,x2); %计算出y2(n) y=filter(num,den,x); %计算出y(n) stem(y1);n=0:40; a=2; b=-3;x1=cos(2*pi*0.1*n); x2=cos(2*pi*0.4*n); x=a*x1+b*x2;num=[2.2403 2.4908 2.2403]; den=[1 -0.4 0.75];y1=filter(num,den,x1); %计算出y1(n) y2=filter(num,den,x2); %计算出y2(n) y=filter(num,den,x); %计算出y(n) stem(y1);>> stem(y2);>> stem(y);２.用MA TLAB程序仿真步骤1给出的系统，对两个不同的输入序列x(n)和x(n-10)，计算并显示相应的输出序列y3(n)和y4(n)。
n=0:40;x1=2*n;num=[2.2403, 2.4908,2.2403];den=[1,-0.4,0.75];ic=[0 0]; %设置零初始条件y3=filter(num,den,x1,ic); %计算输入为x1(n)时的输出y1(n)[y,ny]=seqshift(x1,n,10)y4=filter(num,den,y,ic);subplot(2,1,1)stem(n,y3);ylabel('振幅');title('y3(n)');subplot(2,1,2)stem(ny,y4);ylabel('振幅');title('y4(n)');3．用MA TLAB 程序仿真计算下列两个有限长序列的卷积和并显示图形。
)2(2)1(3)()(1-+-+=n n n n x δδδ )3()()(2--=n u n u n xfunction[y,ny]=convwthn(x,nx,h,nh) nys=nx(1)+nh(1);nyf=nx(end)+nh(end); y=conv(x,h);ny=[nys:nyf];n=0:20;x1=impseq(0,0,20)+3*impseq(1,0,20)+2*impseq(2,0,20) x2=stepseq(0,0,20)-stepseq(3,0,20) subplot(3,1,1) stem(n,x1); subplot(3,1,2) stem(n,x2);[y,ny]=convwthn( x1,n,x2,n); subplot(3,1,3) stem(ny,y);实验四 离散时间信号的DTFT一、实验目的１. 运用MATLAB 计算离散时间系统的频率响应。
２. 运用MATLAB 验证离散时间傅立叶变换的性质。
二、实验原理（一）、计算离散时间系统的DTFT 已知一个离散时间系统∑∑==-=-Nk k Nk kk n x b k n y a)()(，可以用MA TLAB 函数frequz 非常方便地在给定的L 个离散频率点l ωω=处进行计算。
由于)(ωj e H 是ω的连续函数，需要尽可能大地选取L 的值（因为严格说，在MA TLAB 中不使用symbolic 工具箱是不能分析模拟信号的，但是当采样时间间隔充分小的时候，可产生平滑的图形），以使得命令plot 产生的图形和真实离散时间傅立叶变换的图形尽可能一致。
在MATLAB 中，freqz 计算出序列{M b b b ,,,10 }和{N a a a ,,,10 }的L 点离散傅立叶变换，然后对其离散傅立叶变换值相除得到L l e H lj ,,2,1),( =ω。
为了更加方便快速地运算，应将L 的值选为2的幂，如256或者512。
例3.1 运用MA TLAB 画出以下系统的频率响应。
y(n)-0.6y(n-1)=2x(n)+x(n-1) 程序：clf;w=-4*pi:8*pi/511:4*pi; num=[2 1];den=[1 -0.6]; h=freqz(num,den,w); subplot(2,1,1)plot(w/pi,real(h));gridtitle(‘H(e^{j\omega}的实部’)) xlabel(‘\omega/ \pi ’); ylabel(‘振幅’)；subplot(2,1,1)plot(w/pi,imag(h));gridtitle(‘H(e^{j\omega}的虚部’)) xlabel(‘\omega/ \pi ’); ylabel(‘振幅’)；（二）、离散时间傅立叶变换DTFT 的性质。