Stałomomentowa regulacyjna

close all;

VD = [133; 121.6; 117.8; 114; 96.9; 85.5; 76; 62.7; 47.5; 39.9; 36.1];
w1 = [0; 8.93; 14.92; 18.62; 37.34; 48.74; 61.66; 74.58; 91.49; 98.8; 104.5];
P1 = [171; 171; 171; 209; 218.5; 228; 266; 266; 294.5; 313; 313];
P1M = [247; 228; 228; 228; 228; 228; 247; 266; 266; 275.5; 285];
Mst = [0; 0; 0.02; 0.03; 0.05; 0.05; 0.06; 0.06; 0.09; 0.1; 0.1];


[X,Y] = polyfit(VD,w1,3);
Ypoly = polyval(X,VD);

[X1,Y1] = polyfit(VD,P1,3);
Ypoly1 = polyval(X1,VD);

[X2,Y2] = polyfit(Mst,P1M,3);
Ypoly2 = polyval(X2,Mst);


figure(1)
plot(VD,w1,'b*',VD,Ypoly,'b')
grid on;
xlabel('U_{DC}[V]');
ylabel('ω[rad/s]');
title('Charakterystyka biegu jałowego')
subtitle('ω = f(U_{DC})');
legend('ω = f(U_{DC})','aproksymacja')

figure(2)
plot(VD,P1,'b*',VD,Ypoly1,'b')
grid on;
xlabel('U_{DC}[V]');
ylabel('P_{1}[A]');
title('Charakterystyka biegu jałowego')
subtitle('P_{1} = f(U_{DC})');
legend('P_{1} = f(U_{DC})','aproksymacja')

figure(3)
plot(Mst,P1M,'b*',Mst,Ypoly2,'b')
grid on;
xlabel('M_{st}[Nm]');
ylabel('P_{IM}[W]');
title('Charakterystyka biegu jałowego')
subtitle('P_{IM} = f(M_{st})');
legend('P_{IM} = f(M_{st})','aproksymacja')