机电控制系统仿真 实验报告Word文档下载推荐.docx
- 文档编号:17226557
- 上传时间:2022-11-29
- 格式:DOCX
- 页数:16
- 大小:333.50KB
机电控制系统仿真 实验报告Word文档下载推荐.docx
《机电控制系统仿真 实验报告Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《机电控制系统仿真 实验报告Word文档下载推荐.docx(16页珍藏版)》请在冰豆网上搜索。
a=rand(3,2)
a=1;
b=2;
c=3;
x=[5bc;
a*ba+cc/b]
a=linspace(1,10,10)
a=eye(3,2)
b=eye(3)
a=zeros(3,2)
b=zeros(3)
对应于问题②的程序如下:
a=[123;
456;
235];
b=inv(a)
c=det(a)
d=rank(a)
345];
[m,n]=size(a)
length(A)=max(size(A))
456]'
a1=[123;
456].'
b=[1+2i2-7i]'
b2=[1+2i2-7i].'
对应于问题③的程序如下:
xlabel('
time/s'
)
ylabel('
amplitude/1'
axis([0,6.28,-1.5,1.5])
figure
(2)
subplot(2,1,1)
plot(t,y2)
subplot(2,1,2)
plot(t,y3)
t=0:
0.01:
2*pi;
y=sin(t);
y1=sin(t+0.5*pi);
y2=cos(t);
y3=cos(t+0.5*pi);
figure
(1)
subplot(1,2,1)
plot(t,y)
subplot(1,2,2)
plot(t,y1)
对应于问题④的程序如下:
p=[1-12025116];
a=roots(p)
r=[11.7473,2.7028,-1.2251+1.4672i,-1.2251-1.4672i];
p=poly(r)
a=[123];
b=[12];
c=conv(a,b)
对应于问题⑤的程序如下:
A=zeros(9);
fori=1:
1:
9
forj=1:
a(i,j)=i*j;
ifi<
j
a(i,j)=0
end
end
1.3运行结果
对应于问题①,运行程序后得到结果如下:
对应于问题②,运行程序后得到结果如下:
对应于问题③,运行程序后得到结果如下:
对应于问题④,运行程序后得到结果如下:
a=
11.7473
2.7028
-1.2251+1.4672i
-1.2251-1.4672i
p=1.0000-11.9999-0.001525.0012116.0023
c=1476
对应于问题⑤,运行程序后得到结果如下:
第二次上机:
MATLAB的符号运算
2.1问题描述
本次上机共进行了5个问题的调试,分别是:
1命令窗口中符号矩阵的定义;
2命令窗口中符号矩阵的修改;
3计算不定积分
4拉氏变换和反拉氏变换
5方程求解
6用matlab表示传递函数
2.2程序编制
A=sym('
[a,2*b;
3*a,0]'
B=['
[a,2*b]'
;
'
[3*a,0]'
]
A1=subs(A,'
b'
'
c'
)
A1(2,2)='
4*b'
m=int(int('
x*exp(-x*y)'
'
x'
),'
y'
n=laplace(m)
k=ilaplace(n)
f='
a*x^2+b*x+c'
solve(f)%方程f=ax2+bx+c求解
solve(f,'
)%指定求解变量为b
num=[1,1];
den=[1,2,1];
%传递函数分子分母系数(降幂排列)
G=tf(num,den)
2.3运行结果
B=
[a,2*b]
[3*a,0]
A=
[a,2*b]
[3*a,0]
A1=
[a,2*(c)]
[3*a,4*b]
m=
1/y*exp(-x*y)
n=k=
1/y/(s+y)1/y*exp(-y*t)
ans=
1/2/a*(-b+(b^2-4*a*c)^(1/2))
1/2/a*(-b-(b^2-4*a*c)^(1/2))
对应于问题⑥,运行程序后得到结果如下:
Transferfunction:
s+1
-------------
s^2+2s+1
第三次上机:
控制系统的MATLAB分析
3.1问题描述
本次上机共进行了8个问题的调试,分别是:
1系统传递函数的求解;
2带反馈的传递函数系统阶跃响应和单位脉冲响应求解;
3阻尼比变化时的阶跃响应求解
4系数变化时的系统阶跃响应求解
5输入为任意函数的系统响应求解
6判定系统稳定性
7典型环节的波德(Bode)图绘制
8奈奎斯特图的绘制
3.2程序编制
symstsrRcC;
r=1*sym('
heaviside(t)'
)+t*sym('
);
c=t+0.9-0.9*exp(-10*t);
R=laplace(r);
C=laplace(c);
G0=C/R;
G=factor(G0)
num=[20];
den=[1836400];
sys1=tf(num,den);
sys=feedback(sys1,1);
0.1:
10;
y=step(sys,t);
%dc=dcgain(numc,denc)
dc=dcgain(sys)
y=impulse(sys,t);
num=1;
19.9
forbc=0.1:
1
den=[1,2*bc,1];
sys=tf(num,den);
i=i+1;
step(sys,t)
holdon
den=conv(conv([10],[0.51]),[4,1]);
rangek=[1.4,2.3,3.5];
t=linspace(0,20,200)'
forj=1:
3
s1=tf(num*rangek(j),den);
sys=feedback(s1,1);
y(:
j)=step(sys,t);
plot(t,y(:
1:
3)),grid
gtext('
k=1.4'
),gtext('
k=2.3'
k=3.5'
对应于问题⑥的程序如下:
holdon;
plot(t,u,'
:
'
clear
P=[12345];
roots(P)
den=[11];
u=sin(2*t);
lsim(sys,u,t);
%固有频率为1,阻尼比0.1的振荡环节
num=[1];
den=[10.01];
sys2=tf(num,den);
figure(3)
bode(sys2)
%惯性环节
bode(sys1)
%积分环节
den=[10];
bode(sys)
对应于问题⑦的程序如下:
sys=tf(num3,den3);
nyquist(sys);
title('
nyquistplot'
figure
(2)%绘制冲击响应曲线验证判断
[numc,denc]=cloop(num3,den3);
impulse(numc,denc)
num1=[16.70];
den1=conv([0.851],conv([0.251],[0.06251]));
[num2,den2]=cloop(num1,den1);
num3=10*num2;
den3=den2;
[z,p,k]=tf2zp(num3,den3);
%P=num3+den3;
%roots(P)
对应于问题⑧的程序如下:
3.3运行结果
对应于问题⑦,运行程序后得到结果如下:
对应于问题⑧,运行程序后得到结果如下:
第四次上机:
控制系统的Simulink仿真
4.1问题描述
本次上机共进行了1个问题的调试:
带单位响应的阶跃响应模型建立与验证。
4.2程序编制
模型建立如下:
“StartSimulation”之后,双击“Scope”的下图:
“Command”验证程序及结果如下:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 机电控制系统仿真 实验报告 机电 控制系统 仿真 实验 报告