Bety

Caso a
clc
A=1;
d=1/20;
T=1/4;
n=-20:1:20;
x=(A*d/T)*sinc(n*d/T);
Stem(n,x)
Grid;
title('Espectro de frecuencia de la señal rectangular')
-------------------------------------------------

-------------------------------------------------

-------------------------------------------------

Caso b
Clc
A=1;
d=1/20;
T=1/2;
n=-40:1:40;
x=(A*d/T)*sinc(n*d/T);Stem(n,x)
grid;
title('Espectro de frecuencia de la señal rectangular')
Xlabel('W')

-------------------------------------------------

Cambiando de parámetros a T=1
clc
A=1;
d=1/20;
T=1;
n=-80:1:80;
x=(A*d/T)*sinc(n*d/T);
Stem(n,x)
grid;
title('Espectro de frecuencia de la señal rectangular')



-------------------------------------------------

clc;
N=128;a=10;
t=linspace(0,1,N)
x=exp(-a*t);
hold on
plot(t,x,'b+')
X=fft(x);
Z=ifft(X);
plot(t,real(z),'rx');
hold off
grid
title('GRAFICA COMPARATIVA ENTRE LA T DE F Y LA FUNCION EXACTA')
-------------------------------------------------

clc;
N=1024;
a=10;
T=linspace(0,1,N)
x=exp(-a*T);
X=(1/N).*fft(x,N);
Z=abs(X);
F=0:N;
w=2*pi*F;
A=1./sqrt(a*a+w.*w);fprintf('N\T DE F \t \valor exacto\n')
for(q=0:9)
fprintf('%2d \t%0.3f \t\t%0.3f \n',q,Z(q+1),A(q+1))
end

N 0 0.100 0.100
1 0.085 0.085
2 0.063 0.062
3 0.047 0.047
4 0.037 0.037
5 0.030 0.030
6 0.026 0.026
7 0.022 0.022
8 0.020 0.020
9 0.017 0.017
-------------------------------------------------
>>
clc;
N=16;a=10;
T=linspace(0,1,N)
x=exp(-a*T);
X=(1/N).*fft(x,N);
Z=abs(X);
fprintf('N\t fft\n')
for(q=0:N-1)
fprintf('%2d \t%0.3f \n',q,Z(q+1))
end

N fft
0 0.128
1 0.111
2 0.085
3 0.067
4 0.056
5 0.049
6 0.044
7 0.042
8 0.041
9 0.042
10 0.044
11 0.049
12 0.056
13 0.067
14 0.085
15 0.111-------------------------------------------------
>>

clc;
N=16;
t=0:1/N:(1-1/N);
a=10;
x=exp(-a*t);
X=(1/N).*fft(x,N);
m=[0:(N/2)-1];
Z=abs(X);
stem(m,Z(m+1));
grid
title('ESPECTRO POSITIVO DE LA T DE F')
-------------------------------------------------


clc,clf,clear all;
N=16;
t=0:1/N:(1-1/N);
a=10;
x=exp(-a*t);
X=fft(x,N);
z=(1/N)*fftshift(X);
Z=abs(z);m=[(-N/2)+1:1:(N/2)-1];
stem(m,Z(2:N),'o');
grid
title('ESPECTRO SIMETRICO DE LA T DE F')

-------------------------------------------------


clear all,clf,clc;
N=32;
T=1/N;
t=0:T:pi;
f1=5;
f2=10;
w1=2*pi*f1;
w2=2*pi*f2
s1=5*cos(w1*t);
s2=10*cos(w2*t);
x=s1+s2;
subplot(2,1,1)
plot(t,x);
title('GRAFICA DE LA SUMA DE DOS FUNCIONES COSENOIDALES EN EL TIEMPO ')
X=fft(x,N);z=(1/N)*fftshift(X);
f=(-N/2)+1:(N/2)-1;
subplot(2,1,2)
stem(f,real(z(2:N)));
title('ESPECTRO DE LA SUMA DE DOS ONDAS COSENOIDALES')
-------------------------------------------------

clear all,clf,clc;
N=32;
T=1/N;
t=0:T:pi;
f1=5;
f2=10;
w1=2*pi*f1;
w2=2*pi*f2
s1=2*sin(20*w1*t);
s2=5*sin(10*w2*t);
x=s1+s2;
subplot(2,1,1)
plot(t,x);title('GRAFICA DE LA SUMA DE DOS FUNCIONES senoidales EN EL TIEMPO ')
X=fft(x,N);
z=(1/N)*fftshift(X);
f=(-N/2)+1:(N/2)-1;
subplot(2,1,2)
stem(f,real(z(2:N)));
title('ESPECTRO DE LA SUMA DE DOS senoidales ONDAS ')



-------------------------------------------------

clear all,clf,clc;
N=32;
T=1/N;
t=0:T:pi;
f1=5;
f2=10;
w1=2*pi*f1;
w2=2*pi*f2s1=30*sin(10*w1*t);
s2=9*cos(2*w2*t);
x=s1+s2;
subplot(2,1,1)
plot(t,x);
title('GRAFICA DE LA SUMA DE DOS FUNCIONES EN EL TIEMPO ')
X=fft(x,N);
z=(1/N)*fftshift(X);
f=(-N/2)+1:(N/2)-1;
subplot(2,1,2)
stem(f,real(z(2:N)));
title('ESPECTRO DE LA SUMA DE DOS ONDAS ')




-------------------------------------------------





%programa para...