Tarea 7 Matlab

 Universidad Autónoma de Nuevo León
Facultad de Ingeniería Mecánica y Eléctrica

Materia
Programación con Matlab
Tema
Tarea # 7
Catedrático
M.C. César González Cervántes

Alumno
Mario Alberto Aguilar Stackpole

Matricula: 1568092Salón: 7223
Grupo: 010 Hora: M5

Miércoles 30 de Abril de 2014.
Código para Matlab de página 98
Libro: Sistemas de Comunicación Digitales y Analógicos

M=7;
N=2^M;
n=0:1:N-1;
tfin=1;
T=10;
dt=T/N;
t=n*dt;
w=zeros(length(t),1); %X = zeros devuelve el escalar 0.
for(i=1:1:length (w)) %For- END_FOR es una declaración de repetición proporcionando un bucle
%de iteración automática en un intervalo de números u objetos.

if(t(i) <= tfin) %Declara y evalua la condición y actúa de la manera declarada si la condición se cumple

w(i)=1;
end;
end;
W=dt*fft(w);
f=n/T;
pos = index(f,4/tfin); %POS escribe o extrae la parte de la traducción deuna matriz homogénea
plot_pr(2);
plot(t,w)
axis([0 T 0 1.5]);
xlabel('t (sec) -->');
ylabel ('w(t)');
title('Forma de onda de tiempo');
pause;
subplot(311);
plot(f(1:pos),abs(W(1:pos)));
xlaberl('f (Hz) --->');
ylabel('|W(f)|');
title('ESPECTRO DE MAGNITUD hasta el cuarto nulo (NULL');
subplot(312);
plot(f(1:pos),180/pi*angle(W(1:pos)));
xlabel('f(Hz) --->');
ylabel ('thetha (f) (grados)');title('ESPECTRO DE FASE hasta el cuarto nulo (NULL)');
grid;
subplot(313);
plot(f,abs (W));
xlabel('f (Hz) --->');
ylabel ('|W(f)|');
title('ESPECTRO DE MAGNITUD sobre rango completo de frecuencias de la FFT');







Definición de comandos utilizados

M=7;
N=2^M;
n=0:1:N-1;
tfin=1;
T=10;
dt=T/N;
t=n*dt;
w=zeros(length(t),1); %X = zeros devuelve el escalar 0.

for (i=1:1:length (w)) %For- END_FORes una declaración de repetición proporcionando un bucle de iteración automática en un intervalo de números u objetos.

if(t(i) <= tfin) %Declara y evalua la condición y actúa de la manera declarada si la condición se cumple

w(i)=1;
end;
end;
W=dt*fft(w);
f=n/T;

pos = index(f,4/tfin); %POS escribe o extrae la parte de la traducción de una matriz homogénea

plot_pr(2);%plot_pr(m): Este comando informa al usuario de que se generarán"m" gráficas y que se debe presionar la barra espaciadora para pasar a la próxima gráfica.


plot(t,w) %Generar gráfica
axis([0 T 0 1.5]); %Definir el eje
xlabel('t (sec) -->'); %Nombrar el eje x
ylabel ('w(t)'); %Nombrar el eje y
title('Forma de onda de tiempo'); %Insertar un titulo
pause;
subplot(311); %Acomodar gráficas encuadrantes
plot(f(1:pos),abs(W(1:pos)));
xlabel('f (Hz) --->');
ylabel('|W(f)|');
title('ESPECTRO DE MAGNITUD hasta el cuarto nulo (NULL');
subplot(312);
plot(f(1:pos),180/pi*angle(W(1:pos)));
xlabel('f(Hz) --->');
ylabel ('thetha (f) (grados)');
title('ESPECTRO DE FASE hasta el cuarto nulo (NULL)');
grid on; %Generan un efecto de malla en el area de gráficas
subplot(313);
plot(f,abs (W));
xlabel('f(Hz) --->');
ylabel ('|W(f)|');
title('ESPECTRO DE MAGNITUD sobre rango completo de frecuencias de la FFT');


Código para Matlab de página 102
Libro: Sistemas de Comunicación Digitales y Analógicos

M=4;
N=2^M;
fo=10;
wo=2*pi*fo;
n=0:1:N-1;
T=1/fo;
dt=T/N;
t=n*dt;
w=3*sin(wo*t+(pi/180*20));
W=fft(w);
W=W(:);
n1=-N/2:1:N/2;
fn1=n1/T;
fs=1/dt;
cn=1/N*W;
Theta=(180/pi)*angle(cn+0.001);cn=fftshift(cn)';
Theta=fftshift(Theta)';
cn=[cn cn(1)];
Theta=[Theta Theta(1)];
cn=cn(:);
Theta=Theta(:);
plot(t,w);
xlabel('t(sec)-->');
ylabel('w(t)');
pause;
plot(n,abs(W),'o');
for i=1:1:length(n)
line([n(i) n(i)], [0 abs (W(i))]);
end;
xlabel('n');
ylabel('|W(n)|');
title('Puntos de la FFT');
axis([0 16 0 25])
pause;
subplot(211);
plot(fn1,abs(cn),'o')
for i=1:1:length(n1)
line([fn1(i)...