Trabajo Práctico No. 1
Comandos Matlab correspondientes a los Ejemplos de la Sección 3
- Ejemplo 1
t=0:0.0001:0.04;
x=sin(2*pi*50*t')+0.5*sin(2*pi*150*t');
td=0:1/600:0.04;
xn=sin(2*pi*50*td')+0.5*sin(2*pi*150*td');
plot(t,x)
hold on
stem(td,xn)
grid;
xlabel('Tiempo [seg]'); ylabel('x(t)');
- Ejemplo 2
Fs = 500; % frecuencia de muestreo en Hz
t1 = [0:1/(20*Fs):1]'; % vector tiempo "continuo"
F1 = 50; % frecuencia de la señal x1
F2 = F1 + Fs; % frecuencia de la señal x2
x1 = cos(2*pi*F1*t1); % señal x1
x2 = cos(2*pi*F2*t1); % señal x2
t = [0:1/Fs:1]'; % vector tiempo discreto
xm = cos(2*pi*F1*t); % señal muestreada
plot(t1,x1,'b-',t1,x2,'r-',t,xm,'go'); % gráfica
axis([0 0.02 -1 1]); % escalado en los ejes
xlabel('Tiempo [s]'); % rotulado del eje de abscisas
grid; % grilla
title('Señales x1(t), x2(t) y muestreadas'); % título
- Ejemplo 3
n=0:30;
x=exp(j*n/3);
subplot(211), stem(n,real(x));
xlabel('Indice "n"');title('Parte Real');
subplot(212), stem(n,imag(x));
xlabel('Indice "n"');title('Parte Imaginaria');
- Ejemplo 4
[x,y]=meshgrid(-8:0.5:8);
r=sqrt(x.^2+y.^2)+eps;
z=sin(r)./r;
mesh(x,y,z,ones(size(z))); % en colores: mesh(x,y,z)
xlabel('x'); ylabel('y'); zlabel('f(x,y)');
- Ejemplo 5
A = [1,2;3,4];
D = det(A);
if (A(1,1) & D) >0,
disp('A es definida positiva');
else
disp('A no es definida positiva');
end;
- Ejemplo 6
function F = Frobenius(A)
% Esta función calcula la norma de Frobenius de una matriz A
% La norma de Frobenius de una matriz (m x n) A=(a_ij) esta
% definida como:
% F^2= Sum_{i=1,m} Sum_{j=1,n} |a_ij|^2
F = sqrt(sum(sum(abs(A).^2)));
- Ejemplo 7
function [F,N1] = Norma(A)
% Esta función calcula la norma de Frobenius y la norma-1 de la
% matriz A.
% La norma de Frobenius de una matriz (m x n) A=(a_ij) esta
% definida como:
% F^2= Sum_{i=1,m} Sum_{j=1,n} |a_ij|^2
% La norma-1 esta definida como
% N1 = max_j sum_i |a_ij|
F = sqrt(sum(sum(abs(A).^2)));
N1 = max(sum(abs((A))));
Volver a TeSyS