%Zmienne do wszystkich transmitancji: k1=1; k2=2; k3=2.5; k4=4; T1=0.1; T2=1; T3=3.8; T4=5; ksi1=0.5; ksi2=1.5; ksi3=3; Ti1=3.14/2; Ti2=3.14/4; Ti3=3.2; Td1=0.4; Td2=0.9; Td3=1.2; %1 %Licznik i mianownik do I transmitancji: licz1=[0,k1]; licz2=[0,k2]; licz3=[0,k3]; licz4=[0,k4]; mian1=[T1,1]; mian2=[T2,1]; mian3=[T3,1]; %Odpowiedzi dla stałego k i różnego T: t1=tf(licz1,mian1); t2=tf(licz1,mian2); t3=tf(licz1,mian3); step(t1, t2, t3); figure(1); legend('k=1 T=0.1','k=1 T=1','k=1 T=3.8'); figure(2); impulse(t1, t2, t3); legend('k=1 T=0.1','k=1 T=1','k=1 T=3.8'); %Odpowiedzi dla stałego T i róznego k: t1=tf(licz2,mian2); t2=tf(licz3,mian2); t3=tf(licz4,mian2); figure(3); step(t1,t2,t3); legend('k=1 T=1','k=2 T=1','k=4 T=1'); figure(4); impulse(t1,t2,t3); legend('k=2 T=1','k=2.5 T=1','k=4 T=1'); %2 %Licznik i mianownik do II transmitancji: licz1=[0,0,k1]; licz2=[0,0,k2]; licz3=[0,0,k3]; mian1=[T1*T2,T1+T2,1]; mian2=[T1*T3,T1+T3,1]; mian3=[T2*T4,T2+T4,1]; %odpowiedzi dla k=1,T1=0.1,T2=1 i k=2,T1=0.1,T2=0.1; t1=tf(licz1,mian1); t2=tf(licz2,mian2); t3=tf(licz3,mian3); figure(5); step(t1,t2,t3); legend('k=1 T1=0.1 T2=1','k=2 T1=0.1 T2=3.8','k=2.5 T1=1, T2=5'); figure(6); impulse(t1,t2,t3); legend('k=1 T1=0.1 T2=1','k=2 T1=0.1 T2=3.8','k=2.5 T1=1, T2=5'); %3 %Mianowniki do III transmitancji, licznik bez zmian: mian1=[T1^2,2*ksi1*T1,1]; mian2=[T1+0.3^2,2*ksi2*T1,1]; mian3=[T1+0.6^2,2*ksi3*T1,1]; %Odpowiedzi dla ksi>1 i ksi<1: t1=tf(licz1,mian1); t2=tf(licz1,mian2); t3=tf(licz1,mian3); figure(7); step(t1,t2,t3); legend('k=1 T=0.1 ksi=0.5','k=1 T=0.4 ksi=1.5','T=0.7 ksi=3'); figure(8); impulse(t1,t2,t3); legend('k=1 T=0.1 ksi=0.5','k=2 T=0.4 ksi=1.5','T=0.7 ksi=3'); %4 %Licznik i mianownik do IV transmitancji: licz1=[0,k1]; mian1=[T3-0.5*Ti1,Ti1,0]; mian2=[T3*Ti2,Ti2,0]; mian3=[T3+0.5*Ti3,Ti3,0]; %Odpowiedzi i=TT/2, i=TT/4 dla k=1 i T=3.8: t1=tf(licz1,mian1); t2=tf(licz1,mian2); t3=tf(licz1,mian3); figure(9); step(t1,t2,t3); legend('k=1 T=3.3 Ti=TT/2','k=1 T=3.8 Ti=TT/4','T=4.3 Ti=3.2'); figure(10); impulse(t1,t2,t3); legend('k=1 T=3.3 Ti=TT/2','k=1 T=3.8 Ti=TT/4','T=4.3 Ti=3.2'); %5 %Licznik i mianownik do V transmitancji: licz1 = [Td1,0]; licz2 = [Td2,0]; licz3 = [Td3,0]; mian1 = [T2,1]; %Odpowiedzi Td=0.4 i Td=0.9 dla T=1: t1=tf(licz1,mian1); t2=tf(licz2,mian1+0.2); t3=tf(licz3,mian1+0.4); figure(11); step(t1,t2,t3); legend('Td=0.4 T=1','Td=0.9 T=1.2','Td=1.2 T=1.4'); figure(12); impulse(t1,t2,t3); legend('Td=0.4 T=1','Td=0.9 T=1.2','Td=1.2 T=1.4'); %6 %Opóźnienie: n1=1; n2=3; n3=5; theta1=0.5; theta2=2; theta3=2.5; [licz_op1, mian_op1]=pade(theta1, n1); [licz_op2, mian_op2]=pade(theta2, n2); [licz_op3, mian_op3]=pade(theta3, n3); %Licznik i mianownik dla VI transmitancji: licz_inner=[0,k1]; mian_inner=[T2,1]; [licz1,mian1]=series(licz_op1,mian_op1, licz_inner,mian_inner); [licz2,mian2]=series(licz_op2,mian_op2, licz_inner+0.5,mian_inner+1); [licz3,mian3]=series(licz_op3,mian_op3, licz_inner+1.5,mian_inner+0.6); %odpowiedzi z opóźnieniem 0,5s i 2s dla T=1 i k=1; t1=tf(licz1,mian1); t2=tf(licz2,mian2); t3=tf(licz3,mian3); figure(13); step(t1,t2,t3); legend('opoznienie=0.5s n=2 T=1 k=1', 'opoznienie=2s n=3 T=2 k=1.5','opoznienie=2.5s n=5 T=1.6 k=2.5'); figure(14); impulse(t1,t2,t3); legend('opoznienie=0.5s n=2 T=1 k=1','opoznienie=2s n=3 T=2 k=1.5','opoznienie=2.5s n=5 T=1.6 k=2.5');