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