mechaniczne

close all; clear; clc;
P1 = [600; 720; 780; 960; 1080; 1200; 1320; 1560; 1800; 2040; 2220; 2520];
V1 = [10; 15; 30; 40; 46; 55; 61; 68; 75; 80; 85; 90]; %v3
w1 = [73.27; 70.65; 69.35; 68.04; 66.73; 65.42; 64.11; 62.54; 60.19; 58.88; 56.26; 54.95];
Mst1 = [0; 0.45; 1.63; 2.94; 4.07; 5.68; 7.43; 9.25; 11.84;13.93; 16.62; 18.84];
PDC1 = [10; 31.5; 112.5; 200; 271.4; 371.25; 475.8; 578; 712.5; 820; 935; 1035];
spr1 = [1.67; 4.38; 14.43; 20.84; 25.13; 30.94; 36.05; 37.06; 39.59; 40.2; 42.12; 41.08];
cosfi1 = [0.24; 0.28; 0.29; 0.34; 0.37; 0.38; 0.4; 0.43; 0.47; 0.48; 0.49; 0.52];
PIM1 = [480; 540; 720; 900; 1020; 1200; 1440; 1740; 2040; 2400; 2640; 3000];

P2 = [720; 840; 960; 1080; 1320; 1500; 1920; 2160; 2520; 2760; 3060; 3360];
V2 = [11; 25; 35; 45; 55; 65; 75; 85; 95; 98; 105; 110];
w2 = [87.66; 86.35; 85.05; 83.74; 82.43; 81.12; 78.5; 77.2; 73.27; 71.96; 69.35; 68.04];
Mst2 = [0;0.9;1.86;3.01;4.68;6.62;9.46;12.12;15.56; 17.03; 20.07; 22.64];
PDC2 = [16.5;77.5;157.5;252;385;536.25;742.5;935;1140; 1225; 1391.25; 1540];
spr2 = [2.3;9.23;16.41;23.34;29.17;35.75;38.68;43.29;45.24; 44.39; 45.47; 45.84];
cosfi2 = [0.28;0.32;0.34;0.37;0.42;0.45;0.53;0.54;0.57; 0.58; 0.6; 0.61];
PIM2 = [480;600;720;960;1200;1440;1800;2160;2520;2880;3120;3600];


[X,Y] = polyfit(V1,w1,2);
Ypoly = polyval(X,V1);

[X1,Y1] = polyfit(V1,P1,2);
Ypoly1 = polyval(X1,V1);

[X2,Y2] = polyfit(Mst1,PIM1,2);
Ypoly2 = polyval(X2,Mst1);

[X3,Y3] = polyfit(Mst1,PDC1,3);
Ypoly3 = polyval(X3,Mst1);

[X4,Y4] = polyfit(Mst1,spr1,3);
Ypoly4 = polyval(X4,Mst1);

[X5,Y5] = polyfit(Mst1,cosfi1,2);
Ypoly5 = polyval(X5,Mst1);

[X6,Y6] = polyfit(Mst1,w1,2);
Ypoly6 = polyval(X6,Mst1);

[X7,Y7] = polyfit(V2,w2,2);
Ypoly7 = polyval(X7,V2);

[X8,Y8] = polyfit(V2,P2,2);
Ypoly8 = polyval(X8,V2);

[X9,Y9] = polyfit(Mst2,PIM2,2);
Ypoly9 = polyval(X9,Mst2);

[X10,Y10] = polyfit(Mst2,PDC2,3);
Ypoly10 = polyval(X10,Mst2);

[X11,Y11] = polyfit(Mst2,spr2,3);
Ypoly11 = polyval(X11,Mst2);

[X12,Y12] = polyfit(Mst2,cosfi2,2);
Ypoly12 = polyval(X12,Mst2);

[X13,Y13] = polyfit(Mst2,w2,2);
Ypoly13 = polyval(X13,Mst2);

figure(1)
plot(V1,w1,'b*',V1,Ypoly,'b',V2,w2,'r*',V2,Ypoly7,'r');
grid on;
xlabel('U_{DC}[V]');
ylabel('ω[rad/s]');
title('Charakterystyka mechaniczna')
subtitle('ω = f(U_{DC})');
legend('V_{TG}=28','','V_{TG}=33.5')

figure(2)
plot(V1,P1,'b*',V1,Ypoly1,'b',V2,P2,'r*',V2,Ypoly8,'r');
grid on;
xlabel('U_{DC}[V]');
ylabel('P_{1}[W]');
title('Charakterystyka mechaniczna')
subtitle('P_{1} = f(U_{DC})');
legend('V_{TG}=28','','V_{TG}=33.5')

figure(3)
plot(Mst1,PIM1,'b*',Mst1,Ypoly2,'b',Mst2,PIM2,'r*',Mst2,Ypoly9,'r');
grid on;
xlabel('M_{st}[Nm]');
ylabel('P_{IM}[W]');
title('Charakterystyka mechaniczna')
subtitle('P_{IM} = f(M_{st})');
legend('V_{TG}=28','','V_{TG}=33.5')

figure(4)
plot(Mst1,PDC1,'b*',Mst1,Ypoly3,'b',Mst2,PDC2,'r*',Mst2,Ypoly10,'r');
grid on;
xlabel('M_{st}[Nm]');
ylabel('P_{DC}[W]');
title('Charakterystyka mechaniczna')
subtitle('P_{DC} = f(M_{st})');
legend('V_{TG}=28','','V_{TG}=33.5')

figure(5)
plot(Mst1,spr1,'b*',Mst1,Ypoly4,'b',Mst2,spr2,'r*',Mst2,Ypoly11,'r');
grid on;
xlabel('M_{st}[Nm]');
ylabel('η[%]');
title('Charakterystyka mechaniczna')
subtitle('η = f(M_{st})');
legend('V_{TG}=28','','V_{TG}=33.5')

figure(6)
plot(Mst1,cosfi1,'b*',Mst1,Ypoly5,'b',Mst2,cosfi2,'r*',Mst2,Ypoly12,'r');
grid on;
xlabel('M_{st}[Nm]');
ylabel('cosφ');
title('Charakterystyka mechaniczna')
subtitle('cosφ = f(M_{st})');
legend('V_{TG}=28','','V_{TG}=33.5')

figure(7)
plot(Mst1,w1,'b*',Mst1,Ypoly6,'b',Mst2,w2,'r*',Mst2,Ypoly13,'r');
grid on;
xlabel('M_{st}[Nm]');
ylabel('ω[rad/s]');
title('Charakterystyka mechaniczna')
subtitle('ω = f(M_{st})');
legend('V_{TG}=28','','V_{TG}=33.5')