数字信号处理基础实验指导书
《数字信号处理》实验指导书
光电工程学院二○○九年十月
实验一离散时间信号分析
一、实验目的
1．掌握各种常用的序列，理解其数学表达式和波形表示。
2．掌握在计算机中生成及绘制数字信号波形的方法。
3．掌握序列的相加、相乘、移位、反转等基本运算及计算机实现与作用。
4．掌握线性卷积软件实现的方法。
5．掌握计算机的使用方法和常用系统软件及应用软件的使用。
6．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。
二、实验原理
1．序列的基本概念
离散时间信号在数学上可用时间序列来表示，其中代表序列的第n个数字，n代表时间的序列，n的取值范围为的整数，n取其它值没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号进行等间隔采样，采样间隔为T，得到一个有序的数字序列就是离散时间信号，简称序列。
2．常用序列
常用序列有：单位脉冲序列（单位抽样）、单位阶跃序列、矩形序列、实指数序列、复指数序列、正弦型序列等。
3．序列的基本运算
序列的运算包括移位、反转、和、积、标乘、累加、差分运算等。
4．序列的卷积运算
上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。
（1）反褶：先将和的变量换成，变成和，再将以纵轴为对称轴反褶成。
（2）移位：将移位，得。当为正数时，右移位；当为负数时，左
移位。
（3）相乘：将和的对应点值相乘。
（4）求和：将以上所有对应点的乘积累加起来，即得。
三、主要实验仪器及材料
微型计算机、Matlab软件6.5或更高版本。
四、实验内容
1．知识准备
认真复习以上基础理论，理解本实验所用到的实验原理。
2．离散时间信号（序列）的产生
利用MATLAB或C语言编程产生和绘制下列有限长序列：
（1）单位脉冲序列
（2）单位阶跃序列
（3）矩形序列
（4）正弦型序列
（5）任意序列
3．序列的运算
利用MATLAB编程完成上述两序列的移位、反转、加法、乘法等运算，并绘制运算后序列的波形。
4．卷积运算
利用MATLAB编制一个计算两个序列线性卷积的通用程序，计算上述两序列，并绘制卷积后序列的波形。
5．上机调试并打印或记录实验结果。
6．完成实验报告。
五、实验报告要求
1. 简述实验原理及目的。
2. 给出上述序列的实验结果。
3. 列出计算卷积的公式，画出程序框图，并列出实验程序清单
（可略）（包括必要的程序说明）。
4. 记录调试运行情况及所遇问题的解决方法。
5. 给出实验结果，并对结果做出分析。
6. 简要回答思考题。
1 如何产生方波信号序列和锯齿波信号序列?
2 实验中所产生的正弦序列的频率是多少?是否是周期序列?
实验二用FFT进行谱分析
一、实验目的
（1）进一步加深DFT算法原理和基本性质的理解（因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质）。
熟悉FFT程序结构及编程方法。
（2）熟悉应用FFT对确定信号进行谱分析方法，熟悉FFT算法原理和FFT 子程序的应用。
（3）学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应有FFT。
二、实验原理及方法
一个连续信号x a(t)的频谱可以用其傅里叶变换表示
若对x a(t)进行理想取样可得取序列为
对x[k]进行傅里叶变换，可得其频谱为
其中数字频率ω与模拟频率Ω的关系为：
(f s是取样频率)
若x a(t)是带限信号，且在满足取样定理的条件下，X(e jω)是Xa(jΩ)的周期延拓。为了能够在计算机上进行分析计算，常取X(e jω)一个周期（2π）内的样点值X[m]来近似X(e jω)。这样，对于长度为N的有限序列（无限长序列也可用有限长序列来逼近），便可通过DFT求其离散频
谱X[m],即
DFT可以用其快速算法FFT来实现。常用的FFT是以2为基数的，故应取x[k]的长度N=2M（M为正整数）。
有关此部分的详细说明参见教材第3章第3节。
三、实验内容及步骤
1. 复习DFT的定义、性质和用DFT作谱分析的有关内容。
2. 复习FFT算法原理与编程思想，并对照DFT-FFT运算流图和
程序框图，读懂本实验提供的FFT参考程序。
3. 用FFT进行频谱分析，用matlab编制信号产生子程序，产生
以下典型信号供谱分析用：
开始
读入长度N
调用信号产生子程序产生实验信号
调用绘图子程序（函数）绘制时间序列波形图
调用FFT子程序（函数）计算信号的DFT
调用绘图子程序（函数）绘制|Z(k)|曲线
结束
图1 主程序框图
4. 编写主程序
图1给出了主程序框图，供参考。本实验提供FFT子程序和通用绘图子程序。
5. 按实验内容要求，上机实验，并写出实验报告。
四、上机实验内容
1. 运行本实验提供的FFT参考程序，了解程序提供信号的频谱。
2. 用matlab产生信号x1(n)、x2(n)、x3(n)、x4(n)、x5(n)，并编制FFT
函数进行频谱分析。
五、实验报告要求
1.简述实验原理及目的。
2.对实验内容及步骤3中所给出的信号逐个进行谱分析。给出所编制的实验主程序、实验信号序列的时域的频域图形并分析所得图形，说明参数改变对时域和频域信号波形的影响。
3.简要回答以下问题：
①在N=8时，x2(n)和x3(n)的幅频特性会相同吗?为什么?N=16呢?
②FFT在什么条件下也可以用来分析周期信号序列的频谱?如果正弦信号系统sin(2πf0k)，f0=0.1Hz，用16点FFT来做DFT运算，得到的频谱是信号本身的真实谱吗?为什么?
实验三用双线性变换法设计IIR数字滤波器
一、实验目的
1. 掌握用双线性变换法设计IIR DF的原理及具体设计方法，熟悉用双线性变换法设计低通IIR DF的计算机编程。
2. 观察用双线性变换法设计的DF的频响特性，了解双线性变换法的特点。
3. 熟悉用双线性变换法设计BW和CB型DF的全过程。
二、实验原理与方法
滤波器的设计就寻找一个好的系统函数，使其实现对理想特性的最大限度逼近，而且必须是一个因果、稳定和可实现的系统函数。
目前，模拟滤波器（AF）的频域法设计理论比较成熟，借助于AF 的逼近方法，用模拟域到数字域的某种变换，来完成DF的逼近，这是一类简单而行之有效的方法。本实验仅涉及这类方法中的一种——双线性变换法。
双线性变换法的设计准则是使DF的频率响应与参考AF的频率响应相似，它依靠双线性变换式
建立起s平面和z平面的单值映射关系。数字频率Ω与模拟频率ω的关系为
双线性变换法不存在频率混叠问题，但数字域和模拟频率不是线性关系，这种非线性关系使得通带截止频率、过渡带的边缘频率的相对位置都发生了非线性畸变。这种频率的畸变需要通过预畸变来校正。
用双线性变换法设计IIR DF的步骤如下：
1、确定（或给出）DF的通、阻带截频及通带最大允许衰耗和阻带最小衰耗。
2、计算出预畸变后相应模拟低通原型的频率。（通过函数double
*bcg（）实现该计算）。
3、确定低通AF的阶数N。在double *bcg（）函数中
（BW型）
实现语句*n = (int)(fabs(log(c)/log(wp/ws)/2.0) + 0.99999)取整
（CB型）
实现语句*n = (int)(fabs(FNCCH(c/e)/FNCCH(ws/wp))+0.99999)取整
[提示： #define FNCCH (x) log(x+sqrt(x*x-1))]
4、求出低通AF的系统函数H(s)。
对于模方函数令则有：
可求出的极点，即
程序中参数x即为书中的β参数，循环
for(i=0; i<*n; i++)
{
y = (2.0*i+1.0)*PI/(2.0*(*n));
(b+i)->real = -wc*FNSH1(x)*sin(y);
(b+i)->imag = -wc*FNCH1(x)*cos(y);
} 用于计算系统函数的零极点。
[提示： #define FNSSH (x) log(x+sqrt(x*x+1))]
程序中求出了CB滤波器的系数K，其中在第一个等号后的公式，的首项系数为，因而其项的系数为，整个分母多项式的系数为，这也是程序中系数a的由来。
a = pow(wc,(*n))/(e*pow(2.0, (*n)-1));
BW滤波器的传递函数系数a为
a = pow(wc,(double)(*n));
传递函数分母的计算由以下语句完成：
for(i=0; i<*n; i++)
{
for(k=0; k { cs1 = (b1+k)->real - (b1+k+1)->real*(b+i)->real; cs2 = (b1+k)->imag - (b1+k+1)->real*(b+i)->imag; (b2+k+1)->real = cs1 + (b1+k+1)->imag*(b+i)->imag; (b2+k+1)->imag = cs2 - (b1+k+1)->imag*(b+i)->real; } b2->real = -(b1->real*(b+i)->real - b1->imag*(b+i)->imag); b2->imag = -(b1->real*(b+i)->imag + b1->imag*(b+i)->real); (b2+i+1)->real = (b1+i)->real; (b2+i+1)->imag = (b1+i)->imag; for(k=0; k<=i+1; k++) { (b1+k)->real = (b2+k)->real; (b1+k)->imag = (b2+k)->imag; (b2+k)->real = 0.0; (b2+k)->imag = 0.0; } } 由以下语句得到传递函数的表达，完成AF的设计： for(i=0; i<=*n; i++) h[i] = (b1+i)->real/a; 5、由H(s)经双线性变换求低通DF的系统函数，即 6、根据H(z)求出DF的频响，以检验所设计的DF特性是否满足要求。 三、实验内容及步骤 1、读懂所给参考程序，并熟悉程序的整体结构和功能。 2、根据给定指标，设计相应的滤波器（程序IIRDF.EXE）。 1 设计一个CB型低通DF，通带截频f p=3000Hz，衰耗满足A pmax=3dB，阻带截频f T=3400Hz，衰耗A Tmin=31dB，取样频率f s=8000Hz。 2 设计一个BW型低通DF，满足：通带截频f p=100Hz，衰耗满足 A pmax=3dB，阻带截频f T=400Hz，衰耗A Tmin=15dB,取样频率 f s=2000Hz。 3 设计一个BW型高通DF，满足：通带截频f p=400Hz，衰耗满足 A pmax=3dB，阻带截频f T=350Hz，衰耗A Tmin=15dB，取样频率 f s=1000Hz。 4 设计一个CB型带通DF，满足：通带边界频率为100Hz～500Hz，通带衰耗小于3dB，过渡带宽20Hz，阻带衰耗大于15dB，取样频率 f s=2000Hz。 3、通过MATLAB编程完成2中①~③滤波器的设计。 4、人体心电图信号在测量过程中往往受到工业高频干扰，所以必须经过低通滤波处理后，才能作为判断心脏功能的有用信息。下面程序中给出一实际心电图信号采样序列样本x(n)，其中存在高频干扰。在实验 中，以x(n)作为输入序列，滤除其中的干扰成分。我们用双线性变换法设计一个巴特沃斯低通IIR数字滤波器对上述信号进行处理。 直接运行程序，结果输出滤波器幅频特性曲线图，有噪声的心电图采集信号波形图和经过三级二阶滤波器后的心电图信号波形图，可以看出低通滤波器除信号中高频噪声的滤波效果。 四、实验报告要求 1、简述实验目的及原理。 2、给出上述四个滤波器的频响特性，并对结果加以分析。 3、通过MATLAB编程完成2中①~③滤波器的设计。 4、根据对心电信号进行处理的程序，写出用双线性变换法设计的巴特沃斯低通IIR数字滤波器的系统函数。设计指标参数为：在通带内频率低于0.2π时，最大衰减小于1dB；在阻带内[0.3π,π]频率区间上，最小衰减大于15dB。 5、与脉冲响应不变设计法相比较，简述双线性变换设计法的优缺点。 实验四用频率取样法设计FIR数字滤波器 一、实验目的 1、掌握用频率取样法设计线性相位FIR DF的方法，并掌握该方法的计算机编程。 2、熟悉频率取样理论，熟悉内插函数及其应用。 3、了解FIR DF的频率特性和相位特性，观察过渡带取样点对滤波器幅频特性的影响。 二、实验原理方法 FIR DF的系统函数H(z)为 对于理想频率特性，其N个样点值H[m]在一定条件下就可以不失真代表它，即 此时， 将代入求H(z)公式，就可得FIR DF的频响特性，即 式中称为内插函数。 用频率取样法设计FIR DF，就是对理想频率特性进行等间隔频率取样，然后由取样序列H[m]求得系统函数H(z)。 在本实验中，我们要设计的是线性相位FIR DF，因此还需要根据线性相位性质，在四种线性相位FIR DF可以归结为以下六个步骤：（1）确定所设计DF的性能要求和阶数N。 （2）根据要求，合理选择单位脉冲响应的奇、偶对称性，从而确定理想频率响应的幅度和相位。 （3）对理想频率响应在0~2π区间等间隔取样，得到H[m]。 （4）根据内插公式，求出H(z)的幅频特性曲线。 （5）检查幅频特性是否满足性能要求，若不满足，可调整过渡带取样点值，重复第（2）步，直到满足条件为止。 （6）对H[m]作IDFT求得h(k)。 三、实验内容及步骤 1、读懂所给参考程序，根据频率取样法FIR DF的设计步骤，检验程序的各个部分，了解编程思路，然后总结归纳出程序框图。 2、上机调试该程序，并利用该程序设计如下滤波器。 （1）用频率取样法设计一个线性相位LP DF，N=15，0~π之间的幅度取样值为 求出其单位脉冲响应h[k]及幅频和相频特性曲线。 （2）用频率取样法设计一个线性相位FIR LP-DF，N=30，通带频率区间为0~0.6π，阻带频率区间为0.8π~1π，通带内波动不大于0.03dB，阻带衰减不小于30dB，求出h[k]。说明选用了第几类滤波器，并给出其频幅频及相频特性曲线。观察过渡带取样点的设置对滤波器特性的影响。 3、通过MATLAB编程完成上述两个滤波器的设计。 四、实验报告要求 1、简述实验目的及原理 2、根据实验内容要求给出所设计滤波器的系统函数H(z)及相应的幅频和相频特性曲线。检查所设计的滤波器是否符合要求，说明设计过程中应考虑的因素及其对所设计滤波器性能的影响。给出利用MATLAB 设计滤波器的源码。 3、总结用频率取样法设计FIR DF的特点，提出参考程序的改进方案。 4、简要回答以下问题： （1）尽管我们要求设计的是实序列h[k]的FIR DF，可是在调试程序的过程中，有人发现他所得到的h[k]是复序列，即不但实部有值而且虚部也有值。这是什么原因?这是什么原因?问题可能出现在哪里? （2）如果我们希望用频率取样法设计一个带阻滤波器，阻带区间为0.3~0.7π,滤波器长度N是选用奇数好，还是选用偶数好?说明理由。 一、实验目的 1. 通过本次实验回忆并熟悉MATLAB这个软件。 2. 通过本次实验学会如何利用MATLAB进行序列的简单运算。 3. 通过本次实验深刻理解理论课上的数字信号处理的一个常见方法——对时刻n的样本附近的一些样本求平均，产生所需的输出信号。 3. 通过振幅调制信号的产生来理解载波信号与调制信号之间的关系。 二、实验内容 1. 编写程序在MATLAB中实现从被加性噪声污染的信号中移除噪声的算法，本次试验采用三点滑动平均算法，可直接输入程序P1.5。 2. 通过运行程序得出的结果回答习题Q1.31-Q1.33的问题，加深对算法思想的理解。 3. 编写程序在MATLAB中实现振幅调制信号产生的算法，可直接输入程序P1.6。 4. 通过运行程序得出的结果回答习题Q1.34-Q1.35的问题，加深对算法思想的理解。 三、主要算法与程序 1. 三点滑动平均算法的核心程序： %程序P1.5 %通过平均的信号平滑 clf; R=51; d=0.8*(rand(R,1)-0.5);%产生随噪声 m=0:R-1; s=2*m.*(0.9.^m);%产生为污染的信号 x=s+d';%产生被噪音污染的信号 subplot(2,1,1); plot(m,d','r-',m,s,'g--',m,x,'b-.'); xlabel('时间序号n');ylabel('振幅'); legend('d[n]','s[n]','x[n]'); x1=[0 0 x];x2=[0 x 0];x3=[x 0 0]; y=(x1+x2+x3)/3; subplot(2,1,2); plot(m,y(2:R+1),'r-',m,s,'g--'); legend('y[n]','s[n]'); xlabel('时间序号n');ylabel('振幅'); 2. 振幅调制信号的产生核心程序：（由于要几个结果，因此利用subplot函数画图） %程序P1.6 %振幅调制信号的产生 n=0:100; m=0.1;fH=0.1;fL=0.01; m1=0.3;fH1=0.3;fL1=0.03; xH=sin(2*pi*fH*n); xL=sin(2*pi*fL*n); y=(1+m*xL).*xH; xH1=sin(2*pi*fH1*n); xL1=sin(2*pi*fL1*n); y1=(1+m1*xL).*xH; y2=(1+m*xL).*xH1; y3=(1+m*xL1).*xH; subplot(2,2,1); stem(n,y); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.1;fH=0.1;fL=0.01;'); subplot(2,2,2); stem(n,y1); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.1;fL=0.01;'); subplot(2,2,3); stem(n,y2); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.3;fL=0.01;'); subplot(2,2,4); stem(n,y3); grid; 控制工程基础实验指导书 自控原理实验室编印 （内部教材） 实验项目名称: （所属课 程: 院系: 专业班级: 姓名: 学号: 实验日期: 实验地点: 合作者: 指导教师: 本实验项目成绩: 教师签字: 日期: （以下为实验报告正文） 、实验目的 简述本实验要达到的目的。目的要明确，要注明属哪一类实验（验证型、设计型、综合型、创新型）。 二、实验仪器设备 列出本实验要用到的主要仪器、仪表、实验材料等。 三、实验内容 简述要本实验主要内容，包括实验的方案、依据的原理、采用的方法等。 四、实验步骤 简述实验操作的步骤以及操作中特别注意事项。 五、实验结果 给出实验过程中得到的原始实验数据或结果，并根据需要对原始实验数据或结果进行必要的分析、整理或计算，从而得出本实验最后的结论。 六、讨论 分析实验中出现误差、偏差、异常现象甚至实验失败的原因，实验中自己发现了什么问题，产生了哪些疑问或想法，有什么心得或建议等等。 七、参考文献 列举自己在本次准备实验、进行实验和撰写实验报告过程中用到的参考文献资 料。 格式如下 作者，书名（篇名），出版社（期刊名），出版日期（刊期），页码 实验一控制系统典型环节的模拟、实验目的 、掌握比例、积分、实际微分及惯性环节的模拟方法; 、通过实验熟悉各种典型环节的传递函数和动态特性; 、了解典型环节中参数的变化对输出动态特性的影响。 二、实验仪器 、控制理论电子模拟实验箱一台; 、超低频慢扫描数字存储示波器一台; 、数字万用表一只; 、各种长度联接导线。 三、实验原理 运放反馈连接 基于图中点为电位虚地，略去流入运放的电流，则由图 由上式可以求得下列模拟电路组成的典型环节的传递函数及其单位阶跃响应。 、比例环节 实验模拟电路见图所示 U i R i U o 接示波器 以运算放大器为核心元件，由其不同的输入网络和反馈网络组成的各种典型环节，如图所示。图中和为复数阻抗，它们都是构成。 Z2 Z1 Ui ,— U o 接示波器 得: 《数字信号处理》实验指导书 光电工程学院二○○九年十月 实验一离散时间信号分析 一、实验目的 1．掌握各种常用的序列，理解其数学表达式和波形表示。 2．掌握在计算机中生成及绘制数字信号波形的方法。 3．掌握序列的相加、相乘、移位、反转等基本运算及计算机实现与作用。 4．掌握线性卷积软件实现的方法。 5．掌握计算机的使用方法和常用系统软件及应用软件的使用。 6．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列来表示，其中代表序列的第n个数字，n代表时间的序列，n的取值范围为的整数，n取其它值没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号进行等间隔采样，采样间隔为T，得到一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）、单位阶跃序列、矩形序列、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反转、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将和的变量换成，变成和，再将以纵轴为对称轴反褶成。 （2）移位：将移位，得。当为正数时，右移位；当为负数时，左 移位。 （3）相乘：将和的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得。 三、主要实验仪器及材料 微型计算机、Matlab软件6.5或更高版本。 四、实验内容 1．知识准备 认真复习以上基础理论，理解本实验所用到的实验原理。 2．离散时间信号（序列）的产生 利用MATLAB或C语言编程产生和绘制下列有限长序列： （1）单位脉冲序列 （2）单位阶跃序列 （3）矩形序列 （4）正弦型序列 （5）任意序列 3．序列的运算 利用MATLAB编程完成上述两序列的移位、反转、加法、乘法等运算，并绘制运算后序列的波形。 4．卷积运算 利用MATLAB编制一个计算两个序列线性卷积的通用程序，计算上述两序列，并绘制卷积后序列的波形。 5．上机调试并打印或记录实验结果。 6．完成实验报告。 五、实验报告要求 1. 简述实验原理及目的。 2. 给出上述序列的实验结果。 3. 列出计算卷积的公式，画出程序框图，并列出实验程序清单 （可略）（包括必要的程序说明）。 4. 记录调试运行情况及所遇问题的解决方法。 5. 给出实验结果，并对结果做出分析。 6. 简要回答思考题。 1 如何产生方波信号序列和锯齿波信号序列? 2 实验中所产生的正弦序列的频率是多少?是否是周期序列? 实验一 离散时间信号分析 班级 信息131班 学号 201312030103 姓名 陈娇 日期 一、实验目的 掌握两个序列的相加、相乘、移位、反褶、卷积等基本运算。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列)}({n x 来表示，其中)(n x 代表序列的第n 个数字，n 代表时间的序列，n 的取值范围为+∞<<∞-n 的整数，n 取其它值)(n x 没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号)(t x a 进行等间隔采样，采样间隔为T ，得到)}({nT x a 一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）） （n δ、单位阶跃序列)(n u 、矩形序列)(n R N 、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 ∑∞ -∞==-= m n h n x m n h m x n y )(*)()()()( 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将)(n x 和)(n h 的变量n 换成m ，变成)(m x 和)(m h ，再将)(m h 以纵轴为对称轴反褶成)(m h -。 （2）移位：将)(m h -移位n ，得)(m n h -。当n 为正数时，右移n 位；当n 为负数时，左移n 位。 （3）相乘：将)(m n h -和)(m x 的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得)(n y 。 三、主要实验仪器及材料 微型计算机、Matlab6.5 教学版、TC 编程环境。 四、实验内容 （1）用Matlab 或C 语言编制两个序列的相加、相乘、移位、反褶、卷积等的程序； （2）画出两个序列运算以后的图形； （3）对结果进行分析； （4）完成实验报告。 五、实验结果 六、实验总结 实验二二阶系统的瞬态响应分析 一、实验目的 1、熟悉二阶模拟系统的组成。 2、研究二阶系统分别工作在ξ=1，0<ξ<1，和ξ> 1三种状态下的单 位阶跃响应。 3、分析增益K对二阶系统单位阶跃响应的超调量σP、峰值时间tp和调 整时间ts。 4、研究系统在不同K值时对斜坡输入的稳态跟踪误差。 5、学会使用Matlab软件来仿真二阶系统，并观察结果。 二、实验仪器 1、控制理论电子模拟实验箱一台； 2、超低频慢扫描数字存储示波器一台； 3、数字万用表一只； 4、各种长度联接导线。 三、实验原理 图2-1为二阶系统的原理方框图，图2-2为其模拟电路图，它是由惯性环节、积分环节和反号器组成，图中K=R2/R1，T1=R2C1，T2=R3C2。 图2-1 二阶系统原理框图 图2-1 二阶系统的模拟电路 由图2-2求得二阶系统的闭环传递函 12 22 122112 /() (1)()/O i K TT U S K U S TT S T S K S T S K TT ==++++ :而二阶系统标准传递函数为 (1)(2), 对比式和式得 n ωξ== 12 T 0.2 , T 0.5 , n S S ωξ====若令则。调节开环增益K 值，不仅能改变系统无阻尼自然振荡频率ωn 和ξ的值，可以得到过阻尼(ξ>1)、 临界阻尼（ξ=1）和欠阻尼（ξ<1）三种情况下的阶跃响应曲线。 （1）当K >0.625， 0 < ξ < 1，系统处在欠阻尼状态，它的单位阶跃响应表达式为： 图2-3 0 < ξ < 1时的阶跃响应曲线 （2）当K =0.625时，ξ=1，系统处在临界阻尼状态，它的单位阶跃响应表达式为： 如图2-4为二阶系统工作临界阻尼时的单位响应曲线。 (2) +2+=222n n n S S )S (G ωξω ω1 ()1sin( ) (3) 2-3n t o d d u t t tg ξωωωω--=+=式中图为二阶系统在欠阻尼状态下的单位阶跃响应曲线 e t n o n t t u ωω-+-=)1(1)( 实验一 离散傅里叶变换（DFT ）对确定信号进行谱分析 一．实验目的 1．加深对DFT 算法原理和基本性质的理解。 2．熟悉DFT 算法和原理的编程方法。 3．学习用DFT 对信号进行谱分析的方法，了解可能出现的误差及其原因，以便在实际中正确利用。 二．实验原理 一个连续信号)(t x a 的频谱可以用其傅里叶变换表示，即 dt e t x j X t j a a Ω-∞ ∞ -? = Ω)()( 若对)(t x a 进行理想采样可得采样序列 )(|)()(nT x t x n x a nT t a === 对)(n x 进行DTFT ，可得其频谱为： ∑∞ -∞ =-= n n j j e n x e X ωω )()( 其中数字频率ω与模拟频率Ω的关系为： s f T Ω = Ω=ω )(n x 的DFT 为∑∞ -∞ =-= n nk N j e n x k X π 2)()( 若)(t x a 是限带信号，且在满足采样定理的条件下，)(ω j e X 是)(Ωj X a 的周期延拓， )(k X 是)(ωj e X 在单位圆上的等间隔采样值，即k N j e X k X πωω2| )()(= =。 为在计算机上分析计算方便，常用)(k X 来近似)(ω j e X ，这样对于长度为N 的有限 长序列（无限长序列也可用有限长序列来逼近），便可通过DFT 求其离散频谱。 三．实验内容 1．用DFT 对下列序列进行谱分析。 （1）)()04.0sin(3)(100n R n n x π= 1 （2）]0,0,0,0,0,0,0,0,1,1,1,1[)(=n x 2．为了说明高密度频谱和高分辨率频谱之间的区别，考察序列 )52.0cos()48.0cos()(n n n x ππ+= （1）当0≤n ≤10时，确定并画出x(n)的离散傅里叶变换。 （2）当0≤n ≤100时，确定并画出x(n)的离散傅里叶变换。 四．实验结果 1. （1） （2） 实验一 连续时间系统的时域和频域分析相关MATLAB 函数1.设描述连续时间系统的微分方程为：)()()()()()()()(01)1(1)(01)1(1)(t f b t f b t f b t f b t y a t y a t y a t y a m m m m n n n n +'+++=+'+++---- 则可用向量和表示该系统，即 a b ] ,,,,[011a a a a a n n -=],,,,[011b b b b b m m -=注意，向量和的元素一定要以微分方程时间求导的降幂次序排列，且缺项要用0补齐。a b 如微分方程)()()(2)(3)(t f t f t y t y t y +''=+'+''表示该系统的向量为 ]2 3 1[=a ]1 0 1[=b （1）求解冲激响应：impulse()函数impulse()函数有以下四种调用格式： ① impulse(b,a) 该调用格式以默认方式绘制由向量和定义的连续时间系统的冲激响应的时域波形。a b ② impulse(b,a,t)该调用格式绘制由向量和定义的连续时间系统在时间范围内的冲激响应的时a b t ~0域波形。③ impulse(b,a, t1:p:t2)该调用格式绘制由向量和定义的连续时间系统在时间范围内，且以时间间a b 21~t t 隔均匀抽样的冲激响应的时域波形。p ④ y=impulse(b,a,t1:p:t2)该调用格式并不绘制系统冲激响应的波形，而是求出由向量和定义的连续时间系a b 统在时间范围内以时间间隔均匀抽样的系统冲激响应的数值解。21~t t p （2）求解阶跃响应：step()函数 step()函数也有四种调用格式：① step(b,a) ② step(b,a,t) ③ step(b,a, t1:p:t2) ④ y=step(b,a,t1:p:t2) 上述调用格式的功能与impulse()函数完全相同。 实验一 Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容 1、帮助命令 使用help命令，查找sqrt（开方）函数的使用方法； 2、矩阵运算 （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]; B=[1 0 0;0 2 0;0 0 3]; A\B,A/B （3）矩阵的转置及共轭转置 已知A=[5+i,2-i,1;6*i,4,9-i]; 求A.', A' （4）使用冒号选出指定元素 已知：A=[1 2 3;4 5 6;7 8 9]; 求A中第3列前2个元素；A中所有列第2，3行的元素； （5）方括号[] 用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列 3、多项式 （1）求多项式p(x) = x3 - 2x - 4的根 （2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ， 求矩阵A的特征多项式; 求特征多项式中未知数为20时的值； 4、基本绘图命令 （1）绘制余弦曲线y=cos(t)，t∈[0，2π] （2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π] 5、基本绘图控制 绘制[0，4π]区间上的x1=10sint曲线，并要求： （1）线形为点划线、颜色为红色、数据点标记为加号； （2）坐标轴控制：显示范围、刻度线、比例、网络线 （3）标注控制：坐标轴名称、标题、相应文本； 6、基本程序设计 （1）编写命令文件：计算1+2+?+n<2000时的最大n值； （2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。 三、预习要求 利用所学知识，编写实验内容中2到6的相应程序，并写在预习报告上。 数字信号处理实验答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】 实验一熟悉Matlab环境 一、实验目的 1.熟悉MATLAB的主要操作命令。 2.学会简单的矩阵输入和数据读写。 3.掌握简单的绘图命令。 4.用MATLAB编程并学会创建函数。 5.观察离散系统的频率响应。 二、实验内容 认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。在熟悉了MATLAB基本命令的基础上，完成以下实验。 上机实验内容： （1）数组的加、减、乘、除和乘方运算。输入A=[1 2 3 4]，B=[3 4 5 6]，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B并用stem语句画出A、B、C、D、E、F、G。 clear all; a=[1 2 3 4]; b=[3 4 5 6]; c=a+b; d=a-b; e=a.*b; f=a./b; g=a.^b; n=1:4; subplot(4,2,1);stem(n,a); xlabel('n');xlim([0 5]);ylabel('A'); subplot(4,2,2);stem(n,b); xlabel('n');xlim([0 5]);ylabel('B'); subplot(4,2,3);stem(n,c); xlabel('n');xlim([0 5]);ylabel('C'); subplot(4,2,4);stem(n,d); xlabel('n');xlim([0 5]);ylabel('D'); subplot(4,2,5);stem(n,e); xlabel('n');xlim([0 5]);ylabel('E'); subplot(4,2,6);stem(n,f); xlabel('n');xlim([0 5]);ylabel('F'); subplot(4,2,7);stem(n,g); xlabel('n');xlim([0 5]);ylabel('G'); （2）用MATLAB实现下列序列： a) x(n)= 0≤n≤15 b) x(n)=e+3j)n 0≤n≤15 c) x(n)=3cosπn+π)+2sinπn+π) 0≤n≤15 d) 将c)中的x(n)扩展为以16为周期的函数x(n)=x(n+16),绘出四个周期。 武汉工程大学 数字信号处理实验报告 姓名：周权 学号：1204140228 班级：通信工程02 一、实验设备 计算机，MATLAB语言环境。 二、实验基础理论 1.序列的相关概念 2.常见序列 3.序列的基本运算 4.离散傅里叶变换的相关概念 5.Z变换的相关概念 三、实验内容与步骤 1.离散时间信号（序列）的产生 利用MATLAB语言编程产生和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形表示。 四实验目的 认识常用的各种信号，理解其数字表达式和波形表示，掌握在计算机中生成及绘制数字信号波形的方法，掌握序列的简单运算及计算机实现与作用，理解离散时间傅里叶变换，Z变换及它们的性质和信号的频域分 实验一离散时间信号（序列）的产生 代码一 单位样值 x=2; y=1; stem(x,y); title('单位样值 ') 单位阶跃序列 n0=0; n1=-10; n2=10; n=[n1:n2]; x=[(n-n0)>=0]; stem(n,x); xlabel('n'); ylabel('x{n}'); title('单位阶跃序列'); 实指数序列 n=[0:10]; x=(0.5).^n; stem(n,x); xlabel('n'); ylabel('x{n}'); title('实指数序列'); 正弦序列 n=[-100:100]; x=2*sin(0.05*pi*n); stem(n,x); xlabel('n'); ylabel('x{n}'); title('正弦序列'); 随机序列 n=[1:10]; x=rand(1,10); subplot(221); stem(n,x); xlabel('n'); ylabel('x{n}'); title('随机序列'); 实验三：离散LSI 系统的频域分析 一、实验内容 2、求以下各序列的z 变换： 12030() ()sin() ()sin()n an x n na x n n x n e n ωω-=== 程序清单如下： syms w0 n z a; x1=n*a^n;X1=ztrans(x1) x2=sin(w0*n);X2=ztrans(x2) x3= exp(-a*n)*sin(w0*n);X3=ztrans(x3) 程序运行结果如下： X1 =z/(a*(z/a - 1)^2) X2 =(z*sin(w0))/(z^2 - 2*cos(w0)*z + 1) X3 =(z*exp(a)*sin(w0))/(exp(2*a)*z^2 - 2*exp(a)*cos(w0)*z + 1) 3、求下列函数的逆z 变换 0 312342 1 1() () () ()() 1j z z z z X z X z X z X z z a z a z e z ω---= = = = ---- 程序清单如下： syms w0 n z a; X1=z/(z-a);x1=iztrans(X1) X2= z/(a-z)^2;x2=iztrans(X2) X3=z/ z-exp(j*w0);x3=iztrans(X3) X4=(1-z^-3)/(1-z^-1);x4=iztrans(X4) 程序运行结果如下： x1 =a^n x2 =n*a^n/a 课程名称 数字信号 实验成绩 指导教师 实 验 报 告 院系 信息工程学院 班级 学号 姓名 日期 x3 =charfcn[0](n)-iztrans(exp(i*w0),w0,n) x4 =charfcn[2](n)+charfcn[1](n)+charfcn[0](n) 4、求一下系统函数所描述的离散系统的零极点分布图，并判断系统的稳定性 （1） (0.3)()(1)(1) z z H z z j z j -= +-++ z1=[0,0.3]';p1=[-1+j,-1-j]';k=1; [b1,a1]=zp2tf(z1,p1,k); subplot(1,2,1);zplane(z1,p1); title('极点在单位圆外); subplot(1,2,2);impz(b1,a1,20); 由图可见:当极点位于单位圆内，系统的单位序列响应随着频率的增大而收敛；当极点位于单位圆上，系统的单位序列响应为等幅振荡；当极点位于单位圆外，系统的单位序列响应随着频率的增大而发散。由此可知系统为不稳定系统。 -1 -0.5 00.51 -2 -1.5-1-0.500.511.5 2Real Part I m a g i n a r y P a r t 极点在单位圆外 n (samples) A m p l i t u d e Impulse Response 《控制工程基础》实验指导书常熟理工学院机械工程学院 2009.9 目录 1.MATLAB时域分析实验 (2) 2.MATLAB频域分析实验 (4) 3.Matlab校正环节仿真实验 (8) 4.附录：Matlab基础知识 (14) 实验1 MATLAB 时域分析实验 一、实验目的 1． 利用MATLAB 进行时域分析和仿真。 要求：(1)计算连续系统的时域响应（单位脉冲输入，单位阶跃输入，任意输入）。 2．掌握Matlab 系统分析函数impulse 、step 、lsim 、roots 、pzmap 的应用。 二、实验内容 1．已知某高阶系统的传递函数为 ()265432 220501584223309240100 s s G s s s s s s s ++=++++++，试求该系统的单位脉冲响应、单位阶跃响应、单位速度响应和单位加速度响应。 MATLAB 计算程序 num=[2 20 50]; den=[1 15 84 223 309 240 100]; t= (0: 0.1: 20); figure (1); impulse (num,den,t); %Impulse Response figure (2); step(num,den,t);%Step Response figure (3); u1=(t); %Ramp.Input hold on; plot(t,u1); lsim(num,den,u1,t); %Ramp. Response gtext(‘t’); figure (4); u2=(t.*t/2);%Acce.Input u2=(0.5*(t.*t)) hold on; plot(t,u2); lsim(num,den,u2,t);%Acce. Response 实验一 基于Matlab 的数字信号处理基本操作 一、 实验目的：学会运用MA TLAB 表示的常用离散时间信号；学会运用MA TLAB 实现离 散时间信号的基本运算。 二、 实验仪器：电脑一台，MATLAB6.5或更高级版本软件一套。 三、 实验内容： （一） 离散时间信号在MATLAB 中的表示 离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。离散序列通常用)(n x 来表示，自变量必须是整数。 离散时间信号的波形绘制在MATLAB 中一般用stem 函数。stem 函数的基本用法和plot 函数一样，它绘制的波形图的每个样本点上有一个小圆圈，默认是空心的。如果要实心，需使用参数“fill ”、“filled ”，或者参数“.”。由于MATLAB 中矩阵元素的个数有限，所以MA TLAB 只能表示一定时间范围内有限长度的序列；而对于无限序列，也只能在一定时间范围内表示出来。类似于连续时间信号，离散时间信号也有一些典型的离散时间信号。 1. 单位取样序列 单位取样序列)(n δ，也称为单位冲激序列，定义为 ) 0() 0(0 1)(≠=?? ?=n n n δ 要注意，单位冲激序列不是单位冲激函数的简单离散抽样，它在n =0处是取确定的值1。在MATLAB 中，冲激序列可以通过编写以下的impDT .m 文件来实现，即 function y=impDT(n) y=(n==0); %当参数为0时冲激为1,否则为0 调用该函数时n 必须为整数或整数向量。 【实例1-1】 利用MATLAB 的impDT 函数绘出单位冲激序列的波形图。 解：MATLAB 源程序为 >>n=-3:3; >>x=impDT(n); >>stem(n,x,'fill'),xlabel('n'),grid on >>title('单位冲激序列') >>axis([-3 3 -0.1 1.1]) 程序运行结果如图1-1所示。 图1-1 单位冲激序列 数字信号处理实验四 第一题结果： （1）没有增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 %H(3,13) = 0.75;H(5,11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线 （2）增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 H(3) = 0.75;H(13) = 0.75;H(5) = 0.25;H(11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线 第二题结果： 实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验容 1、帮助命令 使用 help 命令，查找 sqrt （开方）函数的使用方法； 2、MATLAB 命令窗口 （1）在MATLAB 命令窗口直接输入命令行计算3 1)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; 机械控制工程基础实验 指导书 Company number【1089WT-1898YT-1W8CB-9UUT-92108】 河南机电高等专科学校《机械控制工程基础》 实验指导书 专业：机械制造与自动化、起重运输机械设计与制造等 机械制造与自动化教研室编 2012年12月 目录 实验任务和要求 一、自动控制理论实验的任务 自动控制理论实验是自动控制理论课程的一部分，它的任务是： 1、通过实验进一步了解和掌握自动控制理论的基本概 念、控制系统的分析方法和设计方法； 2、重点学习如何利用MATLAB工具解决实际工程问题和 计算机实践问题； 3、提高应用计算机的能力及水平。 二、实验设备 1、计算机 2、MATLAB软件 三、对参加实验学生的要求 1、阅读实验指导书，复习与实验有关的理论知识，明确每次实验的目的，了解内容和方法。 2、按实验指导书要求进行操作；在实验中注意观察，记录有关数据和图 像，并由指导教师复查后才能结束实验。 3、实验后关闭电脑，整理实验桌子，恢复到实验前的情况。 4、认真写实验报告，按规定格式做出图表、曲线、并分析实验结果。字迹 要清楚，画曲线要用坐标纸，结论要明确。 5、爱护实验设备，遵守实验室纪律。 实验模块一 MATLAB基础实验 ——MATLAB环境下控制系统数学模型的建立 一、预备知识 的简介 MATLAB为矩阵实验室（Matrix Laboratory）的简称，由美国MathWorks公司出品的商业数学软件。主要用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 来源：20世纪70年代，美国新墨西哥大学计算机科学系主任Cleve Moler 为了减轻学生编程的负担，用FORTRAN编写了最早的MATLAB。1984年由 Little、Moler、Steve Bangert合作成立了的MathWorks公司正式把MATLAB推向市场。到20世纪90年代，MATLAB已成为国际控制界的标准计算软件。 地位：和Mathematica、Maple并称为三大数学软件，在数学类科技应用软件中，在数值计算方面首屈一指。 功能：矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等。 应用范围：工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 图1-1 MATLAB图形处理示例 的工作环境 启动MATLAB，显示的窗口如下图所示。 MATLAB的工作环境包括菜单栏、工具栏以及命令运行窗口区、工作变量区、历史指令区、当前目录窗口和M文件窗口。 (1)菜单栏用于完成基本的文件输入、编辑、显示、MATLAB工作环境交互性设置等操作。 (2)命令运行窗口“Command Window”是用户与MATLAB交互的主窗口。窗口中的符号“》”表示MATLAB已准备好，正等待用户输入命令。用户可以在“》”提示符后面输入命令，实现计算或绘图功能。 说明：用户只要单击窗口分离键，即可独立打开命令窗口，而选中命令窗口中Desktop菜单的“Dock Command Window”子菜单又可让命令窗口返回桌面（MATLAB桌面的其他窗口也具有同样的操作功能）；在命令窗口中，可使用方向 本科生实验报告 实验课程数字信号处理基础 学院名称地球物理学院 专业名称地球物理学 学生姓名 学生学号 指导教师王山山 实验地点5417 实验成绩 二〇一四年十一月二〇一四年十二月 填写说明 1、适用于本科生所有的实验报告（印制实验报告册除外）； 2、专业填写为专业全称，有专业方向的用小括号标明； 3、格式要求： ①用A4纸双面打印（封面双面打印）或在A4大小纸上用蓝黑色水笔书写。 ②打印排版：正文用宋体小四号，1.5倍行距，页边距采取默认形式（上下2.54cm， 左右2.54cm，页眉1.5cm，页脚1.75cm）。字符间距为默认值（缩放100%，间距：标准）；页码用小五号字底端居中。 ③具体要求： 题目（二号黑体居中）； 摘要（“摘要”二字用小二号黑体居中，隔行书写摘要的文字部分，小4号宋体）； 关键词（隔行顶格书写“关键词”三字，提炼3-5个关键词，用分号隔开，小4号黑体)； 正文部分采用三级标题； 第1章××(小二号黑体居中，段前0.5行) 1.1 ×××××小三号黑体×××××（段前、段后0.5行） 1.1.1小四号黑体（段前、段后0.5行） 参考文献（黑体小二号居中，段前0.5行），参考文献用五号宋体，参照《参考文献著录规则（GB/T 7714－2005）》。 实验一生成离散信号并计算其振幅谱 并将信号进行奇偶分解 一、实验原理 单位脉冲响应h(t)=exp(-a*t*t)*sin(2*3.14*f*t)进行离散抽样，分别得到t=0.002s，0.009s，0.011s采样的结果。用Excel软件绘图显示计算结果。并将信号进行奇偶分解，分别得到奇对称信号h(n)-h(-n)与偶对称信号h(n)+h(-n)。用Excel 软件绘图显示计算结果。 二、实验程序代码 （1）离散抽样 double a,t; a=2*f*f*log(m); int i; for(i=0;i 实验1认识实验-MATLAB语言上机操作实践 一、实验目的 ㈠了解MATLAB语言的主要特点、作用。 ㈡学会MATLAB主界面简单的操作使用方法。 ㈢学习简单的数组赋值、运算、绘图、流程控制编程。 二、实验原理 ㈠简单的数组赋值方法 MATLAB中的变量和常量都可以是数组（或矩阵）,且每个元素都可以是复数。 在MATLAB指令窗口输入数组A=[1 2 3；4 5 6；7 8 9]，观察输出结果。然后，键入：A（4，2）= 11 键入：A (5，：) = [-13 -14 -15] 键入：A（4，3）= abs (A(5，1)) 键入：A ([2，5]，：) = [ ] 键入：A/2 键入：A (4，：) = [sqrt(3) (4+5)/6*2 –7] 观察以上各输出结果。将A式中分号改为空格或逗号，情况又如何？请在每式的后面标注其含义。 2．在MATLAB指令窗口输入B=[1+2i，3+4i；5+6i ，7+8i], 观察输出结果。 键入：C=[1，3；5，7]+[2，4；6，8]*i，观察输出结果。 如果C式中i前的*号省略，结果如何？ 键入：D = sqrt (2+3i) 键入：D*D 键入：E = C’, F = conj(C), G = conj(C)’ 观察以上各输出结果, 请在每式的后面标注其含义。 3．在MATLAB指令窗口输入H1=ones(3,2)，H2=zeros(2,3)，H3=eye(4)，观察输出结果。 ㈡、数组的基本运算 1．输入A=[1 3 5]，B= [2 4 6]，求C=A+B，D=A-2，E=B-A 2．求F1=A*3，F2=A.*B，F3=A./B，F4=A.\B, F5=B.\A, F6=B.^A, F7=2./B, F8=B.\2 *3．求B'，Z1=A*B’，Z2=B’*A 观察以上各输出结果,比较各种运算的区别，理解其含义。 ㈢、常用函数及相应的信号波形显示 例1：显示曲线f(t)=2sin(2πt)，（t>0） ⅰ点击空白文档图标（New M-file），打开文本编辑器。 ⅱ键入：t=0:0.01:3; (1) f=2*sin(2*pi*t); (2) plot(t,f); title(‘f（t）-t曲线’)； xlabel（‘t’），ylabel（‘f（t）’）； 前言 《机电系统控制基础》既是一门理论性较强、又紧密联系工程实际的实践性较强的课程，本课程的重点在于培养学生对机电系统进行建模、分析与控制的能力。难点在于如何使机电类专业的学生结合工程实际，特别是结合机械工程实际，从整体分析系统的动态行为，理解和掌握略显深奥、难懂的经典控制理论，并应用经典控制论中的基本概念和基本方法来分析、研究和解决机械工程中的实际问题。 通过实验教学环节使学生验证课堂教学的理论，使学生能够建立机电系统控制的整体概念，加深对经典控制论中基本概念和基本方法的理解，并掌握其在分析、研究和解决实际机械工程控制问题中的应用。通过三方面的实验：原理性仿真实验，面向机电系统中典型物理对象/系统的特性测试与分析实验，和典型机电系统的控制三方面实验。将所学的课程内容融会贯通，培养学生分析和解决问题的能力。 1 机电系统控制基础原理性仿真实验 1.1 实验目的 通过仿真实验，掌握在典型激励作用下典型机电控制系统的时间响应特性，分析系统开环增益、系统阻尼、系统刚度、负载、无阻尼自振频率等机电参数对响应、超调量、峰值时间、调整时间、以及稳态跟踪误差的影响；掌握系统开环传递函数的各参数辨识方法，最后，学会使用matlab软件对机电系统进行仿真，加深理解系统动态响应特性与系统各参数的关系。 1.2系统典型输入的响应实验 1.2.1 实验原理 1．一阶系统的单位脉冲响应 惯性环节（一阶系统）单位脉冲响应simulink实现图，如图1-1所示 (a)可观测到输出曲线 (b)输入、输出曲线均可观测到 图1-1惯性环节（一阶系统）单位脉冲响应simulink实现图 2．一阶系统的单位阶跃响应 一阶系统的单位阶跃响应simulink实现图如图1-2所示。数字信号处理实验一
控制工程基础实验指导书(答案)
数字信号处理基础实验指导书
数字信号处理实验一
控制工程基础实验指导书(答案) 2..
数字信号处理实验
数字信号处理实验1,2,3,4
《控制系统计算机仿真》实验指导书
数字信号处理实验答案完整版
数字信号处理实验报告一
数字信号处理实验三
控制工程-实验指导书-修订版
实验一 基于Matlab的数字信号处理基本
数字信号处理实验4
数字信号处理实验报告(实验1_4)
机械控制工程基础实验指导书
数字信号处理基础实验报告_
数字信号处理实验1认识实验
哈工大_机电系统控制基础实验_指导书