仿真实验报告
系统工程仿真实验报告
实验一：基于VENSIM的系统动力学仿真
一、实验目的
VENSIM是一个建模工具，可以建立动态系统的概念化的，文档化的仿真、分析和优化模型。PLE（个人学习版）是VENSIM的缩减版，主要用来简单化学习动态系统，提供了一种简单富有弹性的方法从常规的循环或储存过程和流程图建立模型。本实验就是运用VENSIM进行系统动力学仿真，进一步加深对系统动力学仿真的理解。
二、实验软件
VENSIM PLE
三、原理
1、在VENSIM中建立系统动力学流图；
2、写出相应的DYNAMO方程；
3、仿真出系统中水准变量随时间的响应趋势；
四、实验内容及要求
某城市国营和集体服务网点的规模可用SD来研究。现给出描述该问题的DYNAMO方程及其变量说明。
L S·K=S·J+DT*NS·JK
N S=90
R NS·KL=SD·K*P·K/(LENGTH-TIME·K)
A SD·K=SE-SP·K
C SE=2
A SP·K=SR·K/P·K
A SR·K=SX+S·K
C SX=60
L P·K=P·J+DT*NP·JK
N P=100
R NP·KL=I*P·K
C I=0.02
其中：LENGTH为仿真终止时间、TIME为当前仿真时刻，均为仿真控制变量；S为个体服务网点数(个)、NS为年新增个体服务网点数(个/年)、SD为实际千人均服务网点与期望差(个/千人)、SE为期望的千人均网点数、SP为的千人均网点数(个/千人)、SX为非个体服务网点数(个)、SR为该城市实际拥有的服务网点数(个)、P为城市人口数(千人)、NP为年新增人口数(千人/年)、I为人口的年自然增长率。
要求：在VENSIM中建立相应的系统动力学流图和DYNAMO方程，进行仿真。
五、实验结果
1、请将VENSIM 中建立的系统动力学流图拷贝如下：
2、画出系统中水准变量随时间的响应趋势。
P
800
600
400
200
0P : Current
S
2,000
1,500
1,000
500
S : Current
实验二：基于EXCEL的随机存储系统仿真
一、实验目的
理解构建仿真模型的步骤；学习利用Excel中的RAND函数产生随机数，利用VLOOKUP 函数找出对应的需求；熟悉应用仿真模型解决库存管理问题；了解利用Excel软件进行仿真的优缺点。
二、实验软件
MICROSOFT Excel
三、实验内容及要求
马克是一家五金商店的店主，他想为电钻制定一种优良的低成本的存货策略。电钻的日需求相对比较低，但同时也比较容易受到变动的影响。在过去的300天里，马克观察到的销售量如表所示。
表1 电钻日需求的概率和随机数区间
当马克通过订货来补充存货时，货物的发送往往会有1-3天的延时。以往接受50笔订单所耗费的天数如表2所示。如果订货提前期为1天，货物是在第二天晚间送达。
表2 再订货的订货提前期概率和随机数区间
马克想要模拟的第一种存货策略是订货量为10，再订货点为5，且考虑在途货物数量，上一订单未到货时不连续订货。假设期初的库存量为10台。电钻的每笔订货成本为10美元，存储电钻每台每天5美分，缺货成本每次6美元。将仿真过程填写在表3中，并计算该存货策略的总成本。
表3 五金商店存货仿真
四、实验步骤
1、定义问题。
2、定义可控和不可控的输入变量。可控变量为订货点和订货量，不可控变量为日需求和订货提前期。
3、建立仿真模型。
4、实现仿真。假设期初的库存量为10台，订货量为10台，再订货点为5台，利用Excel RAND()函数生成随机数，将仿真过程记录于表3中。
5、检验结果。对30天的期间进行仿真，计算订货成本、储存成本、缺货成本以及总的库存成本，并给出更好的订货策略建议。
五、实验结果
1、参数设置的截图：
需求概
下限上线取值
率
0.05 0 0.05 0
0.1 0.05 0.15 1
0.2 0.15 0.35 2
0.4 0.35 0.75 3
0.15 0.75 0.9 4
0.1 0.9 1 5
提前概
下限上线天数
率
0.2 0 0.2 1
0.5 0.2 0.7 2
0.3 0.7 1 3
订货量再订货点
9 5
2、请将仿真界面的截图拷贝如下：
3、更好的订货策略建议。1增加订货量
2降低在订货点
《计算方法》课内实验报告
《计算方法》实验报告 姓名： 班级： 学号： 实验日期: 2011年10月26日
一、实验题目： 数值积分 二、实验目的： 1．熟悉matlab 编写及运行数值计算程序的方法。 2．进一步理解数值积分的基础理论。 3．进一步掌握应用不同的数值积分方法求解给定的积分并给出数据结果及误差分析。 三、实验内容： 1．分别用复合梯形求积公式及复合辛普森求积公式计算积分xdx x ln 10 ? , 要求计算精度达到410-,给出计算结果并比较两种方法的计算节点数. 2．用龙贝格求积方法计算积分dx x x ?+3 021,使误差不超过510-. 3．用3=n 的高斯-勒让德公式计算积分?3 1 sin x e x ,给出计算结果. 4．用辛普森公式(取2==M N ) 计算二重积分.5 .00 5 .00 dydx e x y ? ? - 四、实验结果： 1.（1）复合梯形法: 将区间[a,b]划分为n 等份，分点n k n a b h kh a x k ,2,1,0,,=-=+=在每个区间[1,+k k x x ](k=0,1,2,···n-1)上采用梯形公式，则得 )()]()([2)()(1 11 1 f R x f x f h dx x f dx x f I n n k k k b a n k x x k k ++===∑?∑? -=+-=+ 故)]()(2)([21 1 b f x f a f h T n k k n ++=∑-=称为复合梯形公式 计算步长和划分的区间 Eps=1E-4 h1=sqrt(Eps/abs(-(1-0)/12*1/(2+1))) h1 =0.0600 N1=ceil(1/h1) N1 =17 用复合梯形需要计算17个结点。 复合梯形： function T=trap(f,a,b,n) h=(b-a)/n;
信号与系统仿真实验报告
信号与系统仿真实验报告1．实验目的 了解MATLAB的基本使用方法和编程技术，以及Simulink平台的建模与动态仿真方法，进一步加深对课程内容的理解。 2．实验项目 信号的分解与合成，观察Gibbs现象。 信号与系统的时域分析，即卷积分、卷积和的运算与仿真。 信号的频谱分析，观察信号的频谱波形。 系统函数的形式转换。 用Simulink平台对系统进行建模和动态仿真。 3．实验内容及结果 3.1以周期为T，脉冲宽度为2T1的周期性矩形脉冲为例研究Gibbs现象。 已知周期方波信号的相关参数为：x(t)=∑ak*exp(jkω),ω=2*π/T,a0=2*T1/T,ak=sin(kωT1)/kπ。画出x(t)的波形图（分别取m=1,3,7,19,79,T=4T1）,观察Gibbs现象。 m=1; T1=4; T=4*T1;k=-m:m; w0=2*pi/T; a0=2*T1/T; ak=sin(k*w0*T1)./(k*pi); ak(m+1)=a0; t=0:0.1:40; x=ak*exp(j*k'*w0*t); plot(t,real(x)); 3.2求卷积并画图 (1)已知：x1(t)=u(t-1)-u(t-2), x2(t)=u(t-2)-u(t-3)求：y(t)=x1(t)*x2(t)并画出其波形。 t1=1:0.01:2; f1=ones(size(t1)); f1(1)=0; f1(101)=0; t2=2:0.01:3; f2=ones(size(t2)); f2(1)=0; f2(101)=0; c=conv(f1,f2)/100;
t3=3:0.01:5; subplot(311); plot(t1,f1);axis([0 6 0 2]); subplot(312); plot(t2,f2);axis([0 6 0 2]); subplot(313); plot(t3,c);axis([0 6 0 2]); (2）已知某离散系统的输入和冲击响应分别为：x[n]=[1,4,3,5,1,2,3,5], h[n]=[4,2,4,0,4,2].求系 统的零状态响应，并绘制系统的响应图。 x=[1 4 3 5 1 2 3 5]; nx=-4:3; h=[4 2 4 0 4 2]; nh=-3:2; y=conv(x,h); ny1=nx(1)+nh(1); ny2=nx(length(nx))+nh(length(nh)); ny=[ny1:ny2]; subplot(311); stem(nx,x); axis([-5 4 0 6]); ylabel('输入') subplot(312); stem(nh,h); axis([-4 3 0 5]); ylabel('冲击效应') subplot(313); stem(ny,y); axis([-9 7 0 70]); ylabel('输出'); xlabel('n'); 3.3 求频谱并画图 (1) 门函数脉冲信号x1(t)=u(t+0.5)-u(t-0.5) N=128;T=1; t=linspace(-T,T,N); x=(t>=-0.5)-(t>=0.5); dt=t(2)-t(1); f=1/dt; X=fft(x); F=X(1:N/2+1); f=f*(0:N/2)/N; plot(f,F)
Matlab通信系统仿真实验报告
Matlab通信原理仿真 学号： 2142402 姓名：圣斌
实验一Matlab 基本语法与信号系统分析 一、实验目的： 1、掌握MATLAB的基本绘图方法； 2、实现绘制复指数信号的时域波形。 二、实验设备与软件环境： 1、实验设备：计算机 2、软件环境：MATLAB R2009a 三、实验内容： 1、MATLAB为用户提供了结果可视化功能，只要在命令行窗口输入相应的命令，结果就会用图形直接表示出来。 MATLAB程序如下： x = -pi::pi; y1 = sin(x); y2 = cos(x); %准备绘图数据 figure(1); %打开图形窗口 subplot(2,1,1); %确定第一幅图绘图窗口 plot(x,y1); %以x，y1绘图 title('plot(x,y1)'); %为第一幅图取名为’plot(x,y1)’ grid on; %为第一幅图绘制网格线 subplot(2,1,2) %确定第二幅图绘图窗口 plot(x,y2); %以x，y2绘图 xlabel('time'),ylabel('y') %第二幅图横坐标为’time’,纵坐标为’y’运行结果如下图： 2、上例中的图形使用的是默认的颜色和线型，MATLAB中提供了多种颜色和线型，并且可以绘制出脉冲图、误差条形图等多种形式图： MATLAB程序如下： x=-pi:.1:pi; y1=sin (x); y2=cos (x); figure (1); %subplot (2,1,1); plot (x,y1); title ('plot (x,y1)'); grid on %subplot (2,1,2); plot (x,y2);
太原理工大学数值计算方法实验报告
本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 2016年 4 月 8 日
y = x*x*x + 4 * x*x - 10; return y; } float Calculate(float a,float b) { c = (a + b) / 2; n++; if (GetY(c) == 0 || ((b - a) / 2) < 0.000005) { cout << c <<"为方程的解"<< endl; return 0; } if (GetY(a)*GetY(c) < 0) { return Calculate(a,c); } if (GetY(c)*GetY(b)< 0) { return Calculate(c,b); } } }; int main() { cout << "方程组为：f(x)=x^3+4x^2-10=0" << endl; float a, b; Text text; text.Getab(); a = text.a; b = text.b; text.Calculate(a, b); return 0; } 2.割线法： // 方程求根（割线法）.cpp : 定义控制台应用程序的入口点。// #include "stdafx.h" #include"iostream"
心得体会 使用不同的方法，可以不同程度的求得方程的解，通过二分法计算的程序实现更加了解二分法的特点，二分法过程简单，程序容易实现，但该方法收敛比较慢一般用于求根的初始近似值，不同的方法速度不同。面对一个复杂的问题，要学会简化处理步骤，分步骤一点一点的循序处理，只有这样，才能高效的解决一个复杂问题。
测试装置动态特性仿真实验报告
测试装置动态特性仿真实验 班级：7391 学号：2009301828 姓名：张志鹏 一、实验目的 1、加深对一阶测量装置和二阶测量装置的幅频特性与相频特性的理解； 2、加深理解时间常数变化对一阶系统动态特性影响； 3、加深理解频率比和阻尼比变化对二阶系统动态特性影响； 4、使学生了解允许的测量误差与最优阻尼比的关系。 二、实验原理 1、 一阶测量装置动态特性 一阶测量装置是它的输入和输出关系可用一阶微分方程描述。一阶测量装置的频率响应函数为： 式中：S S 为测量装置的静态灵敏度；τ为测量装置的时间常数。 一阶测量装置的幅频特性和相频特性分别为： 可知，在规定S S =1的条件下，A （ω）就是测量装置的动态灵敏度。 当给定一个一阶测量装置，若时间常数τ确定，如果规定一个允许的幅值误差ε，则允许测量的信号最高频率ωH 也相应地确定。 为了恰当的选择一阶测量装置，必须首先对被测信号的幅值变化范围和频率成分有个初步了解。有根据地选择测量装置的时间常数τ，以保证A （ω）≥1-ε 能够满足。 2、二阶测量装置动态特性 二阶测量装置的幅频特性与相频特性如下： 幅频特性202220)/(4))/(1(/1)(ωωξωωω--=A 相频特性2200))/(1/()/(2()(ωωωωξφ--=arctg w Α（ω）是ξ和ω/0ω的函数，即具有不同的阻尼比ξ的测试装置当输入信??????ωτ+ωτ-ωτ+=ωτ+=ω22s s )(1j ) (11S j 11S )j (H ()()2 11 A ωτ+=ω()ωτ -=ωφarctan
号频率相同时，应具有不同的幅值响应，反之，当不同的频率的简谐信号送入同一测试装置时它们的幅值响应也不相同，同理具有不同的阻尼比ξ的测试装置当输入信号频率相同时，应有不同的相位差。 (1).当ω=0时，Α（ω）=1;(2).当ω→∞，A （ω）=0;(3).当ξ≥0.707时随着输入信号频率的加大，Α（ω）单调的下降, ξ＜0.707时Α（ω）的特性曲线上出现峰值点;(4)如果ξ=0，))/(1/(1))/(1(/1)(202 20ωωωωω-=-=A ，显然，其峰值点出现在ω=0ω处。其值为“∞”，当ξ从0向0.707变化过程中随着的加大其峰值点逐渐左移，并不断减小。 对以上二阶环节的幅频特性的结论论证如下： (1).当ω=0时A(ω)=1 (2).当ω→∞时,A(ω)=0 (3).要想得到A(ω)的峰值就要使202220)/(4))/(1(/1)(A ωωξ-ωω-=ω 中的202220)/(4))/(1(ωωξωω--取最小值。 令:t=20)/(ωω t t t f 224)1()(ξ+-= 对其求导可得t=1-22ξ时,f(t)取最小值.由于t=20)/(ωω≥0,所以1-22ξ≥0, 2ξ必须小于1/2时,f(t)才有最小值,即ξ>2/2时，A(ω)不出现峰值点;当ξ＜2/2时4244)(ξξ-=t f ,f(t)对ξ求导得)21(82ξξ-，可以看出f(t): ξ属于[0， 2/2]时单调递增，于是得A(ω)的峰值点A 为4244/1)(/1ξξ-=t f ; 在ξ属于[0，2/2]递减。 (4).当ξ=0时 A=∞，t=20)/(ωω，ω/0ω=1,即ξ=0时A(ω)的峰值为∞，且必出现在ω/0ω=1时，当ξ=2/2时,t=0→ω=0，A(ω)=1. 还可以看出，在ξ属于[0，2/2]增大时t=1-22ξ就减小，即f(t)的峰值左平移。 （二）阻尼比的优化 在测量系统中，无论是一阶还是二阶系统的幅频特性都不能满足将信号中的所有频率都成比例的放大。于是希望测量装置的幅频特性在一段尽可能宽的范围内最接近于1。根据给定的测量误差，来选择最优的阻尼比。
c 计算器实验报告
简单计算器 姓名: 周吉祥 实验目的：模仿日常生活中所用的计算器，自行设计一个简单的计算器程序，实现简单的计算功能。 实验内容： (1)体系设计： 程序是一个简单的计算器，能正确输入数据，能实现加、减、乘、除等算术运算，运算结果能正确显示，可以清楚数据等。 (2)设计思路： 1)先在Visual C++ 6.0中建立一个MFC工程文件，名为 calculator. 2)在对话框中添加适当的编辑框、按钮、静态文件、复选框和单 选框 3)设计按钮，并修改其相应的ID与Caption. 4)选择和设置各控件的单击鼠标事件。 5)为编辑框添加double类型的关联变量m_edit1. 6)在calculatorDlg.h中添加math.h头文件，然后添加public成 员。 7)打开calculatorDlg.cpp文件，在构造函数中，进行成员初始 化和完善各控件的响应函数代码。 (3)程序清单：
●添加的public成员： double tempvalue; //存储中间变量 double result; //存储显示结果的值 int sort; //判断后面是何种运算：1.加法2.减法3. 乘法 4.除法 int append; //判断后面是否添加数字 ●成员初始化： CCalculatorDlg::CCalculatorDlg(CWnd* pParent /*=NULL*/) : CDialog(CCalculatorDlg::IDD, pParent) { //{{AFX_DATA_INIT(CCalculatorDlg) m_edit1 = 0.0; //}}AFX_DATA_INIT // Note that LoadIcon does not require a subsequent DestroyIcon in Win32 m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME); tempvalue=0; result=0; sort=0; append=0; }
仿真实验报告
大学物理仿真实验报告一一塞曼效应 一、实验简介 塞曼效应是物理学史上一个著名的实验。荷兰物理学家塞曼（Zeeman）在1896年发现把产生光谱的光源置于足够强的磁场中，磁场作用于发光体，使光谱发生变化，一条谱线即会分裂成几条偏振化的谱线，这种现象称为塞曼效应。 塞曼效应是法拉第磁致旋光效应之后发现的又一个磁光效应。这个现象的发现是对光的 电磁理论的有力支持，证实了原子具有磁矩和空间取向量子化，使人们对物质光谱、原子、分子有更多了解。 塞曼效应另一引人注目的发现是由谱线的变化来确定离子的荷质比的大小、符号。根据 洛仑兹（H.A?Lorentz）的电子论，测得光谱的波长，谱线的增宽及外加磁场强度，即可称得离子的荷质比。由塞曼效应和洛仑兹的电子论计算得到的这个结果极为重要，因为它发表在J、 J汤姆逊（J、J ThomSOn）宣布电子发现之前几个月，J、J汤姆逊正是借助于塞曼效应由洛仑 兹的理论算得的荷质比，与他自己所测得的阴极射线的荷质比进行比较具有相同的数量级，从而得到确实的证据，证明电子的存在。 塞曼效应被誉为继X射线之后物理学最重要的发现之一。 1902年，塞曼与洛仑兹因这一发现共同获得了诺贝尔物理学奖（以表彰他们研究磁场对光的效应所作的特殊贡献）。至今，塞曼效应依然是研究原子内部能级结构的重要方法。 本实验通过观察并拍摄Hg（546.1 nm）谱线在磁场中的分裂情况，研究塞曼分裂谱的特征，学习应用塞曼效应测量电子的荷质比和研究原子能级结构的方法。 二、实验目的 1?学习观察塞曼效应的方法观察汞灯发出谱线的塞曼分裂； 2?观察分裂谱线的偏振情况以及裂距与磁场强度的关系； 3?利用塞曼分裂的裂距，计算电子的荷质比 e m e数值。 三、实验原理 1、谱线在磁场中的能级分裂 设原子在无外磁场时的某个能级的能量为E0,相应的总角动量量子数、轨道量子数、 自旋量子数分别为J、L、S。当原子处于磁感应强度为B的外磁场中时，这一原子能级将 分裂为2J 1层。各层能量为 E = E o MgJ B B（1） 其中M为磁量子数，它的取值为J , J -1 ,…，-J共2J 1个；g为朗德因子；J B为 hc 玻尔磁矩（A B= ）；B为磁感应强度。 4兀m 对于L-S耦合
系统仿真实验报告
中南大学系统仿真实验报告 指导老师胡杨 实验者 学号 专业班级 实验日期 2014.6.4 学院信息科学与工程学院
目录 实验一MATLAB中矩阵与多项式的基本运算 (3) 实验二MATLAB绘图命令 (7) 实验三MATLAB程序设计 (9) 实验四MATLAB的符号计算与SIMULINK的使用 (13) 实验五MATLAB在控制系统分析中的应用 (17) 实验六连续系统数字仿真的基本算法 (30)
实验一MATLAB中矩阵与多项式的基本运算 一、实验任务 1．了解MATLAB命令窗口和程序文件的调用。 2．熟悉如下MATLAB的基本运算： ①矩阵的产生、数据的输入、相关元素的显示； ②矩阵的加法、乘法、左除、右除； ③特殊矩阵：单位矩阵、“1”矩阵、“0”矩阵、对角阵、随机矩阵的产生和运算； ④多项式的运算：多项式求根、多项式之间的乘除。 二、基本命令训练 1．eye(m) m=3; eye(m) ans = 1 0 0 0 1 0 0 0 1 2．ones(n)、ones(m,n) n=1;m=2; ones(n) ones(m,n) ans = 1 ans = 1 1
3．zeros(m,n) m=1,n=2; zeros(m,n) m = 1 ans = 0 0 4．rand(m,n) m=1;n=2; rand(m,n) ans = 0.8147 0.9058 5．diag(v) v=[1 2 3]; diag(v) ans = 1 0 0 0 2 0 0 0 3 6．A\B 、A/B、inv(A)*B 、B*inv(A) A=[1 2;3 4];B=[5 6;7 8]; a=A\B b=A/B c=inv(A)*B d=B*inv(A) a = -3 -4 4 5 b = 3.0000 -2.0000 2.0000 -1.0000
MATLAB通信系统仿真实验报告1
MATLAB通信系统仿真实验报告
实验一、MATLAB的基本使用与数学运算 目的：学习MATLAB的基本操作，实现简单的数学运算程序。 内容： 1-1要求在闭区间[0，2π]上产生具有10个等间距采样点的一维数组。试用两种不同的指令实现。 运行代码：x=[0:2*pi/9:2*pi] 运行结果： 1-2用M文件建立大矩阵x x=[0.10.20.30.40.50.60.70.80.9 1.11.21.31.41.51.61.71.81.9 2.12.22.32.42.52.62.72.82.9 3.13.23.33.43.53.63.73.83.9] 代码：x=[0.10.20.30.40.50.60.70.80.9 1.11.21.31.41.51.61.71.81.9 2.12.22.32.42.52.62.72.82.9 3.13.23.33.43.53.63.73.83.9] m_mat 运行结果： 1-3已知A=[5，6;7，8]，B=[9，10;11,12]，试用MATLAB分别计算 A+B,A*B,A.*B,A^3，A.^3，A/B,A\B. 代码：A=[56;78]B=[910;1112]x1=A+B X2=A-B X3=A*B X4=A.*B X5=A^3 X6=A.^3X7=A/B X8=A\B
运行结果： 1-4任意建立矩阵A，然后找出在[10,20]区间的元素位置。 程序代码及运行结果： 代码：A=[1252221417;111024030;552315865]c=A>=10&A<=20运行结果： 1-5总结：实验过程中，因为对软件太过生疏遇到了些许困难，不过最后通过查书与同学交流都解决了。例如第二题中，将文件保存在了D盘，而导致频频出错，最后发现必须保存在MATLAB文件之下才可以。第四题中，逻辑语言运用到了ij，也出现问题，虽然自己纠正了问题，却也不明白错在哪了，在老师的讲解下知道位置定位上不能用ij而应该用具体的整数。总之第一节实验收获颇多。
计算方法实验报告格式
计算方法实验报告格式 小组名称: 组长姓名(班号)： 小组成员姓名(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一个完整的实验，应包括数据准备、理论基础、实验内容及方法，最终对实验结果进行分析，以达到对理论知识的感性认识，进一步加深对相关算法的理解，数值实验以实验报告形式完成，实验报告格式如下： 一、实验名称 实验者可根据报告形式需要适当写出． 二、实验目的及要求 首先要求做实验者明确，为什么要做某个实验，实验目的是什么，做完该实验应达到什么结果，在实验过程中的注意事项，实验方法对结果的影响也可以以实验目的的形式列出． 三、算法描述（实验原理与基础理论） 数值实验本身就是为了加深对基础理论及方法的理解而设置的，所以要求将实验涉及到的理论基础，算法原理详尽列出． 四、实验内容 实验内容主要包括实验的实施方案、步骤、实验数据准备、实验的算法以及可能用到的仪器设备． 五、程序流程图 画出程序实现过程的流程图，以便更好的对程序执行的过程有清楚的认识，在程序调试过程中更容易发现问题． 六、实验结果 实验结果应包括实验的原始数据、中间结果及实验的最终结果，复杂的结果可以用表格
形式列出，较为简单的结果可以与实验结果分析合并出现． 七、实验结果分析 实验结果分析包括对对算法的理解与分析、改进与建议． 数值实验报告范例 为了更好地做好数值实验并写出规范的数值实验报告，下面给出一简单范例供读者参考． 数值实验报告 小组名称: 小组成员(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一、实验名称 误差传播与算法稳定性． 二、实验目的 1．理解数值计算稳定性的概念． 2．了解数值计算方法的必要性． 3．体会数值计算的收敛性与收敛速度． 三、实验内容 计算dx x x I n n ? += 1 10 ，1,2,,10n = ． 四、算法描述 由 dx x x I n n ? += 1 10 ，知 dx x x I n n ?+=--101110，则
物流仿真实验报告
《物流仿真实验》 实验报告书 实验报告题目: 物流仿真实验学院名称: 管理学院 专业: 物流管理 班级: 物流1303 姓名: 孟颖颖 学号: 2 成绩: 2016年7月 实验报告 一、实验名称 物流仿真实验 二、实验要求 ⑴根据模型描述与模型数据对配送中心进行建模;
⑵分析仿真实验结果,进行利润分析,找出利润最大化的策略。 三、实验目的 1、掌握仿真软件Flexsim的操作与应用,熟悉通过软件进行物流仿真建模。 2、记录Flexsim软件仿真模拟的过程,得出仿真的结果。 3、总结Flexsim仿真软件学习过程中的感受与收获。 三、实验设备 (1)硬件及其网络环境 服务器一台:PII400/10、3G/128M以上配置、客户机100台、局域网或广域网。 (2)软件及其运行环境 Flexsim,Windows 2000 Server、SQL Server 7、0以上版本、IIS 5、0、SQL Server 数据库自动配置、IIS 虚拟目录自动配置 四、实验步骤 1 概念模型 1个Sink到操作区,如图:
第二步:连接端口 根据配送流程,对模型进行适宜的连接,所有端口连接均用A连接,如图: 第三步:Source的参数设置 为使Source产生实体不影响后面Processor的生产,尽可能的将时间间隔设置尽可能的小,并对三个Source做出同样的设定。 打开Source参数设置窗口,将时间到达间隔设置为常数1,同时为对三个实体进行区别,进行设置产品颜色,点击触发器,打开离开触发的下拉菜单,点击设置临时实体类型,设置不同实体类型,颜色自然发生变化。并对另外两个Source 进行同样的设置,如图:
交通运输系统仿真实验报告
一、系统描述 1.1.系统背景 本系统将基于下面的卫星屏幕快照创建一个模型。当前道路网区域的两条道路均为双向，每个运动方向包含一条车道。Tapiolavagen路边有一个巴士站，Menninkaisentie路边有一个带五个停车位的小型停车场。 1.2.系统描述 （1）仿真十字路口以及三个方向的道路，巴士站，停车点；添加小汽车、公交车的三维动画，添加红绿灯以及道路网络描述符； （2）创建仿真模型的汽车流程图，三个方向产生小汽车，仿真十字路口交通运行情况。添加滑条对仿真系统中的红绿灯时间进行实时调节。添加分析函数，统计系统内汽车滞留时间，用直方图进行实时展示。 二、仿真目标 1、timeInSystem值：在流程图的结尾模块用函数统计每辆汽车从产生到丢弃的，在系统中留存的时间。 2、p_SN为十字路口SN方向道路的绿灯时间，p_EW为十字路口EW方向道路的绿灯时间。 3、Arrival rate：各方向道路出现车辆的速率（peer hour）。
三、系统仿真概念分析 此交通仿真系统为低抽象层级的物理层模型，采用离散事件建模方法进行建模，利用过程流图构建离散事件模型。 此十字路口交通仿真系统中，实体为小汽车和公交车，可以源源不断地产生；资源为道路网络、红绿灯时间、停车点停车位和巴士站，需要实施分配。系统中小汽车（car）与公共汽车（bus）均为智能体，可设置其产生频率参数，行驶速度，停车点停留时间等。 四、建立系统流程 4.1.绘制道路 使用Road Traffic Library中的Road模块在卫星云图上勾画出所有的道路，绘制交叉口，并在交叉口处确保道路连通。 4.2.建立智能体对象 使用Road Traffic Library中的Car type模快建立小汽车（car）以及公共汽车（bus）的智能体对象。 4.3.建立逻辑 使用Road Traffic Library中的Car source、Car Move To、Car Dispose、
仿真实验报告经典案例概述
XXXXX 实验报告 学院（部）XX学院 课程名称生产系统仿真实验 学生姓名 学号 专业 2012年9月10日
《生产系统仿真》实验报告 年月日 学院年级、专业、班实验时间9月10日成绩 课程名称生产系统仿真 实训项目 名称 系统仿真软件的基础应 用 指导 教师 一、实验目的 通过对Flesim软件进一步的学习，建立模型，运用Flesim软件仿真该系统，观察并分析运行结果，找出所建模型的问题并进行改进，再次运行循环往复，直到找出构建该系统更为合理的模型。 二、实验内容 １、建立生产模型。 该模型生产三种产品，产品到达速率服从均值为20、方差为2的正态分布；暂存器的最大容量为25个；检测器的检测时间服从均值为30的指数分布，预制时间为10s；传送带的传送速率为1m/s，带上可容纳的最大货件数为10个。 ２、运行生产模型。 ３、对运行结果进行分析，提出改进方案在运行，直到找到更为合理的模型。 三、实验报告主要内容 １、根据已有数据建立生产模型。 将生产系统中所需实体按组装流程进行有序的排列，并进行连接如图１所示
图1 ２、分别对发生器、暂存器、检验台和传送带进行参数设置。 （１）发生器的产品到达速率服从均值为20、方差为2的正态分布。如图2所示。 （２）暂存器的最大容量设置为25件。如图3所示。 （3）设置检验台的检测时间服从均值为30s的指数分布，预制时间为10s.如图4所示。 （4）传送带的传送速率为1m/s，最大容量为10件。如图5所示 图2 图3 图4 图5 3、对发生器及暂存器进一步设置。 （1）发生器在生成产品时设置三种不同类型的产品，通过颜色区分。如图6所示。 （2）暂存器在输出端口通过设置特定函数以使不同颜色的产品在不同的检验台检验。如图7所示。
控制系统仿真实验报告
哈尔滨理工大学实验报告 控制系统仿真 专业：自动化12-1 学号：1230130101 姓名：
一．分析系统性能 课程名称控制系统仿真实验名称分析系统性能时间8.29 地点3# 姓名蔡庆刚学号1230130101 班级自动化12-1 一．实验目的及内容： 1. 熟悉MATLAB软件的操作过程； 2. 熟悉闭环系统稳定性的判断方法； 3. 熟悉闭环系统阶跃响应性能指标的求取。 二．实验用设备仪器及材料： PC, Matlab 软件平台 三、实验步骤 1. 编写MATLAB程序代码； 2. 在MATLAT中输入程序代码，运行程序； 3.分析结果。 四．实验结果分析： 1.程序截图
得到阶跃响应曲线 得到响应指标截图如下
2.求取零极点程序截图 得到零极点分布图 3.分析系统稳定性 根据稳定的充分必要条件判别线性系统的稳定性最简单的方法是求出系统所有极点,并观察是否含有实部大于0的极点，如果有系统不稳定。有零极点分布图可知系统稳定。
二．单容过程的阶跃响应 一、实验目的 1. 熟悉MATLAB软件的操作过程 2. 了解自衡单容过程的阶跃响应过程 3. 得出自衡单容过程的单位阶跃响应曲线 二、实验内容 已知两个单容过程的模型分别为 1 () 0.5 G s s =和5 1 () 51 s G s e s - = + ，试在 Simulink中建立模型，并求单位阶跃响应曲线。 三、实验步骤 1. 在Simulink中建立模型，得出实验原理图。 2. 运行模型后，双击Scope，得到的单位阶跃响应曲线。 四、实验结果 1．建立系统Simulink仿真模型图，其仿真模型为
通信工程系统仿真实验报告
通信原理课程设计 实验报告 专业：通信工程 届别：07 B班 学号：0715232022 姓名：吴林桂 指导老师：陈东华
数字通信系统设计 一、 实验要求： 信源书记先经过平方根升余弦基带成型滤波，成型滤波器参数自选，再经BPSK ，QPSK 或QAM 调制（调制方式任选），发射信号经AWGN 信道后解调匹配滤波后接收，信道编码可选（不做硬性要求），要求给出基带成型前后的时域波形和眼图，画出接收端匹配滤波后时域型号的波形，并在时间轴标出最佳采样点时刻。对传输系统进行误码率分析。 二、系统框图 三、实验原理： QAM 调制原理：在通信传渝领域中，为了使有限的带宽有更高的信息传输速率，负载更多的用户必须采用先进的调制技术，提高频谱利用率。QAM 就是一种频率利用率很高的调制技术。 t B t A t Y m m 00sin cos )(ωω+= 0≤t ≤Tb 式中 Tb 为码元宽度t 0cos ω为 同相信号或者I 信号； t 0s i n ω 为正交信号或者Q 信号； m m B A ,为分别为载波t 0cos ω,t 0sin ω的离散振幅； m 为 m A 和m B 的电平数，取值1 , 2 , . . . , M 。 m A = Dm*A ；m B = Em*A ； 式中A 是固定的振幅，与信号的平均功率有关，（dm ，em ）表示调制信号矢量点在信号空
间上的坐标，有输入数据决定。 m A 和m B 确定QAM 信号在信号空间的坐标点。称这种抑制载波的双边带调制方式为 正交幅度调制。 图3.3.2 正交调幅法原理图 Pav=（A*A/M ）*∑(dm*dm+em*em) m=(1,M) QAM 信号的解调可以采用相干解调,其原理图如图3.3.5所示。 图3.3.5 QAM 相干解调原理图 四、设计方案： (1)、生成一个随机二进制信号 (2)、二进制信号经过卷积编码后再产生格雷码映射的星座图 (3)、二进制转换成十进制后的信号 (4)、对该信号进行16-QAM 调制 (5)、通过升余弦脉冲成形滤波器滤波，同时产生传输信号 (6)、增加加性高斯白噪声，通过匹配滤波器对接受的信号滤波 (7)、对该信号进行16-QAM 解调 五、实验内容跟实验结果：
计算方法实验报告 拟合
南京信息工程大学实验（实习）报告 一、实验目的: 用最小二乘法将给定的十个点拟合成三次多项式。 二、实验步骤： 用matlab编制以函数为基的多项式最小二乘拟合程序，并用于对下列数据作三次多项式最小二乘拟合(取权函数wi=1) x -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 y -2.30 -1 -0.14 -0.25 0.61 1.03 1.75 2.75 4.42 6.94 给定直线方程为：y=1/4*x3+1/2*x2+x+1 三、实验结论： 最小二乘法：通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。 一般地。当测量数据的散布图无明显的规律时，习惯上取n次代数多项式。 程序运行结果为： a = 0.9731 1.1023 0.4862 0.2238 即拟合的三次方程为：y=0.9731+1.1023x+0.4862*x2+0.2238*x3
-2.5 -2-1.5-1-0.5 00.51 1.52 2.5 -4-20246 81012 x 轴 y 轴 拟合图 离散点 y=a(1)+a(2)*x+a(3)*x.2+a(4)*x.3 结论： 一般情况下，拟合函数使得所有的残差为零是不可能的。由图形可以看出最小二乘解决了残差的正负相互抵消的问题，使得拟合函数更加密合实验数据。 优点：曲线拟合是使拟合函数和一系列的离散点与观测值的偏差平方和达到最小。 缺点：由于计算方法简单，若要保证数据的精确度，需要大量的数据代入计算。
计算机仿真实验报告实验
《计算机仿真》上机实验报告 姓名： 学号： 2012104021 专业：测控 班级： 12级
实验一常微分方程的求解及系统数学模型的转换一．实验目的 通过实验熟悉计算机仿真中常用到的Matlab指令的使用方法，掌握常微分方程求解指令和模型表示及转换指令，为进一步从事有关仿真设计和研究工作打下基础。 二. 实验设备 个人计算机，Matlab软件。 三. 实验准备 预习本实验有关内容（如教材第2、3、5章中的相应指令说明和例题），编写本次仿真练习题的相应程序。 四. 实验内容 1. Matlab中常微分方程求解指令的使用 题目一：请用MATLAB的ODE45算法分别求解下列二个方程。要求:1.编写出Matlab 仿真程序；2.画出方程解的图形并对图形进行简要分析；3.分析下列二个方程的关系。 1．2． 1.function fun=funl(t,x) fun=-x^2;
[t,x]=ode45('fun1',[0,20],[1]); figure(1);plot(t,x); grid 2.function fun=fun2(t,x) fun=x^2; [t,x]=ode45('fun2',[0,20],[-1]); figure(2);plot(t,x); grid
题目二：下面方程组用在人口动力学中，可以表达为单一化的捕食者-被捕食者模式（例如，狐狸和兔子）。其中1x 表示被捕食者， 2x 表示捕食者。如果被捕食者有无限的食物，并且不会出现捕食者。于是有1'1x x ，则这个式子是以指数形式增长的。大量的被捕食者将会使捕食者的数量增长；同样，越来越少的捕食者会使被捕食者的数量增长。而且，人口数量也会增长。请分别调用ODE45、ODE23算法求解下面方程组。要求编写出Matlab 仿真程序、画出方程组解的图形并对图形进行分析和比较。 1.ODE45
OFDM系统仿真实验报告
无线通信——OFDM系统仿真
一、实验目的 1、了解OFDM 技术的实现原理 2、利用MATLAB 软件对OFDM 的传输性能进行仿真并对结论进行分析。 二、实验原理与方法 1 OFDM 调制基本原理 正交频分复用(OFDM)是多载波调制(MCM)技术的一种。MCM 的基本思想是把数据流串并变换为N 路速率较低的子数据流，用它们分别去调制N 路子载波后再并行传输。因子数据流的速率是原来的1/N ，即符号周期扩大为原来的N 倍，远大于信道的最大延迟扩展，这样MCM 就把一个宽带频率选择性信道划分成N 个窄带平坦衰落信道，从而“先天”具有很强的抗多径衰落和抗脉冲干扰的能力，特别适合于高速无线数据传输。OFDM 是一种子载波相互混叠的MCM ，因此它除了具有上述毗M 的优势外，还具有更高的频谱利用率。OFDM 选择时域相互正交的子载波，创门虽然在频域相互混叠，却仍能在接收端被分离出来。 2 OFDM 系统的实现模型 利用离散反傅里叶变换( IDFT) 或快速反傅里叶变换( IFFT) 实现的OFDM 系统如图1 所示。输入已经过调制(符号匹配) 的复信号经过串P 并变换后,进行IDFT 或IFFT 和并/串变换,然后插入保护间隔,再经过数/模变换后形成OFDM 调制后的信号s (t ) 。该信号经过信道后,接收到的信号r ( t ) 经过模P 数变换,去掉保护间隔以恢复子载波之间的正交性,再经过串/并变换和DFT 或FFT 后,恢复出OFDM 的调制信号,再经过并P 串变换后还原出输入的符号。 图1 OFDM 系统的实现框图 从OFDM 系统的基本结构可看出, 一对离散傅里叶变换是它的核心,它使各子载波相互正交。设OFDM 信号发射周期为[0,T],在这个周期内并行传输的N 个符号为001010(,...,)N C C C -，,其中ni C 为一般复数, 并对应调制星座图中的某一矢量。比如00(0)(0),(0)(0)C a j b a b =+?和分别为所要传输的并行信号, 若将
四旋翼飞行器仿真-实验报告
动态系统建模仿真实验报告（2） 四旋翼飞行器仿真 姓名： 学号： 指导教师： 院系： 2014.12.28
1实验容 基于Simulink建立四旋翼飞行器的悬停控制回路，实现飞行器的悬停控制； 建立GUI界面，能够输入参数并绘制运动轨迹； 基于VR Toolbox建立3D动画场景，能够模拟飞行器的运动轨迹。 2实验目的 通过在 Matlab 环境中对四旋翼飞行器进行系统建模，使掌握以下容： 四旋翼飞行器的建模和控制方法 在Matlab下快速建立虚拟可视化环境的方法。 3实验器材 硬件：PC机。 工具软件：操作系统：Windows系列；软件工具：MATLAB及simulink。 4实验原理 4.1四旋翼飞行器 四旋翼飞行器通过四个螺旋桨产生的升力实现飞行，原理与直升机类似。四个旋翼位于一个几何对称的十字支架前，后，左，右四端，如图 1 所示。旋翼由电机控制；整个飞行器依靠改变每个电机的转速来实现飞行姿态控制。 图1四旋翼飞行器旋转方向示意图
在图 1 中， 前端旋翼 1 和后端旋翼 3 逆时针旋转， 而左端旋翼 2 和右端的旋翼 4 顺时针旋转， 以平衡旋翼旋转所产生的反扭转矩。 由此可知， 悬停时， 四只旋翼的转速应该相等，以相互抵消反扭力矩；同时等量地增大或减小四只旋翼的转速，会引起上升或下降运动；增大某一只旋翼的转速，同时等量地减小同组另一只旋翼的转速，则产生俯仰、横滚运动；增大某一组旋翼的转速，同时等量减小另一组旋翼的转速，将产生偏航运动。 4.2建模分析 四旋翼飞行器受力分析,如图 2 所示 图2四旋翼飞行器受力分析示意图 旋翼机体所受外力和力矩为： 重力mg , 机体受到重力沿w z -方向； 四个旋翼旋转所产生的升力i F (i= 1 , 2 , 3 , 4)，旋翼升力沿b z 方向； 旋翼旋转会产生扭转力矩i M (i= 1 , 2 , 3 , 4)。i M 垂直于叶片的旋翼平面，与旋转矢量相反。 力模型为：2i F i F k ω= ，旋翼通过螺旋桨产生升力。F k 是电机转动力系数， 可取826.1110/N rpm -?，i ω为电机转速。旋翼旋转产生旋转力矩Mi(i=1,2,3,4)，
通信系统仿真实验报告(DOC)
通信系统实验报告——基于SystemView的仿真实验 班级： 学号： 姓名： 时间：
目录 实验一、模拟调制系统设计分析 -------------------------3 一、实验内容-------------------------------------------3 二、实验要求-------------------------------------------3 三、实验原理-------------------------------------------3 四、实验步骤与结果-------------------------------------4 五、实验心得------------------------------------------10 实验二、模拟信号的数字传输系统设计分析------------11 一、实验内容------------------------------------------11 二、实验要求------------------------------------------11 三、实验原理------------------------------------------11 四、实验步骤与结果------------------------------------12 五、实验心得------------------------------------------16 实验三、数字载波通信系统设计分析------------------17 一、实验内容------------------------------------------17 二、实验要求------------------------------------------17 三、实验原理------------------------------------------17 四、实验步骤与结果------------------------------------18 五、实验心得------------------------------------------27