motor_3_xgfx.m

来自「用matlab实现motor运动控制的控制器设计仿真」· M 代码 · 共 24 行

%Correlation analysis method
t=0:0.1:12.6;           %Get the plus response
num=0.01;
den=[0.005 0.06 0.1001];
sys=tf(num,den);
y=impulse(sys,t);
y=y';
plot(t,y,'b');grid on
hold on
a=1;%Simulate
Ts=0.1;
T=8;
L=127;
for k=1:T
    u((k-1)*L+1:k*L)=idinput(L,'prbs');
end
t1=0:0.1:(8*L-1)/10;
y=lsim(sys,u,t1);
Y=prbs(a,Ts,L,u,y,T);
plot(t,Y,'r');
xlabel('t');
ylabel('y');
legend('Actual response',2,'Theory response',1);