Laboratorio 1.docx

UNIVERSIDAD NACIONAL DE SAN AGUSTÍN FACULTAD DE PRODUCCIÓN Y SERVICIOS

ESCUELA PROFESIONAL DE INGENIERÍA EN TELECOMUNICACIONES

ASIGNATURA: PROCESAMIENTO DIGITAL DE SEÑALES

TEMA: INTRODUCCIÓN A SCILAB

ALUMNO: ROMERO ZAVALA EDWARD AARON

GRUPO: A DOCENTE: ING. RUBEN MATHEUS HERRERA

AREQUIPA – PERÚ

2019

Ejercicio 1 Construya la matriz MODELO y luego escriba los comandos que responden a las siguientes consignas: A=[1:20:181;2:20:182;3:20:183;4:20:184;5:20:185;6:20:186;7:20:187 ;8:20:188;9:20:189; 10:20:190;11:20:191;12:20:192;13:20:193;14:20:194;15:20:195;16:2 0:196;17:20:197;18:20:198;19:20:199;20:20:200]

Asignar a la variable a, el valor de la matriz MODELO que se encuentra en la 3er fila y 2da columna. Comando a=A(3,2) Respuesta

Asignarles al vector b la 3er fila de la matriz MODELO. Comando b=A(3,:) Respuesta

Asignarles al vector c la 2da columna de la matriz MODELO. Comando c=A(:,2) Respuesta

A partir de la matriz MODELO crear otra matriz D que contenga solo los elementos de filas pares y las columnas impares. Comando D=A(2:2:20,1:2:10) Respuesta

Construya una matriz E igual a la matriz MODELO pero donde se hayan intercambiado la 3ra y la última columna. Comando E=A; E(:,3)=A(:,10); E(:,10)=A(:,3)

Respuesta

Calcule la matriz F que tenga el valor: A’*A. ¿Es lo mismo que hacer A*A’? Comando F=A'*A F1=A*A' Respuesta de A'*A

Respuesta de A*A'

A’*A es diferente a A*A’

Elimine de la matriz MODELO la última fila Comando A(20,:)=[] Respuesta

Ejercicio 2 Graficar las siguientes funciones f ( x ) = x 2 - 1 En el dominio x = [-3.3, -1.3] x = -3.3:0.1:-1.3 y = sqrt((x.^2-1)) plot(x,y) Grafico

1 f ( x) = x 2 seno( ) x En el dominio x = [ -2, 2]

x=linspace(-2,2,100); y=(x.^2).*sin(1./x); plot(x,y) Gráfico

f ( x, y ) =

seno( x 2 + y 2 ) x 2 + y 2 + 0.1

En el dominio

x = [ -10,10] e y = [ -10,10]

x=linspace(-2,-1e-10,100); y=(x.^2).*sin(1./x); plot(x,y) x=linspace(1e-10,2,100); y=(x.^2).*sin(1./x); plot(x,y) plot(0,0,'o') Gráfico

Ejercicio 3 Crear una function de Scilab que realice el producto matricial entre dos matrices. Comparar el resultado con el producto * del software. function P=Pmatriz(A, B) [m q]=size(A); [q n]=size(B); for i=1:m for j=1:n suma=0; for p=1:q suma=suma + A(i,p)*B(p,j); end P(i,j)=suma; end end endfunction

Ejercicio 4 Utilizando la ecuación de Van der Waals:

(P +

a = 3.592 lt2 .atm/mol2

a )(V - b) = RT V2

b = 0.04267 lt/mol R = 0.082054 lt.atm/mol.ºK

a) Construir una function de Scilab tal que P=vanderwaals(V,T). b) Graficar la dependencia de la presión en función del volumen a las siguientes temperaturas: T1 = 100 °K T2 = 300 °K T3 = 600 °K y en todos los casos comparar la solución con la obtenida usando la ley de los Gases Ideales. Nota: Incluir todo en un mismo gráfico

function P=vanderwaals(v, T) R=0.082054; b = 0.04267; a=3.592; P = R*T/(v-b) - a/v^2; endfunction v=linspace(1,35,100); T1=100; T2=300; T3=600; for i=1:length(v); p1(i) = vanderwaals(v(i),T1); pideal1(i) = 0.082054*T1/(v(i)); p2(i) = vanderwaals(v(i),T2); pideal2(i) = 0.082054*T2/(v(i)); p3(i) = vanderwaals(v(i),T3); pideal3(i) = 0.082054*T3/(v(i)); end plot(v,p1,'r') plot(v,p2,'b') plot(v,p3,'y') plot(v,pideal1,'r--') plot(v,pideal2,'b--') plot(v,pideal3,'y--') xlabel('Volumen molar (l/mol)') ylabel('Presión (atm)') legend('Presión real a 100 K','Presión real a 300 K','Presión real a600 K','Presión ideal a 100 K','Presión ideal a 300 K','Presión ideal a 600 K')

Gráfico