Podwojna inercja

z opoznieniem

te=ScopeData.time;

ue=ScopeData.signals(1).values;

he=ScopeData.signals(2).values;

% subplot(211)

% plot(te,ue);grid

% subplot(212)

% plot(te,he);grid F5

n=min(find(ue>0.3));

t=te(n:end)-te(n);

u=ue(n:end)-ue(n-1);

h=he(n:end)-he(n);

plot(t,h);grid procenty

hu=h(end);

U=u(end);

k=hu/U F5

dh=diff(h)./diff(t);

dt=t(1:end-1);

plot(t,h,dt,dh);grid

legend('h','dh'); F5

[mdh n]=max(dh);

hp=h(n);

factor=hp/hu

tp=t(n);

plot(t,h,tp,hp,'*');grid F5

n10=max(find(h<=0.1*hu));

t10=t(n10);

h10=h(n10);

n90=max(find(h<=0.9*hu));

t90=t(n90);

h90=h(n90);

T=(t90-t10)/3.3

tau=t10-0.53*T F5

[Lp Mp]=pade(tau,12);

L=k*Lp;

M=[T 1];

M1=conv(conv(M,M),Mp);

hm=lsim(L,M1,u,t);

plot(t,h,t,hm);grid F5

kp=0.68*T/(k*tau)

Ti=2*T

Td=Ti/4

tr=7.2*tau

P=kp

I=kp/Ti

D=kp*Td

N=8/Td F5

Inercja

te=ScopeData.time;

ue=ScopeData.signals(1).values;

he=ScopeData.signals(2).values;

% subplot(211)

% plot(te,ue);grid

% subplot(212)

% plot(te,he);grid F5

n=min(find(ue>0.03));

t=te(n:end)-te(n);

u=ue(n:end)-ue(n-1);

h=he(n:end)-he(n);

figure

plot(t,h);grid -procenty

hu=h(end);

U=u(end);

k=hu/U F5

n10=max(find(h<=0.1*hu));

h10=h(n10);

t10=t(n10);

n90=max(find(h<=0.9*hu));

h90=h(n90);

t90=t(n90);

T=(t90-t10)/2.2

L=k;

M=[T 1];

hm=lsim(L,M,u,t);

figure

plot(t,h,t,hm);grid F5

tr=60*10

kp=4*T/(k*tr)

Ti=T

P=kp

I=kp/Ti F5

Inercja z opoznieniem

te=ScopeData.time;

ue=ScopeData.signals(1).values;

he=ScopeData.signals(2).values;

% subplot(211)

% plot(te,ue);grid

% subplot(212)

% plot(te,he);grid F5

n=min(find(ue<=0.3));

t=te(n:end)-te(n);

u=ue(n:end)-ue(n-1);

h=he(n:end)-he(n);

plot(t,h);grid -procenty

hu=h(end);

U=u(end);

k=hu/U F5

n10=max(find(h<=0.1*hu));

h10=h(n10);

t10=t(n10);

n90=max(find(h<=0.9*hu));

h90=h(n90);

t90=t(n90);

T=(t90-t10)/2.2

tau=t10-0.1*T F5

[Lp Mp]=pade(tau,12);

L=k*Lp;

M=[T 1];

M1=conv(M,Mp);

hm=lsim(L,M1,u,t);

plot(t,h,t,hm);grid F5

kp=0.34*T/(k*tau)

Ti=T

Tr=7.2*tau

P=kp

I=kp/Ti F5