控制系统仿真实验报告0717014636Word格式文档下载.docx
《控制系统仿真实验报告0717014636Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《控制系统仿真实验报告0717014636Word格式文档下载.docx（13页珍藏版）》请在冰豆网上搜索。
3、调试M文件，验证KCL、KVL；
4、掌握用simulink模块搭建电路并且进行仿真。
二、实验内容
电路如图1所示，该电路是一个分压电路，已知R13，R27，VS20V。
试求恒
压源的电流I和电压V1、V2。
I
R1V1
+
VS
_
R2V2
图1
三、列写电路方程
（1）用欧姆定律求出电流和电压
（2）通过KCL和KVL求解电流和电压
（1）I=Vs/（R1+R2）=2A,V1=I*R1=6V,V2=I*R2=14V
（2）I*R1+I*R2-Vs=0,V1=I*R1,V2=I*R2,
=>
I=2A,V1=6V,V2=14V.
四、编写M文件进行电路求解
（1）M文件源程序
（2）M文件求解结果
（1）M文件源程序
R1=3;
R2=7;
Vs=20;
I=Vs/（R1+R2）
V1=I*R1
V2=Vs-V1
（2）M文件求解结果
I=2
V1=6
V2=14
五、用simulink进行仿真建模
（1）给出simulink下的电路建模图
（2）给出simulink仿真的波形和数值
电流I波形I=2A
电压U1波形,U1=6V
电压U2波形,U2=14V
六、结果比较与分析
根据M文件编程输入到matlab中，实验结果与理论计算结果一致。
实验二数值算法编程实现
掌握各种计算方法的基本原理，在计算机上利用MATLAB完成算法程序的编写拉
格朗日插值算法程序，利用编写的算法程序进行实例的运算。
二、实验说明
1．给出拉格朗日插值法计算数据表；
2．利用拉格朗日插值公式，编写编程算法流程，画出程序框图，作为下述编程的依
据；
3．根据MATLAB软件特点和算法流程框图，利用MATLAB软件进行上机编程；
4．调试和完善MATLAB程序；
5．由编写的程序根据实验要求得到实验计算的结果。
三、实验原始数据
上机编写拉格朗日插值算法的程序，并以下面给出的函数表为数据基础，在整个插
值区间上采用拉格朗日插值法计算
f（0.6），写出程序源代码，输出计算结果：
x
-2.15
-1.00
0.01
1.02
2.03
3.25
f（x）
17.03
7.24
1.05
17.06
23.05
四、拉格朗日插值算法公式及流程框图
N=5,则使用六点插值公式：
（x
x1）（x
x2）（x
x3）（x
x4）（x
x5）
L5（x）
y0（x0x1）（x0
x2）（x0
x3）（x0x4）（x0x5）
x0）（x
y1（x1
x0）（x1
x2）（x1
x3）（x1
x4）（x1
x1）（x
y2（x2
x0）（x2
x1）（x2
x3）（x2
x4）（x2
y3（x3
x0）（x3
x1）（x3
x2）（x3
x4）（x3
y4（x4
x0）（x4
x1）（x4
x2）（x4
x3）（x4
x4）
y5（x5
x0）（x5
x1）（x5
x2）（x5
x3）（x5
x0
2.15
x1
1.00
x2
x3
x4
x5
3.25x0.6
y0
y1
y2
y3
y4
y5
输入X,Y、n、x
y0
k=1,⋯n
p1
j=1⋯n
k=j?
pp*（xxj）/（xkxj）
yyp*yk
输出x，y
五、程序代码
functiony=lagrange（x0,y0,x）;
x0=[-2.15
3.25];
y0=[17.037.24
23.05];
x=[0.6]
n=length（x0）;
m=length（x）;
fori=1:
m
z=x（i）;
s=0.0;
fork=1:
n
p=1.0;
forj=1:
ifj~=k
p=p*（z-x0（j））/（x0（k）-x0（j））;
end
s=p*y0（k）+s;
y（i）=s;
六、计算结果
f（0.6）=0.0201
实验三动态电路的建模及仿真
1．了解动态电路的理论，掌握动态电路建模的基本原理；
2．熟悉MATLAB的Simulink模块，并掌握使用模块搭建过程。
电力系统是一个大规模、时变的复杂系统，主要由发电、变电、输电、配电和用电
等环节组成，在国民经济中占有非常重要的作用。
动态过程是电力系统中的存在的常态
结构形式，为了更好的理解动态过程的理论，掌握动态过程的物理本质，本实验利用
MATLAB搭建一个包含RLC元件的简单动态系统，采用编程和数值模型仿真的方法分别进行分析计算，得出计算结果。
三、实验内容
电路如图2所示，该电路是一个
RLC电路，已知
200
6
VS
，
L
410（H）
V
C4106（F），R1.5
。
试求电感的电流IL和电容的电压UC。
IL
R
C
UC
图2
由基尔霍夫定律分析得：
dUC
1UC
1IL
dt
RC
dIL
1
US
functionsecord=zy3（t,y）
US=200;
L=0.000004;
C=0.000004;
R=1.5;
secord=[（-1/（R*C））*y
（1）+（1/C）*y
（2）;
（-1/L）*y
（1）+（US/L）];
然后在commandandwindow上输入：
clear;
tspan=[0,10];
y0=[0;
0];
tspan=[0,0.0001];
[t,yy]=ode45（'
zy3'
tspan,y0）;
plot（t,yy）;
xlabel（'
t'
）;
在使用MATLAB解决正弦稳态电路的时候一般可以采用欧拉公式法和相量法两种思路进行理论分析计算。
掌握SIMULINK建模的方法，使用SIMULINK对正弦稳态电路进行建模分析。
实验四正弦稳态电路的设计及仿真
1、掌握正弦稳态电路分析的方法；
2、并掌握使用模块搭建过程。
电力系统在正弦稳态电源的激励下会在整个系统中产生正弦稳态的响应，一般在对
正弦稳态电路进行分析的时候采用相量法进行分析。
在使用MATLAB解决正弦稳态电
路的时候一般可以采用欧拉公式法和相量法两种思路进行理论分析计算。
掌握
SIMULINK建模的方法，使用SIMULINK对正弦稳态电路进行建模分析。
电路如图3所示，该电路是一个正弦稳态电路，已知R1R2R35，XL2，
XC3，US1150（V），US250（V），IS30（A）。
试求各支路电流，并且绘
制各支路电流的时间曲线。
is
i2i3
R2
i5
i1
i4
R3
R1
US1_
US2
图3
R1=5;
R2=5;
R3=5;
XL=2;
XC=3;
US1=15;
US2=5;
IS=3;
R=[11100;
00-111;
R1-j*XL000;
R10-R2-R30;
000-R3-j*XC];
U=[-IS;
IS;
US1;
US2;
US2];
I=R\U
UR=I
（1）*R1
UC=-j*XC*I（5）
w=314;
t=0:
0.001:
0.1;
I1=abs（I
（1））*sin（w*t+angle（I
（1）））;
I2=abs（I
（2））*sin（w
升级会员 