matlab zad1 RC1

%PARAMETRY SINUSOIDY
A=1
f=0.2
phi=0

%PARAMTETRY SYGNAŁU
s_r=10 %SAMPLE RATE
t=10

x=0:1/s_r:t
y=A*sin(2*pi*f*x+phi)

%KWANTYZACJA
y1=uencode(y,4)

%WYKRES 1
%SYGNAŁ SINUSOIDALNY I JEGO KWANTYZACJA

figure
plot(x,y,'red')
hold on
h=stem(x,y1,'filled')
hold off
xlabel('t')
ylabel('sin(2\pi ft+\phi)')
title('Sygał sinusoidalny i jego kwantyzacja')

%ZAMIANA NA CIĄG BITÓW
b_mx=de2bi(y1)
b_mx1 = transpose(b_mx)
b_vec=reshape(b_mx1,1,numel(b_mx))

%Wykres 2
%Skwantowany sygnałi jego reprezentacja bitowa
figure
stairs(x,y1)
grid on
hold on
b_x=0:t/(length(b_vec)-1):t
stairs(b_x,7+b_vec,'LineWidth',2,'Marker','d','MarkerFaceColor','c')
  %HISTOGRAM BYŁBY CZYTELNIEJSZY
hold off

%GENEROWANIE ZAKŁÓCEŃ
n=29 %ZAKŁUCANY JEST CO N_TY BIT
dis=b_vec
for i=1:n:length(b_vec)
    % dis(i)=1
    dis(i)=~dis(i)
end

%ZAMIANA ZAKŁÓCONEGO CIAGU BITÓW NA SKWANTOWANY SYGNAŁ
TEMP=reshape(dis,4,numel(dis)/4)
TEMP2=transpose(TEMP)
dec=bi2de(TEMP2)
dec=transpose(dec)

%WYKRES 3
%ZAKŁOCANY SYGNAŁ
figure
stairs(x,dec)
grid on

%ZAMIANA CIAGU BITÓW NA SKWANTOWANY SYGNAŁ
D=reshape(b_vec,4,numel(b_vec)/4)
D2=transpose(D)
dec=bi2de(D2)
dec=transpose(dec)