Ejercicios en matlab

Solo disponible en BuenasTareas
  • Páginas : 9 (2224 palabras )
  • Descarga(s) : 0
  • Publicado : 3 de marzo de 2011
Leer documento completo
Vista previa del texto
CARLOS OMAR RAMOS L. 2007 EJERCICIOS DE MATLAB x=(-2*pi:pi/10:2*pi); y=2*sin(x); z=5*cos(1.5*x); suma=y+z; prod=y.*z; subplot(2,2,1), plot(x,y), grid on, axis tight axis([0 2*pi -2 2]) xlabel('x'), ylabel('2sen(x)'),title('función seno') subplot(2,2,2), plot(x,z), grid on, axis tight axis([0 2*pi -5 5]) xlabel('x'), ylabel('5cos(1.5x)'),title('función coseno') subplot(2,2,3), plot(x,suma), gridon, axis tight xlabel('x'), ylabel('2sen(x)+5cos(1.5x)'),title('función suma') subplot(2,2,4), plot(x,prod), grid on, axis tight xlabel('x'), ylabel('2sen(x)·5cos(1.5x)'),title('función producto') figure(2) plotyy(x,y,x,z), grid on, axis tight xlabel('x'), ylabel('y,z'),title('dos funciones'), figure(3) plot(x,y./(1+z)), grid on, axis tight >> x=0:1/500:2; vm=4*cos(2*pi*2*x); r=(5+vm);vc=cos(50*pi*x); modam=vc.*r; plot(x,modam) axis([0 2 -10 10]) xlabel('Tiempo') ylabel('Amplitud') title('AM con índice de mod. m = 0,8') % Este programa mostrará cómo obtener una onda cuadrada % a partir de los coeficientes de Fourier y luego su espectro de amplitud % s = 30; % valor del rango de los coeficientes N = 3; % número de periodos que se mostrarán n = -s:s; % rango de la sumatoria de loscoeficientes A = 8; % valor de la amplitud máxima B = -2; % valor de la amplitud mínima f = 1000; % valor de la frecuencia T = 1/f; % Valor del periodo tau = T/4; % ancho del pulso respecto del periodo d = tau/T; % ciclo útil Cn = (A-B)*sin(pi*n*d)./(pi*n); % fórmula para hallar el valor % del coeficiente n. Observar el punto . antes de dividir %que indica que es un escalar que se divide entre la matriz % %Matlab inicia el valor de los vectores indexándolos, comenzado por 1. % Para la fórmula indicada anteriormente, se generaría un error, de tal manera que % se tomará el valor de Co, que corresponde al índice s+1 Cn(s+1)=(A-B)*d + B; % valor de Co, observe cómo se agrega el offset % mediante la adición del valor mínimo B t = 0:T/100:N*T; % se plicarán 100 puntos por periodo en resolución y semuestran 3 periodos % % t es el vector de la variable tiempo vt = zeros(1,size(t)); %inicializa el vector con sólo ceros del tamaño del vector de tiempo t % Ahora se calcularán los puntos para conformar una onda cuadrada for k = 1:2*s+1 vt = Cn(k)*exp(i*2*pi*n(k)*f*t)+ vt; % Aquí se calculan los puntos que han de conformar end % la onda cuadrada (por eso se suma el punto anterior vt); en resumen es lasumatoria subplot(3,1,1), plot(t,vt); % se programa una ventana para tres gráficas, y se dispone grid % la primera, usando cuadrícula xlabel('tiempo') ylabel('Amplitud') 1

n=-6:20 cerosizq=zeros(1,6) m=1:18 coseno=cos(pi*m/8) cerosder=zeros(1,3) tabla3=[cerosizq,coseno,cerosder] stem(n,tabla3) axis([-6 20 -1.2 1.2]) title('cos(pi*n/8)(u[n]-u[n-17])') gris

CARLOS OMAR RAMOS L. 2007 % Ahorase mostrarán las magnitudes de los coeficientes % en el dominio de la frecuencia Amplitud = abs(Cn); Kn = n*f; subplot(3,1,2), stem(Kn,Amplitud);% se grafica el espectro de amplitud grid xlabel('frecuencia') ylabel('Amplitud') % Ahora se mostrará la potencia de los coeficientes % en el dominio de la frecuencia Potencia = Cn.^2; Rn = n*f; subplot(3,1,3), stem(Rn,Potencia);% se grafica la densidadespectral de potencia grid xlabel('frecuencia') ylabel('Potencia') t0=0:2; t1=2:3; t2=3:5; t3=5:6; f0=3/2*t0; f1=9-t1*3; f2=3*t2/2-9/2; f3=18-3*t3; plot(t0,f0,'b',t1,f1,'b',t2,f2,'b',t3,f3,'b') grid b=m%define una matriz de 1x5 llena de ceros tabla2=10*ones(1,11)%define una matriz de 1x9 llena de unos x=9:-1:1 c=x%define una matriz de 1x7 llena de ceros tabla3=[b,tabla2,c]%define una matriz con losdatos de las anteriores subplot(2,1,1) stem(n,tabla3,'r') grid dy=diff(tabla3)./diff(n) subplot(2,1,2) stem(dy) grid dy=sym('dy')

tabla1=ones(1,5) tabla2=zeros(1,4) tabla3=[tabla1,tabla2] stem(tabla3)

n=-5:10 x=6:10 tabla1=zeros(1,5) tabla2=ones(1,11) tabla3=[tabla1,tabla2] stem(n,tabla3) axis([-5 10 0 2])

n=1:20 %define el rango de las muestras Ln=length(n) %toma el valor de la...
tracking img