matlab扰动观测法,扰动观测器DOB matlab仿真

【实例简介】

【实例截图】

【核心代码】

tol=400*10^(-6);

[np,dp]=pade(tol,6);

delay=tf(np,dp);

delta=tf(np,dp)-1;

sys=1/delta;

figure(1);

bode(1/delta,'r',{5,10^5});

grid on;

tol1=0.00035;

Q1=tf([3*tol1,1],[tol1^4,4*tol1^3,6*tol1^2,4*tol1,1]);

hold on;

bode(Q1,'y');

tol2=0.00125;

Q2=tf([3*tol2,1],[tol2^4,4*tol2^3,6*tol2^2,4*tol2,1]);

hold on;

bode(Q2,'b');

tol3=0.00425;

Q3=tf([3*tol3,1],[tol3^4,4*tol3^3,6*tol3^2,4*tol3,1]);

hold on;

bode(Q1,'g');

Gp=tf([56870],[0.00053,1,112,94790]);

Gn=tf([56865],[0.00050,1,111,94785]);

Q1=tf([3*tol1,1],[tol1^4,4*tol1^3,6*tol1^2,4*tol1,1]);

bode(Q1);

QGn=Q1/Gn;

[num,den]=tfdata(QGn,'v');

Y=sim('temp');