f=40;%Częstotliwość bo okres wynosi 0.025 sekundy, także 1/0.025=40
A=1;%Amplituda sygnału
CzasTrawnia=0.100;%Czas trwania sygnału, 100 milisekund
Fp=10000;%Częstotliwość próbkowania
Tp=1/Fp;%Okres próbkowania
t=0:Tp:CzasTrawnia;% Ciag liczb odpowiadajacych momentom probkowania
N=length(t);%Liczba próbwekk
y1=A*sin(2*pi*f*t); %Generowanie sygnalu sinusoidualnego
y2=A*sin(2*pi*f*t);

for k=1:length(y1); %Przekształcenie funkcji sinusoidalnej na protokątną
if (y1(k) >=0) y1(k)=1;
else y1(k)=-1;
end
end

for k=1:length(y2); %Przekształcenie funkcji sinusoidalnej na protokątną
if (y2(k) >=0) y2(k)=1;
else y2(k)=0;
end
end

subplot(211),plot(t,y1), title('Sygnał prostokątny unipolarny'), stem (t,y1),grid on, xlabel('Czas[s]'),ylabel('Amplituda'),axis([0 0.1 -1.2 1.2])
subplot(212),plot(t,y2), title('Sygnał prostokątny unipolarny'),stem (t,y1),grid on, xlabel('Czas[s]'),ylabel('Amplituda'),axis([0 0.1 0 1.2])