Imię Nazwisko: Maciej Woźniak Termin: środa, 17.05

Nr albumu: 218467 Data lab.: 13.01.2016

Karta sprawozdania nr. 12 i 13: Człony dynamiki: dziedzina czasu i częstotliwości

Dla wybranej transmitancji z miniprojektu dwoma sposobami:

I Schemat Simulink

II Funkcja

Wykresy są identyczne dla obu sposobów.

Dla tej transmitancji utworzono także charakterystykę Bodego. Naniesiono na wykres 3 asymptoty. Z nich wynika, że wykres jest poprawny dla podanej transmitancji.

Skrypt

clear;

cp=1000;

rop=1.2;

qk=20000;

Tkz=30;

Tzew=-20;

Tw1=20;

Tw2=15;

fp=qk/(cp*rop*Tkz);

Ks=(qk-cp*rop*fp*Tw2)/(Tw1+Tw2-2*Tzew);

K0=(qk-cp*rop*fp*Tw1-((qk-cp*rop*fp*Tw2)/(Tw1+Tw2-2*Tzew))*(Tw1-Tzew))/(Tw1-Tw2);

V1=200; %m3

V2=150; %m3

Cv1=cp*rop*V1;

Cv2=cp*rop*V2;

%transmitancja

X=[(cp*rop*fp)^2+cp*rop*fp*K0];

Z=[Cv1*Cv2 (cp*rop*fp+K0+Ks)*Cv2+(cp*rop*fp+K0+Ks)*Cv1 (cp*rop*fp)^2+Ks^2+2*cp*rop*fp*Ks+cp*rop*fp*K0+2*K0*Ks];

%obiekty w dziedzinie czasu

sim('trans');

hold on;

grid on;

figure(1);

plot(t,simout);

K=(cp*rop*fp)^2+cp*rop*fp*K0;

s=tf('s');

M=(Cv1*Cv2*s^2+((cp*rop*fp+K0+Ks)*Cv2+(cp*rop*fp+K0+Ks)*Cv1)*s+(cp*rop*fp)^2+Ks^2+2*cp*rop*fp*Ks+cp*rop*fp*K0+2*K0*Ks);

H=K/M;

figure(2);

step(H);

%obiekty w dziedzinie czestotliwosci

figure(3);

bode(H);

---