备注：（1）、按照要求独立完成实验内容。
（2）、实验结束后，把电子版实验报告按要求格式改名（例：09号_张三_实验七.doc）后，实验室统一刻盘留档。
实验三零极点分布对系统频率响应的影响
一、实验目的
学习用分析零极点分布的几何方法分析研究信号和系统频率响应。
二、实验原理
如果知道信号的Z变换以及系统的系统函数H(z)，可以得到它们的零极点分布，由零极点分布可以很方便地对它们的频率响应进行定性分析。信号的幅度特性由零点矢量长度之积除以极点矢量的长度之积，当频率ω从0变化到2π时，观察零点矢量长度和极点矢量长度的变化，重点观察那些矢量长度较短的情况。另外, 由分析知道, 极点主要影响频率响应的峰值，极点愈靠近单位圆，峰值愈尖锐；零点主要影响频率特性的谷值，零点愈靠近单位圆，谷值愈深，如果零点在单位圆上，那么频率特性为零。根据这些规律可以定性画出频率响应的幅度特性。
峰值频率和谷值频率可以近似用响应的极点和零点的相角表示，例如极点
z1=0.9ejπ/4，峰值频率近似为π/4，极点愈靠近单位圆，估计法结果愈准确。
本实验借助计算机分析信号和系统的频率响应，目的是掌握用极、零点分布的几何分析法分析频率响应，实验时需要将z=ejω代入信号的Z变换和系统函数中，再在0~2π之间，等间隔选择若干点，并计算它的频率响应。
三、实验内容（包括代码与产生的图形）
要求：不仅打印幅度特性曲线，而且要有系统频率特性的文字分析。
1. 假设系统用下面差分方程描述：
y(n)=x(n)+ay(n-1)
假设a=0.7, 0.8, 0.9 ，分别在三种情况下分析系统的频率特性，并打印幅度特性曲线。a=0.7
代码：B=1;a=0.7
A=[1,-a];
subplot(3,1,3);zplane(B,A);
xlabel('êμ2?Re');
ylabel('Dé2?Im');
title('y(n)=x(n)-ay(n-1)′?ê?oˉêyá??￠??μ?·?2?');
grid on
[H,w]=freqz(B,A,'whole');
subplot(3,1,2);
plot(w/pi,abs(H),'linewidth',2);
grid on;
xlabel('\omega/\pi');
ylabel('|H(e^j^\omega)|');
title('·ù?μ?ìó|ì?D?');
axis([0,2,0,6]);
subplot(3,1,1);
plot(w/pi,angle(H),'linewidth',2);
grid on;
axis([-0.1,2.1,-3,3]);
xlabel('\omega/\pi');
ylabel('\phi(\omega)');
title('?à?μ?ìó|ì?D?');
图像：
-5
05
-101
实部Re
虚部I m
y(n)=x(n)-ay(n-1)传输函数零、极点分布
00.20.4
0.60.81 1.2 1.4 1.6 1.82
5
ω/π
|H (e j ω)|
幅频响应特性
0.2
0.4
0.6
0.8
1 1.
2 1.4
1.6
1.8
2
-2
2
ω/π
φ(ω)
a=0.8
代码：B=1;a=0.8
A=[1,-a];
subplot(3,1,3);zplane(B,A); xlabel('êμ2?Re'); ylabel('Dé2?Im');
title('y(n)=x(n)-ay(n-1)′?ê?oˉêyá??￠??μ?·?2?'); grid on
[H,w]=freqz(B,A,'whole'); subplot(3,1,2);
plot(w/pi,abs(H),'linewidth',2); grid on ;
xlabel('\omega/\pi'); ylabel('|H(e^j^\omega)|'); title('·ù?μ?ìó|ì?D?'); axis([0,2,0,6]);
subplot(3,1,1);
plot(w/pi,angle(H),'linewidth',2); grid on ;
axis([-0.1,2.1,-3,3]);
xlabel('\omega/\pi'); ylabel('\phi(\omega)'); title('?à?μ?ìó|ì?D?');
图像：
-6
-4
-2
02
4
6
实部Re
虚部I m
y(n)=x(n)-ay(n-1)传输函数零、极点分布
00.20.4
0.60.81 1.2 1.4 1.6 1.8
2
5
ω/π
|H (e j ω)|
幅频响应特性
0.2
0.4
0.6
0.8
1 1.
2 1.4
1.6
1.8
2
-2
02
ω/π
φ(ω)
相频响应特性
a=0.9
代码：B=1;a=0.9
A=[1,-a];
subplot(3,1,3);zplane(B,A); xlabel('êμ2?Re'); ylabel('Dé2?Im');
title('y(n)=x(n)-ay(n-1)′?ê?oˉêyá??￠??μ?·?2?'); grid on
[H,w]=freqz(B,A,'whole'); subplot(3,1,2);
plot(w/pi,abs(H),'linewidth',2); grid on ;
xlabel('\omega/\pi'); ylabel('|H(e^j^\omega)|'); title('·ù?μ?ìó|ì?D?'); axis([0,2,0,6]);
subplot(3,1,1);
plot(w/pi,angle(H),'linewidth',2); grid on ;
axis([-0.1,2.1,-3,3]); xlabel('\omega/\pi'); ylabel('\phi(\omega)'); title('?à?μ?ìó|ì?D?');
图像：
-5
05
-101
实部Re
虚部I m
y(n)=x(n)-ay(n-1)传输函数零、极点分布
00.20.4
0.60.81 1.2 1.4 1.6 1.82
5
ω/π
|H (e j ω)|
幅频响应特性
0.2
0.4
0.6
0.8
1 1.
2 1.4
1.6
1.8
2
-2
2
ω/π
φ(ω)
相频响应特性
分析：由y (n )=x (n )+ay (n -1)可知:H[z]=B[z]/A[z]=1/(1-az^(-1))
系统极点z=a ，零点z=0，当B 点从w=0逆时针旋转时，在w=0点，由于极点向量长度最短，形成波峰,并且当a 越大，极点越接近单位圆，峰值愈高愈尖锐；在w=pi 点形成波谷；z=0处零点不影响幅频响应。
2. 假设系统用下面差分方程描述： y (n ) = x (n ) +ax (n -1)
假设a =0.7, 0.8, 0.9 ， 分别在三种情况下分析系统的频率特性，并打印幅度特性曲线。
a =0.7 代码：A=1;a=0.7
B=[1,a];
subplot(3,1,3);zplane(B,A); xlabel('êμ2?Re');
ylabel('Dé2?Im');
title('y(n)=x(n)-ax(n-1)′?ê?oˉêyá??￠??μ?·?2?'); grid on
[H,w]=freqz(B,A,'whole'); subplot(3,1,2);
plot(w/pi,abs(H),'linewidth',2); grid on ;
xlabel('\omega/\pi'); ylabel('|H(e^j^\omega)|'); title('·ù?μ?ìó|ì?D?'); axis([0,2,0,6]);
subplot(3,1,1);
plot(w/pi,angle(H),'linewidth',2); grid on ;
axis([-0.1,2.1,-3,3]); xlabel('\omega/\pi'); ylabel('\phi(\omega)'); title('?à?μ?ìó|ì?D?');
图像：
-6
-4
-2
02
4
6
实部Re
虚部I m
y(n)=x(n)-ax(n-1)传输函数零、极点分布
00.20.4
0.60.81 1.2 1.4 1.6
1.82
5
ω/π
|H (e j ω)|
幅频响应特性
0.2
0.4
0.6
0.8
1 1.
2 1.4
1.6
1.8
2
-2
02
ω/π
φ(ω)
相频响应特性
a=0.8
代码：A=1;a=0.8
B=[1,a];
subplot(3,1,3);zplane(B,A); xlabel('êμ2?Re'); ylabel('Dé2?Im');
title('y(n)=x(n)-ax(n-1)′?ê?oˉêyá??￠??μ?·?2?'); grid on
[H,w]=freqz(B,A,'whole'); subplot(3,1,2);
plot(w/pi,abs(H),'linewidth',2); grid on ;
xlabel('\omega/\pi'); ylabel('|H(e^j^\omega)|'); title('·ù?μ?ìó|ì?D?'); axis([0,2,0,6]);
subplot(3,1,1);
plot(w/pi,angle(H),'linewidth',2); grid on ;
axis([-0.1,2.1,-3,3]); xlabel('\omega/\pi'); ylabel('\phi(\omega)'); title('?à?μ?ìó|ì?D?'); 图像：
-6
-4
-2
02
4
6
实部Re
虚部I m
y(n)=x(n)-ax(n-1)传输函数零、极点分布
00.20.4
0.60.81 1.2 1.4 1.6 1.8
2
5
ω/π
|H (e j ω)|
幅频响应特性
0.2
0.4
0.6
0.8
1 1.
2 1.4
1.6
1.8
2
-2
02
ω/π
φ(ω)
相频响应特性
a =0.9
代码：A=1;a=0.9
B=[1,a];
subplot(3,1,3);zplane(B,A); xlabel('êμ2?Re'); ylabel('Dé2?Im');
title('y(n)=x(n)-ax(n-1)′?ê?oˉêyá??￠??μ?·?2?'); grid on
[H,w]=freqz(B,A,'whole'); subplot(3,1,2);
plot(w/pi,abs(H),'linewidth',2); grid on ;
xlabel('\omega/\pi'); ylabel('|H(e^j^\omega)|'); title('·ù?μ?ìó|ì?D?'); axis([0,2,0,6]);
subplot(3,1,1);
plot(w/pi,angle(H),'linewidth',2); grid on ;
axis([-0.1,2.1,-3,3]); xlabel('\omega/\pi'); ylabel('\phi(\omega)'); title('?à?μ?ìó|ì?D?'); 图像：
-6
-4
-2
02
4
6
实部Re
虚部I m
y(n)=x(n)-ax(n-1)传输函数零、极点分布
00.20.4
0.60.81 1.2 1.4 1.6
1.82
5
ω/π
|H (e j ω)|
幅频响应特性
0.2
0.4
0.6
0.8
1 1.
2 1.4
1.6
1.8
2
-2
02
ω/π
φ(ω)
相频响应特性
分析：由y (n )=x (n )+ax (n -1)可知:H[z]=B[z]/A[z]= (1-az^(-1)) /1
系统极点z=0，零点z=a，当B点从w=0逆时针旋转时，在w=0点，由于零点向量长度最长，形成波峰：在w=pi点形成波谷；z=a处极点不影响相频响应。
3. 假设系统函数用下式描述：
y(n)=1.273y(n-1)-0.81y(n-2)+x(n)+x(n-1)
试分析它的频率特性，要求打印其幅度特性曲线，并求出峰值频率和谷值频率。
代码：A=[1,-1.273,0.81];
B=[1,1];
subplot(3,1,3);zplane(B,A);
xlabel('êμ2?Re');
ylabel('Dé2?Im');
title('y(n)=1.273y(n-1)-0.81y(n-2)+x(n)+x(n-1)′?ê?oˉêyá??￠??μ?·?2?'); grid on
[H,w]=freqz(B,A,'whole');
subplot(3,1,2);
plot(w/pi,abs(H),'linewidth',2);
grid on;
xlabel('\omega/\pi');
ylabel('|H(e^j^\omega)|');
title('·ù?μ?ìó|ì?D?');
axis([0,2,0,20]);
subplot(3,1,1);
plot(w/pi,angle(H),'linewidth',2);
grid on;
axis([-0.1,2.1,-3,3]);
xlabel('\omega/\pi');
ylabel('\phi(\omega)');
title('?à?μ?ìó|ì?D?');
图像：
-6
-4
-2
02
4
6
实部Re
虚部I m
y(n)=1.273y(n-1)-0.81y(n-2)+x(n)+x(n-1)传输函数零、极点分布0
0.2
0.4
0.6
0.8
1 1.
2 1.4 1.6
1.82
01020
ω/π
|H (e j ω)|
幅频响应特性
0.2
0.4
0.6
0.8
1 1.
2 1.4
1.6
1.8
2
-2
02
ω/π
φ(ω)
分析：由y (n )=1.273y (n -1)-0.81y (n -2)+x (n )+x (n -1)可知:
H[z]=B[z]/A[z]=(1+z^(-1))/(1-1.273 z^(-1)+0.81 z^(-2))
系统极点z1=0.79+j0.62*1.62^(-2)，z2=0.79-j0.62*1.62^(-2)零点z1=-1，z2=0当B 点从w=0逆时针旋转时，当旋转到接近极点z1=0.79+j0.62*1.62^(-2)是极点向量长度最短，幅度特性出现峰值。当转到w=pi 点形成波谷；z=a 处零点不影响幅频响应。当旋转到接近极点z2=0.79-j0.62*1.62^(-2)是极点向量长度再次最短，幅度特性再次出现峰值。 四、总结
1、将实验结果与理论结果进行对比。
从图形可得出零极点分布，根据相频响应幅度图和幅频响应幅度图可以得出零极点对幅
度特性的影响。实验结果与理论结果符合。
2、总结零、极点分布对频率响应的影响。
当旋转点转到极点附近时，幅度特性出现峰值，并且极点越靠近单位圆，峰值越尖锐。如果极点在单位圆上，系统不稳定。当旋转点转到零点附近时，幅度特性出现谷值，并且零点越靠近单位圆，谷值越接近零，零点在单位圆上时，谷值为零。所以：极点位置主要影响频响的峰值位置及尖锐程度，零点主要影响频响的谷值位置及形状
—完—
实验五：FIR数字滤波器设计与软件实现 信息学院 10电本2班王楚炘 2010304224 10.5.1 实验指导 1．实验目的 （1）掌握用窗函数法设计FIR数字滤波器的原理和方法。 （2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。 （3）掌握FIR滤波器的快速卷积实现原理。 （4）学会调用MATLAB函数设计与实现FIR滤波器。 2．实验内容及步骤 （1）认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理； （2）调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图10.5.1所示； 图10.5.1 具有加性噪声的信号x(t)及其频谱如图（3）请设计低通滤波器，从高频噪声中提取xt中的单频调幅信号，要求信号幅频失真小于0.1dB，将噪声频谱衰减60dB。先观察xt的频谱，确定滤波器指标参数。 （4）根据滤波器指标选择合适的窗函数，计算窗函数的长度N，
调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序，调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 （4）重复（3），滤波器指标不变，但改用等波纹最佳逼近法，调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。 提示：MATLAB函数fir1和fftfilt的功能及其调用格式请查阅本书 第7章和第？章； 采样频率Fs=1000Hz，采样周期T=1/Fs； 根据图10.6.1(b)和实验要求，可选择滤波器指标参数：通带截止频率fp=120Hz，阻带截至频率fs=150Hz，换算成数字频率，通带截止频率，通带最大衰为0.1dB，阻带截至频率，阻带最小衰为60dB。]实验程序框图如图10.5.2所示，供读者参考。 Fs=1000，T=1/Fs xt=xtg 产生信号xt, 并显示xt及其频谱 用窗函数法或等波纹最佳逼近法 设计FIR滤波器hn 对信号xt滤波：yt=fftfilt(hn,xt) 1、计算并绘图显示滤波器损耗函数 2、绘图显示滤波器输出信号yt End 图10.5.2 实验程序框图 4.思考题 (1)如果给定通带截止频率和阻带截止频率以及阻带最小衰减,如何用窗函数法设计线性相位低通滤波器?请写出设计步骤. 答：用窗函数法设计线性相位低通滤波器的设计步骤: a.根据对阻带衰减及过渡带的指标要求，选择窗函数的类型，并估计窗口的长度N； b.构造希望逼近的频率响应函数； c.计算h d(n)； d.加窗得到设计结果h(n)=h d(n)w(n)。 (2)如果要求用窗函数法设计带通滤波器,且给定通带上、下截止频率为和，阻带上、下截止频率为和，试求理想带通滤波器的截止频率。 答：希望逼近的理想带通滤波器的截止频率分别为：
实验报告 ( 2016-2017年度第二学期) 名称：《现代控制理论基础》 题目：状态空间模型分析 院系：控制科学与工程学院 班级： ___ 学号： __ 学生姓名： ______ 指导教师： _______ 成绩： 日期： 2017年 4月 15日
线控实验报告 一、实验目的： l.加强对现代控制理论相关知识的理解； 2.掌握用 matlab 进行系统李雅普诺夫稳定性分析、能控能观性分析； 二、实验内容 1 第一题：已知某系统的传递函数为G (s) S23S2 求解下列问题： （1）用 matlab 表示系统传递函数 num=[1]; den=[1 3 2]; sys=tf(num,den); sys1=zpk([],[-1 -2],1); 结果： sys = 1 ------------- s^2 + 3 s + 2 sys1 = 1 ----------- (s+1) (s+2) （2）求该系统状态空间表达式： [A1,B1,C1,D1]=tf2ss(num,den); A = -3-2 10 B = 1 C = 0 1
第二题：已知某系统的状态空间表达式为： 321 A ,B,C 01：10 求解下列问题： （1）求该系统的传递函数矩阵： （2）该系统的能观性和能空性： （3）求该系统的对角标准型： （4）求该系统能控标准型： （5）求该系统能观标准型： （6）求该系统的单位阶跃状态响应以及零输入响应：解题过程： 程序： A=[-3 -2;1 0];B=[1 0]';C=[0 1];D=0; [num,den]=ss2tf(A,B,C,D); co=ctrb(A,B); t1=rank(co); ob=obsv(A,C); t2=rank(ob); [At,Bt,Ct,Dt,T]=canon(A,B,C,D, 'modal' ); [Ac,Bc,Cc,Dc,Tc]=canon(A,B,C,D, 'companion' ); Ao=Ac'; Bo=Cc'; Co=Bc'; 结果： （1） num = 0 01 den = 1 32 （2）能控判别矩阵为： co = 1-3 0 1 能控判别矩阵的秩为： t1 = 2 故系统能控。 （3）能观判别矩阵为： ob = 0 1
实验一 离散时间信号分析 班级 信息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. 通过本次实验回忆并熟悉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;
数字信号处理实验实验一离散时间信号与系统及MA TLAB实现 1.单位冲激信号： n = -5:5; x = (n==0); subplot(122); stem(n, x); 2.单位阶跃信号： x=zeros(1,11); n0=0; n1=-5; n2=5; n = n1:n2; x(:,n+6) = ((n-n0)>=0); stem(n,x); 3.正弦序列： n = 0:1/3200:1/100; x=3*sin(200*pi*n+1.2); stem(n,x); 4.指数序列 n = 0:1/2:10; x1= 3*(0.7.^n); x2=3*exp((0.7+j*314)*n); subplot(221); stem(n,x1); subplot(222); stem(n,x2); 5.信号延迟 n=0:20; Y1=sin(100*n); Y2=sin(100*(n-3)); subplot(221); stem(n,Y1); subplot(222); stem(n,Y2);
6.信号相加 X1=[2 0.5 0.9 1 0 0 0 0]; X2=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7]; X=X1+X2; stem(X); 7.信号翻转 X1=[2 0.5 0.9 1]; n=1:4; X2=X1(5-n); subplot(221); stem(n,X1); subplot(222); stem(n,X2); 8.用MATLAB计算序列{-2 0 1 –1 3}和序列{1 2 0 -1}的离散卷积。a=[-2 0 1 -1 3]; b=[1 2 0 -1]; c=conv(a,b); M=length(c)-1; n=0:1:M; stem(n,c); xlabel('n'); ylabel('幅度'); 9．用MA TLAB计算差分方程 当输入序列为时的输出结果。 N=41; a=[0.8 -0.44 0.36 0.22]; b=[1 0.7 -0.45 -0.6]; x=[1 zeros(1,N-1)]; k=0:1:N-1; y=filter(a,b,x); stem(k,y) xlabel('n'); ylabel('幅度') 10.冲激响应impz N=64; a=[0.8 -0.44 0.36 0.22];
实验一 离散傅里叶变换（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）
实验一 IIR数字滤波器设计及软件实现 一、实验目的 （1）熟悉用双线性变换法设计IIR数字滤波器的原理与方法； （2）学会调用MATLAB信号处理工具箱中滤波器设计函数（或滤波器设计分析工具fdatool）设计各种IIR数字滤波器，学会根据滤波需求确定滤波器指标参数。 （3）掌握IIR数字滤波器的MATLAB实现方法。 （3）通过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的概念。 二、实验原理 设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法。基本设计过程是：①先将给定的数字滤波器的指标转换成过渡模拟滤波器的指标；②设计过渡模拟滤波器；③将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。MATLAB信号处理工具箱中的各种IIR数字滤波器设计函数都是采用双线性变换法。第六章介绍的滤波器设计函数butter、cheby1 、cheby2 和ellip可以分别被调用来直接设计巴特沃斯、切比雪夫1、切比雪夫2和椭圆模拟和数字滤波器。本实验要求同学调用如上函数直接设计IIR数字滤波器。 本实验的数字滤波器的MATLAB实现是指调用MATLAB信号处理工具箱函数filter对给定的输入信号x(n)进行滤波，得到滤波后的输出信号y(n）。 三、实验内容及步骤 （1）调用信号产生函数mstg，产生由三路抑制载波调幅信号相加构成的复合信号st，该函数还会自动绘图显示st的时域波形和幅频特性曲线，如图1-1所示。由图可见，三路信号时域混叠无法在时域分离。但频域是分离的，所以可以通过滤波的方法在频域分离，这就是本实验的目的。
图1-1 三路调幅信号st 的时域波形和幅频特性曲线 （2）要求将st 中三路调幅信号分离，通过观察st 的幅频特性曲线，分别确定可以分离st 中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率。要求滤波器的通带最大衰减为0.1dB,阻带最小衰减为60dB 。 提示：抑制载波单频调幅信号的数学表示式为 0001()cos(2)cos(2)[cos(2())cos(2())]2 c c c s t f t f t f f t f f t ππππ==-++ 其中，cos(2)c f t π称为载波，f c 为载波频率，0cos(2)f t π称为单频调制信号，f 0为调制正弦波信号频率，且满足0c f f >。由上式可见，所谓抑制载波单频调幅信号，就是2个正弦信号相乘，它有2个频率成分：和频0c f f +和差频0c f f -，这2个频率成分关于载波频率f c 对称。所以，1路抑制载波单频调幅信号的频谱图是关于载波频率f c 对称的2根谱线，其中没有载频成分，故取名为抑制载波单频调幅信号。容易看出，图1-1中三路调幅信号的载波频率分别为250Hz 、500Hz 、1000Hz 。如果调制信号m(t)具有带限连续频谱，无直流成分，则()()cos(2)c s t m t f t π=就是一般的抑制载波调幅信号。其频谱图是关于载波频率f c 对称的2个边带（上下边带），在专业课通信原理中称为双边带抑制载波 (DSB-SC) 调幅信号,简称双边带 (DSB) 信号。如果调制信号m(t)有直流成分，则()()cos(2)c s t m t f t π=就是一般的双边带调幅信号。其频谱图是关于载波频率f c 对称的2个边带（上下边带），并包含载频成分。
《现代控制理论》实验报告 状态空间极点配置控制实验 一、实验原理 经典控制理论的研究对象主要是单输入单输出的系统，控制器设计时一般需要有关被控对象的较精确模型，现代控制理论主要是依据现代数学工具，将经典控制理论的概念扩展到多输入多输出系统。极点配置法通过设计状态反馈控制器将多变量系统的闭环系统极点配置在期望的位置上，从而使系统满足瞬态和稳态性能指标。 １．状态空间分析 对于控制系统X = AX + Bu 选择控制信号为：u = ?KX 式中：X 为状态向量（ n 维）u 控制向量（纯量） A n × n维常数矩阵 B n ×1维常数矩阵 求解上式，得到 x(t) = (A ? BK)x(t) 方程的解为： x(t) = e( A?BK )t x(0) 状态反馈闭环控制原理图如下所示： 从图中可以看出，如果系统状态完全可控，K 选择适当，对于任意的初始状态，当t趋于无穷时，都可以使x(t)趋于0。 ２．极点配置的设计步骤 1) 检验系统的可控性条件。 2) 从矩阵 A 的特征多项式 来确定 a1, a2,……,an的值。 3) 确定使状态方程变为可控标准型的变换矩阵 T：T = MW 其中 M 为可控性矩阵， 4) 利用所期望的特征值，写出期望的多项式 5) 需要的状态反馈增益矩阵 K 由以下方程确定： 二、实验内容 针对直线型一级倒立摆系统应用极点配置法设计控制器，进行极点配置并用Matlab进行仿真实验。 三、实验步骤及结果 1.根据直线一级倒立摆的状态空间模型，以小车加速度作为输 入的系统状态方程为： 可以取1 l 。则得到系统的状态方程为： 于是有：
直线一级倒立摆的极点配置转化为： 对于如上所述的系统，设计控制器，要求系统具有较短的调整时间（约 3 秒）和合适的阻尼（阻尼比? = 0.5）。 2.采用四种不同的方法计算反馈矩阵 K。 方法一：按极点配置步骤进行计算。 1) 检验系统可控性，由系统可控性分析可以得到，系统的状态完全可控性矩阵的秩等于系统的状态维数(4)，系统的输出完全可控性矩阵的秩等于系统输出向量y 的维数(2)，所以系统可控。 倒立摆极点配置原理图 2) 计算特征值 根据要求，并留有一定的裕量（设调整时间为 2 秒），我们选取期望的闭环极点s =μi (i = 1,2,3,4) ，其中： 其中，μ 3,μ 4 使一对具有的主导闭环极点，μ 1 ,μ 2 位于 主导闭环极点的左边，因此其影响较小，因此期望的特征方程为： 因此可以得到： 由系统的特征方程： 因此有 系统的反馈增益矩阵为： 3) 确定使状态方程变为可控标准型的变换矩阵 T：T = MW 式中： M = 0 1.0000 0 0 1.0000 0 0 0 0 0.7500 0 5.5125 0.7500 0 5.5125 0 W = 0 -7.3500 -0.0000 1.0000 -7.3500 -0.0000 1.0000 0 -0.0000 1.0000 0 0 1.0000 0 0 0 于是可以得到： T = -7.3500 -0.0000 1.0000 0 0 -7.3500 -0.0000 1.0000 0 -0.0000 0.7500 0 -0.0000 0 -0.0000 0.7500 T’= -7.3500 0 0 -0.0000 -0.0000 -7.3500 -0.0000 0 1.0000 -0.0000 0.7500 -0.0000 0 1.0000 0 0.7500
数字信号处理实验答案 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),绘出四个周期。
实验三 离散付里叶变换（DFT ） 一、实验目的： 1． 掌握离散付里叶级数 2． 掌握DFT 变换。 3． 掌握DFT 特性。 4． 掌握利用DFT 计算线性卷积。 5． 掌握快速付里叶变换(FFT)。 二、实验原理： 1．离散付里叶级数（DFS ） )(~n x 为周期序列，其频率为基本频率（N /2π ）的倍数（或谐波）。其离散付叶级 数（DFS ）为：；∑ -=-±== 1 2,,1,0,)(~)(~N n kn N j k e n x k X π IDFS 为：∑ -== 1 2)(~ 1)(~N k kn N j e k X N n x π 2．离散付里叶变换（DFT ） )(n x 为长度N 的有限长序列，其DFT 为：21 1()()N j kn N k x n X k e N π--== ∑ IDFT 为：21 ()(),0,1,N j kn N n X k x n e k π-== =±∑ 3.DFT 的特性： （1） 线性性：)]([)]([)]()([2121n x bDFT n x aDFT n bx n ax DFT +=+ （2） 循环折叠（圆周对称）性：?? ? -≤≤-==-1 1) (0) 0())((N n n N x n x n x N （3） 共轭性: N k X n x DFT ))(()]([* *-= （4） 实序列的对称性（圆周共轭对称性）：N k X k X ))(()(* -= （5） 序列的圆周移位：N m n x m n x ))(()(~ -=-
（6） 频域中的圆周移位：)())(()]([|ln k R l k X n x W DFT N N N -=- （7） 时域循环卷积：)()()]()([2121k X k X n x n x DFT =? （8） 频域循环卷积（乘法性）：)()(1)]()([2121k X k X N n x n x DFT ?= （9） 帕塞瓦尔（Parseval ）定理：∑ ∑ -=-== = 1 2 1 2 ) (1)(N k N n x k X N n x E 4.用DFT 计算线性卷积: 设)(1n x 为1N 点序列，)(2n x 为2N 点序列，)(3n x 为)(1n x 和)(2n x 的线性卷积，其为121-+N N 点序列，)(4n x 为)(1n x 和)(2n x 的圆卷积，其长度为 N ，当121-+=N N N 时，)()(43n x n x =。实际中，采用分段卷积法，即重 叠保留法和重叠相加法。需要对数据流进行分块处理，这时直接采用DFT 计算线性卷积会产生一些问题，而应该将)(n x 通过重复前M-1个取样进行分块，这样可得到正确结果。 5.快速付里叶变换（FFT ）： 掌握基2-时域抽取FFT(DIT-FFT)和基2-频域抽取FFT(DIF-FFT)。MATLAB 提供fft 函数来计算x 的DFT 。fft 函数是用机器语言写的，采用混合基法，其调用形式为：),(N x fft X =。如N 为2的幂，则得到高速的基2-FFT 算法；若N 不是2的乘方，则将N 分解成质数，得到较慢的混合基FFT 算法；最后，若N 为质数，则fft 函数采用的是原始的DFT 算法。 三、实验步骤： 1.离散付里叶级数（DFS ） （1）自已动手:编写实现离散付里叶级数和逆离散付里叶级数的函数。 （2）已知周期性序列如下所示：}3,2,1,0,3,2,1,0,3,2,1,0{)(~ ↑ =n x 求其离散付里叶级数。 2.离散付里叶变换（DFT ） (1) 编写实现DFT 和IDFT 的函数。 (2) 已知)(n x 是一个六点序列，如下所示： ?? ?≤≤=e ls e n n x 0 501 )( 要求计算该序列的离散时间的付里叶变换和离散付里叶变换，并绘出它们的幅度和相
数字信号处理软件实验 MATLAB 仿真 2015年12月16日
实验一：数字信号的 FFT 分析 ● 实验目的 通过本次实验，应该掌握： (a) 用傅立叶变换进行信号分析时基本参数的选择。 (b) 经过离散时间傅立叶变换（DTFT ）和有限长度离散傅立叶变换（DFT ）后信号频谱上的区别，前者 DTFT 时间域是离散信号，频率域还是连续的，而 DFT 在两个域中都是离散的。 (c) 离散傅立叶变换的基本原理、特性，以及经典的快速算法（基2时间抽选法），体会快速算法的效率。 (d) 获得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨率的概念，为将来进一步进行时频分析（例如小波）的学习和研究打下基础。 (e) 建立 DFT 从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念的一个典型应用是数字音频压缩中的分析滤波器，例如 DVD AC3 和MPEG Audio 。 ● 实验内容及要求 ? 离散信号的频谱分析 设信号 此信号的0.3pi 和 0.302pi 两根谱线相距很近，谱线 0.45pi 的幅度很小，请选择合适的序列长度 N 和窗函数，用 DFT 分析其频谱，要求得到清楚的三根谱线。 ? DTMF 信号频谱分析 用计算机声卡采用一段通信系统中电话双音多频（DTMF ）拨号数字 0～9的数据，采用快速傅立叶变换（FFT ）分析这10个号码DTMF 拨号时的频谱。 00010450303024().*cos(.)sin(.)cos(.)x n n n n ππππ=+--
●MATLAB代码及结果 ?离散信号的频谱分析 clf; close all; N=1000; n=1:1:N; x=0.001*cos(0.45*n*pi)+sin(0.3*n*pi)-cos(0.302*n*pi-pi/4); y=fft(x,N); mag=abs(y); w=2*pi/N*[0:1:N-1]; stem(w/pi,mag); axis([0.25 0.5 0 2]); xlabel('频率'); ylabel('X(k)'); grid on;
现代控制理论课 程设计报告 题目打印机皮带驱动系统能控能观和稳定性分析 项目成员史旭东童振梁沈晓楠 专业班级自动化112 指导教师何小其 分院信息分院 完成日期 2014-5-28
目录 1. 课程设计目的 (3) 2.课程设计题目描述和要求 (3) 3.课程设计报告内容 (4) 3.1 原理图 (4) 3.2 系统参数取值情况 (4) 3.3 打印机皮带驱动系统的状态空间方程 (5) 4. 系统分析 (7) 4.1 能控性分析 (7) 4.2 能观性分析 (8) 4.3 稳定性分析 (8) 5. 总结 (10)
项目组成员具体分工 打印机皮带驱动系统能控能观和稳定性 分析 课程设计的内容如下： 1.课程设计目的 综合运用自控现代理论分析皮带驱动系统的能控性、能观性以及稳定性，融会贯通并扩展有关方面的知识。加强大家对专业理论知识的理解和实际运用。培养学生熟练运用有关的仿真软件及分析，解决实际问题的能力，学会使用标准、手册、查阅有关技术资料。加强了大家的自学能力，为大家以后做毕业设计做很好的铺垫。 2.课程设计题目描述和要求 (1)环节项目名称：能控能观判据及稳定性判据 (2)环节目的： ①利用MATLAB分析线性定常系统的可控性和客观性。 ②利用MATLAB进行线性定常系统的李雅普诺夫稳定性判据。 (3)环节形式：课后上机仿真 (4)环节考核方式： 根据提交的仿真结果及分析报告确定成绩。 (5)环节内容、方法： ①给定系统状态空间方程，对系统进行可控性、可观性分析。 ②已知系统状态空间方程，判断其稳定性，并绘制出时间响应曲线验
证上述判断。 3.课程设计报告内容 3.1 原理图 在计算机外围设备中，常用的低价位喷墨式或针式打印机都配有皮带驱动器。它用于驱动打印头沿打印页面横向移动。图1给出了一个装有直流电机的皮带驱动式打印机的例子。其光传感器用来测定打印头的位置，皮带张力的变化用于调节皮带的实际弹性状态。 图1 打印机皮带驱动系统 3.2 系统参数取值情况 表1打印装置的参数
实验三：离散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
DSP运行实验报告 一、实验目的 熟悉CCS软件仿真下，DSP程序的下载和运行；熟悉借助单片机的DSP程序下载和运行； 熟悉借助仿真器的DSP程序下载和运行；熟悉与DSP程序下载运行相关的CCS编程环境。 二、实验原理 CCS软件仿真下，借用计算机的资源仿真DSP的内部结构，可以模拟DSP程序的下载和运行。 如果要让程序在实验板的DSP中运行、调试和仿真，可以用仿真器进行DSP程序下载和运行。初学者也可以不用仿真器来使用这款实验板，只是不能进行程序调试和仿真。 在本实验板的作用中，单片机既是串口下载程序的载体，又是充当DSP 的片外存储器（相对于FLASH），用于固化程序。 三、实验设备、仪器及材料 安装有WINDOWS XP操作系统和CCS3.3的计算机。 四、实验步骤（按照实际操作过程） 1、CCS软件仿真下，DSP程序的下载和运行。 第一步：安装CCS，如果不使用仿真器，CCS 的运行环境要设置成一个模拟仿真器（软仿真）。
第二步：运行CCS，进入CCS 开发环境。 第三步：打开一个工程。 将实验目录下的EXP01目录拷到D:\shiyan下（目录路径不能有中文），用[Project]\[Open]菜单打开工程，在“Project Open”对话框中选 EXP01\CPUtimer\CpuTimer.pjt，选“打开”， 第四步：编译工程。 在[Project]菜单中选“Rebuild All”，生成CpuTimer.out文件。 第五步：装载程序。 用[File]\[Load Program]菜单装载第四步生成CpuTimer.out文件，在当前工程目录中的Debug 文件夹中找到CpuTimer.out文件，选中，鼠标左键单击“打开”。
计算机控制理论与设计作业 题目：基于极点配置方法的直流调速系统的控制器设计
摘要 本文目的是用极点配置方法对连续的被控对象设计控制器。基本思路是对连续系统进行数学建模，将连续模型进行离散化，针对离散的被控对象，用极点配置的方法分别在用状态方程和传递函数两种描述方法下设计前馈和反馈控制器，并用MATLAB仿真。文中具体以直流调速系统作为研究对象，对直流调速系统的组成和结构进行了分析，把各个部分进行数学建模，求出其传递函数，组成系统结构框图，利用自控原理的知识对结构图化简，求出被控对象的传递函数和状态方程，进一步得将其离散化。第一种是通过极点配置设计方法的原理，用状态方程设计被控对象的控制律，因为直流调速系统存在噪声，实际状态不可测，故选择了全阶的观测器，又因为采样时间小于计算延时，所以选择了预报观测器。利用所学知识对此闭环系统设计前馈和反馈控制器[1]。第二种利用传统的离散传递函数，从代数多项式的角度进行复合控制器的设计，在保证系统稳定的情况下，分析系统的可实现性，稳定性，静态指标，动态指标，抗干扰等方面性能研究前馈反馈相结合控制器设计。重点是保证被控对象的不稳定的零极点不能被抵消。最后利用MATLAB的Simulink进行仿真，观察系统的输出的y和u和收敛性，并加入扰动看其抗干扰性能，得出结论。 经研究分析，对于直流调速系统，基于极点配置设计的前馈反馈相结合的控制器，具有良好的稳定性能和抗干扰性能。运行结果符合实际情况。 关键词：极点配置；状态方程；直流调速系统；代数多项式；Matlab；
1绪论 1.1论文的背景及意义 在工业生产和日常生活中，自动控制系统分为确定性系统和不确定性系统两类，确定性系统是指系统的结构和参数是确定的，确定的输入下，输出也确定的一类系统。确定性系统相对于不确定性系统而言的。在确定的系统中所用的变量都可用确切的函数关系来描述，系统的运动特性可以完全确定。以确定性系统为研究对象的控制理论称为确定性控制理论。本文以直流调速系统为研究对象，利用极点配置的设计方法，包括利用状态空间模型和传递函数模型分别描述线性系统，采用闭环极点为指标的控制器设计的理论和方法，设计出前馈和反馈控制器，组建闭环控制系统，用Matlab进行仿真可以逼真地还原出实际系统。 1.2 论文的主要内容 本文直流电机的调速系统的模型作为研究对象，利用线性系统极点配置的设计方法，设计前馈反馈控制器。论文研究的主要内容: (1)阅读学习国内外期刊文献，研究了极点配置的基本原理和Matlab的实现方法。 (2)系统的说明直流电机的系统结构和工作原理并分析，建立直流调速系统的数学模型，将其进行离散化，并讨论其传递函数与状态方程之间的关系。 (3)分析极点配置控制器的设计原理，利用状态方程设计控制器。 （4）将被控对象的传递函数离散化，利用传递函数模型设计控制器。 (4)在MATLAB中建立闭环直流调速系统的模型，根据闭环极点配置的设计步骤编写程序，用Simulink搭建仿真系统，对闭环直流调速系统的输出进行仿真分析。 (5)对仿真结果分析。将仿真结果与实际直流调速系统的阶跃响应的各项参数相比较，得出结论。
数字信号处理实验四 第一题结果： （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))'); 单位脉冲响应曲线 幅频和相频特性曲线 第二题结果：
子程序： function myplot(B,A) %myplot(B,A) %时域离散系统损耗函数绘图 %B为系统函数分子多项式系数向量 %A为系统函数分母多项式系数向量 [H,W]=freqz(B,A,1000); m=abs(H); plot(W/pi,20*log10(m/max(m)));grid on; xlabel('\omega/\pi');ylabel('幅度(dB)') axis([0,1,-80,5]);title('损耗函数曲线'); function tplot(xn,T,yn) %时域序列连续曲线绘图函数 % xn:信号数据序列，yn:绘图信号的纵坐标名称（字符串） % T为采样间隔 n=0:length(xn)-1;t=n*T; plot(t,xn); xlabel('t/s');ylabel(yn); axis([0,t(end),min(xn),1.2*max(xn)]) 程序： %实验4程序exp4.m % IIR数字滤波器设计及软件实现 clear all;close all Fs=10000;T=1/Fs; %采样频率 %调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st st=mstg; %低通滤波器设计与实现========================================= fp=280;fs=450; wp=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 y1t=filter(B,A,st); %滤波器软件实现 % 低通滤波器设计与实现绘图部分 figure(2);subplot(3,1,1); myplot(B,A); %调用绘图函数myplot绘制损耗函数曲线 yt='y_1(t)'; subplot(3,1,2);tplot(y1t,T,yt); %调用绘图函数tplot绘制滤波器输出波形 %带通滤波器设计与实现==================================================== 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;