%% Cviceni 13032026

clc
clear 
close all
% A) Vytvoreni signalu

f = 50; %Hz
omega = 2*pi*f; %rad*s^-1
delka_zaznamu = 1; %s
fvz = 2000; %Hz , Vzorkovaci frekvence
t = 0:1/fvz:delka_zaznamu;  %Delta casu
t_pul = delka_zaznamu/2 * fvz;

%%stav S1
Urms_S1 = [230; 230; 230];
fi_S1 = deg2rad([0; 120; -120]);

%%stav S2
Urms_S2 = [250; 200; 210];
fi_S2 = deg2rad([10; -110; 130]);

%%Funkce uS1 a uS2
u_S1 = sqrt(2).*Urms_S1 .* sin(omega* t(1:t_pul) + fi_S1);
u_S2 = sqrt(2) .* Urms_S2 .* sin(omega * t(t_pul+1:end) + fi_S2);

u = [u_S1 u_S2]; %%Sjednoceni 

figure
plot(t, u, '*-') %%Graf us1 a us2 v case t

%%B)

K = fvz/f;
Uc = zeros(size(u)); %%cos
Us = Uc;   %%sin

for m = 1:length(u) - K+1
    for k = m:m+K-1
        Uc(:, m) = Uc(:, m) + sqrt(2)/K * u(:, k) * cos(-k*2*pi/K);
        Us(:, m) = Us(:, m) + sqrt(2)/K * u(:, k) * sin(-k*2*pi/K);
    end
end

U = Uc + 1i * Us

U_abs = abs(U);
U_angle = angle(U);

figure
subplot(2, 1, 1)
plot(t, U_abs)
subplot(2, 1, 2)
plot(t, rad2deg(U_angle))

figure
compass(U(:, 1))

% Ulozeni dat
U_angle = rad2deg(U_angle);
U_save = [U_abs; U_angle];

save data_cvic_13032026.mat U_save