no1.m

电力牵引电机控制

% 方框图初始化
R1=1000;R2=1000;R3=1000;R4=1000;
Km=10;J=2;Ra=1;La=1;
b=0.5;Kb=0.1;Kt=0.1;
n1=(R1+R2)*R4;
d1=R1*(R3+R4);
n2=540;d2=1;
n3=Km;d3=Ra+La;
n4=1;d4=[J b];
n5=Kb;d5=1;
n6=Kt;d6=1;
nblocks=6;                                      % 方框个数
blkbuild;                                       % 建立方框图
sys=ss(a,b,c,d);                                % 建立方框图的对角非连接形式的状态空间模型
q=[1 0 0;2 1 -6;3 2 -5;4 3 0;5 4 0;6 4 0];
sysc=connect(sys,q,1,6);                        %求解方框图的状态空间模型
[A,B,C,D]=ssdata(sysc);                         % 数据扩展
[num,den]=ss2tf(A,B,C,D,1);                     %将状态空间模型转换为传递函数模型
% 显示结果
disp('State-Space Model of the block-diagram is:');
A
B
C
D
disp('Transfer Function of the block-diagram is:');
num
den
printsys(num,den)
sys2=tf(num,den);
t=[0:0.001:1];
step(sys2,t);
axis([0 1 0 1.2]);
title('单位阶跃响应')
xlabel('时间')
ylabel('振幅')
grid on;
M=ctrb(A,B);
rank(M)
N=obsv(A,C);
rank(N)