matlab sprawko6

clear all;
close all;

T1=1;
dT=0.001;

K=1;

%==== 1 ====%

M1=[T1, 1];

T2=0;
M11=[T2, 1];

T2_1=T1/10;
M21=[T2_1, 1];

T2_2=T1/2;
M31=[T2_2, 1];

T2_3=T1 + dT;
M41=[T2_3, 1];
%sim("schemat");
%plot([0:10], ans.iner);

%==== 2 ====%
T1=2;
M2=[T1, 0];

T2=0;
M12=[T2, 1];

T2_1=T1/100;
M22=[T2_1, 1];

T2_2=T1/10;
M32=[T2_2, 1];

T2_3=T1 + dT;
M42=[T2_3, 1];

T2_4=10*T1;
M52=[T2_4, 1];

%==== 3 ====%

Td=1.5;
T2a_0=0.0001;
T2a_1=Td/100;
T2a_2=Td/10;
T2a_3=Td + dT;
T2a_4=10*Td;



T_1=2.14;
T_2=2.43;
T_3=3.29;
T_4=4;

sim("schemat");
%{
figure (1);
subplot(321);
plot(ans.czas, ans.rozn1);
xlabel("czas [s]");
ylabel("wartosc");
title("Człon Różniczkowy T_2 = 0.0001, K = 1, Td = 1.5");
subplot(322);
plot(ans.czas, ans.rozn2);
xlabel("czas [s]");
ylabel("wartosc");
title("Człon Różniczkowy T_2 = Td/100, K = 1, Td = 1.5");
subplot(323);

plot(ans.czas, ans.rozn3);
xlabel("czas [s]");
ylabel("wartosc");
title("Człon Różniczkowy T_2 = Td/10, K = 1, Td = 1.5");

subplot(324);
plot(ans.czas, ans.rozn4);
xlabel("czas [s]");
ylabel("wartosc");
title("Człon Różniczkowy T_2 = Td + dT, K = 1, Td = 1.5");

subplot(313);
plot(ans.czas, ans.rozn5);
xlabel("czas [s]");
ylabel("wartosc");
title("Człon Różniczkowy T_2 = 10*Td, K = 1, Td = 1.5");


%=================%
figure (2);
subplot(221);
plot(ans.czas, ans.iner1);
hold on;
plot(ans.czas, ans.iner5);
xlabel("czas [s]");
ylabel("wartosc");
legend("Model Inercyjny: T_2 = 0, K = 1, T_1 = 1", "Model Kupfmullera: T_0 = 1, T = 2.14");


subplot(222);
plot(ans.czas, ans.iner2);
hold on;
plot(ans.czas, ans.iner6);
xlabel("czas [s]");
ylabel("wartosc");
legend("Model Inercyjny: T_2 = T_1/10, K = 1, T_1 = 1", "Model Kupfmullera: T_0 = 1.07, T = 2.43");


subplot(223);
plot(ans.czas, ans.iner3);
hold on;
plot(ans.czas, ans.iner7);
xlabel("czas [s]");
ylabel("wartosc");
legend("Model Inercyjny: T_2 = T_1/2, K = 1, T_1 = 1", "Model Kupfmullera: T_0 = 1.21, T = 3.29");

subplot(224);
plot(ans.czas, ans.iner4);
hold on;
plot(ans.czas, ans.iner8);
xlabel("czas [s]");
ylabel("wartosc");
legend("Model Inercyjny: T_2 \approx T_1, K = 1, T_1 = 1", "Model Kupfmullera: T_0 = 1.29, T = 4");





%
sim("schemat");
figure (1);
plot(ans.czas, ans.iner);
xlabel("czas [s]");
ylabel("wartosc");
title("Człon inercyjny");
legend("T_2 = 0", "T_2 = T_1/10", "T_2 = T_1/2", "T_2 \approx T_1");


figure (2);
plot(ans.czas, ans.calk);
title("Człon calkujacy");
xlabel("czas [s]");
ylabel("wartosc");
legend("T_2 = 0", "T_2 = T_i/10", "T_2 = T_i/2", "T_2 \approx T_i", "T_2 = 10T_i");

figure (3);
plot(ans.czas, ans.rozn);
title("Człon różniczkujący");
xlabel("czas [s]");
ylabel("wartosc");
legend("T_2 = 0", "T_2 = T_i/10", "T_2 = T_i/2", "T_2 \approx T_i", "T_2 = 10T_i");
%}