电气控制CAD.doc

3．绘制下面的图形。
（1）sin(1/t),-1<t<1
t=-1:0.01:1
y=sin(1./t)
plot(t,y)
grid on
xlabel('t')
ylabel('y')
gtext('sin(1/t)')

y=1-cos(7.*t).^3

t=-1:0.001:1
y=1-cos(7.*t).^3
plot(t,y)
grid on
xlabel('t')
ylabel('t')
gtext('y=1-cos(7.*t).^3')

微分方程
num=[10]
den=[1,8,36,40,10]
step(num,den)
grid on
xlabel('time')
ylabel('Amplitude')
gtext('step response')

2.控制工具箱
num=[0 0 0 0 10];
den=[1 8 36 40 10];
[A B C D]=tf2ss(num,den)
A =
-8 -36 -40 -10
1 0 0 0
0 1 0 0
0 0 1 0
B =
1
0
0
0
C = 0 0 0 10
D = 0
function dx=wffc(t,x)
u=1;
dx=[-8*x(1)-36*x(2)-40*x(3)-10*x(4)+u;x(1);x(2);x(3)]
[t,x]=ode45('wffc',[0 8],[0;0;0;0]);
y=10*x(:,4);
plot(t,y);
grid on

Simulink：
simulink分析

输入为t
6．某小功率随动系统动态结构如图所示，已知：
若系统输入分别为，适用simulink分析系统的输出分别如何？

输入为t

输入为1