EE354HW01

function HW1(problem_no)

    switch(problem_no)
        case(1)

            % 1.
            T = 2;
            f = 1/T;
            w0 = 2*pi*f; % w0 better equal pi
            Amp = 1.5;
            t = linspace(0,3*T,1+2^7)';
            % x = Amp*(2/pi)*asin(sin(w0*t));
            x = Amp*(2*(floor(t)/2==floor(floor(t)/2))-1);
            plot(t,x);

            A = [];
            % A0 = integral(x(t)*e^0),t,0,T);
            % A0 = int(x(t),t,0,T)
            % A0 = 0
            % A1 = int(x(t)*e^(-j*k*w0*t),t,0,T);
            % A1 = int(x(t)*e^(-j*w0*t),t,0,T)
            % A1 = int1 + int2
            % int1 = int(1.5*e^(-j*w0*t),t,0,1)
            % int2 = int(-1.5*e^(-j*w0*t),t,1,2)
            % int1 = 1.5*(e^(-j*w0*1)-e^0)
            % int1 = 1.5*(-1 - 1) = -3
            % int2 = -1.5*(e^(-j*w0*2)-e^(-j*w0*1))
            % int2 = -1.5*(1 - -1) = -3
            % A1 = -3 + -3
            % A1 = -6
            % An1 =
    case(2)
    case(3)


    % 3.

    % H(j*w) = (w<=50);

    T = pi/6;
    f = 1/T;
    w0 = 2*pi*f;

    % what are values kf = [k ; k+1] where A(kf(1)) = 1 & A(kf(2)) = 0
    % (lowpass)

    kf = [0;0];
    kf(2) = 0;
    kf(1) = w0;
    k = 0;
    while (kf(2) <= 50)
    k = k+1;
    kf(1) = kf(2);
    kf(2) = w0*k;
    end;

    otherwise
    disp('no problem selected');
    end
end