实验五 信号与系统的复频域分析
王靖
08通信 12号
实验目的
（1）掌握利用MA TLAB 进行连续时间信号与系统的复频域分析。
（2）掌握利用MA TLAB 进行离散系统的复频域分析。
实验环境
安装MATLAB7.0以上版本的计算机
实验内容
1. 利用help 命令了解以下命令的基本用法
residue ，roots ，pzmap ，cart2pol ，residuez ，tf2zp ，zplane
2. 部分分式展开的MATLAB 实现
用部分分式展开法求X(s)的反变换。
2321
()452s X s s s s +=+++
步骤一：建立新的m 文件，保存并命名为program1.m 。
步骤二：输入以下命令，理解每条命令的含义。 %program1，部分分式展开法求反变换
[10
1];[1452];[,,](,)
n u m d en r p k resid u e n u m d en ===
步骤三：保存程序并运行，记录得到的结果。
如右图所示
步骤四：由得到的结果可以直接获得X(s)展开表示式
25
4
2
()21(1)X s s s s =-++++：
步骤五：由此可得到X(s)反变换的原函数，记录。
X(t)=(5exp(-2*t)-4exp(-t)+2texp(-t))
思考：将其转换成极坐标形式，应该如何使用cart2pol 命令？离散系统的部分分式展开，如何使用命
令residuez ，得到的结果如何利用？
将笛卡尔坐标转化为极坐标用 [angle,mag]=cart2pol(real(r),imag(r))
[r,p,k] = residuez(nun,,den)
3. H(s)的零极点与系统特性的MATLAB 计算
假设系统函数为：
321
()221H s s s s =+++
试作出零极点分布图，求系统的单位冲激响应h(t)和系统的幅度响应|H(jw)|，并判定系统的稳定性。 步骤一：建立m 文件，保存并命名为program2.m 。
步骤二：输入以下命令，理解每条命令的含义。
%program2
[1];
[1221];
(,);
(1);();
0:0.02:10;
0:0.02:5;
(,,);
(2);(,);
('Im R e ')
(,,);
(3);(,());
('n u m d en sys tf n u m d en fig u re p zm a p sys t w h im p u lse n u m d en t fig u re p lo t t h title p u lse sp o n e H freq s n u m d en w fig u re p lo t w a b s H title M a g n i ======-=R e ')tu d e sp o n e -
步骤三：保存、运行。记录输出波形。判断系统的稳定性，观察冲激响应与幅度响应的波形。
4. H(z)的零极点计算
求解因果离散时间LTI 系统的零极点：
1231230.1453(133)()10.16280.34030.0149z
z z H z z z z -------+-=+++
步骤一：建立m 文件，保存并命名为program3.m 。
步骤二：将系统函数改写为z 的正幂形式，利用tf2zp 函数求解系统的零极点。
32320.1453(331)
()0.16280.34030.0149z z z H z z z z -+-=+++
步骤三：输入以下命令，理解每条命令的含义。
%program3
0.1453*[1331];
[10.16280.3403
0.0149];[,,]2(,)
n u m d en z p k tf zp n u m d en =---== 步骤四：保存、运行，记录运行结果。由结果写出系统的零极点。(如右图)
步骤五：使用命令：zplane(num,den)，可直接在z 平面画出单位圆、零点和极点 。比较理论值和作出图形是否一致。
相关练习：
考虑如何求解因果离散时间LTI 系统的单位脉冲响应和系统的幅度响应。
答：利用h=impz （num,den,t ）求解脉冲响应利用stem(t,h)画图，用[H,w]=freqz(num,den)
求解幅度响应，利用plot(w/pi,abs(H))画出幅度响应
西华大学实验报告（理工类） 开课学院及实验室： 实验时间 ： 年 月 日 一、实验目的 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 。毒=：
实验五 信号与系统的复频域分析 王靖 08通信 12号 实验目的 （1）掌握利用MA TLAB 进行连续时间信号与系统的复频域分析。 （2）掌握利用MA TLAB 进行离散系统的复频域分析。 实验环境 安装MATLAB7.0以上版本的计算机 实验内容 1. 利用help 命令了解以下命令的基本用法 residue ，roots ，pzmap ，cart2pol ，residuez ，tf2zp ，zplane 2. 部分分式展开的MATLAB 实现 用部分分式展开法求X(s)的反变换。 2321 ()452s X s s s s +=+++ 步骤一：建立新的m 文件，保存并命名为program1.m 。 步骤二：输入以下命令，理解每条命令的含义。 %program1，部分分式展开法求反变换 [10 1];[1452];[,,](,) n u m d en r p k resid u e n u m d en === 步骤三：保存程序并运行，记录得到的结果。 如右图所示 步骤四：由得到的结果可以直接获得X(s)展开表示式 25 4 2 ()21(1)X s s s s =-++++： 步骤五：由此可得到X(s)反变换的原函数，记录。 X(t)=(5exp(-2*t)-4exp(-t)+2texp(-t)) 思考：将其转换成极坐标形式，应该如何使用cart2pol 命令？离散系统的部分分式展开，如何使用命 令residuez ，得到的结果如何利用？ 将笛卡尔坐标转化为极坐标用 [angle,mag]=cart2pol(real(r),imag(r)) [r,p,k] = residuez(nun,,den)
实验一 离散时间系统的时域分析 一、实验目的 １. 运用MATLAB 仿真一些简单的离散时间系统，并研究它们的时域特性。 ２. 运用MATLAB 中的卷积运算计算系统的输出序列，加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。 二、实验原理 离散时间系统其输入、输出关系可用以下差分方程描述： ∑=∑=-=-M k k N k k k n x p k n y d 00] [][ 当输入信号为冲激信号时，系统的输出记为系统单位冲激响应 ][][n h n →δ，则系统响应为如下的卷积计算式： ∑∞ -∞=-= *=m m n h m x n h n x n y ][][][][][ 当h[n]是有限长度的（n ：[0，M]）时，称系统为FIR 系统；反之，称系统为IIR 系统。在MA TLAB 中，可以用函数y=Filter(p,d,x) 求解差分方程，也可以用函数 y=Conv(x,h)计算卷积。 例1 clf; n=0:40; a=1;b=2; x1= 0.1*n; x2=sin(2*pi*n); x=a*x1+b*x2; num=[1, 0.5,3]; den=[2 -3 0.1]; ic=[0 0]; %设置零初始条件 y1=filter(num,den,x1,ic); %计算输入为x1(n)时的输出y1(n) y2=filter(num,den,x2,ic); %计算输入为x2(n)时的输出y2(n) y=filter(num,den,x,ic); %计算输入为x (n)时的输出y(n) yt= a*y1+b*y2; %画出输出信号 subplot(2,1,1) stem(n,y); ylabel(‘振幅’)； title(‘加权输入a*x1+b*x2的输出’)；
学生实验报告 （理工类） 课程名称：信号与线性系统专业班级：M11通信工程 学生学号：1121413017 学生姓名：王金龙 所属院部：龙蟠学院指导教师：杨娟
20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制 实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求
实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。
实验项目名称：常用连续信号的表示 实验学时： 2学时 同组学生姓名： 无 实验地点： A207 实验日期： 11.12.6 实验成绩： 批改教师： 杨娟 批改时间： 一、实验目的和要求 熟悉MATLAB 软件；利用MATLAB 软件，绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机，装有MATLAB7.0软件 三、实验过程 1. 绘制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?=； 2. 绘制指数信号at Ae t (f =），其中A=1，0.4a -=； 3. 绘制矩形脉冲信号，脉冲宽度为2； 4. 绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5； 5. 对上题三角波脉冲信号进行尺度变换，分别得出)2t (f ，)2t 2(f -； 6. 绘制抽样函数Sa （t ）,t 取值在-3π到+3π之间； 7. 绘制周期矩形脉冲信号，参数自定； 8. 绘制周期三角脉冲信号，参数自定。 四、实验结果与分析 1．制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?= 实验代码： A=1;
管理信息系统实验报告 实验3 系统分析 课程名称：管理信息系统 指导教师：王玮 班级：信管1401 学号： 姓名：唐赛赛 时间： 2016.04.06 地点： 3 号机房
一、实验目的 1．了解开发Visio解决方案的基本概念和关于Visio工具的一些基本的操作和应用； 2．掌握系统分析阶段数据流程图的画法； 二、实验步骤和实验结果： 使用Visio中提供的“组织结构图”模具，绘制下面例题的组织结构图，附在图后。 2、使用Visio绘制“业务流程图模具”和“数据流程图模具”（1）创建“业务流程图模具” 先在“框图”－〉“基本形状”中找到圆角矩形，右击选择“添加到我的形状”－〉“添加到新模具”。之后出现“另存为”对话框，把新模具命名为“业务流程图”，把圆角矩形形添加到了新模具“业务程图”中。用同样的思路，先在“框图”－〉“基本形状”中找到圆形，右击选择“添加到我的形状”－〉“添加到模具“业务程图”中；在“框图”－〉“基本形状”找到矩形，在“流程图”中的“IDEFO图表形状”找到动态连接线，在“流程图”中的“SDL图表形状”中找到文档，多文档，添加到模具“业务程图”中。可以通过设置“动态连接线”属性来改变其形状。如下图：
添加完成后，我们就可以在画业务流程图时打开该模具，业务流程图所有的元素都会在一个模具中显示出来。（2）创建“数据流程图模具”先在“框图”－〉“基本形状”中找到圆形（或是“流程图”中的“混合流程图形状”中找到外部实体2 ），右击选择“添加到我的形状”－〉“添加到新模具”（注，使用外部实体2来表示外部实体的时候，请将之旋转180度使用）。之后出现“另存为”对话框，把新模具命名为“数据流程图”，这样我们就把圆形形添加
北京理工大学信号与系统实验报告5-连续时间系统的复频域分析
实验5连续时间系统的复频域分析 （综合型实验） 一、实验目的 1）掌握拉普拉斯变换及其反变换的定义并掌握MATLAB 实现方法。 2）学习和掌握连续时间系统函数的定义及复频域分析方法。 3）掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号x(t)的拉普拉斯变换定义为 (s)(t)e st X x dt +∞ --∞ = ? （1） 拉普拉斯反变换为1 (t)(s)e 2j st j x X ds j σσπ+∞ -∞ =? （2） MATLAB 中相应函数如下： (F) L laplace = 符号表达式F 拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 (F,t)L laplace =用t 替换结果中的变量s 。 () F ilaplace L =以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。 (,) F ilaplace L x =用x 替换结果中的变量t 。
的连续时间系统，其系统函数为s 的有理函数 110 110 ...(s)...M M M M N N N N b s b s b H a s a s a ----+++= +++ （7） 3.连续时间系统的零极点分析 系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s 平面上，零点用O 表示，极点用?表示，这样得到的图形为零极点分布图。可以通过利用MATLAB 中的求多项式根的roots 函数来实现对（7）分子分母根的求解，调用格式如下： r=roots(c)，c 为多项式的系数向量，返回值r 为多项式的根向量。 求取零极点以及绘制系统函数的零极点分布图可以采用pzmap 函数，调用格式如下： pzmap(sys)绘出由系统模型sys 描述的系统的零极点分布图。 [p,z]=pzmap(sys)这种调用方式返回极点与零点，不绘出零极点分布图。 还有两个专用函数tf2zp 和zp2tf 可实现系统的传递函数模型和零极点增益模型的转换。调用格
实验五：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)如果要求用窗函数法设计带通滤波器,且给定通带上、下截止频率为和，阻带上、下截止频率为和，试求理想带通滤波器的截止频率。 答：希望逼近的理想带通滤波器的截止频率分别为：
实验5 连续时间系统的复频域分析 一、实验目的 1.掌握拉普拉斯变换及其反变换的定义，并掌握MATLAB 实现方法。 2.学习和掌握连续时间系统系统函数的定义及复频域分析方法。 3.掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号)(t x 的拉普拉斯变换定义为 )1.....(..........)()(dt e t x s X st ? +∞ ∞ --= 拉普拉斯反变换定义为 )2....(..........)(21)(ds e s X j t x j j st ?∞ +∞ -=σσπ 在MATLAB 中，可以采用符号数学工具箱的laplace 函数和ilaplace 函数进行拉氏变换和反拉氏变换。 L=laplace(F)符号表达式F 的拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 L=laplace(F,t)用t 替换结果中的变量s 。 F=ilaplace(L)以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。 F=ilaplace(L,x)用x 替换结果中的变量t 。 除了上述ilaplace 函数，还可以采用部分分式法，求解拉普拉斯逆变换，具体原理如下： 当 X (s )为有理分式时，它可以表示为两个多项式之比： )3.(..........)()()(0 110 11a s a s a b s b s b s D s N s X N N N N M M M M +?+++?++==---- 式（3）可以用部分分式法展成一下形式 )4.....(.............)(2211N N p s r p s r p s r s X -++-+-= 通过查常用拉普拉斯变换对，可以由式（1-2）求得拉普拉斯逆变换。 利用 MATLAB 的residue 函数可以将 X (s )展成式（1-2）所示的部分分式展开式，该 函数的调用格式为：[r,p,k] = residue(b,a) 其中b 、a 为分子和分母多项式系数向量，r 、p 、k 分别为上述展开式中的部分分式系数、极点和直项多项式系数。 2.连续时间系统的系统函数
天津职业技术师范大学课程设计大学学籍管理系统的设计与开发 专业：软件工程 班级学号：软件1002-17 学生姓名：靳利强 指导教师：龚良波老师 二〇一三年七月
一．需求分析 1.课程名称：大学教务信息系统的设计与开发 2.设计目的： 为方便学校做好学生学籍管理工作，设计一个学生学籍管理系统，在设计过程中作了系统分析和总体设计，软件设计采取模块化的设计思路。 3.需求概述 该学生学籍管理系统主要对学生学籍信息、成绩信息进行管理，提供一个平台，供学籍管理人员增删改查学生信息、学生成绩信息。系统分为学生信息管理、学生成绩管理、信息查询等几个模块。学籍管理人员登录成功后可以对学生信息管理、学生成绩管理、信息查询等模块进行操作，如学生信息添加、修改、删除和查询；学生成绩登记、修改、删除和查询；查询信息等。 4功能需求： 1)功能齐全：界面操作灵活方便，设计包括以下基本功能： 2)学生信息管理、教师信息管理、财务信息管理、班级信息管理、课 程信息管理、成绩信息管理、打印信息管理、教室信息管理、综合信息查询、系统管理等，至少实现其中的三个功能，且每个功能至少包括两个子功能。 3)按照软件工程的要求进行分析、设计和开发。 4)界面友好：界面友好、输入有提示、尽量展示人性化。 5)可读性强：源程序代码清晰、有层次、主要程序段有注释。
6)健壮性好：用户输入非法数据时，系统应及时给出警告信息。 二.概要设计 1.功能模块： 2数据流图： （1）学生端
（2）管理员端
学生端功能： A 登录,学生登录后，验证成功，进入其信息展示页。 管理员端功能： B 登录，管理员登录后，验证成功，进入学生信息列表，可以对学生信息进行修改，删除，按班级查询，按学号查询，按名字查询。上传图片，更新图片等操作。 三．详细设计及实现 数据库设计： 学生表： 教师表：
实验六信号与系统复频域分析 一、实验目的 1.学会用MATLAB进行部分分式展开； 2.学会用MATLAB分析LTI系统的特性； 3.学会用MATLAB进行Laplace正、反变换。 4.学会用MATLAB画离散系统零极点图； 5.学会用MATLAB分析离散系统的频率特性； 二、实验原理及内容 1．用MATLAB进行部分分式展开 用MATLAB函数residue可以得到复杂有理分式F(s)的部分分式展开式，其调用格式为 其中，num,den分别为F(s)的分子和分母多项式的系数向量，r为部分分式的系数，p为极点，k为F(s)中整式部分的系数，若F(s)为有理真分式，则k为零。 例6-1 用部分分式展开法求F(s)的反变换 解:其MATLAB程序为 format rat; num=[1,2]; den=[1,4,3,0]; [r,p]=residue(num,den) 程序中format rat是将结果数据以分数形式显示
F(s)可展开为 210.536()13 F s s s s --=++++ 所以，F(s)的反变换为 3211()()326t t f t e e u t --??=--???? 2．用MATLAB 分析LTI 系统的特性 系统函数H （s ）通常是一个有理分式，其分子和分母均为多项式。计算H （s ）的零极点可以应用MATLAB 中的roots 函数，求出分子和分母多项式的根，然后用plot 命令画图。 在MATLAB 中还有一种更简便的方法画系统函数H （s ）的零极点分布图，即用pzmap 函数画图。其调用格式为 pzmap(sys) sys 表示LTI 系统的模型，要借助tf 函数获得，其调用格式为 sys=tf(b,a) 式中，b 和a 分别为系统函数H （s ）的分子和分母多项式的系数向量。 如果已知系统函数H （s ），求系统的单位冲激响应h(t)和频 率响应H ω（j ）可以用以前介绍过的impulse 和freqs 函数。 例6-2 已知系统函数为 321221 s s s +++H(s)= 试画出其零极点分布图，求系统的单位冲激响应h(t)和频率响应H ω（j ），并判断系统是否稳定。 解：其MATLAB 程序如下： num=[1];
王锋 实验一：利用FFT 作快速相关估计 一、实验目的 a.掌握信号处理的一般方法，了解相关估计在信号分析与处理中的作用。 b.熟悉FFT算法程序；熟练掌握用FFT作快速相关估计的算法。 c.了解快速相关估计的谱分布的情况。 二、实验内容 a.读入实验数据[1]。 b.编写一利用FFT作相关估计的程序[2]。 c.将计算结果表示成图形的形式，给出相关谱的分布情况图。 注[1]：实验数据文件名为“Qjt.dat”。 实验数据来源：三峡前期工程 “覃家沱大桥” 实测桥梁振动数据。 实验数据采样频率：50Hz。 可从数据文件中任意截取几段数据进行分析，数据长度N 自定。 注[2]：采用Matlab 编程。 三、算法讨论及分析 算法为有偏估计，利用FFT计算相关函数 Step 1: 对原序列补N个零，得新序列x2N(n) Step2: 作FFT[x2N(n)]得到X2N(k) Step 3: 取X2N(k)的共轭，得 Step 4: 作 Step 5: 调整与的错位。 四、实验结果分析 1. 该信号可以近似为平稳信号么？ 可以近似为平稳信号，随机过程的统计特性不随样本的采样时刻而发生变化。取N=8192，分别取间隔m=500,m=700,m=1000,所得到的均值均为0.5366，方差为47369，与时间无关。
图1-1 自相关函数图 (上图表示的R0，下图为调整后的R0) 2. 该信号是否具有周期性，信噪比如何？ >> load Qjt.dat; %加载数据 N=32768; %数据长度 i=1:1:N; %提取数据 plot(i,Qjt(i)); 抛去几个极值点，从图1-2可以看出，数据具有一定的周期性，杂音比较少，说明信噪比较高。 图1-2 数据图
实验八系统的复频域 分析
一、实验目的 1、掌握系统的复频域分析方法。 2、掌握测试系统的频率响应的方法。 二、预习内容 1、系统频响的方法。（见第四章波特图的介绍） 三、实验原理 1. N 阶系统系统的传递函数 用微分方程描述的N 阶系统为： 根据零状态响应（起始状态为零），则对其进行拉氏变换有： 则系统传递函数可表达为： 用差分方程描述的N 阶系统为： 根据零状态响应（起始状态为零），则对其进行拉氏变换有： 则系统传递函数可表达为： 2.根据系统传递函数的零极点图分析系统 零点：传递函数分子多项式的根。 极点：传递函数分母多项式的根。 根据零极点图的不同分布分析系统。 3.涉及到的Matlab 函数 (1)freqz 函数：实验六中出现过，可用来求单位圆上的有理z 变换的值。调用格式：同实验六 (2)zplane 函数：得到有理z 变换的零极点图。 调用格式：zplane(num,den)
其中，num和 den是按z ?1 的升幂排列的、z 变换分子分母多项式系数的行向量。 (3)roots 函数：求多项式的根。 调用格式：r=roots(c), c 为多项式系数向量；r 为根向量。 四、实验内容 1.系统零极点的求解 (1)求解系统和的零极点，验 证下面程序的运行结果，根据系统零极点图分析系统性质。 b=[1,0,-1]; a=[1,2,3,2]; zr=roots(b); pr=roots(a); plot(real(zr),imag(zr),'go',real(pr),imag(pr),'mx','markersize',12,'linewidth',2); grid; legend('零点','极点'); figure; zplane(b,a); (2)参考上述程序，绘制系统和 的零极点图，并分析系统性质。与用zplane 函数直接绘制系统零极点图（注:圆心的圆圈并非系统的零点）做比较。
本科实验报告 课程名称：系统分析与设计 实验项目：《》实验实验地点： 专业班级：学号： 学生姓名： 指导教师： 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)：文本
王锋 实验五：多种功率谱估计的比较 一、实验目的 a.了解功率谱估计在信号分析中的作用； b.掌握随机信号分析的基础理论，掌握参数模型描述形式下的随机信 号的功率谱的计算方法； c.掌握在计算机上产生随机信号的方法； d.了解不同的功率谱估计方法的优缺点。 二、实验准备 有三个信号源，分别代表三种随机信号（序列）。 信号源1： 123()2cos(2)2cos(2)2cos(2)()x n f n f n f n z n πππ=+++ 其中，1230.08,=0.38,0.40f f f == z(n)是一个一阶 AR 过程，满足方程： ()(1)(1)()z n a z n e n =--+ (1)0.823321a =- e(n)是一高斯分布的实白噪声序列，方差20.1σ= 信号源2和信号源3： 都是4阶的AR 过程，它们分别是一个宽带和一个窄带过程，满足方程： ()(1)(1)(2)(2)(3)(3)(4)(4)()x n a x n a x n a x n a x n e n =--------+ e(n)是一高斯分布的实白噪声序列，方差2σ，参数如下： 三、实验内容 a. 描绘出这三个实验信号的真实功率谱波形。 b. 在计算机上分别产生这个三个信号，令所得到的数据长度 N= 256 。 注意：产生信号的时候注意避开起始瞬态点。例如，可以产生长度为512 的信号序列，然后取后面256 个点作为实验数据。 c. 分别用如下的谱估计方法，对三个信号序列进行谱估计。 1、经典谱估计 周期图法 自相关法 平均周期图法（Bartlett 法）
Welch法（可选每段64 点，重叠32 点，用Hamming 窗）2、现代谱估计 Yule - Walker方程（自相关法） 最小二乘法 注：阶次p可在3－20之间，由自己给定。 四、实验结果分析 生成的信号源
实验四：连续系统的复频域分析 一、实验目的： 1、掌握连续与离散时间系统的正反复频域与Z域变换 2、掌握利用MATLAB进行零极点分析，进一步了解零极点对整个系统的影响 3、掌握simulink环境下系统建模与仿真以及系统求解。 二、实验内容： 1、已知某连续系统的系统函数为： （1）利用[r, p, k]=residue(num, den),求H（s）的极零点以及多项式系数; （2）画出系统的零极点分布图，判断系统得稳定性。 （3）求h(t)，判断系统得稳定性。 2、已知某离散系统的系统函数为：， （1）利用[r, p, k]=residuez(num, den)求H（z）的极零点以及多项式系数; （2）画出零极点分布图，判断系统得稳定性。 （3）求单位函数响应用impz(b, a)，判断系统是否稳定； 3、已知线性时不变微分方程 在Simulink环境下搭建起系统的仿真模型，并查看仿真结果曲线。（1）写出传递函数H（s），绘出系统模拟框图； （2）当f(t)分别为，，的零状态响应；且当与课本P81的结果进行比较（3）方程的初值为, ,求全响应； 4、已知某信号，n(t)为正态噪声干扰且服从N(0,0.22)分布，对此信号进行采样，采样间隔为0.001s，之后对此信号进行Botterworth低通滤波，从信号中过滤10HZ的输出信号，试对系统进行建模与仿真。 三、实验数据处理与结果分析： 第一题：题1_1：
>> num=[2,5]; den=[1,1,3,2]; [r,p,k]=residue(num,den) r = -0.5750 - 0.7979i -0.5750 + 0.7979i 1.1499 p =-0.1424 + 1.6661i -0.1424 - 1.6661i -0.7152 k =[]
《管理信息系统》 实验二 题目：系统分析 专业：信息管理与信息系统 班级：1106班 姓名 ************************* 指导教师：贺玉珍老师 完成日期：2014.4.28
运城学院超市管理系统设计分析说明书 一、系统目标：随着小超市规模的发展不断扩大，商品数量急剧增加，有关商品的各种信息量也成倍增长。超市时时刻刻都需要对商品各种信息进行统计分析。而大型的超市管理系统功能过于强大而造成操作繁琐降低了小超市的工作效率。 超市管理系统是市场上最流行的超市上常用的系统之一，它主要包含以下几个模块：系统权限的设定、原始数据录入、数据的汇总及查询等。从而，实现对进货、销售及员工信息等实现全面、动态、及时的管理。 本文系统的分析了软件开发的背景以过程；首先介绍了软件的开发环境，其次介绍了本软件的详细设计过程：数据库的设计、各个模块的设计和实现，以及具体界面的设计和功能。 二、系统的初步调查 通过实地参观和学习，对超市的整体情况进行调研。了解超市的组织机构划分，充分了解超市进销存的流程的整体情况，对开发新系统的态度等。通过召开座谈会和个人访谈方法了解各个部门的主要职能及具体运作方式、过程等。 进行初步调研的具体内容为： （1）员工的规模：大约有多少员工，有多少是稳定的，有多少是浮动的； （2）员工管理人员的数量； （3）超市的商品销售状况 （4）客户编码方式； 三、可行性分析： 1.技术可行性研究，在IT行业中从业的工作人员一般都要求掌握计算机技术，具有一定的软硬件基础，会使用各种管理软件，熟悉IT产品。因为，有的超市对员工的素质要求比较高，从管理层到下面的销售人员，都要求具有一定的计算机基础，所以在新系统投入使用时，只要对员工进行少量的培训，系统的功能和使用方法就基本上能够是系统顺利运行。 2经济可行性研究，因为通过网络传递销售信息可以不受距离的限制，因此可以借阅许多的人力和物力，方便管理，由此可以减少不必要的开支，同时该系统可以提高超市的销售效率，即提高了超市的经济效益，所以从经济上完全是可行的，(1)超市有能力承担系统开发费用，(2)新系统将为企业带来经济效益3操作可行性研究，本系统采用基于Windows的图形用户界面，而该系统是大家熟悉的操作系统，对于那些有一般的计算机知识的人员就可以轻松上手。而整个超市管理系统采用最友好的交互界面，简介明了，不需要对数据库进行深入的
实验一 LabVIEW中的信号分析与处理 一、实验目的： 1、熟悉各类频谱分析VI的操作方法； 2、熟悉数字滤波器的使用方法； 3、熟悉谐波失真分析VI的使用方法。 二、实验原理： 1、信号的频谱分析是指用独立的频率分量来表示信号；将时域信号变换到频域，以显示在时域无法观察到的信号特征，主要是信号的频率成分以及各频率成分幅值和相位的大小，LabVIEW中的信号都是数字信号，对其进行频谱分析主要使用快速傅立叶变换(FFT)算法：·“FFT Spectrum(Mag-Phase).vi”主要用于分析波形信号的幅频特性和相频特性，其输出为单边幅频图和相频图。 ·“FFT.vi”以一维数组的形式返回时间信号的快速傅里叶运算结果，其输出为双边频谱图，在使用时注意设置FFT Size为2的幂。 ·“Amplitude and Phase Spectrum .vi”也输出单边频谱，主要用于对一维数组进行频谱分析，需要注意的是，需要设置其dt（输入信号的采样周期）端口的数据。 2、数字滤波器的作用是对信号进行滤波，只允许特定频率成份的信号通过。滤波器的主要类型分为低通、高通、带通、带阻等，在使用LabVIEW中的数字滤波器时，需要正确设置滤波器的截止频率（注意区分模拟频率和数字频率）和阶数。 3、“Harmonic Distortion Analyzer .vi”用于分析输入的波形数据的谐波失真度（THD），该vi还可分析出被测波形的基波频率和各阶次谐波的电平值。 三、实验容： (1) 时域信号的频谱分析 设计一个VI，使用4个Sine Waveform.vi（正弦波形）生成频率分别为10Hz、30Hz、50Hz、100Hz，幅值分别为1V、2V、3V、4V的4个正弦信号（采样频率都设置为1kHz，采样点数都设置为1000点），将这4个正弦信号相加并观察其时域波形，然后使用FFT Spectrum(Mag-Phase).vi对这4个正弦信号相加得出的信号进行FFT频谱分析，观察其幅频和相频图，并截图保存。
实验报告 课程：管理信息系统 一、实验目的 验证有关概念和理论，加深对概念和知识的理解和认识；熟悉和掌握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、档案查询 功能：对档案卡片的查询功能，在这里可以查到符合程序要求的任何信息。
信号分析与处理实验实验一时间信号的产生 班级：自动化1101班 姓名：陈宝平 学号： 成绩：
1. 实验目的 数字信号处理系统中的信号都是以离散时间形态存在的。研究离散时间信号，首先 需要产生出各种离散时间信号。使用MATLAB 软件可以方便的产生各种常见的离散时间信号，还具有强大的绘图功能，便于用户直观地处理输出结果。 通过本实验，学习用MATLAB 产生一些常见的离散时间信号，并通过MATLAB 绘图工具对产生的信号进行观察加深对常见离散信号和信号卷积和运算的理解。 2. 实验原理 离散时间信号用x(n)来表示，自变量n 必须是整数；连续时间信号用x(t)来表示。常见的时间信号如下： (1) 单位冲激序列δ(n)=???≠=0,00,1n n ; 如图(a); 单位冲激信号δ(t)=? ??≠=0,00 ,1t t ;如图(b): (a) (b) 如果δ(n)在时间轴上延迟了k 个单位，得到δ(n-k)=? ??≠=k n k n ,0,1。 (2) 单位阶跃序列u(n)=???<≥0,00,1n n ; 如图(c)：单位阶跃信号u(t)=???<≥0 ,00 ,1t t ; 如图(d): (c) (d) 如果u(n)在时间轴上延迟了k 个单位,得到u(n-k)=???<≥k n k n ,0,1。
(3) 矩形序列R N (n)=???≥<-≤≤),0(,0) 10(,1N n n N n ,矩形序列有一个重要的参数，就是序列宽度 N 。R N (n)与u(n)之间的关系为R N (n)= u(n)- u(n-N)。如图(e): 单位矩形信号R (t)=? ??≥<≤≤),0(,0) 0(,1T t t T t , R(t)=u(t)-u(t-T),如图(f): (e) (f) （4）正弦序列x(n)=Acos(ω0n+?)。只有当 2ωπ 为有理数时，正弦序列具有周期性, 如图(g): 正弦信号x(t)=Acos(?ω+t 0),如图(h); (g) (h) (5)单边实指数序列x(n)=a n u(n)，当a>0时，该序列均取正值，当a<0时，序列在正负摆动。如图分别为：x(n)=1.2n , x(n)=(-1.2)n , x(n)=0.8n , x(n)=(-0.8)n 的图。
信号与系统 实验报告 实验五连续系统的复频域分析 实验五连续系统的复频域分析 一、实验目的 1. 深刻理解拉普拉斯变换、逆变换的定义，掌握用MATLAB实现拉普拉斯变换、逆变换的方法。 2会求几种基本信号的拉氏变换。 3 掌握用MATLAB绘制连续系统零、极点的方法。 4 求解系统函数H(s)。 二
1已知连续时间信号f（t）=sin（t）u（t）、求出该信号的拉普拉斯变换，并用MATLAB 绘制拉普拉斯变换的曲面图。 syms t; ft=sin(t)*heaviside(t); Fs=Laplace(ft); a=-0.5:0.08:0.5; b=-2:0.08:2; [a,b]=meshgrid(a,b); c=a+i*b; d=ones(size(a)); c=c.*c; c=c+d; c=1./c; c=abs(c); mesh(a,b,c); surf(a,b,c) axis([-0.5,0.5,-2,2,0,10]) colormap(hsv
) 2求[(1-e^(-at))]/t的拉氏变换。 syms t s a f1=(1-exp(-a*t))/t; F=laplace(f1,t,s) F = log(s+a)-log(s) 3求F（s）=-log（s）+ log（s+a）的拉氏逆变换syms t s a F =log(s+a)-log(s); f1=ilaplace(F,s,t) f1 = (1-exp(-a*t))/t
4已知某连续系统的系统函数为： H（s）=(s^2+3s+2)/(8s^4+2s^3+3s^2+5)试用MATLAB求出该系统的零极点，画出零极点分布图。 b=[1 3 2]; a=[8 2 3 0 5]; zs=roots(b); ps=roots(a); hold on plot(real(zs),imag(zs),'o'); plot(real(ps),imag(ps),'x'); grid axis([-2.5,1,-1,1]) 5已知H(s)=（s+1）/（s^2+s+1）,绘制阶跃响应图形，冲激响应图形，频率激响应图形。 syms t s H=(s+1)/(s^2+s+1); f1=ilaplace(H,s,t); f2=heaviside(t);