Построение функции с кусочными условиями

clc; close all; clear all;
t=linspace(0,0.5,1000);
G1 = [];
G2 = [];
fc = 50;
wc = 2*pi*fc;
a1=0.00375;
a2=0.006;
g1=0.00375;
g2=0.004;
for k=1:100;
    if (((2*k+1)*pi+a1+g1)/wc)<t<(((2*k+2)*pi+a1)/wc)
        G1 = 0;
    else G1 = 1;
    end
    hold on
    plot3(t,ones(1,length(t))*k, ones(1,length(t))*G1); grid on
    ylim([0 100])
    zlim([0 1])
    xlabel('time')
    ylabel('k(iteration)')
    zlabel('logical level')
    drawnow; pause(.1); view(140,50)
end

clear all, close all
t=linspace(0,0.5,1000);
g1 = zeros(1, length(t));
g2 = zeros(1, length(t));
fc=50;% частота сети
%k=10;
alp1=0.00375;
alp2=0.006;
gam1=0.00375;
gam2=0.004;
for k = 1:length(t);
    g1(t<((2*k*pi+alp1)/(2*pi*fc)))=0;
    g2(t<((2*k*pi+alp1)/(2*pi*fc)))=1;
    g1(((2*k*pi+alp1)/(2*pi*fc))<t&t<((2*k*pi+alp1+gam1)/(2*pi*fc)))=1;
    g2(((2*k*pi+alp1)/(2*pi*fc))<t&t<((2*k*pi+alp1+gam1)/(2*pi*fc)))=1;
    g1(((2*k*pi+alp1+gam1)/(2*pi*fc))<t&t<(((2*k+1)*pi+alp2)/(2*pi*fc)))=1;
    g2(((2*k*pi+alp1+gam1)/(2*pi*fc))<t&t<(((2*k+1)*pi+alp2)/(2*pi*fc)))=0;
    g1((((2*k+1)*pi+alp2)/(2*pi*fc))<t&t<(((2*k+1)*pi+alp2+gam2)/(2*pi*fc)))=1;
    g2((((2*k+1)*pi+alp2)/(2*pi*fc))<t&t<(((2*k+1)*pi+alp2+gam2)/(2*pi*fc)))=1;
    g1((((2*k+1)*pi+alp2+gam2)/(2*pi*fc))<t&t<(((2*k+2)*pi+alp2)/(2*pi*fc)))=0;
    g1(t>(((2*k+2)*pi+alp2)/(2*pi*fc)))=0;
    g2((((2*k+1)*pi+alp2+gam2)/(2*pi*fc))<t&t<(((2*k+2)*pi+alp2)/(2*pi*fc)))=1;
    g2(t>(((2*k+2)*pi+alp2)/(2*pi*fc)))=1;
    
    hold all
    subplot(211); plot(t,g1); ylim([-0.1 1.1]); grid on;
    hold all
    subplot(212); plot(t, g2); ylim([-0.1 1.1]); grid on;
end