clear all;
close all;
M1=[5	211
4.6	215
4.5	217
4.21	193
4	181
3.75	164
3.5	149
3.25	135
3	122
2.5	97
2	75
1.5	55
1	40
];
M2=[5.02	212
4.8	210
4.5	208
4.26	190
4	173
3.72	156
3.5	144
3.26	131
3.01	119
2.75	106
2.5	84
2	74
1.5	55
];
M3=[5	211
4.75	208
4.5	207
4.25	189
4	173
3.75	159
3.5	144
3.25	131
3	118
2.75	106
2.5	95
2.25	84
2	74
1.5	56
];
M4=[5	207
4.75	203
4.49	200
4.25	185
4	170
3.75	155
3.51	144
3.25	129
3.01	117
2.75	107
2.5	97
2.25	86
2	76
1.5	61
];
M5=[5	214
4.75	212
4.5	210
4.25	193
4	177
3.75	162
3.5	149
3.25	136
3	123
2.75	111
2.5	99
2.25	88
2	79
1.75	70
1.5	65
];

%Coefficients M1:
  M1_p1 = 48.991;
  M1_p2 = -18.714;
%Coefficients M2:
  M2_p1 = 49.574;
  M2_p2 = -27.685;
%Coefficients M3:
  M3_p1 = 48.812;
  M3_p2 = -24.225;
%Coefficients M4:
  M4_p1 = 45.888;
  M4_p2 = -15.729;
%Coefficients M5:
  M5_p1 = 47.6;
  M5_p2 = -15.5;
  
  p1=(M1_p1+M2_p1+M3_p1+M4_p1+M5_p1)/5
  p2=(M1_p2+M2_p2+M3_p2+M4_p2+M5_p2)/5
  U_steu=(1.5:0.05:4.5)';
  I=p1*M1(:,1)+p2;
  
%figure(1)
plot(M1(:,1),M1(:,2),'+k',M2(:,1),M2(:,2),'*k',M3(:,1),M3(:,2),'ok',M4(:,1),M4(:,2),'xk',M5(:,1),M5(:,2),'dk',M1(:,1),I,'k')
xlabel('U_{steu} [V]')
ylabel('I [mA]')
grid on
legend('p_{vor}= 14 mbar','p_{vor}= 140 mbar','p_{vor}= 250 mbar','p_{vor}= 370 mbar','p_{vor}= 496 mbar','approx. Kennlinie',4)
P=I*24/1000
figure(2)
[ax,h1,h2]=plotyy(M1(:,1),I,M1(:,1),P)
set(ax(1),'YColor','k'); 
axes(ax(1)); ylabel('I [mA]');
set(h1,'Color','k','LineStyle','-','LineWidth',1.25)
set(ax(2),'YColor','k');
axes(ax(2)); ylabel('P [W]');
set(h2,'Color','k','LineStyle',':','LineWidth',1.25);
xlabel('U_{steu} [V]')
%legend([h1 h2], 'Stromaufnahme', 'Leistungsaufnahme', 'Line1', 'Line2');
grid on
hold(ax(1),'on') 
hold(ax(2),'on')
plot(ax(1),M1(:,1),M1(:,2),'+k',M2(:,1),M2(:,2),'*k',M3(:,1),M3(:,2),'ok',M4(:,1),M4(:,2),'xk',M5(:,1),M5(:,2),'dk');
legend(ax(1),'Stromaufnahme','p_{vor}= 14 mbar','p_{vor}= 140 mbar','p_{vor}= 250 mbar','p_{vor}= 370 mbar','p_{vor}= 496 mbar')
legend(ax(2),'Leistungsaufnahme')
%hold(ax(1),'off') 
%hold(ax(2),'off')