计算机仿真习题及答案
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
计算机仿真试题
1．编写一个函数，使其能够产生如下的分段函数：错误！未找到引用源。
并调用此函数，绘制x=[0,+2]范围内的f(x)*f(x+2) 。（10分）function y=f(x)
if x<=2
y=0.5*x;
else if x>6
y=0.5;
else y=1.5-0.25*x;
end
end
x=0:0.05:2;
y= f(x)’*f(x+2));
plot(x,y)
图 1-1
2．已知4阶龙格-库塔算法如下：
试利用该算法求解以下微分方程：
（15分）本题可以调用MATLAB函数中龙格-库塔算法函数ode45，首先编写m文件：
function dy=func(x,y)
dy=-y+1;
end
再在主窗口调用此文件：
[x,y]=ode45('func',[0,5],0)%这里的[0,5]为任取区间，表示方程在此范围的解。运行结果如下:
x =
0.0001
0.0001
0.0002
0.0002
0.0005
0.0007
0.0010
0.0012
0.0025
0.0037
0.0050
0.0062
0.0125
0.0188
0.0251
0.0313 0.0627
0.0941
0.1255
0.1569
0.2819
0.4069
0.5319
0.6569
0.7819
0.9069
1.0319
1.1569
1.2819
1.4069
1.5319
1.6569
1.7819
1.9069
2.0319
2.1569
2.2819
2.4069
2.5319
2.6569
2.7819
2.9069
3.0319
3.1569
3.2819
3.4069
3.5319
3.6569
3.7819
3.9069
4.0319
4.1569
4.2819
4.4069
4.5319
4.6569
4.7427
4.8285
4.9142
5.0000
y =
0.0001
0.0001
0.0002
0.0002
0.0005
0.0007
0.0010
0.0012
0.0025 0.0037 0.0050 0.0062 0.0124 0.0186 0.0248 0.0309 0.0608 0.0898 0.1180 0.1452 0.2457 0.3343
0.4125
0.4816
0.5425
0.5963
0.6437
0.6855
0.7225
0.7551
0.7839
0.8093
0.8317
0.8515
0.8689
0.8843
0.8979
0.9099
0.9205
0.9298
0.9381
0.9454
0.9518
0.9574
0.9624
0.9669
0.9708
0.9742
0.9772
0.9799
0.9823
0.9843
0.9862
0.9878
0.9892
0.9905
0.9913
0.9920
0.9927
0.9933
为只管起见，我们使用函数命令画出x-y（plot(x,y)）的关系如下图：
图1-2
3．用matlab语言求下列系统的状态方程、传递函数、零极点增益、和部分分式形式的模型参数，并分别写出其相应的数学模型表达式：（15分）
（1）G（s）=
32
432
72424
10355024
s s s
s s s s
+++
++++
（2）
.
X=
2.25 -5 -1.25 -0.54
2.25 -4.25 -1.25 -0.252
0.25 -0.5 -1.25 -12
1.25 -1.75 -0.25 -0.75 0
X
⎡⎤⎡⎤
⎢⎥⎢⎥
⎢⎥⎢⎥
+⎢⎥⎢⎥
⎢⎥⎢⎥
⎣⎦⎣⎦u
y= [0 2 0 2] X
解：（1）
a)求对应状态方程参数：
num=[1 07 24 24]; den=[1 10 35 50 24]; [A,B,C,D]=tf2ss(num,den) 运行结果：
A =
-10 -35 -50 -24
1 0 0 0
0 1 0 0
0 0 1 0
B =
1
C =
1 7 24 24
D =