实验3 离散系统的变换域分析
一、实验目的：
加深对离散系统的频率响应分析和零、极点分布的概念理解。
二、实验原理：
离散系统的时域方程为
∑∑==-=-M
m m N k n m n x b k n y a
00)()( 其变换域分析方法如下：
X(z)H(z)Y(z) )()()()()(=?-=
*=∑∞-∞=m m n h m x n h n x n y
系统函数为 N
N M
M z a z a a z b z b b z X z Y z H ----++++++==......)()()(110110 分解因式 ∏∏∑∑=-=-=-=---==k k
M m m k k k M m m m z d z c K z a z b z H 111100)1()
1()( ， 其中 m c 和 k d 称为零、极点。
在MATLAB 中，可以用函数[z ，p ，K]=tf2zp （num ，den ）求得有理分式形式的系统函数的零、极点，用函数zplane （z ，p ）绘出零、极点分布图；也可以用函数zplane （num ，den ）直接绘出有理分式形式的系统函数的零、极点分布图。使h=freqz(num,den,w)函数可求系统的频率响应，w 是频率的计算点，如w=0:pi/255:pi, h 是复数，abs(h)为幅度响应，angle(h)为相位响应。另外，在MATLAB 中，可以用函数 [r ，p ，k]=residuez （num ，den ）完成部分分式展开计算；可以用函数sos=zp2sos （z ，p ，K ）完成将高阶系统分解为2阶系统的串联。
（在实验报告中对这几种函数的使用方法及参数含义做出说明，这一部分手写）
三、实验内容
例1 求下列直接型系统函数的零、极点，并将它转换成二阶节形式
解 用MATLAB 计算程序如下：
num=[1 -0.1 -0.3 -0.3 -0.2];
den=[1 0.1 0.2 0.2 0.5];
[z,p,k]=tf2zp(num,den);
disp('零点');disp(z);
disp('极点');disp(p);
disp('增益系数');disp(k); sos=zp2sos(z,p,k);
disp('二阶节');disp(real(sos));
zplane(num,den)
输入到“num”和“den”的分别为分子和分母多项式的系数。计算求得零、极点增益系数和二阶节的系数：
零点
0.9615
-0.5730
-0.1443 + 0.5850i
-0.1443 - 0.5850i
极点
0.5276 + 0.6997i
0.5276 - 0.6997i
-0.5776 + 0.5635i
-0.5776 - 0.5635i
增益系数
1
二阶节
1.0000 -0.3885 -0.5509 1.0000 1.1552 0.6511
1.0000 0.2885 0.3630 1.0000 -1.0552 0.7679
系统函数的二阶节形式为：
极点图如右图。
例2 差分方程
)
3(02.0)2(36.0)1(44.0)(8.0 )
3(6.0)2(45.0)1(7.0)(-+-+--=-----+n x n x n x n x n y n y n y n y 所对应的系统的频率响应。
解：差分方程所对应的系统函数为
3213
216.045.07.0102.036.044.08.0)(--------+++-=z
z z z z z z H 用MATLAB 计算的程序如下：
k=256;
num=[0.8 -0.44 0.36 0.02];
den=[1 0.7 -0.45 -0.6];
w=0:pi/k:pi;
h=freqz(num,den,w);
subplot(2,2,1);
plot(w/pi,real(h));grid
title('实部')
xlabel('\omega/\pi');ylabel('幅度')
subplot(2,2,2);
plot(w/pi,imag(h));grid
title('虚部')
xlabel('\omega/\pi');ylabel('Amplitude')
subplot(2,2,3);
plot(w/pi,abs(h));grid
title('幅度谱')
xlabel('\omega/\pi');ylabel('幅值')
subplot(2,2,4);
plot(w/pi,angle(h));grid
title('相位谱') xlabel('\omega/\pi');ylabel('弧度')
练习1.求系统
54321543212336.09537.08801.14947.28107.110528.0797.01295.01295.00797.00528.0)(-----------+-+-+++++=z z z z z z z z z z z H 的零、极点和幅度频率响应和相位响应。
要求：绘出零、极点分布图，幅度频率响应和相位响应曲线。 代码如下：
clear;
num=[0.0528 0.0797 0.1295 0.1295 0.797 0.0528]; den=[1 -1.8107 2.4947 -1.8801 0.9537 0.2336];
[z,p,k]=tf2zp(num,den);
disp('零点');disp(z);
disp('极点');disp(p);
disp('增益系数');disp(k);
subplot(2,2,1);
zplane(num,den);
title('零极点分布图');
k=256;
w=0:pi/k:pi;
h=freqz(num,den,w);
subplot(2,2,2);
plot(w/pi,abs(h));grid
title('幅度谱')
xlabel('\omega/\pi');ylabel('幅值')
subplot(2,2,3);
plot(w/pi,angle(h));grid
title('相位谱')
xlabel('\omega/\pi');ylabel('弧度')
输出结果：
练习2.求一因果线性移不变系统()0.81(2)()(2)
y n y n x n x n =-+--的单位抽样响应()h n ，单位阶跃响应()g n ，并绘出()jw H e
的幅频和相频特性。
代码：
clear;
w=[0:1:500]*pi/500;
a=[1,0,-0.81];
b=[1,0,-1];
x=[1 zeros(1,25)];
h=filter(b,a,x);
n=0:1:25;
subplot(2,2,1);
stem(n,h,'filled');
title('单位抽样相应');
step=ones(1,25);
subplot(2,2,2);
g=conv(h,step);
N1=0:49;
stem(n1,g,'filled');
title('单位阶跃响应');
[r,p,k]=residuez(b,a);
h=freqz(b,a,w);
magH=abs(h);
phah=angle(h);
subplot(2,2,3);
plot(w/pi,magH,'b'); grid
title('幅度特性曲线');
subplot(2,2,4);
plot(w/pi,phah/pi,'b'); grid
title('相位特性曲线')
输出结果：
练习3.利用MATLAB 计算12318()1834F z z z z
---=
+-- 的部分分式展开式。 代码：
clear;
a=[18 3 -4 -1];
b=[18];
[r,p,k]=residuez(b,a);
实验心得：
西华大学实验报告（理工类） 开课学院及实验室： 实验时间 ： 年 月 日 一、实验目的 1)初步掌握电力系统物理模拟实验的基本方法。 2)加深理解功率极限的概念，在实验中体会各种提高功率极限措施的作用。 3)通过对实验中各种现象的观察，结合所学的理论知识，培养理论结合实际及分析问题的能力。 二、实验原理 所谓简单电力系统，一般是指发电机通过变压器、输电线路与无限大容量母线联接而且不计各元件的电阻和导纳的输电系统。 对于简单系统，如发电机至系统d 轴和g 轴总电抗分别为d X ∑和q X ∑，则发电机的功率特性为 当发电机装有励磁调节器时，发电机电势q E 随运行情况而变化，根据一般励磁调节器的性能，可认为保持发电机'q E （或' E ）恒定。这时发电机的功率特性可表示成 或 这时功率极限为 随着电力系统的发展和扩大，电力系统的稳定性问题更加突出，而提高电力系统稳定性和输送能力的最重要手段之一，就是尽可能提高电力系统的功率极限。从简单电力系统功率极限的表达式看，要提高功率极限，可以通过发电机装设性能良好的励磁调节器，以提高发电机电势、增加并联运行线路回路数；或通过串联电容补偿等手段，以减少系统电抗，使受端系统维持较高的运行电压水平；或输电线采用中继同步调相机、中继电力系统等手段以稳定系统中继点电压。 (3)实验内容 1)无调节励磁时，功率特性和功率极隈的测定 ①网络结构变化对系统静态稳定的影响（改变戈）： 在相同的运行条件下（即系统电压U-、发电机电势E 。保持不变．罚芳赆裁Ll=E 。），分别 测定输电线单回线和双回线运行时，发电机的功一角特性曲线，＆豆甍辜授冁蝮和达到功率极 限时的功角值。同时观察并记录系统中其他运行参数（如发电极端毫玉萼蔫交化。将两种 情况下的结果加以比较和分析。 实验步骤如下： a)输电线路为单回线； b)发电机与系统并列后，调节发电机，使其输出的有功和无ZZ 蔓专零： c)功率角指示器调零； d)逐步增加发电机输出的有功功率，而发电机不调节震磁： e)观察并记录系统中运行参数的变化，填入表1.3中： f)输电线路为双回线，重复上述步骤，将运行参数填入表l 。毒=：
实验一 离散时间信号分析 班级 信息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 离散系统的时域及变换域分析 一、实验目的： 1.加深对离散系统的差分方程、单位抽样响应和卷积分析方法的理解。 2.加深对离散系统的频率响应分析和零、极点分布的概念理解。 二、实验原理： 1.时域 离散系统 其输入、输出关系可用以下差分方程描述： ∑∑==-=-M m m N k n m n x b k n y a )()( 输入信号分解为冲激信号， ∑∞ -∞ =-= m m n m x n x )()()(δ 系统单位抽样序列h （n ）， 则系统响应为如下的卷积计算式： ∑∞ -∞ =-= *=m m n h m x n h n x n y )()()()()( 当0 0≠a N k a k ,...2,1,0==时，h(n)是有限长度的（n ：[0，M]），称系统为FIR 系统；反之，称系统为IIR 系统。 在MATLAB 中，可以用函数y=filter(b,a,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积。 2.变换域 离散系统的时域方程为 ∑∑==-=-M m m N k n m n x b k n y a )()(
其变换域分析方法如下： X(z)H(z) Y(z) )()()()()(=?-= *=∑∞ -∞ =m m n h m x n h n x n y 系统函数为 N N M M z a z a a z b z b b z X z Y z H ----++++++= =......)()()(110110 分解因式 ∏∏∑∑=-=-=-=---== N k k M m m N k k k M m m m z d z c K z a z b z H 1 1 11 ) 1() 1()( ， 其中 m c 和 k d 称为零、极点。 在MATLAB 中，可以用函数[z ，p ，K]=tf2zp （num ，den ）求得有理分式形式的系统函数的零、极点，用函数zplane （z ，p ）绘出零、极点分布图；也可以用函数zplane （num ，den ）直接绘出有理分式形式的系统函数的零、极点分布图。使用h=freqz(num,den,w)函数可求系统的频率响应，w 是频率的计算点，如w=0:pi/255:pi, h 是复数，abs(h)为幅度响应，angle(h)为相位响应。另外，在MATLAB 中，可以用函数 [r ，p ，k]=residuez （num ，den ）完成部分分式展开计算；可以用函数sos=zp2sos （z ，p ，K ）完成将高阶系统分解为2阶系统的串联。 三 、实验内容 1.时域 （1.）编制程序求解下列系统的单位抽样响应，并绘出其图形。 )1()()2(125.0)1(75.0)(--=-+-+n x n x n y n y n y 解 用MATLAB 计算程序如下： N=15; n=0:N-1; b=[1,-1]; a=[1,0.75,0.125]; x=[n==0]; y=filter(b,a,x); subplot(3,2,1); stem(n,y,'.'); axis([0,N,-1,2]); ylabel('y(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;
管理信息系统实验报告 实验3 系统分析 课程名称：管理信息系统 指导教师：王玮 班级：信管1401 学号： 姓名：唐赛赛 时间： 2016.04.06 地点： 3 号机房
一、实验目的 1．了解开发Visio解决方案的基本概念和关于Visio工具的一些基本的操作和应用； 2．掌握系统分析阶段数据流程图的画法； 二、实验步骤和实验结果： 使用Visio中提供的“组织结构图”模具，绘制下面例题的组织结构图，附在图后。 2、使用Visio绘制“业务流程图模具”和“数据流程图模具”（1）创建“业务流程图模具” 先在“框图”－〉“基本形状”中找到圆角矩形，右击选择“添加到我的形状”－〉“添加到新模具”。之后出现“另存为”对话框，把新模具命名为“业务流程图”，把圆角矩形形添加到了新模具“业务程图”中。用同样的思路，先在“框图”－〉“基本形状”中找到圆形，右击选择“添加到我的形状”－〉“添加到模具“业务程图”中；在“框图”－〉“基本形状”找到矩形，在“流程图”中的“IDEFO图表形状”找到动态连接线，在“流程图”中的“SDL图表形状”中找到文档，多文档，添加到模具“业务程图”中。可以通过设置“动态连接线”属性来改变其形状。如下图：
添加完成后，我们就可以在画业务流程图时打开该模具，业务流程图所有的元素都会在一个模具中显示出来。（2）创建“数据流程图模具”先在“框图”－〉“基本形状”中找到圆形（或是“流程图”中的“混合流程图形状”中找到外部实体2 ），右击选择“添加到我的形状”－〉“添加到新模具”（注，使用外部实体2来表示外部实体的时候，请将之旋转180度使用）。之后出现“另存为”对话框，把新模具命名为“数据流程图”，这样我们就把圆形形添加
实验一 离散傅里叶变换（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）
实验6 离散时间系统的z 域分析 一、实验目的 1.掌握z 变换及其反变换的定义，并掌握MATLAB 实现方法。 2.学习和掌握离散时间系统系统函数的定义及z 域分析方法。 3.掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理 1. Z 变换 序列x(n)的z 变换定义为 ()()n n X z x n z +∞ -=-∞ = ∑ Z 反变换定义为 1 1 ()()2n r x n X z z dz j π-= ? 在MATLAB 中，可以采用符号数学工具箱的ztrans 函数和iztrans 函数计算z 变换和z 反变换： Z=ztrans(F) 求符号表达式F 的z 变换。 F=ilaplace(Z) 求符号表达式Z 的z 反变换。 2.离散时间系统的系统函数 离散时间系统的系统函数H(z)定义为单位抽样响应h(n)的z 变换 ()()n n H z h n z +∞ -=-∞ = ∑ 此外，连续时间系统的系统函数还可以由系统输入和输出信号的z 变换之比得到 ()()/()H z Y z X z =
由上式描述的离散时间系统的系统函数可以表示为 101101()M M N N b b z b z H z a a z a z ----+++= +++…… 3.离散时间系统的零极点分析 离散时间系统的零点和极点分别指使系统函数分子多项式和分母多项式为零的点。在MATLAB 中可以通过函数roots 来求系统函数分子多项式和分母多项式的根，从而得到系统的零极点。 此外，还可以利用MATLAB 的zplane 函数来求解和绘制离散系统的零极点分布图，zplane 函数调用格式为： zplane(b,a) b,a 为系统函数的分子、分母多项式的系数向量（行向量）。 zplane(z,p) z,p 为零极点序列（列向量）。 系统函数是描述系统的重要物理量，研究系统函数的零极点分布不仅可以了解系统单位抽样响应的变化，还可以了解系统的频率特性响应以及判断系统的稳定性： ①系统函数的极点位置决定了系统单位抽样响应h(n)的波形，系统函数零点位置只影响冲激响应的幅度和相位，不影响波形。 ②系统的频率响应取决于系统的零极点，根据系统的零极点分布情况，可以通过向量分析系统的频率响应。 ③因果的离散时间系统稳定的充要条件是H(z)的全部极点都位于单位圆内。 三、实验内容 (1)已知因果离散时间系统的系统函数分别为： ①23221()0.50.0050.3 z z H z z z z ++=--+
天津职业技术师范大学课程设计大学学籍管理系统的设计与开发 专业：软件工程 班级学号：软件1002-17 学生姓名：靳利强 指导教师：龚良波老师 二〇一三年七月
一．需求分析 1.课程名称：大学教务信息系统的设计与开发 2.设计目的： 为方便学校做好学生学籍管理工作，设计一个学生学籍管理系统，在设计过程中作了系统分析和总体设计，软件设计采取模块化的设计思路。 3.需求概述 该学生学籍管理系统主要对学生学籍信息、成绩信息进行管理，提供一个平台，供学籍管理人员增删改查学生信息、学生成绩信息。系统分为学生信息管理、学生成绩管理、信息查询等几个模块。学籍管理人员登录成功后可以对学生信息管理、学生成绩管理、信息查询等模块进行操作，如学生信息添加、修改、删除和查询；学生成绩登记、修改、删除和查询；查询信息等。 4功能需求： 1)功能齐全：界面操作灵活方便，设计包括以下基本功能： 2)学生信息管理、教师信息管理、财务信息管理、班级信息管理、课 程信息管理、成绩信息管理、打印信息管理、教室信息管理、综合信息查询、系统管理等，至少实现其中的三个功能，且每个功能至少包括两个子功能。 3)按照软件工程的要求进行分析、设计和开发。 4)界面友好：界面友好、输入有提示、尽量展示人性化。 5)可读性强：源程序代码清晰、有层次、主要程序段有注释。
6)健壮性好：用户输入非法数据时，系统应及时给出警告信息。 二.概要设计 1.功能模块： 2数据流图： （1）学生端
（2）管理员端
学生端功能： A 登录,学生登录后，验证成功，进入其信息展示页。 管理员端功能： B 登录，管理员登录后，验证成功，进入学生信息列表，可以对学生信息进行修改，删除，按班级查询，按学号查询，按名字查询。上传图片，更新图片等操作。 三．详细设计及实现 数据库设计： 学生表： 教师表：
数字信号处理实验答案 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 )( 要求计算该序列的离散时间的付里叶变换和离散付里叶变换，并绘出它们的幅度和相
本科实验报告 课程名称：系统分析与设计 实验项目：《》实验实验地点： 专业班级：学号： 学生姓名： 指导教师： 2016年11月日
一、实验目的 通过《系统分析与设计》实验，使学生在实际的案例中完成系统分析与系统设计中的主要步骤，并熟悉信息系统开发的有关应用软件，加深对信息系统分析与设计课程基础理论、基本知识的理解，提高分析和解决实际问题的能力，使学生在实践中熟悉信息系统分析与设计的规范，为后继的学习打下良好的基础。 二、实验要求 学生以个人为单位完成，自选题目，班内题目不重复，使用UML进行系统分析与设计，并完成实验报告。实验报告（A4纸+电子版)在最后一次上课时提交（10周）。 三、实验主要设备：台式或笔记本计算机 四、实验内容 1 选题及项目背景 学生填写自选题目 2 定义 学生填写（对自选项目系统进行描述200-400字） 3 参考资料 学生填写 4 系统分析与设计 4.1需求分析 4.1.1识别参与者 学生填写 4.1.2 对需求进行捕获与描述 学生填写时删除以下括号内容 （内容要求1：对每个用例进行概要说明，参考以下格式： 用例名称：删除借阅者信息执行者：管理员 目的：完成一次删除借阅者信息的完整过程。） （内容要求2：选择其中一个用例（如下订单）给出其用例描述。格式参考下表
） 4.1.3 用例图 通过已掌握的需求，初步了解系统所要完成的功能。下面给出用例图。 4.1.4 分析与讨论 1)建模用例图的步骤、方法？ 2)如何识别系统的参与者？应该如何划分用例，应注意哪些问题？ 3)心得 4.2 建立对象模型 4.2.1 候选类的数据字典 学生填写 4.2.2定义类 （内容以“书籍信息”类为例列出该类的属性和操作如下： “书籍信息”类 ?属性 国际标准书号(ISBN)：文本(String) 书名(name)：文本
实验三：离散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
实验2 离散系统的变换域分析 一、实验目的 1、熟悉对离散系统的频率响应分析方法； 2、加深对零、极点分布的概念理解。 二、实验原理 离散系统的时域方程为 其变换域分析方法如下： 频域: 系统的频率响应为: Z域: 系统的转移函数为: 分解因式: ， 其中 和 称为零、极点。 三、预习要求 1. 在MATLAB中，熟悉函数tf2zp、zplane、freqz、residuez、 zp2sos的使用，其中：[z，p，K]=tf2zp（num，den）求得有 理分式形式的系统转移函数的零、极点；zplane（z，p）绘 制零、极点分布图；h=freqz(num,den,w)求系统的单位频率 响应；[r，p，k]=residuez（num，den）完成部分分式展开 计算；sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系 统的串联。 2. 阅读文中的范例，学习频率分析法在MATLAB中的实现； 3. 编程实现系统参数输入，绘出幅度频率响应和相位响应曲线
和零、极点分布图。 四、实验内容 求系统 的零、极点和幅度频率响应和相位响应。 五、范例 求下列直接型系统函数的零、极点，并将它转换成二阶节形式 解:用MATLAB计算程序如下： num=[1 -0.1 -0.3 -0.3 -0.2]; den=[1 0.1 0.2 0.2 0.5]; [z,p,k]=tf2zp(num,den); m=abs(p); disp('零点');disp(z); disp('极点');disp(p); disp('增益系数');disp(k); sos=zp2sos(z,p,k); disp('二阶节');disp(real(sos)); zplane(num,den) 输入到“num”和“den”的分别为分子和分母多项式的系数。计算求得零、极点增益系数和二阶节的系数： 零点: 0.9615 -0.5730 -0.1443 + 0.5850i -0.1443 - 0.5850i 极点: 0.5276+0.6997i 0.5276-0.6997i -0.5776+0.5635i -0.5776-0.5635i 增益系数: 1 二阶节: 1.0000 -0.3885 -0.5509 1.0000 1.15520 0.6511
《管理信息系统》 实验二 题目：系统分析 专业：信息管理与信息系统 班级：1106班 姓名 ************************* 指导教师：贺玉珍老师 完成日期：2014.4.28
运城学院超市管理系统设计分析说明书 一、系统目标：随着小超市规模的发展不断扩大，商品数量急剧增加，有关商品的各种信息量也成倍增长。超市时时刻刻都需要对商品各种信息进行统计分析。而大型的超市管理系统功能过于强大而造成操作繁琐降低了小超市的工作效率。 超市管理系统是市场上最流行的超市上常用的系统之一，它主要包含以下几个模块：系统权限的设定、原始数据录入、数据的汇总及查询等。从而，实现对进货、销售及员工信息等实现全面、动态、及时的管理。 本文系统的分析了软件开发的背景以过程；首先介绍了软件的开发环境，其次介绍了本软件的详细设计过程：数据库的设计、各个模块的设计和实现，以及具体界面的设计和功能。 二、系统的初步调查 通过实地参观和学习，对超市的整体情况进行调研。了解超市的组织机构划分，充分了解超市进销存的流程的整体情况，对开发新系统的态度等。通过召开座谈会和个人访谈方法了解各个部门的主要职能及具体运作方式、过程等。 进行初步调研的具体内容为： （1）员工的规模：大约有多少员工，有多少是稳定的，有多少是浮动的； （2）员工管理人员的数量； （3）超市的商品销售状况 （4）客户编码方式； 三、可行性分析： 1.技术可行性研究，在IT行业中从业的工作人员一般都要求掌握计算机技术，具有一定的软硬件基础，会使用各种管理软件，熟悉IT产品。因为，有的超市对员工的素质要求比较高，从管理层到下面的销售人员，都要求具有一定的计算机基础，所以在新系统投入使用时，只要对员工进行少量的培训，系统的功能和使用方法就基本上能够是系统顺利运行。 2经济可行性研究，因为通过网络传递销售信息可以不受距离的限制，因此可以借阅许多的人力和物力，方便管理，由此可以减少不必要的开支，同时该系统可以提高超市的销售效率，即提高了超市的经济效益，所以从经济上完全是可行的，(1)超市有能力承担系统开发费用，(2)新系统将为企业带来经济效益3操作可行性研究，本系统采用基于Windows的图形用户界面，而该系统是大家熟悉的操作系统，对于那些有一般的计算机知识的人员就可以轻松上手。而整个超市管理系统采用最友好的交互界面，简介明了，不需要对数据库进行深入的
数字信号处理实验四 第一题结果： （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))'); 单位脉冲响应曲线 幅频和相频特性曲线 第二题结果：
实验一 时域离散信号与系统变换域分析 一、实验目的 1．了解时域离散信号的产生及基本运算实现。 2．掌握离散时间傅里叶变换实现及系统分析方法。 3. 熟悉离散时间傅里叶变换性质。 4. 掌握系统Z 域分析方法。 5. 培养学生运用软件分析、处理数字信号的能力。 二、实验设备 1、计算机 2、Matlab7.0以上版本 三、实验内容 1、对于给定的时域离散信号会进行频谱分析，即序列的傅里叶变换及其性质分析。 2、对于离散系统会进行频域分析及Z 域分析。包括频谱特性、零极点画图、稳定性分析。 3、对于差分方程会用程序求解，包括求单位冲击序列响应，零输入响应、零状态响应、全响应，求其系统函数，及其分析。 4、信号时域采样及其频谱分析，序列恢复。 5、扩展部分主要是关于语音信号的读取及其播放。 四、实验原理 1、序列的产生及运算 在Matlab 中自带了cos 、sin 、exp （指数）等函数，利用这些函数可以产生实验所需序列。 序列的运算包括序列的加法、乘法，序列)(n x 的移位)(0n n x -，翻褶)(n x -等。序列的加法或乘法指同序号的序列值逐项对应相加或相乘，但Matlab 中“+”“.*”运算是对序列的值直接进行加或乘，不考虑两序列的序号是否相同，因此编程时考虑其序号的对应。 2、序列的傅里叶变换及其性质 序列的傅里叶变换定义：)(|)(|)()(ω?ωωω j j n n j j e e X e n x e X ==∑∞-∞=-，其幅度特性为|)(|ωj e X ， 在Matlab 中采用abs 函数；相位特性为)(ω?，在Matlab 中采用angle 函数。 序列傅里叶变换的性质：
实验报告 课程：管理信息系统 一、实验目的 验证有关概念和理论，加深对概念和知识的理解和认识；熟悉和掌握Visual Basic 6.0 软件的使用方法；初步具备信息管理知识和制作数据字典、系统数据流程图的能力。运用课程讲授的管理信息系统的系统分析方法、模块化系统设计方法以及系统的调试方法进行人事档案管理信息系统的分析、设计、开发、实现与调试。 二、实验方法 面向对象法 三、实验环境及开发工具 1.硬件环境 在最低配置的情况下，系统的性能往往不尽如人意，但现在的硬件性能已经相当的出色，而且价格便宜，因此通常给服务器的配置高性能的硬件。 处理器：Interl Pentium II 266 MX 或更高 内存：64M 硬盘空间：2 GB 显卡：SVGA 显示适配器 显示器：液晶17寸 2.软件环境 操作系统：Windows/98/ME/2000/XP或更高版本 数据库：Microsoft Access 2000 3.实验开发工具：Visual Bisic 6.0程序系统 四、实验内容
（一）、系统分析 1、系统数据流程图 2、数据字典 3、系统中所有实体（包括实体的属性）以及实体之间的联系类型分析 人员的个人资料经过专业的处理部门的处理形成个人档案。档案包括自然情况，工作情况，简历，政治情况等各方面信息，内容比较庞大复杂。将档案信息传送到人员信息库。同时还综合考虑档案管理工作的性质，总结归纳出所需实现
的功能。为人事档案进行服务，对人事的变动、人事资料、以及人事资料的查询，统计等功能。总体上说具有编辑，查询，用户管理，图表统计等功能。然后将最终结果提交到人力资源管理部门，由人力资源管理人员进行审查，以便于对职工的调配。 4、典型处理的表达 档案完整添加用户档案到档案库 个人信息成功添加到档案库 修改用户档案信息 失败退回用户档案 退回用户档案 （二）、系统设计 1、子系统划分（或功能划分或模块划分） 功能划分 1、用户管理 功能：设置使用人事管理系统的用户及其使用权限。整个人事管理系统由多个功能模块组成，不同的模块完成不同的功能，所以可以为不同的职工分配不同的功能，使其具有不同的权限，完成其权限所对应的功能，从而很好地管理好整个系统。 2、辅助表管理 功能：通过它的这个功能可以有效的对本单位人事部门的扩充进行及时的计算机管理。只要管理员进行简单的数据字段添加即可。辅助表管理功能是高级管理员及中级管理员拥有的权限，它的功能是对数据库进行新表的添加。 3、档案编辑 功能：档案编辑模块中有4个子模块。他们是档案卡片、个人简历、家庭成员、历史档案等功能。这些功能因管理员的权限不同所表示出的功能使用也不同，普通管理员没有数据修改及删除的权利。在这些功能里详细的记录了所有单位员工的资料。 4、档案查询 功能：对档案卡片的查询功能，在这里可以查到符合程序要求的任何信息。
实验3 离散系统的变换域分析 一、实验目的： 加深对离散系统的频率响应分析和零、极点分布的概念理解。 二、实验原理： 离散系统的时域方程为 ∑∑==-=-M m m N k n m n x b k n y a 00)()( 其变换域分析方法如下： X(z)H(z)Y(z) )()()()()(=?-= *=∑∞-∞=m m n h m x n h n x n y 系统函数为 N N M M z a z a a z b z b b z X z Y z H ----++++++==......)()()(110110 分解因式 ∏∏∑∑=-=-=-=---==k k M m m k k k M m m m z d z c K z a z b z H 111100)1() 1()( ， 其中 m c 和 k d 称为零、极点。 在MATLAB 中，可以用函数[z ，p ，K]=tf2zp （num ，den ）求得有理分式形式的系统函数的零、极点，用函数zplane （z ，p ）绘出零、极点分布图；也可以用函数zplane （num ，den ）直接绘出有理分式形式的系统函数的零、极点分布图。使h=freqz(num,den,w)函数可求系统的频率响应，w 是频率的计算点，如w=0:pi/255:pi, h 是复数，abs(h)为幅度响应，angle(h)为相位响应。另外，在MATLAB 中，可以用函数 [r ，p ，k]=residuez （num ，den ）完成部分分式展开计算；可以用函数sos=zp2sos （z ，p ，K ）完成将高阶系统分解为2阶系统的串联。 （在实验报告中对这几种函数的使用方法及参数含义做出说明，这一部分手写） 三、实验内容 例1 求下列直接型系统函数的零、极点，并将它转换成二阶节形式 解 用MATLAB 计算程序如下： num=[1 -0.1 -0.3 -0.3 -0.2]; den=[1 0.1 0.2 0.2 0.5]; [z,p,k]=tf2zp(num,den); disp('零点');disp(z); disp('极点');disp(p); disp('增益系数');disp(k); sos=zp2sos(z,p,k);