lab4

Tn = 0;
Tv = 2*pi;
Tdelta = 0.01;
t=Tn:Tdelta:Tv;
x = 2+3*cos(2*t+5)-4*sin(8*t+10)+0.04*cos(0.5*t+10);
fun=@(t) 2+3*cos(2*t+5)-4*sin(8*t+10)+0.04*cos(0.5*t+10);
a0 = 1/Tv*integral(fun,0,2*pi);
A0 = a0;
w0 = 2*pi/Tv;
for n=1:1:21
    funa =@(t) (2+3*cos(2*t+5)-4*sin(8*t+10)+0.04*cos(0.5*t+10)).*cos((n)*t*w0); 
    funb =@(t) (2+3*cos(2*t+5)-4*sin(8*t+10)+0.04*cos(0.5*t+10)).*sin((n)*t*w0); 
    a(n) = 2/Tv.*integral(funa,0,2*pi);
    b(n) = 2/Tv.*integral(funb,0,2*pi);
    A(n) = sqrt(a(n).*a(n)+b(n).*b(n));
    w(n) = n*w0;
    fi(n) = atan(-b(n)/a(n));
    while abs(fi(n)) >= 2*pi
        fi(n) = abs(fi(n))-2*pi;
    end
end
figure(1);
stem(w,A,'k');
hold on;
grid on;
figure(2);
stem(w,fi,'k');
hold on;
grid on;
for n=1:1:21
    f(n) = A0+
end