Your report must contain necessary block diagrams and characteristics of dynamic links in question

K=1;

T=1;

w=0:0.05:3;

C=w.^2;

D=T.^2;

E=1+C*D;

B=sqrt(E)

Gfreq=K/B;

plot(w, Gfreq)

title('Gain Frequency Response')

xlabel('Omega')

ylabel('Gain')

grid on

Matlab script:

w=0:0.1:6;

t3=4;

t4=2;

x1=w*t3;

x2= w*t4;

phi=-atan(x1)-atan(x2)

y=phi*180/pi;

plot(w,y)

title('Phase Frequency Response')

xlabel('Omega')

ylabel('Phase angle')

grid on

w=0.6;

t=2;

a=w^2*t^2;

a = 1.44

w=0;

t=2;

zeta=0.7;

a=(2*zeta*w*t)/(1-w.^2*t.^2);

alpha=-atan(a)

y=alpha*180/pi;

a,alpha,y

Response:

a =

alpha =

y =

w=3;

t=2;

zeta=0.7;

a=(2*zeta*w*t)/(1-w.^2*t.^2);

alpha=-atan(a)

y=alpha*180/pi;

a,alpha,y

Response:

a =

-0.2400

alpha =

0.2355

y =

13.4957

w=6;

t=2;

zeta=0.7;

a=(2*zeta*w*t)/(1-w.^2*t.^2);

alpha=-atan(a)

y=alpha*180/pi;

a,alpha,y

Response:

a =

-0.1175

alpha =

0.1169

y =

6.7005

Checkpoint:

w=3;

t=4;

zeta=0.7;

a=(2*zeta*w*t)/(1-w.^2*t.^2);

alpha=atan(a)

y=alpha*180/pi;

a,alpha,y

Response:

alpha =

-0.1169

a =

-0.1175

alpha =

-0.1169

y =

-6.7005

w=1;

t=2;

zeta=0.7;

a=(2*zeta*w*t)/(1-w.^2*t.^2);

alpha=abs(atan(a))

y=-alpha*180/pi;

a,alpha,y

alpha =

0.7509

a =

-0.9333

alpha =

0.7509

y =

-43.0251

w=0.5;

t=2;

zeta=0.7;

a=(2*zeta*w*t)/(1-w.^2*t.^2);

alpha=abs(atan(a))

y=-alpha*180/pi;

a,alpha,y

zeta=0.7;

w=0:0.1:1

t=1;

a=(2*zeta*w*t)/(1-w.^2*t.^2);

alpha=-atan(a)

y=alpha*180/pi;

a,alpha,y

plot(w,y,'--rs')

title('Phase Frequency Response')

xlabel('Omega')

ylabel('Phase angle')

grid on

hold on

hold off

w=0:0.1:2;

for k=1:4

t=[1 2 3 4];

zeta=[0.2 0.5 0.7 0.9];

a=(2*zeta(k)*w*t(k))/(1-w.^2*t(k).^2);

alpha=-atan(a)

y=alpha*180/pi;

a,alpha,y

plot(w,y)

title('Phase Frequency Response')

xlabel('Omega')

ylabel('Phase angle')

grid on

hold on

end

hold off

Поиск по сайту: