clear;
clc;
%DANE DO ZADANIA DLA STUDENT O NUMERZE n = 2
f1 = 20;
A = 17;
phi = 0.2 * pi;
T = 1/f1;
gamma = 50;
%SYGNAŁ ANALOGOWY
t = 0:0.0001:T;
y = A * square(2*pi*f1*t+phi,gamma);
figure;
subplot(3,1,1);
plot(t,y)
title("y(t)")
xlabel("t")
ylabel("y")
%SYGNAŁ DYSKRETNY y1
N1 = 20;
k1=N/20;
t1=t(1:k1:N); %na podstawie funkcji alias.m
y1=y(1:k1:N);
subplot(3,1,2)
stem(k1,y1)
title("y1[k]")
xlabel("k")
ylabel("y1")
%SYGNAŁ DYSKRETNY y2
N2 = 100
k2=N/100;
t2=t(1:k2:N); %próbkowanie czasu 100 razy na okres
y2=y(1:k2:N); %próbkowanie wartości y 100 razy na okres
subplot(3,1,3)
stem(k2,y2)
title("y2[k]")
xlabel("k")
ylabel("y2")
%WIDMO AMPLITUDOWE I FAZOWE y1
n1 = length(y1)
cm1 = (1/n1)*fft(y1)
c_abs1 = abs(cm1);
phi1 = angle(cm1);
m1=1:1:n1;
figure;
subplot(2,1,1);
stem(m1,c_abs1)
title('Widmo amplitudowe dla y1')
xlabel('m [-]')
ylabel('c_m [-]')
subplot(2,1,2);
stem(m1,phi1)
title('Widmo fazowe dla y1')
xlabel('degree')
ylabel('phi')
%WIDMO AMPLITUDOWE I FAZOWE y2
n2 = length(y2);
cm2 = (1/n2)*fft(y2)
c_abs2 = abs(cm2);
phi2 = angle(cm2);
m2=1:1:n2;
figure;
subplot(2,1,1);
stem(m2,c_abs2)
title('Widmo amplitudowe dla y2')
xlabel('m [-]')
ylabel('c_m [-]')
subplot(2,1,2);
stem(m2,phi2)
title('Widmo fazowe dla y2')
xlabel('degree')
ylabel('phi')