unidad 1-modo consola

Upload: leo

Post on 05-Jul-2018

220 views

Category:

Documents


0 download

TRANSCRIPT

  • 8/16/2019 Unidad 1-Modo Consola

    1/23

     

    MODO CONSOLA

    ASPECTOS GENERALES

    Modo consola

    En el modo básico de funcionamicientífica, escribiendo las expresiinmediata. Este es el llamado modo

    El modo consola de Scilab está spuede desplazar libremente por lapuede escribir a continuación del dde cursor derecha e izquierda permcorrecciones, mientras que las teclpara utilizarlas nuevamente tal cualya hecho.

    Nombre de variables

    En Scilab los nombres de las variab

    1- Pueden estar formados pocorrespondientes al alfabetNombres válidos: velocid 

    Nombres no validos: diáme

     2- Los nombres deben comen

    Nombres válidos: fuerza1  

    Nombres no válidos: 1pres 3- Scilab hace diferencia entr

    distintas.

    4- No se pueden utilizar espanombre radio 1, pero sí s

    Se recomienda no utilizar la diferede esta manera es preferible trab

    radio y Radio. Esto se debe a q

    producir resultados incorrectos de

    Operaciones y precedencia

    En Scilab se pueden realizar todas l

    OPotRadMulDiviAdiSus

     1 Como se explicará más adelante, el expon

    Ing. Juan E. Núñez Mc Leod

    nto de Scilab, se lo puede utilizar como si snes matemáticas y obteniendo los resultadoconsola.

    portado a través de un editor de línea, esto santalla, como por ejemplo cuando se utiliza unnominado “prompt” que está compuesto por loite el desplazamiento sobre la línea que está sis arriba y abajo permiten acceder a las líneasestán o previa modificación, algo muy útil cuan

    les deben responder a las siguientes reglas:

    r caracteres alfanuméricos (letras y números),inglés (no se pueden utilizar ñ, acentos, etc.).

    ad , altura.

    tro, año.

    ar con un carácter alfabético; es decir, deben c, giro2.

    ion, 2volumen.

    e mayúsculas y minúsculas, por lo que las va

    ios, pero sí es posible usar el guión bajo _. Ase puede usar radio_1.

    cia entre mayúsculas y minúsculas para repre jar con las variables r_menor  y r_mayor, qu

    e un error de escritura; por ejemplo, escribir ra

    uy difícil corrección.

    as operaciones matemáticas habituales, tales c

    peración Simbología Ejemplonciación ^ 2^3 = 8

    icación ^ 16^(1/4)=2

    iplicación * 3*4 = 12

    sión / 12/6 = 2

    ión + 4+6=10

    tracción - 8-5=3

    nte debe estar entre paréntesis, ya que la potenciación tien

    Página 1 de 23

    e tratara de una calculadoras de las mismas en forma

    ignifica que el usuario no serocesador de textos. Sólo secaracteres “-->”. Las teclas

    ndo escrita para poder hacerpreviamente escritas, ya seao se quiere repetir un cálculo

    pero las letras deben ser las

    menzar con una letra.

    riables radio  y Radio  son

    í por ejemplo, no es válido el

    entar dos variables distintas,e con las variables llamadas

    dio en lugar de Radio, va a

    mo:

    e prioridad sobre el cociente.

  • 8/16/2019 Unidad 1-Modo Consola

    2/23

     

    Las expresiones matemáticas se rPrimero se calculan los paréntesis,sumas y restas. Si hay operaciones

    Los paréntesis son utilizados paraevaluados. Así por ejemplo, la ex

    resultado 2 . Debido a que la radicaser colocado entre paréntesis, de

    calcula primero el exponente, que r

    Si no fuera por la presencia del pa

    luego el resultado se dividiría por 2  

    En un cálculo como -1^2 Scilab i

    prioridad que la potencia, dando

    paréntesis, escribiendo (-1)^2, en

    Como se verá más adelante, Scila

    calcular raíces de índice par y radicy teniendo en cuenta la discusión pr

    -->(-1)^(1/2)ans =

    6.123D-17 + i

    Se observa que el resultado obteprácticamente nula, comparada codel valor correcto.

    Estos errores de cálculo, que da

    computadora y en general no consque permiten realizar ciertos cálcultiene la función llamada sqrt, por l

    -->sqrt(-1)ans =

    i

    Obteniendo como resultado el valor

    Como conclusión se puede decir qque no solo se obtiene una mejoprocesador.

    Se puede mencionar como otro ejScilab incorpora una función espec

    puede escribir la siguiente instrucci -->exp(2)ans =

    7.3890561

    Se insiste en la necesidad de util

    incorpore, con el objeto de optimiza

    Ing. Juan E. Núñez Mc Leod

    suelven con las mismas prioridades utilizadasluego las potencias y raíces, luego los productode igual prioridad se resuelven de izquierda a d

      ambiar las prioridades de cálculo, ya que los mipresión 10/2+3   da como resultado  8 , mientr

    ión se resuelve como potencia de exponente frsta manera es calculado antes que la potenci

    esulta igual a 0.5 , y luego se evalúa la potencia

    éntesis, quedaría 9^1/2, y en tal caso se calc

    ando finalmente 4.5 .

    nterpreta el signo “-“ como símbolo de sustrac

    omo resultado final el valor -1. Para evitar e

    cuyo caso se obtiene el valor 1 como respuesta.

      b maneja tanto el álgebra como el análisis co

    ndos negativos. Así por ejemplo, si se desea cevia, se puede escribir:

    ido no es el esperado, pero analizándolo sela parte imaginaria, por lo que se puede consi

    n resultados aproximados y no los correctos,

    ituyen un problema. Sí se debe tener en cuent  s de manera más precisa; así por ejemplo, para

    que para obtener la raíz cuadrada de -1 se pu

    correcto i.

    e si Scilab tiene una función de cálculo específir precisión en los resultados, sino también se

    mplo el caso del cálculo de exponenciales delífica llamada exp. De esta manera si se desea

    n:

    izar las instrucciones específicas de cálculo,

    los tiempos de cálculo y la precisión de los mis

    Página 2 de 23

    en una calculadora científica.y cocientes, y finalmente las

    recha.

    smos son los primeros en seras que 10/(2+3)  da como

    ccionario, el exponente debe. Por ejemplo: 9^(1/2), se

    9^0.5 que resulta igual a 3 .

    la primero la potencia 9^1 y

    ión, por lo tanto tiene menor

    ste problema se deben usar

    plejo, por lo que es posible

    lcular la raíz cuadrada de -1,

    aprecia que la parte real eserar como una aproximación

    son comunes en cualquier

    que Scilab ofrece funcionescalcular raíces cuadradas sede escribir:

    co, ésta debe ser utilizada yaconsume menos tiempo de

    l tipo , para las cualescalcular el valor de  se

    n la medida que Scilab las

    os.

  • 8/16/2019 Unidad 1-Modo Consola

    3/23

     

    Constantes predefinidas

    Scilab define las siguientes constan

    Nombre2 

    %pi

    %e

    %i

    %inf

    %t

    %f

    %nan

    Así por ejemplo si se escribe %pi y

    -->%pi%pi =

    3.1415927

    La primera línea está compuesta p

    texto %pi  para luego presionar lasiguientes, que constituyen el result

    Tipos de datos 

    Los tipos de datos básicos que man

    - Datos numéricos: que pued- Cadenas de caracteres: qu

    comillas dobles o simples (- Datos booleanos: son todo

    Por ejemplo: ¿la variable x

    Notación científica 

    La representación de valores en nen ingeniería. Es de gran utilidadformato muy reducido. Así por ejem

     

    Descripción Not

    Masa del electrón

    Masa de la Tierra

    2 Es importante notar que todos los nombres

    3 Scilab no solo maneja el álgebra de númer

    calcular funciones con argumento complejo.4  El apóstrofe no se debe confundir con l

    inmediatamente a la derecha del 0 (cero) jun

    Ing. Juan E. Núñez Mc Leod

    es matemáticas:

    Significado

    Valor de la contante .

    Valor del número e, base de los logaritmos natur

    Valor de la unidad imaginaria  √ 1.

    Representa el infinito matemático ∞.

    Representa el valor lógico verdadero (true).

    Representa el valor lógico false (false).

    Indica que la expresión matemática no es calcul(n ot a  n umber)

    se presiona la tecla Enter, Scilab responde con l

    or el prompt -->, el cual no debe ser escrito, s

    tecla Enter, momento en el cual Scilab responado de la operación.

    eja Scilab son los siguientes:

    en ser tanto valores reales como complejos3.

    corresponden a todo tipo de textos, los cualespóstrofes)

    4, para que Scilab no las confunda co

    s aquellos datos que se pueden expresar en tés igual a 5?

    tación científica, también llamada notación detanto para representar números muy grandesplo:

    ación científica Notación de pu

    9.11 x 10-31

     kg 0.00000000000000000000

    5.9722 x 1024

     kg 597220000000

    de constantes de Scilab están precedidas por el signo %.

    o complejos, sino también el análisis matemático con dicho

    s acentos graves y agudos. El apóstrofe o comilla simplto con el signo de cierre de interrogación “?”.

    Página 3 de 23

    ales.

    ble

    a siguiente información:

    olamente se debe ingresar el

    de mostrando las dos líneas

    deben estar encerrados entrenombres de variables.

    rminos de verdadero  o falso .

    unto flotante, es muy comúncomo muy pequeños, en un

    to fijo

    000000000911 kg

    000000000000 kg

    números. De esta manera se pueden

    se ubica generalmente en la tecla

  • 8/16/2019 Unidad 1-Modo Consola

    4/23

     

    Las partes constitutivas de la notatiene el valor 1.58 x 10

    -5 sus distinta

    - Mantisa: 1.58

    - Base: 10- Exponente: -5

    Tal como se observa a continuación

     Si bien en ingeniería la base 10 espara representar números en notaci

    Para representar esta notación en ucual representa la base 10 y se pdicho sin dejar espacios intermedisigno, el cual puede ser omitido si e

    De esta manera en el ejemplo ant

    representa como 9.11d-31 y la m

    En lugar de utilizar notación científsiguiente manera:

    -->5.9722*10^24ans =

    5.972D+24 

    Si bien el resultado es el correctexpresión matemática para obtenerque este método es ineficiente y po

    Especificación del formato numé

    Por defecto Scilab escribe las cantutiliza un máximo de 10 posiciones

    Si esta forma de representación no

    cual debe mostrar los valores numé

    Esta función tiene dos argumentos:notación decimal, representado po

    por “e”) y el segundo indica la ca

    distintos ejemplos:

    5 Sólo es representado visualmente el signo

    6 Un lugar se reserva para el signo del núme

    Ing. Juan E. Núñez Mc Leod

    ión científica son: la mantisa, la base y el expos partes son:

    :

    el valor por defecto, en informática los microprón científica, aunque esto es totalmente transpa

    na computadora, se escribe primero la mantisa,ede escribir en mayúscula o minúscula) y fina

    os y cada elemento, ya sea mantisa o expons positivo.

    erior se debe escribir: 1.58d-5. Del mismo m

      sa de la Tierra es 5.9722d24.

    ica, se puede caer en la tentación de escribir

    , no es conveniente utilizar este método yael resultado, desperdiciando tiempo de cómputlo tanto no es recomendable.

    ico en pantalla

    idades numéricas en campos de 10 caracterespara representar dígitos, signos

    5 y punto decimal

      es adecuada, la función format permite indic

    ricos en pantalla.

    el primero indica si se debe utilizar notación de“v”) o punto flotante (en álgebra llamada not

    tidad máxima de dígitos a escribir6. En la sigui

    negativo. En caso del signo positivo Scilab deja un espacio e

    ro y otro para el punto decimal.

    Página 4 de 23

    nente. Así por ejemplo, si se

    ocesadores utilizan la base 2ente al usuario.

    luego la letra d o la letra e (lalmente el exponente, todo lonte, con su correspondiente

    odo, la masa del electrón se

    l valor correspondiente de la

    que Scilab debe resolver lao del microprocesador, por lo

    de tamaño; es decir, que sel.

    rle a Scilab el formato con el

    punto fijo (en álgebra llamadaación científica, representado

    nte tabla se pueden analizar

    n blanco. 

  • 8/16/2019 Unidad 1-Modo Consola

    5/23

     

    Instrformat(

    format(

    format(format(

     Se debe tener presente que la fun

    no a la precisión con que dichos vaprecisión disponible.

    Precisión numérica de una comp 

    En la matemática formal existen vconstituyen ejemplos clásicos de nno tienen una representación decim

    Las cantidades numéricas son almlo que es imposible almacenar estede dígitos, que por lo general soninstrucción format(“v”,25), la cadmitido por esta función), lo quedecimal y el signo.

    Si a continuación se escribe 1/3, el

    -->1/3ans =

    0.3333333333333333148

     

    Como se observa, los primero 1resultado artificial, al exigirle a la cintrínseca (la precisión del procesad

    Si se desea sumar 1 + 10-16

     el resul

    -->1+1d-16ans =

    1.

    Analizando el resultado se observaque le exige al microprocesador traserá:

    -->1+1d-15ans =

    1.00000000000000111022

    Se observa que si bien el resultadoprovenientes de las limitaciones nu Si bien, en este texto, se ha hechose advierte al lector que este temadeterminadas circunstancias se pumodelos matemáticos correctos y

    limitaciones tecnológicas de las co 

    Ing. Juan E. Núñez Mc Leod

    cción Cálculo Resultado en pantav”,8) 1/3 0.33333

    v”,5) -1/8 -0.12

    e”,10) 4^10 1.049D+06e”,15) -5.2^10 -1.44555106D+0

      ión format afecta a la representación de los v

    lores son almacenados, ya que Scilab los calcu

    tadora

    alores numéricos que tienen infinitos dígitos, tmeros irracionales y otros como 1/3 que si bienal finita, son los llamados números periódicos.

    cenadas en formato decimal y la memoria de l  tipo de números, por lo tanto las computadoras

    entre 14 y 16. Esto se puede comprobar fácilal solicita 25 espacios para la representación

    equivale a 23 dígitos un vez restados los luga

    esultado que se verá en pantalla será:

    96

    dígitos son correctos y después aparecenomputadora que muestre una cantidad de dígitor numérico del microprocesador).

    ado obtenido será el siguiente:

    que el cálculo está por fuera de la precisión nu  bajar con 17 dígitos. Por otro lado, si se hace el

    es correcto hasta el decimal 15, luego sigue uéricas del microprocesador.

    una discusión muy somera sobre la precisión nes mucho más amplio, y es algo que el ingeniden producir problemas de precisión numéricalos datos son correctos, pero los resultados

    putadoras.

    Página 5 de 23

    lla

    7

    alores numéricos en pantalla,

    la y almacena con la máxima

    l es el caso de √ 2,,  queson número racionales, pero

    s computadoras es finita, porlmacenan una cantidad finita

    mente en Scilab mediante ladecimal (es el máximo valor

    res reservados para el punto

    alores que responden a unos que exceden su precisión

    érica de la computadora, yacálculo 1 + 10

    -15, el resultado

    a serie de valores residuales

    mérica en una computadora,ero debe conocer, ya que enraves; es decir, se tienen loson incorrectos debido a las

  • 8/16/2019 Unidad 1-Modo Consola

    6/23

     

    VECTORES Y MATRICES 

    La información numérica manejadaparticular de estas últimas los vecto

    Declaración

    La primera operación que se debeecuaciones lineales como el que se

    Desde el punto de vista matricial, etérminos independientes B y la mat

    En donde cada matriz está definida

      La matriz A se puede definir en Scil

    1- Escribir A=[2 -3 2;-2 6

    2- Escribir  A=[2,-3,2;-2,6,

     3- Escribir

     A=[2 -3 2  presio

    -2 6 -2  presiona

    1 -1 -1]  y presio

    En todos los casos la matriz comieprimeros casos se observa como ucada fila por espacios (primer caso)filas de la matriz. En el tercer casescribir el punto y coma. En este

    terminó el ingreso de la matriz.

    Sin importar cuál de los métodos semuestra el resultado para poder revi

     A =2. - 3. 2.

    - 2. 6. - 2.1. - 1. - 1.

    Si el resultado no se correspondeescribiendo todo de nuevo o utilizan

    7 Se dan tres posibles formas, el lector debe

    Ing. Juan E. Núñez Mc Leod

    por Scilab es almacenada en forma de escalarres fila y columna.

    prender es cómo definir o ingresar matrices. Sda a continuación:

    2 3 2 42 6 2 10 11 ste sistema está definido por la matriz de los c

      iz de las incógnitas X, tal cual se observa a conti

     .  de la siguiente manera:

    2 3 22 6 21 1 1 ; 41011 ;

    ;ab de acuerdo a las instrucciones dadas a contin

     -2;1 -1 -1] y presionar la tecla Enter.

    -2;1,-1,-1] y presionar la tecla Enter.

    ar la tecla Enter y luego escribir,

    r Enter y terminar escribiendo,

    nar Enter nuevamente.

    za con un corchete izquierdo y termina con unna matriz se puede ingresar en una sola línea,o con comas (segundo caso). El punto y coma s, cada fila de la matriz se escribe en una nuevúltimo caso cuando se escribe el corchete d

      utilice, en el momento de terminar el ingreso yisar que todo esté en orden, tal como se muestr

    on el esperado, se debe proceder a ingresardo las flechas del cursor arriba y abajo para corr

    probarlas y adoptar la que le parezca más sencilla o con la

    Página 6 de 23

    s y matrices, siendo un caso

    a por ejemplo, un sistema de

    eficientes A, la matriz de losinuación:

    uación7:

    corchete derecho. En los dosseparando los elementos deirve para separar las distintasa línea, acá no es necesariorecho, Scilab interpreta que

    resionar la tecla Enter, Scilaba continuación:

    uevamente la matriz, ya seaegir el error cometido.

    ue se sienta más cómodo.

  • 8/16/2019 Unidad 1-Modo Consola

    7/23

     

    Siguiendo el primer método, para

    definido un vector columna (tres fil

    fila (una fila y tres columnas). Desddebe utilizar la que corresponda en

    Para resolver el sistema de ecuaciinversa de A por la matriz B. Parala matriz inversa del argumento, el

    -->X=inv(A)*B

    X =

    - 4.2.5.

    Además de lo visto, hay tres formas

    1- Indicando el valor inicial, elvayan de 0 a 10, yendo de

    resultado que se obtiene es

    -->x=[0:0.5:10]x =

    column 1 t0. 0.5 1.

    column 12 t5.5 6. 6.5

    La forma particular en que

    ventana de la consola decolumna 1 a la 11 y a contin

    2- La otra posibilidad es indiccaso se debe utilizar unatermine en 36, pero que co

    es resultado es el siguiente:

    --> x=linspace(0,36,x =

    0. 4. 8.

    Estos valores constituyen ucontiguos cualesquiera es c

     3- La última opción es similar

    potencias de base 10, de eejemplo:

    -->x=logspace(2,6,5)x =

    100. 1000.

    Para observar mejor el ressiguientes instrucciones:

    Ing. Juan E. Núñez Mc Leod

    ingresar la matriz B se escribe B. Obsérvese

    s y una columna), pero si se escribiera B=[-4

    el punto de vista del álgebra de matrices una ecada caso; es decir, no son intercambiables.

    ones se debe hacer el cálculo   . , queacer esto en Scilab se escribe X=inv(A)*B, d

    ual debe ser una matriz cuadrada. El resultado

    adicionales para definir vectores fila:

    l incremento y el valor final. Así si se desea d0.5 en 0.5, se puede escribir x=[0:0.5:10]

    el siguiente:

    111.5 2. 2.5 3. 3.5

    217. 7.5 8. 8.5 9.

    se muestra la información en este caso se de

    cilab para mostrar la fila completa, por lo queuación de la 12 a la 21.

    r el valor inicial, el valor final y la cantidad de v  función específica. Por ejemplo, se desea un

    ntenga 10 valores, para lo cual se debe escribi

    :

    0)

    12. 16. 20. 24. 28.

    na distribución que se denomina lineal, ya que lonstante.

    a la anterior, sin embargo los valores iniciales yta manera se logra una distribución logarítmica,

    10000. 100000. 1000000. 

    ultado, se expresarán los valores en notación c

    Página 7 de 23

    que de esta manera queda

    10 -11] se define un vector

    s la traspuesta de la otra y se

    significa multiplicar la matriznde la función inv() calcula

    btenido es el siguiente:

    finir un vector cuyos valoresy presionar la tecla Enter. El

    . 4.5 5.

    .5 10.

    e a la falta de espacio en la

    se observa primero desde la

    lores deseados, pero en estevector que comience en 0 yr x=linspace(0,36,10), y

    32. 36. 

    diferencia entre dos valores

    inales son los exponentes detal cual se ve en el siguiente

    ientífica, según muestran las

  • 8/16/2019 Unidad 1-Modo Consola

    8/23

     

    -->format("e",9)

    -->x=logspace(2,6,5)

    x =

    1.00D+02 1.00

    Esta distribución se denoconsecutivos cualesquiera

    Uso de las distribuciones lineale

    Como se discutió anteriormente,linealmente, ya que la diferencia eserie de valores dicha diferencia es

    0, 10, 20, 30, 40, ...

    Sin embargo, en ingeniería es mutrabajar con concentraciones de hidamplio de valores; así en el casohabiendo una diferencia de 14 órdede Reynolds el mismo puede ir desen algunos casos pueden superar ese habla de un conjunto de valoresvalores consecutivos es constante.

    Sea la siguiente serie de valores:

    10

    2

    , 10

    4

    , 10

    6

    , 10

    8

    , ...

    Los logaritmos decimales de cada u

    2, 4, 6, 8, ...

    Siendo la diferencia entre dos vallogarítmico.

    En un curso de análisis matemáticnormal graficar las mismas con esvisualmente los extremos relativosituación es distinta, ya que se des

    de la variable dependiente conocien

    Para comprender la necesidad de u

    Con el objeto de calcular la caídafricción denominado Fanning, el culisos.

    En las siguientes figuras se ha reprmillón. En la primera imagen el ejede ordenadas es lineal.

    8 El orden de magnitud es el valor del expon

    difieren en tres órdenes de magnitud. 

    Ing. Juan E. Núñez Mc Leod

    +03 1.00D+04 1.00D+05 1.00

      ina logarítmica ya que la diferencia entre loe mantiene constante.

    y logarítmicas

    la función linspace  permite obtener un conjntre dos valores consecutivos es constante, asde 10:

    común la necesidad de utilizar espaciados lo  rogeniones, números de Reynolds, etc. En estos

      e hidrogeniones se puede ir de concentracionnes de magnitud entre los dos valores

    8, mientra

    e valores muy bajos como algunas decenas y lll millón, habiendo en este caso 5 órdenes de mcon espaciado logarítmico significa que la difere

    no de dichos valores son:

    res consecutivos igual a 2 en todos los casos

    o, donde se estudian las características y propcalas lineales, ya que se busca representar sós, periodicidad, pendiente, positividad, negatian graficar expresiones funcionales, con el obje

    do el valor de la independiente.

    ilizar ejes logarítmicos se verá el siguiente ejem

      de presión en una tubería, es necesario conocal depende del número de Reynolds exclusiva

    sentado el valor del Fanning vs Re para un Ree abscisas es lineal y en la segunda es logarít

    nte al representar la cantidad en notación científica; así por

    Página 8 de 23

    +06

    s logaritmos de dos valores

    unto de valores espaciadosí por ejemplo en la siguiente

    garítmicos, por ejemplo paracasos se tiene un rango muy

    es de 10-14

      hasta 1 mol/litro,que para el caso del número

    egar a valores muy altos, quegnitud de diferencia. Cuandoncia de los logaritmos de dos

    , lo que indica un espaciado

    iedades de las funciones, eslo una porción, para analizaridad, etc. En ingeniería, lao de permitir conocer el valor

    plo:

    er el valor del coeficiente demente para el caso de tubos

    nolds que va desde 100 a unico, y en ambos casos el eje

    ejemplo, el valor 2 y el valor 2000

  • 8/16/2019 Unidad 1-Modo Consola

    9/23

     

    Fig.1 . Representación d

    Fig. 2. Representación del Fanning vs Rey

    Ing. Juan E. Núñez Mc Leod

    l Fanning vs Reynolds para tubos lisos con escalas lineales

    nolds para tubos lisos con escala logarítmica en el eje de ab

    Página 9 de 23

    en ambos ejes.

    cisas y lineal en el de ordenadas.

  • 8/16/2019 Unidad 1-Modo Consola

    10/23

     

    Si se desea conocer el valor del Faalrededor de 0.0025 sin importar ellineal da un valor que puede ir desleer un valor de Fanning cercano

    imposible determinar el Fanning cdando un valor alrededor de 0.014.

    Otra forma de ver la bondad del esde hidrogeniones desde 10

    -14 hasta

    Scilab selecciona automáticamente

    -->x=linspace(1d-14,1,15)x =

    column 1 to 5

    1.000D-14 0.0714286

    column 6 to 110.3571429 0.4285714

    column 12 to 150.7857143 0.8571429

    Como se observa el resultado no ren su lugar, la mayoría de los valor

    Utilizando la función logspace, sevalores más homogénea en el rang

    -->x=logspace(-14,0,15)

    x =column 1 to 5

    1.000D-14 1.000D-13

    column 6 to 111.000D-09 1.000D-08

    column 12 to 150.001 0.01 0.1

    Analizando el resultado es clara lalos argumentos de logspace son la

    Evitar el eco a pantalla

    Como se ha visto, cada vez que spantalla el resultado correspondientresultados a medida que se va catedioso: imagínese el caso de gene

    En caso que no se desee que seinstrucción correspondiente, comoinstrucción, no envía información a l

    -->x=linspace(1,50,1000);

    -->

    9 La incertidumbre se debe a la elevada pen

    10 Teniendo en cuenta que se representan di

    Ing. Juan E. Núñez Mc Leod

    ning para un Reynolds de un millón, no es difícigráfico utilizado. Sin embargo, si el valor del Ree 0.008 hasta 0.0114

    9, mientras que en la repr

    0.008. Siguiendo con la lectura de valores, p

    n la escala lineal; sin embargo, la lectura es f

    aciado logarítmico es para el caso de querer re1 mol/l. Al utilizar un espaciado lineal el resultala notación científica al tratarse de números muy

    0.1428571 0.2142857 0.2857

    0.5 0.5714286 0.6428571

    0.9285714 1.

    presenta lo que se busca, una distribución hos (en términos relativos) se aproximan a 1.

    logra un espaciado logarítmico el cual se tr  de concentraciones requerida.

    1.000D-12 1.000D-11 1.000D

      1.000D-07 0.000001 0.00001

    1.

    distribución pareja de valores. Para lograr estopotencia inicial, la potencia final

    10 y la cantidad d

    e escribe alguna instrucción y se presiona la te, esto es lo que se llama un eco. Este fenómenlculando; sin embargo, en otros casos puedear un vector con la función linspace de 1000

    roduzca un eco en pantalla, se debe escribir ue ve a continuación, donde una vez que se pra pantalla y muestra nuevamente el prompt para

    iente de la curva en las cercanías del origen de coordenada

      ichos valores en notación científica.

    Página 10 de 23

    l encontrar que el mismo valenolds es de 10000, el gráficosentación logarítmica es fácilra un Reynolds de 1000, es

    ácil en la escala logarítmica,

    resentar las concentracioneso es el siguiente (nótese quepequeños o muy grandes):

    43

    0.7142857

    ogénea desde 10-14

     hasta 1,

    duce en una distribución de

    10

    0.0001

    se debe tener en cuenta quee valores.

    cla Enter, Scilab muestra enes muy útil para conocer los

    esultar innecesario o inclusolementos.

    n punto y coma al final de lasiona Enter Scilab ejecuta lala siguiente instrucción.

    s.

  • 8/16/2019 Unidad 1-Modo Consola

    11/23

     

    Si aparece el prompt únicamente,éxito, ya que de lo contrario Scilabcálculos.

    -->x=0;

    -->y=1/x;!--error 27

    División por cero...-->

    Acceso a elementos individuales

    Sean las matrices A y B declaradascolumna de 3 x 1, tal cual se reprod

    Para acceder a los elementos indivi

    fila y j la columna del elemento descolocar el número uno.

    Así por ejemplo, para conocer el ele

    -->A(1,1)ans =

    2.

    Mientras que si se desean trabajar

    -->B(2,1)ans =

    10.

    En el caso particular de vectores filescribir únicamente la dimensión qescribir:

    -->B(2)ans =

    10.

    Acceso a submatrices 

    Se verá a continuación la forma de

    Se puede acceder a una columna c

    como resultado la segunda column

    Ing. Juan E. Núñez Mc Leod

    in ninguna información adicional, significa que lmostrará un mensaje de error, como se observa

    de una matriz 

    anteriormente, en donde A es una matriz cuadruce a continuación:

      2 3 22 6 21 1 1 ; 41011duales de una matriz se escribe por ejemplo A(

      ado, obsérvese que en el caso de la matriz B,

    mento de la primera fila y la primera columna, s

    on el elemento de la segunda fila de la matriz B,

    a o columna, al ser una de las dimensiones fija,e cambia, sin importar si es la fila o la columna,

    cceder a varios elementos a la vez. Sea la matr

      2 3 22 6 21 1 1 

    ompleta de la misma mediante una instrucción

    completa:

    Página 11 de 23

    operación fue realizada conen la siguiente secuencia de

    ada de 3 x 3 y B es un vector

    ,j) o B(i,1), donde i es laomo columna sólo es posible

    escribe:

    se escribe:

    Scilab nos permite obviarla yen el caso anterior se podría

    iz A definida anteriormente:

    omo la siguiente, obteniendo

  • 8/16/2019 Unidad 1-Modo Consola

    12/23

     

    -->A(:,2)ans =

    - 3.6.

    - 1.

    De la misma manera se puede obte

    -->A(3,:)ans =

    1. - 1. - 1.

    Finalmente se podría escribir laelementos ubicados en la primera y

    -->A(1:2,2:3)ans =

    - 3. 2.6. - 2.

    Matrices características

    Existe una serie de matrices queanteriormente, se dispone de funcio

    Si se desea obtener una matriz uni

    siguiente:

     A =1. 0. 0.0. 1. 0.0. 0. 1.0. 0. 0.

    Si se desea obtener una matriz nul

    siguiente:

     A =0. 0. 0.0. 0. 0.0. 0. 0.

    0. 0. 0.

    Si se desea obtener una matriz cu

    obteniendo como resultado lo siguie 

     A =1. 1. 1.

    1. 1. 1.1. 1. 1.1. 1. 1.

    Existen otras funciones para realiza

    Ing. Juan E. Núñez Mc Leod

    ner la tercera fila completa escribiendo la siguien

    siguiente instrucción, que nos devuelve la ssegunda filas, y en la segunda y tercera column

      son características, que si bien se pueden dnes específicas para facilitar esta tarea.

    taria de 4 x 4 se puede escribir A=eye(4,4),

    0.0.0.1.

    la de 4 x 4 se puede escribir A=zeros(4,4),

    0.0.0.

    0.

    yos elementos sean todos uno, de 4 x 4, se p

    nte:

    1.

    1.1.1.

    r operaciones con matrices que se detallan en la

    Página 12 de 23

    te instrucción:

    ubmatriz compuesta por loss.

    finir con las técnicas vistas

    bteniendo como resultado lo

    bteniendo como resultado lo

    uede escribir A=ones(4,4),

    siguiente tabla.

  • 8/16/2019 Unidad 1-Modo Consola

    13/23

     

    Función

     A’11  Devuel

    diag(A) Devuel

    tril(A) Devuel

    triu(A) Devuel

     prod(A)Devueleleme

    sum(A)Devueleleme

    det(A) Calcul

    inv(A) Calcul

    norm(x) Calcul

    size(A) Devuel

    length(A) Devuel

    FUNCIONES MATEMÁTICAS

    En la siguiente tabla se presenta ulas características de cada una deargumentos de las funciones puedScilab diferencia entre mayúsculaexpresión correcta es sin(x), por

    Funciónsin(x) Seno d

    sind(x) Seno d

    sinh(x) Seno hi

    asin(x) Arcose

    asind(x) Arcose

    asinh(x) Funció

    cos(x) Cosen

    cosd(x) Cosen

    cosh(x) Cosen

    acos(x) Arcoco

    acosd(x) Arcocoacosh(x) Funció

    tan(x) Tangen

    tand(x) Tangen

    tanh(x) Tangen

    atan(x) Arcotan

    atand(x) Arcotan

    atanh(x) Funció

    sec(x) Secant

    secd(x) Secant

    sech(x) Secant

    11 El símbolo a utilizar es el apóstrofe, no se

    Ing. Juan E. Núñez Mc Leod

    Descripción

    ve la matriz traspuesta de A.

    ve un vector con la diagonal principal de la matri

    ve la matriz triangular inferior de A.

    ve la matriz triangular superior de A.

    ve un valor escalar que equivale al producto de ttos de la matriz A.ve un valor escalar que equivale a la suma de totos de la matriz A.

    el valor del determinante de la matriz A.

    la matriz inversa de A.

    el módulo del vector x (norma euclideana).

    ve un vector con la cantidad de filas y columnas

    ve la cantidad total de elementos de la matriz A.

    a lista de las funciones matemáticas más comellas. Se aclara que no es una lista exhaustiv

    en ser matrices, no solo escalares. También siy minúsculas, por lo que si se desea calcu

    l contrario Sin(x) o SIN(x), dará el error “Var

    Descripciónl ángulo x expresado en radianes.

    l ángulo x expresado en grados sexagesimales.

      perbólico de x.

    o de x, el resultado estará en radianes.

    o de x, el resultado estará en grados sexagesi

     inversa del seno hiperbólico de x.

    del ángulo x expresado en radianes.

    del ángulo x expresado en grados sexagesimal

     hiperbólico de x.

    eno de x, el resultado estará en radianes.

    eno de x, el resultado estará en grados sexage inversa del coseno hiperbólico de x.

    te del ángulo x expresado en radianes.

    te del ángulo x expresado en grados sexagesim

    te hiperbólica de x.

    gente de x, el resultado estará en radianes.

    gente de x, el resultado estará en grados sexag

     inversa de la tangente hiperbólica de x.

    del ángulo x expresado en radianes.

    del ángulo x expresado en grados sexagesimal

     hiperbólica de x.

    debe confundir con los acentos graves o agudos.

    Página 13 de 23

    z A.

    odos los

    dos los

    de la matriz A.

    nes, haciéndose mención dey se debe recordar que los

    empre se debe recordar quelar el seno de un ángulo laiable indefinida”.

    ales.

    s.

    imales.

    les.

    simales.

    es.

  • 8/16/2019 Unidad 1-Modo Consola

    14/23

     

    Funciónasec(x) Arcose

    asecd(x) Arcose

    asech(x) Funció

    csc(x) Coseca

    cscd(x) Coseca

    csch(x) Coseca

    acsc(x) Arcoco

    acscd(x) Arcoco

    acsch(x) Funció

    cotg(x) Cotang

    cotd(x) Cotang

    coth(x) Cotang

    acot(x) Arcocot

    acotd(x) Arcocotacoth(x) Funció

    log(x) Logarit

    log10(x) Logarit

    exp(x) Calcula

    sqrt(x) Calcula

    abs(x) Calcula

    FUNCIONES DE MANEJO DE NÚ

    Se detalla a continuación un conju

    actual curso no se estima utilizar elos números complejos, junto con s

    Función

    complex(x,y)Cresernú

    real(z) Devimag(z) Dev

    isreal(z)Devlo c

    conj(z) Dev

    Es importante mencionar que todasde los números complejos. De estalogaritmos negativos, etc. Por otroejemplo se pueden ver los siguiente

    -->sqrt(-4)ans =

    2.i

    -->log(-1)ans =

    3.1415927i

    -->z=1+%i

    Ing. Juan E. Núñez Mc Leod

    Descripciónante de x, el resultado estará en radianes.

    ante de x, el resultado estará en grados sexage

     inversa de la secante hiperbólica de x.

    nte del ángulo x expresado en radianes.

    nte del ángulo x expresado en grados sexagesi

    nte hiperbólica de x.

    ecante de x, el resultado estará en radianes.

    ecante de x, el resultado estará en grados sexa

     inversa de la cosecante hiperbólica de x.

    ente del ángulo x expresado en radianes.

    ente del ángulo x expresado en grados sexagesi

    ente hiperbólica de x.

    angente de x, el resultado estará en radianes.

    angente de x, el resultado estará en grados sex inversa de la cotangente hiperbólica de x.

    o natural de x.

    o en base 10 de x.

    el antilogaritmo natural de x.

    la raíz cuadrada de x.

    el valor absoluto de x.

    EROS COMPLEJOS

    to de funciones específicas para aplicar a núm

    ste conjunto numérico; sin embargo, es necesas funciones, tienen gran aplicación en el campo

    Descripcióna el número complejo x + i y, donde los argummatrices, tienen que tener las mismas dimeeros reales.uelve la parte real del número complejo z.uelve la parte imaginaria del número complejo z.uelve verdadero (%t) si el argumento z es un nontrario devuelve falso (%f).uelve el complejo conjugado de z

    las funciones matemáticas mencionadas anterimanera se pueden calcular raíces de índiceslado estas funciones también admiten argum

    s cálculos:

    Página 14 de 23

    simales.

    ales.

    esimales.

    males.

    gesimales.

    eros complejos. Si bien en el

    rio indicar que el conjunto dede la ingeniería.

    ntos x e y, densiones y ser

    úmero real, de

    rmente trabajan en el campoares y radicandos negativos,ntos complejos. A modo de

  • 8/16/2019 Unidad 1-Modo Consola

    15/23

     

    z =1. + i

    -->sin(z)ans =

    1.2984576 + 0.6349639i

    POLINOMIOS

    Declaración

    Desde el punto de vista algebraico,sus raíces, de esta manera uncontinuación:

    !" " 

    Siendo !# los coeficientes del polinutilizar la forma factorizada, la misen el caso del ejemplo es !".

    Sea por ejemplo el siguiente polino 

    $% Para definir este polinomio en Scila 

    1- Haciendo uso de los coefici

     p=poly([6 1 -4 1],”x

     Nótese que los coeficientecomenzando por el coeficitercer orden, se escribe la lfinalmente la expresión “c

    polinomio. Una vez presicomprobar que el resultado

     p =2 3

    6 + x - 4x + x

    2- Haciendo uso de las raíces

     p=poly([-1 2 3],”x”,

     En este caso se indica con

    raíces del polinomio. Scilab

     p =2 3

    6 + x - 4x + x

    Como se observa en ambos casospolinomios que sean necesarios, y

    12 Se puede utilizar el apóstrofe en lugar de l

    Ing. Juan E. Núñez Mc Leod

    un polinomio puede ser expresado ya sea a trolinomio de tercer grado tiene las dos repr

    !  !  !&  !"% '(% '(% '"( 

    omio y '#  las raíces del mismo. Es importantea debe estar multiplicada por el coeficiente del

    io:

    (  " 4 6  % 1(% 2(% 3( 

    , se puede recurrir a cualquiera de los métodos

    entes se puede escribir la instrucción:

    ,”coef”)12 

    se escriben como un vector (los valores estánnte independiente y avanzando en orden crec

    etra que representará a la variable del polinomief”, indicando que los datos ingresados corres

    nada la tecla Enter, Scilab responde con laconcuerde con el polinomio original:

    se puede escribir la instrucción:

    roots”) 

    la expresión “roots” que la información dada

    responde con el siguiente polinomio:

    se obtiene el mismo polinomio. De esta maneracada uno se le asigna un nombre para poder di

      as comillas, en este caso se podría escribir también p=pol

    Página 15 de 23

    vés de sus coeficientes o deesentaciones que se dan a

    destacar que en el caso detérmino de mayor grado, que

    siguientes:

    encerrados entre corchetes)iente hasta el coeficiente de, encerrada entre comillas, y

    ponden a los coeficientes del

    siguiente información, para

    n el vector corresponde a las

    se pueden declarar todos losferenciarlos.

    y([6 1 -4 1],'x','coef') 

  • 8/16/2019 Unidad 1-Modo Consola

    16/23

     

    Operaciones con polinomios

    Una vez declarados los polinomios,división y potenciación.

    Sean por ejemplo los polinomios $%las siguientes instrucciones: p=poly([6 1 -4 1],”x”,”coe

     q=poly([-4 1 -2],”x”,”coef

     En la siguiente tabla se muestraresultados devueltos por Scilab:

    Operació

     p+q

     p-q

     p*q

    q^2

    Para el caso especial de la divisiónvalores, el cociente y el resto de la

    q , se puede escribir la siguiente inst

     [r,c]=pdiv(p,q)

    En donde la expresión [r,c] hac

    obtendrá el resto y en la segundaresultados:

    c =

    1.75 - 0.5xr =

    13 - 2.75x

    Para verificar que estos resultados

    como resultado tiene que dar el divi

    -->c*q+r

    ans =2 3

    6 + x - 4x + x

    Evaluación de un polinomio para

    Para obtener el valor de un polinom

    Por ejemplo, si se desea obtener el

    se hace uso de la instrucción horn

     

    Ing. Juan E. Núñez Mc Leod

    se pueden realizar las operaciones básicas co

    ( " 4  6  y )%(  2 4,”)

    las distintas operaciones matemáticas anteri

    Resultado

    2 32 + 2x - 6x + x 2 3

    10 - 2x + x

    2 3 4 5- 24 + 2x + 5x - 10x + 9x - 2x

    2 3 416 - 8x + 17x - 4x + 4x

    de polinomios, se debe recurrir a una función esoperación. Si lo que se desea hacer es dividir el

    rucción:

    referencia a un vector de dos dimensiones, e

    el cociente de la división. De esta manera S

    son correctos, se puede multiplicar el cociente

    endo, tal como se observa a continuación:

    un valor determinado de la variable

    io, cuando la variable adopta valores específicos

    valor del polinomio p anteriormente visto, cuandr(p,5), la cual devuelve el valor 36, como se v

    Página 16 de 23

    o suma, resta, multiplicación,

    definidos en Scilab mediante

    iormente mencionadas y los

    pecífica, la cual devuelve dospolinomio p por el polinomio

    la primera de las cuales se

    cilab devuelve los siguientes

    or el divisor y sumar el resto,

    , se utiliza la función horner.

    la variable x toma el valor 5,

    e a continuación:

  • 8/16/2019 Unidad 1-Modo Consola

    17/23

     

    -->horner(p,5)ans =

    36.

    Derivadas de un polinomio

    Supóngase que el polinomio p(x),determinado sistema de referencia

    1

    móvil y la derivada segunda represderivat que permite obtener la de

    v=derivat(p)a=derivat(v)

    Scilab devuelve los siguientes resul 

    v = 21 - 8x + 3x

    a =- 8 + 6x

    Si se desea conocer el valor de launidades de tiempo, se puede escri

    Dato a

    Posición

    Velocidad

    Aceleración

    Otra forma más rápida de hacermatemático, por lo que se puede de

    x=[0 10]

    Para posteriormente obtener los res

    CPosición iniVelocidad iAceleració

    Determinación de las raíces de u

    Dado un polinomio se pueden obte

    los siguientes ejemplos, que se bas

    13 En este caso la variable x representa el ti

     

    Ing. Juan E. Núñez Mc Leod

    anteriormente declarado, indica la posición3, entonces la derivada primera de este polinomnta la aceleración. Para obtener estas expresiorivada de un polinomio, según las siguientes inst

      ados:

    posición, velocidad y aceleración del móvil enbir lo siguiente:

    calcular  Instrucción ResultInicial horner(p,0) 6.

    Final horner(p,10) 616.

    Inicial horner(v,0) 1.

    Final horner(v,10) 221.

     Inicial horner(a,0) -8.

    Final horner(a,10) 52.

      lo mismo es recordando que Scilab tiene unfinir un vector con los dos valores a calcular, de l

    ultados buscados según la siguiente tabla:

    lculo Instrucción Result  icial y final horner(p,x) 6. 6

      nicial y final horner(v,x) 1.  2 inicial y final horner(a,x) -8.

    polinomio

    er las raíces del mismo mediante la instrucción

    an en los polinomios p(x) y q(x) vistos anteriorm

      iempo.

    Página 17 de 23

    e un móvil respecto de unio representa la velocidad deles, se hace uso de la funciónrucciones:

    el momento inicial y a las 10

    do

    a base matricial de manejola siguiente manera:

    do16.

    1.

    52.

    roots, como se observa en

    nte:

  • 8/16/2019 Unidad 1-Modo Consola

    18/23

  • 8/16/2019 Unidad 1-Modo Consola

    19/23

     

    PRECAUCIONES AL ESCRIBIR E

    Siempre se debe recordar que Scilocasionar situaciones que, en el p

    mensajes de error que alerten al us Funciones con matrices como ar Cualquier función matemática adelemento a elemento, devolviendoAsí se pueden realizar los siguiente

    -->x=[90 45 60;0 180 270]x =

    90. 45. 60.0. 180. 270.

    -->y=sind(x)y =

    1. 0.7071068 0.80. 0. - 1.

    -->x=[1 4 9;36 100 144]

    x =

    1. 4. 9.36. 100. 144.

    -->y=sqrt(x)y =

    1. 2. 3.6. 10. 12.

    De esta manera se puede decir quuna matriz del mismo tamaño, en laal elemento original.

    En un lenguaje simbólico matemátic

    15 Nunca se debe olvidar que el responsable

    computadora está entregando resultados vál16 Se debe recordar que la función sind(x) c

    17 Se debe recordar que en lenguaje simbóli

    Ing. Juan E. Núñez Mc Leod

    PRESIONES MATEMÁTICAS

    ab tiene una estructura de cálculo matemáticoor de los casos, originen resultados erróneos

    15

      ario de que algo anda mal.

    umentos

    ite matrices como argumentos, en tales cauna matriz de las mismas dimensiones con loss cálculos

    16:

    60254

    al aplicar una función matemática a una matrizcual cada elemento de la nueva matriz será igu

    o, esto se escribiría de la siguiente manera17

    :

      * +, - , * + 

    / % (  %( * %+(, - , %( * %+( 

    de los resultados, sobre los que se toman decisiones, es el i

    idos.lcula el seno de ángulo expresados en sistema sexagesimal

      o las matrices se representan con letras mayúsculas y los e

    Página 19 de 23

    matricial. Olvidar esto puede, y en el mejor de los casos,

    os Scilab aplica la funcióndistintos valores funcionales.

    , se obtendrá como resultadol a la aplicación de la función

    ingeniero, quien debe verificar que la

    . lementos con minúsculas. 

  • 8/16/2019 Unidad 1-Modo Consola

    20/23

     

    Problemas con los argumentos m A continuación se verán tres casosel lector debe estar prevenido que

    matemáticas con Scilab, y evaluar s

    Multiplicación y división

    Si se desea evaluar la expresióninstrucciones:

    -->x=2x =

    2.

    -->y=x*sin(x)y =

    1.8185949

    Siendo este último un resultado cor

    Difícilmente se utilice Scilab paracalculadora científica. Un caso másDe esta manera se piensa en gencomo Scilab trabaja en forma matri

    Para calcular diez valores de la fun

    -->x=linspace(0,2*%pi,10)x =

    column 1 to 60. 0.6981317 1.3

    column 7 to 10

    4.1887902 4.8869219

    -->y=x*sin(x)!--error 10

     Multiplicación inconsisten

     Sorprendentemente, en lugar deindicando que la multiplicación ematemático matricial, por lo que tod

    El operador *  es, por lo tanto, remultiplicación de matrices, tal que sm x s y el resultado será una matriser igual a la cantidad de filas de la

    En el caso visto anteriormente, se h

    primera tiene diez columnas mieninconsistente, no se puede realizar.

    Si se piensa en qué es lo que sematricial, sino de una lista de valorese denomina cálculo elemento a el

    elemento de sin(x), luego hacer l

    Ing. Juan E. Núñez Mc Leod

    atriciales

    en los cuales el uso de argumentos matricialesstos son solo ejemplos, y debe recordarlos cad

    i debe tomar o no precauciones adicionales.

     %(  para x=2, entonces se pueden es

    ecto (suponiendo que x esté expresado en radia

      realizar un cálculo de esta naturaleza, ya querazonable sería tener que armar una tabla derar un vector con una serie de valores de x, yial debería devolver otro vector con todos los re

      ión entre 0 y 2π, se podrían escribir las siguient

    62634 2.0943951 2.7925268

    5.5850536 6.2831853

    e.

    btener los diez resultados esperados, sólo ainconsistente. Esto sucede justamente porq

    as la operaciones las interpreta en este contexto

      presentativo del producto matricial, en el cual

    i se desea multiplicar una matriz n x m, sólo esn x s. En otras palabras, la cantidad de colum

    segunda.

    a multiplicado una matriz x de 1 x 10 por otra s

      tras que la segunda tiene una única fila. Cla

    pretende hacer con el cálculo, es claro que ns a la cual se debe aplicar la función de maneraemento; es decir, tomar el primer elemento de

    s mismo con los segundos elementos y así suc

    Página 20 de 23

    ocasiona problemas. Aunquevez que escriba expresiones

    ribir en Scilab las siguientes

    nes).

    seguramente se usaría unaalores para hacer un gráfico.evaluar la expresión anterior,ultados.

    s instrucciones:

    .4906585

    parece un mensaje de errore Scilab es un procesador

    .

    se debe cumplir la regla de

    osible hacerlo por otra matrizas de la primera matriz debe

    n(x), también de 1 x 10. La

    ramente la multiplicación es

    o se trata de una operaciónindependiente, esto es lo que

    y multiplicarlo por el primer

    sivamente hasta terminar.

  • 8/16/2019 Unidad 1-Modo Consola

    21/23

     

    Para realizar multiplicaciones de m

    seguido inmediatamente del asterisresultado obtenido es el deseado, el

    -->y=x.*sin(x)y =column 1 to 6

    0. 0.4487504 1.3

    column 7 to 10- 3.6275987 - 4.8126786

    El mismo problema se tiene con el

    una división elemento a elemento,

    intermedios).

    Para analizar esta última operació

    responde a la conocida ley . 5Suponiendo que la constante valgaa 10 atmósferas, se podrían escribir -->k=100;

    -->p=linspace(1,10,5) p =

    1. 3.25 5.5 7

    -->v=k/p

    v =

    0.49535601.60990712.72445823.83900934.9535604

    Analizando el resultado obtenido sconstituidos por un vector fila. Esteanda mal.

    Por otro lado, haciendo un análisi

    operación 10071 (la constantevolumen se calcula como 1 0Revisando los valores del vector v 

    cálculo no es correcto; sin embargo

    Lo que Scilab ha hecho es un cálcucomo resultado 100 (el valor de la c

    18 Para que esta multiplicación se pueda re19

     Este valor debería ser 0 (cero), sin embartrata de un valor muy pequeño.

    Ing. Juan E. Núñez Mc Leod

    atrices, pero elemento a elemento18

    , se utiliza

    o, sin dejar espacio intermedios). De esta manl cual se muestra a continuación:

    5051 1.8137994 0.9551004 - 1

    - 3.5900033 - 1.539D-1519 

    operador /, utilizado para realizar divisiones mse debe utilizar el operador ./  (un punto segu

    se verá el caso de una transformación isotér

    , donde k es una constante.

    100, para evaluar el volumen de un gas ideal culas siguientes instrucciones en Scilab:

    75 10.

    observa que el mismo es un vector columna,cambio en la naturaleza de los vectores es un

    más puntual para calcular el primer valor del

    dividida 1 atmósfera), la cual da como resultad710  (la constante k dividida 10 atmósferas),o se encuentra ninguno de estos valores citado

    , extrañamente Scilab no ha dado ningún error.

    lo matricial, dando por resultado un vector tal qonstante), según se ve a continuación:

    lizar las dos matrices deben tener las mismas dimensiones.

    o por problemas de precisión numérica, discutidos anterior

    Página 21 de 23

    el operador .*  (es un puntora, al utilizar este operador el

    .1938755

    triciales. Si se desea realizar

    ido de la barra, sin espacios

    ica de un gas ideal, la cual

    ando la presión aumenta de 1

    mientras que los datos estánaviso importante de que algo

    volumen se debe realizar la

    o 100, mientras que el últimoque da como resultado 10., dando un claro aviso que el

    e multiplicando al primero da

    ente, el valor no es nulo aunque se

  • 8/16/2019 Unidad 1-Modo Consola

    22/23

     

    -->p*vans =

    100.

    Para poder calcular lo que realmenobserva en la siguiente serie de inst -->k=100;

    -->p=linspace(1,10,5) p =

    1. 3.25 5.5 7

    -->v=k./pv =

    100. 30.769231 1

    Donde no se observa cambio en lacon los datos y se obtiene como rcalculados la discusión previa previ Para hacer el cálculo anterior se d

    directamente con el valor de la cons

    -->v=100./p Warning: "100./ ..." is iwise operation

    v =

    0.49535601.60990712.72445823.83900934.9535604

    Analizando el resultado se obtieneparte del número, por lo que la divdebe utilizar un espacio despuésdivisión como una operación eleme

    -->v=100 ./pv =

    100. 30.769231 1

    Esta discusión también es válida pa

    Potenciación

    Si bien el operador ^ permite calcul

    debe utilizar la notación de element

    Así por ejemplo, si se quiere tenevalores de x desde -5 a 5, se puede

    Ing. Juan E. Núñez Mc Leod

    te se desea, es necesario realizar un cálculo elrucciones:

    75 10.

    .181818 12.903226 10.

    naturaleza de los vectores de datos y resultad  sultado otro vector fila), además los valores in

    mente.

    finió la variable k con el valor 100; sin embar

    tante, tal cual se ve a continuación:

    terpreted as "100.0/ ...". Use "10

    un mensaje de error que indica que el puntoisión se calcula como una operación matricial.el número, para que el punto sea interpretadoto a elemento, tal cual se observa a continuació

      .181818 12.903226 10.

    ra el caso de la multiplicación elemento a eleme

    ar la potencia de escalares, cuando la base est

    a elemento; es decir, el operador .^ en lugar d

    una tabla de 10 valores de la función 3   escribir lo siguiente:

    Página 22 de 23

    mento a elemento, según se

    os (se parte de un vector filaicial y final coinciden con los

    o, se podría hacer el cálculo

    0 ./ ..." for element

    ecimal es interpretado comoPara evitar este problema seen conjunto con la barra de

    n:

    to vista anteriormente.

    constituida por un vector, se

    el primero mencionado.

    8 1  que abarquen los

  • 8/16/2019 Unidad 1-Modo Consola

    23/23

     

    x=linspace(-5,5,10)

    y=3*x.^2-5*x+1

    Obteniendo por resultado lo que se

    x =column 1 to

    - 5. - 3.8888889column 7 t

    1.6666667 2.7

    y =column 1 to

    101. 65.81481column 7 t

    1. 10.259259

    Si por el contrario se escribiera

    -->y=3*x^2-5*x+1

    ! Warning: Syntax "vector ^Use "vector .^ scalar" ins 

    y =column 1 to

    101. 65.81481column 5 to

    4.7037037 - 0.8column 10

    51.

    Como se aprecia se obtienen los rindicando que la posibilidad de utili

    en una versión próxima del program Más allá del tema de sintaxis; enmatricial, ya que no es necesariohacer con una calculadora, sino qurepresenta a la variable independiecon muy pocas instrucciones.

    AYUDA 

    En modo consola se puede solicit

    función requerida. Así si se desea l

    -->help sin

    Y se abrirá una nueva ventana denla función.

    expresa a continuación:

    6- 2.7777778 - 1.6666667 - 0.5555510

    77778 3.8888889 5.

    638.037037 17.666667 4.703

    1026.925926 51.

    calar" is obsolete. It will be remoead.

    438.037037 17.666667

    9

    18519 1. 10.259259 26.9259

    sultados correctos; sin embargo, Scilab muestrar el operador ^ para cálculos donde las bases

    a.

    este ejemplo, se observa la ventaja que ofrecealizar el cálculo de la expresión matemáticae Scilab evalúa la función para cada uno de lonte. De esta manera, fácilmente se pueden gen

    r ayuda sobre cualquier función de Scilab, es

    er información sobre la función seno, se debe e

      ominada Navegador de Ayuda, donde se despl

    6 0.5555556

    037 - 0.8518519

    ed in Scilab 6.0.

    6

    un mensaje de advertencia,ean vectores, será eliminada

    e Scilab al trabajar en modoiez veces, como se debería

    s elementos de la matriz queerar grandes tablas de datos

    ribiendo help  seguido de la

    cribir:

    gará la información relativa a