第一章习题
3.请指出以下的变量名（函数名、M文件名）中，哪些是合法的？Abc 2004x lil-1 wu_2004 a&b qst.u _xyz 解：合法的变量名有：Abc wu_2004
4．指令窗操作
（1）求[12+2×（7-4）]÷32的运算结果
解：>> [12+2*(7-4)]/3^2
ans =
2
（2）输入矩阵A=[1，2，3；4，5，6；7，8，9]，观察输出。解：>> A=[1,2,3;4,5,6;7,8,9]
A =
1 2 3
4 5 6
7 8 9
（3）输入以下指令，观察运算结果；
clear;x=-8:0.5:8;
y=x';
X=ones(size(y))*x;
Y=y*ones(size(x));
R=sqrt(X.^2+Y.^2)+eps;
Z=sin(R)./R;
mesh(X,Y,Z);
colormap(hot)
xlabel('x'),ylabel('y'),zlabel('z')
解：
7．指令行编辑
（1）依次键入以下字符并运行：y1=2*sin(0.3*pi)/(1+sqrt(5))
解：>>y1=2*sin(0.3*pi)/(1+sqrt(5))
y1 =
0.5000
(2)通过反复按键盘的箭头键，实现指令回调和编辑，进行新的计算；y2=2*cos(0.3*pi)/(1+sqrt(5))
解：>>y2=2*cos(0.3*pi)/(1+sqrt(5))
y2 =
0.3633
11.编写题4中（3）的M脚本文件，并运行之。解：
第二章习题
1.在指令窗中键入x=1:0.2:2和y=2:0.2:1,观察所生成的数组。 解:>> x=1:0.2:2 x =
1.0000 1.2000 1.4000 1.6000 1.8000
2.0000 >> y=2:0.2:1 y =
Empty matrix: 1-by-0
2．要求在[0，2π]上产生50个等距采样数据的一维数组，试用两种不同的指令实现。
解： y1=0:2*pi/49:2*pi
y2=linspace(0,2*pi,50)
3.计算e -2t sint ，其中t 为[0，2π]上生成的10个等距采样的数组。 解：>> t=linspace(0,2*pi,10); x=exp(-2*t).*sin(t) x =
0 0.1591 0.0603 0.0131 0.0013 -0.0003 -0.0002 -0.0001 -0.0000 -0.0000
4.已知A=??????4321 , B=???
???8765,计算矩阵A 、B 乘积和点乘. 解：>> A=[1,2;3,4]; B=[5,6;7,8]; x=A*B x =
19 22 43 50 >> x=A.*B x =
5 12
21 32
5.已知A=??????05314320，B=???
???05314320，计算A&B, A|B, ~A, A==B, A>B. 解：>> A=[0,2,3,4;1,3,5,0]; B=[1,0,5,3;1,5,0,5]; a1=A&B
a2=A|B a3=~A a4=(A==B) a5=(A>B) a1 =
0 0 1 1 1 1 0 0 a2 =
1 1 1 1 1 1 1 1 a3 =
1 0 0 0 0 0 0 1 a4 =
0 0 0 0 1 0 0 0 a5 =
0 1 0 1 0 0 1 0
7.将题5中的A 阵用串转换函数转换为串B ，再size 指令查看A 、B 的结构，有何不同？
解：>> A=[0,2,3,4;1,3,5,0] B=num2str(A)
size(A)
size(B)
A =
0 2 3 4
1 3 5 0
B =
0 2 3 4
1 3 5 0
ans =
2 4
ans =
2 10
第三章习题
1.已知系统的响应函数为
)sin(1
1)(θββ
ε+-
=-t e t y t ,其中
???
?
?
?-=-=εεθεβ2
2
1arctan ,1 ,要求用不同线型或颜色,在同一图上绘制ε取值分别为0.2、0.4、0.6、0.8时,系统在t ∈[0,18] 区间的响应曲线,并要求用ε=0.2和 ε=0.8对他们相应的两条曲线进行文字标志。
解： clc
close all clear all t=0:0.02:18;
xi=[0.2,0.4,0.6,0.8]'; sxi=sqrt(1-xi.^2); sita=atan(sxi./xi);
y=1-exp(-xi*t).*sin(sxi*t+sita*ones(1,901))./(sxi*ones(1,901))
plot(t,y(1), 'r-', t,y(2), ' b*', t,y(3), ' g+', t,y(4), ' k.')
text(4.2,1.4,'\xi =0.2') text(3.8,0.9,'\xi=0.8')
2.用plot3、mesh 、surf 指令绘制
()
()
2
2
2
2
111
y
x y x z +++
+-=
三维图(x,y 围自定)。 解：
clc;close all ;clear all ; x=-5:0.1:5;y=-5:0.1:5; [X,Y]=meshgrid(x,y); a=sqrt((1-X).^2+Y.^2); b=sqrt((1+X).^2+Y.^2); Z=1./(a+b);
a1=sqrt((1-x).^2+y.^2); b1=sqrt((1+x).^2+y.^2); z=1./(a1+b1);
subplot(1,3,1),plot3(x,y,z),xlabel('x'),ylabel('y'),zlabel('z');box on ; subplot(1,3,2),surf(X,Y,Z),xlabel('x'),ylabel('y'),zlabel('z');box on ; subplot(1,3,3),mesh(X,Y,Z),xlabel('x'),ylabel('y'),zlabel('z');box on ;
3.对向量t进行以下运算可以构成三个坐标的值向量:x=sin(t),y=cos(t),z=t.利用指令plot3,并选用绿色的实线绘制相应的三维曲线.
解：
t=(0:0.01:2)*pi;
x=sin(t);
y=cos(t);
z=t;
plot3(x,y,z,'b-');box on
第四章习题
1.请分别用for 和while 循环语句计算K=∑=63
02i i 的程序，再写出一种避免循环的
计算程序。(提示：可考虑利用MATLAB 的sum(X,n)函数，实现沿数组X 的第n 维求和。) 解： 1）K=0; for i=0:63; K=K+2^i; end K
K =1.8447e+019 2）i=0;K=0; while i<=63; K=K+2^i; i=i+1; end; K
K =1.8447e+019 3）i=0; X=0:63; for i=0:63;
X(i+1)=2^i; end
sum(X,2)
ans =1.8447e+019
第五章习题
1.将下列系统的传递函数模型用MATLAB 语言表达出来。
(1))170046842541254289()
1700109329135()(23452341+++++++++=s s s s s s s s s s G
解：
num=[1,35,291,1093,1700]; den=[1,289,254,2541,4684,1700]; sys=tf(num,den)
(2))
15).(5).(1()
3(15)(2++++=
s s s s s G 解： z=-3; p=[-1,-5,-15]; k=15; sys=zpk(z,p,k)
(3))252).(1).(1()
23.()2.(.100)(23223+++-++++=s s s s s s s s s s G
解： z=[0,-2,-2]; p=[-1,1]; k=100; sys1=zpk(z,p,k);
num=[1,3,2];
den=[1,2,5,2];
sys2=tf(num,den);
sys=series(sys1,sys2)
4.求题3中的系统模型的等效传递函数模型和零极点模型。解：
A=[3,2,1;0,4,6;0,-3,-5];
B=[1,2,3]' ;
C=[1,2,5];
D=0;
sys=ss(A,B,C,D);
systf=tf(sys)
syszpk=zpk(sys)
Transfer function:
20 s^2 - 83 s + 138
---------------------
s^3 - 2 s^2 - 5 s + 6
Zero/pole/gain:
20 (s^2 - 4.15s + 6.9)
-----------------------
(s-3) (s-1) (s+2)
5.已知系统的动力学方程如下，试用MATLAB 语言写出它们的传递函数。 (1))(2)()(500)(50)(15)(.
.....)3(t r t r t y t y t y t y +=+++ 解： num=[1,2,0]; den=[1,15,50,500]; sys=tf(num,den) Transfer function: s^2 + 2 s
------------------------- s^3 + 15 s^2 + 50 s + 500
(2) )(4)(4)(6)(3)(.
..
t r dt t y t y t y t y =+++? 解： num=[4,0]; den=[1,3,6,4]; sys=tf(num,den) Transfer function: 4 s
--------------------- s^3 + 3 s^2 + 6 s + 4
6.试用MATLAB 语言表示图5-13所示系统。当分别以y =x 2和f 为系统输出、输
入时的传递函数模型和状态空间模型(图中k=7N/m,c1=0.5N/m.s-1, c2=0.2N/m.s-1,m1=3.5kg, m2=5.6kg)。
f
解：)(t
k=7;
c1=0.5;
c2=0.2;
m1=3.5;
m2=5.6;
num=[m1,c1,k];
den=[m1*m2,c1*m1+c2*m1+c1*m2,c1*c2+m2*k,c1*k+c2*k,0];
sys=tf(num,den)
Transfer function:
3.5 s^2 + 0.5 s + 7
--------------------------------------
19.6 s^4 + 5.25 s^3 + 39.3 s^2 + 4.9 s
7.试用MATLAB语言分别表示图5-14所示系统质量m1,m2的位移x1,x2对输入f
的传递函数X2(s)/F(s)和X1(s)/F(s),其中m1=12kg, m2=38kg,k=1000N/m, c=0.1N/m.s-1。
解：
m1=12;
m2=38;
k=1000;
c=0.1;
num=[c,k];
den=[m1*m2,m1*c+m2*c,m1*k+m2*k,0,0];
sys1=tf(num,den)
num=[m1,c,k];
den=[m1*m2,m1*c+m2*c,m1*k+m2*k,0,0];
sys2=tf(num,den)
Transfer function:
0.1 s + 1000
---------------------------
456 s^4 + 5 s^3 + 50000 s^2
Transfer function:
12 s^2 + 0.1 s + 1000
---------------------------
456 s^4 + 5 s^3 + 50000 s^2
补充题
求图示传递函数
sys1=tf([1,2],[1,3,4]);
sys2=tf([1,4,5] ,[1,6,7,8]);
sys3=tf([1,0],[1,2]);
sys4=tf([1],[1,3]);
sys5=parallel(sys3,sys4);
sys=feedback(sys1*sys2*sys5,1,-1)
结果
s^5 + 10 s^4 + 39 s^3 + 74 s^2 + 66 s + 20
----------------------------------------------------------------- s^7 + 14 s^6 + 81 s^5 + 262 s^4 + 530 s^3 + 684 s^2 + 538 s + 212
第六章习题
2.将例6-2中的微分方程改写为以下形式：
1
)0(,0)0(0.)1.(.
.
2
..
===+--y y y y y y μ
求μ分别为1、2时，在时间区间t=[0,20]微分方程的解。 解： M 函数文件
function dx=wffc(t,x,flag,ps) dx=zeros(2,1); dx(1)=x(2);
dx(2)=ps*(1-x(1)^2)*x(2)-x(1); 调用程序
clc;close all ;clear all ; tspan=[0,20]; x0=[0,1]; ps=1;
[T1,X1]=ode45('wffc',tspan,x0,odeset,ps); ps=2;
[T2,X2]=ode45('wffc',tspan,x0,odeset,ps); plot(T1,X1(:,1),'r',T2,X2(:,1),'b-.') X1(:,1)
X2(:,1)
3.对图6-18所示反馈系统进行单位阶跃响应和方波响应(方波周期为30s)仿真。要求：
(1)利用MATLAB模型连接函数求出系统闭环传递函数。
(2)利用step函数求单位阶跃响应。
(3)利用gensig函数产生方波信号，利用lsim函数求方波响应。
解：
clc;close all;clear all;
% (1)
sys1=tf([1,0.5],[1,0.1]);
sys2=ZPK([],[0,-2,-10],20);
sys3=series(sys1,sys2);
sys4=feedback(sys3,1,-1);
% (2)
第1章建模与仿真的基本概念 参照P8例子，列举一个你相对熟悉的简单实际系统为例，采用非形式描述出来。 第2章建模方法论 1、什么是数学建模形式化的表示？试列举一例说明形式化表示与非形式化表示的区别。 模型的非形式描述是说明实际系统的本质，但不是详尽描述。是对模型进行深入研究的基础。主要由模型的实体、包括参变量的描述变量、实体间的相互关系及有必要阐述的假设组成。模型的非形式描述主要说明实体、描述变量、实体间的相互关系及假设等。 例子：环形罗宾服务模型的非形式描述： 实体 CPU，USR1，…，USR5 描述变量 CPU:Who,Now(现在是谁)----范围{1,2,…,5}; Who.Now=i表示USRi由CPU服务。 USR：Completion.State（完成情况）----范围[0，1]；它表示USR完成整个程序任务的比例。参变量 X-----范围[0，1]；它表示USRi每次完成程序的比率。 i 实体相互关系 （1）CPU 以固定速度依次为用户服务，即Who.Now为1，2，3，4，5，1，2…..循环运行。 X工作。假设：CPU对USR的服务时间固定，不（2）当Who.Now=I,CPU完成USRi余下的 i X决定。 依赖于USR的程序；USRi的进程是由各自的参变量 i 2、何谓“黑盒”“白盒”“灰盒”系统？ “黑盒”系统是指系统内部结构和特性不清楚的系统。对于“黑盒”系统，如果允许直接进行实验测量并通过实验对假设模型加以验证和修正。对属于黑盒但又不允许直接实验观测的系统，则采用数据收集和统计归纳的方法来假设模型。 对于内部结构和特性清楚的系统，即白盒系统，可以利用已知的一些基本定律，经过分析和演绎导出系统模型。 3、模型有效性和模型可信性相同吗？有何不同？ 模型的有效性可用实际系统数据和模型产生的数据之间的符合程度来度量。它分三个不同级别的模型有效：复制有效、预测有效和结构有效。不同级别的模型有效，存在不同的行为水平、状态结构水平和分解结构水平的系统描述。 模型的可信度指模型的真实程度。一个模型的可信度可分为： 在行为水平上的可信性，即模型是否重现真实系统的行为。 在状态结构水平上可信性，即模型能否与真实系统在状态上互相对应，通过这样的模型可以对未来的行为进行唯一的预测。 在分解结构水平上的可信性，即模型能否表示出真实系统内部的工作情况，而且是惟一表示出来。 不论对于哪一个可信性水平，可信性的考虑贯穿在整个建模阶段及以后各阶段，必须考虑以下几个方面： 1在演绎中的可信性。2在归纳中的可信性。3在目的方面的可信性。 4、基于计算机建模方法论与一般建模方法论有何不同？（P32） 经典的建模与仿真的主要研究思路，首先界定研究对象-实际系统的边界和建模目标，利用已有的数学建模工具和成果，建立相应的数学模型，并用计算装置进行仿真。这种经典的建
MATLAB 仿真实验报告 课题名称：MATLAB 仿真——图像处理 学院：机电与信息工程学院 专业：电子信息科学与技术 年级班级：2012级电子二班 一、实验目的 1、掌握MATLAB处理图像的相关操作，熟悉相关的函数以及基本的MATLAB语句。 2、掌握对多维图像处理的相关技能，理解多维图像的相关性质 3、熟悉Help 命令的使用，掌握对相关函数的查找，了解Demos下的MATLAB自带的原函数文件。 4、熟练掌握部分绘图函数的应用，能够处理多维图像。 二、实验条件
MATLAB调试环境以及相关图像处理的基本MATLAB语句，会使用Help命令进行相关函数查找 三、实验内容 1、nddemo.m函数文件的相关介绍 Manipulating Multidimensional Arrays MATLAB supports arrays with more than two dimensions. Multidimensional arrays can be numeric, character, cell, or structure arrays. Multidimensional arrays can be used to represent multivariate data. MATLAB provides a number of functions that directly support multidimensional arrays. Contents ： ●Creating multi-dimensional arrays 创建多维数组 ●Finding the dimensions寻找尺寸 ●Accessing elements 访问元素 ●Manipulating multi-dimensional arrays操纵多维数组 ●Selecting 2D matrices from multi-dimensional arrays从多维数组中选择二维矩 阵 (1)、Creating multi-dimensional arrays Multidimensional arrays in MATLAB are created the same way as two-dimensional arrays. For example, first define the 3 by 3 matrix, and then add a third dimension. The CAT function is a useful tool for building multidimensional arrays. B = cat(DIM,A1,A2,...) builds a multidimensional array by concatenating（联系起来）A1, A2 ... along the dimension DIM. Calls to CAT can be nested（嵌套）. (2)、Finding the dimensions SIZE and NDIMS return the size and number of dimensions of matrices. (3)、Accessing elements To access a single element of a multidimensional array, use integer subscripts（整数下标）. (4)、Manipulating multi-dimensional arrays
M/M/1排队系统实验报告 一、实验目的 本次实验要求实现M/M/1单窗口无限排队系统的系统仿真，利用事件调度法实现离散事件系统仿真，并统计平均队列长度以及平均等待时间等值，以与理论分析结果进行对比。 二、实验原理 根据排队论的知识我们知道，排队系统的分类是根据该系统中的顾客到达模式、服务模式、服务员数量以及服务规则等因素决定的。 1、 顾客到达模式 设到达过程是一个参数为λ的Poisson 过程，则长度为t 的时间内到达k 个呼叫的概 率 服从Poisson 分布，即e t k k k t t p λλ-=!)()(，?????????=,2,1,0k ，其中λ>0为一常数，表示了 平均到达率或Poisson 呼叫流的强度。 2、 服务模式 设每个呼叫的持续时间为i τ，服从参数为μ的负指数分布，即其分布函数为{}1,0t P X t e t μ-<=-≥ 3、 服务规则 先进先服务的规则（FIFO ） 4、 理论分析结果 在该M/M/1系统中，设 λρμ=，则稳态时的平均等待队长为1Q ρλρ=-，顾客的平均等待时间为T ρ μλ=-。 三、实验内容 M/M/1排队系统：实现了当顾客到达分布服从负指数分布，系统服务时间也服从负指数分布，单服务台系统，单队排队，按FIFO （先入先出队列）方式服务。 四、采用的语言 MatLab 语言 源代码： clear; clc;
%M/M/1排队系统仿真 SimTotal=input('请输入仿真顾客总数SimTotal='); %仿真顾客总数；Lambda=0.4; %到达率Lambda； Mu=0.9; %服务率Mu； t_Arrive=zeros(1,SimTotal); t_Leave=zeros(1,SimTotal); ArriveNum=zeros(1,SimTotal); LeaveNum=zeros(1,SimTotal); Interval_Arrive=-log(rand(1,SimTotal))/Lambda;%到达时间间隔Interval_Serve=-log(rand(1,SimTotal))/Mu;%服务时间 t_Arrive(1)=Interval_Arrive(1);%顾客到达时间 ArriveNum(1)=1; for i=2:SimTotal t_Arrive(i)=t_Arrive(i-1)+Interval_Arrive(i); ArriveNum(i)=i; end t_Leave(1)=t_Arrive(1)+Interval_Serve(1);%顾客离开时间LeaveNum(1)=1; for i=2:SimTotal if t_Leave(i-1) 基于Matlab/Simulink 的BOOST电路仿真 姓名： 学号： 班级： 时间：2010年12月7日 1引言 BOOST 电路又称为升压型电路, 是一种直流- 直流变换电路, 其电路结构如图1 所示。此电路在开关电源领域内占有非常重要的地位, 长期以来广泛的应用于各种电源设备的设计中。对它工作过程的理解掌握关系到对整个开关电源领域各种电路工作过程的理解, 然而现有的书本上仅仅给出电路在理想情况下稳态工作过程的分析, 而没有提及电路从启动到稳定之间暂态的工作过程, 不利于读者理解电路的整个工作过程和升压原理。采用matlab仿真分析方法, 可直观、详细的描述BOOST 电路由启动到达稳态的工作过程, 并对其中各种现象进行细致深入的分析, 便于我们真正掌握BOO ST 电路的工作特性。 图1BOO ST 电路的结构 2电路的工作状态 BOO ST 电路的工作模式分为电感电流连续工作模式和电感电流断续工作模式。其中电流连续模式的电路工作状态如图2 (a) 和图2 (b) 所示, 电流断续模式的电路工作状态如图2 (a)、(b)、(c) 所示, 两种工作模式的前两个工作状态相同, 电流断续型模式比电流连续型模式多出一个电感电流为零的工作状态。 (a) 开关状态1 (S 闭合) (b) 开关状态2 (S 关断) (c) 开关状态3 (电感电流为零) 图2BOO ST 电路的工作状态 3matlab仿真分析 matlab 是一种功能强大的仿真软件, 它可以进行各种各样的模拟电路和数字电路仿真，并给出波形输出和数据输出, 无论对哪种器件和哪种电路进行仿真, 均可以得到精确的仿真结果。本文应用基于matlab软件对BOO ST 电路仿真, 仿真图如图3 所示,其中IGBT作为开关, 以脉冲发生器脉冲周期T=0.2ms,脉冲宽度为50%的通断来仿真图2 中开关S的通断过程。 图3BOO ST 电路的PSp ice 模型 3.1电路工作原理 在电路中IGBT导通时，电流由E经升压电感L和V形成回路，电感L储能；当IGBT关断时，电感产生的反电动势和直流电源电压方向相同互相叠加，从而在负载侧得到高于电源的电压，二极管的作用是阻断IGBT导通是，电容的放电回路。调节开关器件V的通断周期，可以调整负载侧输出电流和电压的大小。负载侧输出电压的平均值为: （3-1） 式（3-1）中T为开关周期, 为导通时间，为关断时间。 实验一 MATLAB 及仿真实验（控制系统的时域分析） 一、实验目的 学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点 1、 系统的典型响应有哪些 2、 如何判断系统稳定性 3、 系统的动态性能指标有哪些 三、实验方法 （一） 四种典型响应 1、 阶跃响应： 阶跃响应常用格式： 1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。 2、),(Tn sys step ；表示时间范围0---Tn 。 3、),(T sys step ；表示时间范围向量T 指定。 4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。 2、 脉冲响应： 脉冲函数在数学上的精确定义：0 ,0)(1)(0 ?==?∞ t x f dx x f 其拉氏变换为：) ()()()(1)(s G s f s G s Y s f === 所以脉冲响应即为传函的反拉氏变换。 脉冲响应函数常用格式： ① )(sys impulse ； ② ); ,();,(T sys impulse Tn sys impulse ③ ),(T sys impulse Y = （二） 分析系统稳定性 有以下三种方法： 1、 利用pzmap 绘制连续系统的零极点图； 2、 利用tf2zp 求出系统零极点； 3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析 Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim. Wireless Network Experiment Three: Queuing Theory ABSTRACT This experiment is designed to learn the fundamentals of the queuing theory. Mainly about the M/M/S and M/M/n/n queuing MODELS. KEY WORDS: queuing theory, M/M/s, M/M/n/n, Erlang B, Erlang C. INTRODUCTION A queue is a waiting line and queueing theory is the mathematical theory ofwaiting lines.More generally, queueing theory is concerned with the mathematical modeling and analysisof systems that provide service to random demands. Incommunication networks, queues are encountered everywhere. For example, theincoming data packets are randomly arrived and buffered, waiting for the routerto deliver. Such situation is considered as a queue. A queueing model is an abstract description of such a system. Typically, a queueing model represents (1) thesystem's physical configuration, by specifying the number and arrangement of theservers, and (2) the stochastic nature of the demands, by specifying the variabilityin the arrival process and in the service process. The essence of queueing theory is that it takes into account the randomness ofthe arrival process and the randomness of the service process. The most commonassumption about the arrival process is that the customer arrivals follow a Poisson process, where the times between arrivals are exponentially distributed. Theprobability of the exponential distribution function matlab 单服务台排队系统实验报告 一、实验目的 本次实验要求实现M/M/1单窗口无限排队系统的系统仿真，利用事件调度法实现离散事件系统仿真，并统计平均队列长度以及平均等待时间等值，以与理论分析结果进行对比。 二、实验原理 根据排队论的知识我们知道，排队系统的分类是根据该系统中的顾客到达模式、服务模式、服务员数量以及服务规则等因素决定的。 1、 顾客到达模式 设到达过程是一个参数为λ的Poisson 过程，则长度为t 的时间内到达k 个呼 叫的概率 服从Poisson 分布，即 e t k k k t t p λλ-= !)()(，?????????=,2,1,0k ，其中λ>0为一 常数，表示了平均到达率或Poisson 呼叫流的强度。 2、 服务模式 设每个呼叫的持续时间为i τ，服从参数为μ的负指数分布，即其分布函数为 {}1,0t P X t e t μ-<=-≥ 3、 服务规则 先进先服务的规则（FIFO ） 4、 理论分析结果 在该M/M/1系统中，设λρμ= ，则稳态时的平均等待队长为1Q ρλ ρ= -，顾客 的平均等待时间为 T ρμλ= -。 三、实验内容 M/M/1排队系统：实现了当顾客到达分布服从负指数分布，系统服务时间也服 从负指数分布，单服务台系统，单队排队，按FIFO 方式服务。 四、采用的语言 MatLab 语言 源代码： clear; clc; %M/M/1排队系统仿真 SimTotal=input('请输入仿真顾客总数SimTotal='); %仿真顾客总数；Lambda=0.4; %到达率Lambda； Mu=0.9; %服务率Mu； t_Arrive=zeros(1,SimTotal); t_Leave=zeros(1,SimTotal); ArriveNum=zeros(1,SimTotal); LeaveNum=zeros(1,SimTotal); Interval_Arrive=-log(rand(1,SimTotal))/Lambda;%到达时间间隔Interval_Serve=-log(rand(1,SimTotal))/Mu;%服务时间 t_Arrive(1)=Interval_Arrive(1);%顾客到达时间 ArriveNum(1)=1; for i=2:SimTotal t_Arrive(i)=t_Arrive(i-1)+Interval_Arrive(i); ArriveNum(i)=i; end t_Leave(1)=t_Arrive(1)+Interval_Serve(1);%顾客离开时间LeaveNum(1)=1; for i=2:SimTotal if t_Leave(i-1) 基于matlab的电路仿真 杨泽辉51130215 %基于matlab的电路仿真 %关键词: RC电路仿真, matlab, GUI设计 % 基于matlab的电路仿真 %功能：产生根据输入波形与电路的选择产生输出波形 close all;clear;clc; %清空 figure('position',[189 89 714 485]); %创建图形窗口，坐标（189，89），宽714，高485；Na=['输入波形[请选择]|输入波形:正弦波|',... '输入波形:方形波|输入波形:脉冲波'];%波形选择名称数组； Ns={'sin','square','pulse'}; %波形选择名称数组； R=2; % default parameters: resistance 电阻值 C=2; % default parameters: capacitance电容值 f=10; % default parameters: frequency 波形频率 TAU=R*C; tff=10; % length of time ts=1/f; % sampling length sys1=tf([1],[1,1]); % systems for integral circuit %传递函数； sys2=tf([1,0],[1,1]); % systems for differential circuit a1=axes('position',[0.1,0.6,0.3,0.3]); %创建坐标轴并获得句柄； po1=uicontrol(gcf,'style','popupmenu',... %在第一个界面的上方创建一个下拉菜单'unit','normalized','position',[0.15,0.9,0.2,0.08],... %位置 'string',Na,'fontsize',12,'callback',[]); %弹出菜单上的字符为数组Na，字体大小为12， set(po1,'callback',['KK=get(po1,''Value'');if KK>1;',... 'st=char(Ns(KK-1));[U,T]=gensig(st,R*C,tff,1/f);',... 'axes(a1);plot(T,U);ylim([min(U)-0.5,max(U)+0.5]);',... 'end;']); %pol触发事件：KK获取激发位置，st为当前触发位置的字符串，即所选择的波形类型； %[U,T]，gensing,产生信号，类型为st的值，周期为R*C，持续时间为tff， %采样周期为1/f，U为所产生的信号，T为时间； %创建坐标轴al；以T为x轴，U为y轴画波形，y轴范围。。。 Ma=['电路类型[请选择]|电路类型:积分型|电路类型:微分型']; %窗口2电路类型的选择数组； a2=axes('position',[0.5,0.6,0.3,0.3]);box on; %创建坐标轴2； set(gca,'xtick',[]);set(gca,'ytick',[]); %去掉坐标轴的刻度 po2=uicontrol(gcf,'style','popupmenu',... %在第二个窗口的位置创建一个下拉菜单，同1 'unit','normalized','position',[0.55,0.9,0.2,0.08],... 'string',Ma,'fontsize',12,'callback',[]); set(po2,'callback',['KQ=get(po2,''Value'');axes(a2);',... %po2属性设置，KQ为选择的电路类型，'if KQ==1;cla;elseif KQ==2;',... %1则清除坐标轴，2画积分电路，3画微分电路 'plot(0.14+0.8i+0.02*exp(i*[0:.02:8]),''k'');hold on;',... 'plot(0.14+0.2i+0.02*exp(i*[0:.02:8]),''k'');',... 'plot(0.84+0.2i+0.02*exp(i*[0:.02:8]),''k'');',... 'plot(0.84+0.8i+0.02*exp(i*[0:.02:8]),''k'');',... 'plot([0.16,0.82],[0.2,0.2],''k'');',... 'plot([0.16,0.3],[0.8,0.8],''k'');',... 'plot([3,4,4,3,3]/10,[76,76,84,84,76]/100,''k'');',... 'plot([0.4,0.82],[0.8,0.8],''k'');',... 'plot([0.6,0.6],[0.8,0.53],''k'');',... 'plot([0.6,0.6],[0.2,0.48],''k'');',... 'plot([0.55,0.65],[0.53,0.53],''k'');',... 'plot([0.55,0.65],[0.48,0.48],''k'');',... 'text(0.33,0.7,''R'');',... 增量调制MATLAB仿真实验 增量调制(DM)实验 一、实验目的 （1）进一步掌握MATLAB的应用。 （2）进一步掌握计算机仿真方法。 （3）学会用MATLAB软件进行增量调制(DM)仿真实验。 二、实验原理 增量调制是由PCM发展而来的模拟信号数字化的一种编码方式，它是PCM的一种特例。增量调制编码基本原理是指用一位编码，这一位码不是表示信号抽样值的大小，而是表示抽样幅度的增量特性，即采用一位二进制数码“1”或“0”来表示信号在抽样时刻的值相对于前一个抽样时刻的值是增大还是减小，增大则输出“1”码，减小则输出“0”码。输出的“1”，“0”只是表示信号相对于前一个时刻的增减，不表示信号的绝对值。 增量调制最主要的特点就是它所产生的二进制代码表示模拟信号前后两个抽样值的差别(增加、还是减少)而不是代表抽样值本身的大小，因此把它称为增量调制。在增量调制系统的发端调制后的二进制代码1和0只表示信号这一个抽样时刻相对于前一个抽样时刻是增加(用1码)还是减少(用0码)。收端译码器每收到一个1码，译码器的输出相对于前一个时刻的值上升一个量化阶，而收到一个0码，译码器的输出相对于前一个时刻的值下降一个量化阶。 增量调制(DM)是DPCM的一种简化形式。在增量调制方式下,采用1比特量化器,即用1位二进制码传输样值的增量信息,预测器是 一个单位延迟器,延迟一个采样时间间隔。预测滤波器的分子系数向量是[0,1],分母系数为1。当前样值与预测器输出的前一样值相比较,如果其差值大于零,则发1码,如果小于零则发0码。 三、实验内容 增量调制系统框图如图一所示,其中量化器是一个零值比较器,根据输入的电平极性,输出为 δ,预测器是一个单位延迟器,其输出为前一个采样时刻的解码样值,编码器也是一个零值比较器,若其输入为负值,则编码输出为0,否则输出为1。解码器将输入1,0符号转换为 δ,然后与预测值相加后得出解码样值输出,同时也作为预测器的输入 输入样值 e n e n =δsgn(e n ) 传输 n ) n n-1+δsgn(e n ) x n + - + + 预测输出 + n-1 + 预测输出 解码样值输出 x n-1 预测输入x n =x n-1+δsgn(e n ) 图一 增量调制原理框图 设输入信号为: x(t)=sin2π50t+0.5sin 2π150t 增量调制的采样间隔为1ms,量化阶距δ=0.4,单位延迟器初始值为0。建立仿真模型并求出前20个采样点使客商的编码输出序列以 解码 编码 二电平量化 单位延迟 单位 延迟 《电力系统设计》报告 题目: 基于MATLAB的电力系统仿学院：电子信息与电气工程学院 班级: 13级电气 1 班 姓名:田震 学号: 20131090124 日期：2015年12月6日 基于MATLAB的电力系统仿真 摘要：目前，随着科学技术的发展和电能需求量的日益增长，电力系统规模越来越庞大，超高压远距离输电、大容量发电机组、各种新型控制装置得到了广泛的应用，这对于合理利用能源，充分挖掘现有的输电潜力和保护环境都有重要意义。另一方面，随着国民经济的高速发展，以城市为中心的区域性用电增长越来越快，大电网负荷中心的用电容量越来越大，长距离重负荷输电的情况日益普遍，电力系统在人们的生活和工作中担任重要角色，电力系统的稳定运行直接影响着人们的日常生活。从技术和安全上考虑直接进行电力试验可能性很小，因此迫切要求运用电力仿真来解决这些问题。 电力系统仿真是将电力系统的模型化、数学化来模拟实际的电力系统的运行，可以帮助人们通过计算机手段分析实际电力系统的各种运行情况，从而有效的了解电力系统概况。本文根据电力系统的特点，利用MATLAB的动态仿真软件Simulink搭建了无穷大电源的系统仿真模型，得到了在该系统主供电线路电源端发生三相短路接地故障并由故障器自动跳闸隔离故障的仿真结果，并分析了这一暂态过程。通过仿真结果说明MATLAB电力系统工具箱是分析电力系统的有效工具。 关键词：电力系统；三相短路；故障分析；MATLAB仿真 目录 一．前言 (4) 二．无穷大功率电源供电系统仿真模型构建 (5) 1.总电路图的设计 (5) 2.各个元件的参数设定 (6) 2.1供电模块的参数设定 (6) 2.2变压器模块的参数设置 (6) 2.3输电线路模块的参数设置 (7) 2.4三相电压电流测量模块 (8) 2.5三相线路故障模块参数设置 (8) 2.6三相并联RLC负荷模块参数设置 (9) 3.仿真结果 (9) 高频电子线路Matlab 仿真实验要求 1. 仿真题目 （1） 线性频谱搬移电路仿真 根据线性频谱搬移原理，仿真普通调幅波。 基本要求：载波频率为8kHz ，调制信号频率为400Hz ，调幅度为0.3；画出调制信号、载波信号、已调信号波形，以及对应的频谱图。 扩展要求1：根据你的学号更改相应参数和代码完成仿真上述仿真；载波频率改为学号的后5位，调制信号改为学号后3位，调幅度设为最后1位/10。（学号中为0的全部替换为1，例如学号2010101014，则载波为11114Hz ，调制信号频率为114，调幅度为0.4）。 扩展要求2：根据扩展要求1的条件，仿真设计相应滤波器，并获取DSB-SC 和SSB 的信号和频谱。 （2） 调频信号仿真 根据调频原理，仿真调频波。 基本要求：载波频率为30KHz ，调制信号为1KHz ，调频灵敏度32310f k π=??，仿真调制信号，瞬时角频率，瞬时相位偏移的波形。 扩展要求：调制信号改为1KHz 的方波，其它条件不变，完成上述仿真。 2. 说明 （1） 仿真的基本要求每位同学都要完成，并且记入实验基本成绩。 （2） 扩展要求可以选择完成。 1.0 >> ma = 0.3; >> omega_c = 2 * pi * 8000; >> omega = 2 * pi * 400; >> t = 0 : 5 / 400 / 1000 : 5 / 400; >> u_cm = 1; >> fc = cos(omega_c * t); >> fa = cos(omega * t); >> u_am = u_cm * (1 + fa).* fc; >> U_c =fft(fc,1024); >> U_o =fft(fa,1024); >> U_am =fft(u_am, 1024); >> figure(1); >> subplot(321);plot(t, fa, 'k');title('调制信号');grid;axis([0 2/400 -1.5 1.5]); >> subplot(323);plot(t, fc, 'k');title('高频载波');grid;axis([0 2/400 -1.5 1.5]); >> subplot(325);plot(t, u_am, 'k');title('已调信号');grid;axis([0 2/400 -3 3]); >> fs = 5000; >> w1 = (0:511)/512*(fs/2)/1000; >> subplot(322);plot(w1, abs([U_am(1:512)']),'k');title('调制信号频谱');grid;axis([0 0.7 0 500]); >> subplot(324);plot(w1, abs([U_c(1:512)']),'k');title('高频载波频谱');grid;axis([0 0.7 0 500]); >> subplot(326);plot(w1, abs([U_am(1:512)']),'k');title('已调信号频谱');grid;axis([0 0.7 0 500]); 1.1 >> ma = 0.8; >> omega_c = 2 * pi * 11138; >> omega = 2 * pi * 138; >> t = 0 : 5 / 400 / 1000 : 5 / 400; >> u_cm = 1; >> fc = cos(omega_c * t); 《现代机械工程基础实验》之机械工程控制基础综合实验报告 姓名 学号 班级 山东建筑大学机电工程学院 2012.06.04~06 第一部分 Matlab 编程方法及仿真实验 实验1. 三维曲面的绘制（略） 实验2. 系统零极点绘制例：求部分分式展开式和)(t g 一个线性定常系统的传递函数是 1 5422 3)(2 3 ++++= s s s s s G （1） 使用MATLAB 建立传递函数，并确定它的极点和零点，写出)(s G 的部分分式展开式并绘制 系统的脉冲响应。 实验结果：零点-0.6667 极点-0.8796 + 1.1414i -0.8796 - 1.1414i -0.2408 实验3. 系统的阶跃响应 例. )(s G 的阶跃响应 对例2中由（1）式给出的传递函数)(s G ，增加一个0=s 处的极点，使用impulse 命令绘制其拉普拉斯反变换式曲线，得到阶跃响应图。将该响应与对)(s G 使用step 命令所得到的响应比较，确定系统的DC 增益。利用初值定理和终值定理来校验结果。 实验结果：DC 增益= 2 实验4. 双输入反馈系统单位阶跃响应 考虑一个如图1所示的反馈系统，它既有参考输入也有干扰输入，其中对象和传感器的传递函数是 )12)(15.0(4)(++=s s s G p ，105.01 )(+=s s H 控制器是一个增益为80，有一个在3-=s 处的零点，极点/零点比15=α超前控制器。推导 两个独立的MATLAB 模型，其中一个模型的输入为)(s R ，另一个输入为)(s D 。使用这些模型确定闭环零点和极点，并在同一坐标系内绘制它们的阶跃响应。 D (s ) 图1 具有参考和干扰输入的反馈系统方框图 实验结果： 参考输入的CL 极点：-49.3658 -7.3336 + 7.9786i -7.3336 - 7.9786i -3.4670 参考输入的DC 增益：320 干扰输入的CL 零点：-45 干扰输入的CL 极点：-49.3658 -7.3336 + 7.9786i -7.3336 - 7.9786i -3.4670 干扰输入的DC 增益：4 -20 MATLAB的实验仿真 目录 实验一MATLAB在控制系统模型建立与仿真中的应用 (1) 实验二典型系统的时域响应分析 (13) 实验三线性控制系统的根轨迹与频域分析 (17) 实验四线性系统的校正 (22) 附录一 MATLAB6.5 控制系统工具箱函数和结构化的控制语句 (30) 附录二 SIMULINK 基本模块介绍 (34) 实验一 MATLAB 在控制系统模型建立与仿真中的应用 一、 MATLAB 基本操作与使用 1. 实验目的 1) 掌握MATLAB 仿真软件的安装及启动，熟悉 MATLAB 工作环境平台。 2) MATLAB 命令窗口，包括工具条以及菜单选项的使用；MATLAB 语言的基本规定，包括数值的表示、变量命名规定、基本运算符、预定义变量以及表达式等。 3) MATLAB 图形绘制功能、M 文件程序设计和线性控制系统传递函数模型的建立等。 2. 实验仪器 PC 计算机一台，MATLAB 软件1套 3. 实验内容 1) MATLAB 的启动 这里介绍MATLAB 装入硬盘后，如何创建MATLAB 的工作环境。 方法一 MATLAB 的工作环境由matlab.exe 创建，该程序驻留在文件夹matlab\bin\ 中。它的图标是 matlab 。只要从<我的电脑>或<资源管理器>中去找这个程序，然后双击此图标，就会自动创建如图1所示的MATLAB6.5 版的工作平台 。 图1 在英文Windows 平台上的MATLAB6.5 MATLAB 工作平台 方法二 假如经常使用MATLAB ，则可以在Windows 桌面上创建一个MATLAB 快捷方式图标。具体办法为： 把<我的电脑>中的 matlab 图标用鼠标点亮，然后直接把此图标拖到Windows 桌面上即可。此后，直接双击Windows 桌面上的matlab 图标，就可建立图1所示的 MATLAB 工作平台。 2) MATLAB 工作环境平台 桌面平台是各桌面组件的展示平台，默认设置情况下的桌面平台包括 6 个窗口，具体如下： ① MATLAB 窗口 Command Window Wireless Network Experiment Three: Queuing Theory ABSTRACT This experiment is designed to learn the fundamentals of the queuing theory. Mainly about the M/M/S and M/M/n/n queuing models. KEY WORDS: queuing theory, M/M/s, M/M/n/n, Erlang B, Erlang C. INTRODUCTION A queue is a waiting line and queueing theory is the mathematical theory of waiting lines. More generally, queueing theory is concerned with the mathematical modeling and analysis of systems that provide service to random demands. In communication networks, queues are encountered everywhere. For example, the incoming data packets are randomly arrived and buffered, waiting for the router to deliver. Such situation is considered as a queue. A queueing model is an abstract description of such a system. Typically, a queueing model represents (1) the system's physical configuration, by specifying the number and arrangement of the servers, and (2) the stochastic nature of the demands, by specifying the variability in the arrival process and in the service process. The essence of queueing theory is that it takes into account the randomness of the arrival process and the randomness of the service process. The most common assumption about the arrival process is that the customer arrivals follow a开关电源《基于MatlabSimulink的BOOST电路仿真》
MATLAB仿真实验全部
queuing modeling排队论的matlab仿真(包括仿真代码)
matlab单服务台排队系统实验报告
基于matlab的电路仿真
增量调制MATLAB仿真实验
基于MATLAB的电力系统仿真
高频电子线路Matlab仿真实验
Matlab 编程方法及仿真实验
Matlab仿真实验教程
排队论地matlab仿真(包括仿真代码)
相关主题
文本预览