function [ txt ]=fazis_script(t,U,I,f) T=1/f; x=1; k=floor(t(length(t))/T); while t(x)<=k*T x=x+1; end P_t=U(1:x).*I(1:x); int_t=t(1:x); P=(1/(k*T))*trapz(int_t,P_t); U_cs=max(U(1:x)); I_cs=max(I(1:x)); cosfi=(P*2)/(U_cs*I_cs); fi=radtodeg(acos(cosfi)); P_str='P hatásos: '; disp(P_str); disp(P); fi_str='fi fáziskülönbség: '; disp(fi_str); disp(fi); ft_str='fázistényező (cos fi): '; disp(ft_str); disp(cosfi); txt=zeros(3,2); txt=num2cell(txt); txt{1,1}=P_str; txt{1,2}=num2str(P); txt{2,1}=fi_str; txt{2,2}=num2str(fi); txt{3,1}=ft_str; txt{3,2}=num2str(cosfi); figure(3); plot(t(1:x),P_t, 'g'); xlabel('Time (s)'), ylabel('P(t)') xlim([0 (4*T)]) grid on end