Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
0% encontró este documento útil (0 votos)
18 vistas4 páginas

Ejercicios Matlab

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1/ 4

EJERCICIOS MATLAB

1 Sistema de ecuaciones lineales. Regla de Cramer

Sea un sistema lineal de tres ecuaciones con tres incógnitas. Se aplica la regla de Cramer

( )( ) ( )
a11 a12 a13 x1 b1
a21 a22 a23 x 2 = b2 Ec. 1
a31 a32 a33 x 3 b3

| | | | | |
b 1 a 12 a13 a11 b1 a13 a11 a12 b1
b 2 a 22 a23 a21 b2 a23 a21 a22 b2
b 3 a 32 a33 a31 b3 a33 a31 a32 b 3
x 1= x 2= x 1=

| | | | | |
, , Ec. 2
a 11 a12 a13 a 11 a12 a13 a 11 a12 a13
a21 a22 a23 a21 a22 a 23 a21 a22 a23
a31 a32 a33 a31 a32 a33 a31 a32 a33

A ∙ X=B Ec. 3

Donde A es la matriz de los coeficientes, B es el vector de los términos independientes y X es


el vector de las incógnitas.

Se elabora un script para resolver un sistema de ecuaciones lineales que es compatible y


determinado y se aplica al siguiente ejemplo.

{
4 x 2−2 x3 −2 x 4 =−4
x1 +2 x 2+ 4 x 3−3 x 4=5
Ec. 4
−3 x1−3 x 2 +8 x 3−2 x 4=7
−x1 + x 2 +6 x 3−3 x 4=7

( )( ) ( )
0 4 −2 −2 x 1 −4
1 2 4 −3 x 2 = 5
Ec. 5
−3 −3 8 −2 x3 7
−1 1 6 −3 x 4 7

Se determina el vector de las incógnitas y se comprueba que el producto A ∙ x es el vector B de


los términos independientes

2 Rellenar un área de color

Se utilizan los vectores en combinación con los operadores lógicos para producir ciertos efectos
que pueden ser útiles en las representaciones gráficas.

El comando fill rellena un recinto cerrado del color especificado. El recinto está descrito por dos
vectores xx e yy que contienen las abscisas y ordendas de los puntos del recinto. Estos
vectores están formados por elementos y por porciones de otros vectores extraídos mediante
operadores lógicos.

Primero, se extrae un vector de otro mediante operadores lógicos.

Se extraen aquellos elementos cuyo índice se corresponde con el valor uno.

S dibuja la función −x 2+ 3 x + 4 y se define la región comprendida entre la curva, el eje X y las


rectas x=1 y x=3 que se quiere colorear con el comando fill. La definición de la región es:

Punto (1 , 0) , seguido de (1 , f (1)) , todos los puntos de abscisa 1< x <3 y sus
correspondientes ordenadas, el punto (3 , f (3)) y finalmente, el punto (3 , 0) . En el vector xx
se guardan las abscisas y en el vector yy , las ordenadas. Se llama al comando fill y se le pasa los
dos vectores y el color de relleno.

El resultado del programa debe ser:

3 Máximos de una función

Se representa la función suma de cuatro armónicos de frecuencias angulares ¿ 1 ,3 ,3.5 , 4 y 6


rad / s en el intervalo (0.1 , 10)

x (t)=cos(t)+ 0.5 ⋅cos (3t )+ 0.4 ⋅cos (3.5 t)+0.7 ⋅cos (4 t)+0.2 ⋅cos (6 t) Ec. 6

La gráfica de la función tiene la siguiente forma:


Analizando el primer máximo que está en el intervalo (1 , 2) . La función crece hasta t=1.7 y
luego decrece. Se elabora la siguiente tabla utilizando el comando diff que calcula la diferencia

x i+1−x i , i=1 , 2 ,3 , 4. . . Ec. 7

dx ( dx >0 )=1
i t x dx=diff ( x ) d 2 x =diff =dx
dx ( dx≤0 )=−1

10 1.0 −0.5948 0.2258 1 0


11 1.1 −0.3690 0.2698 1 0
12 1.2 −0.0992 0.2778 1 0
13 1.3 0.1786 0.2599 1 0
14 1.4 0.4385 0.2216 1 0
15 1.5 0.6601 0.1630 1 0
16 1.6 0.8231 0.0808 1 −2
17 1.7 0.9039 −0.0268 −1 0
18 1.8 0.8771 −0.1550 −1 0
19 1.9 0.7220 −0.2896 −1 0
20 2.0 0.4324 −0.4078 −1 0
El primer máximo corresponde al índice i=17 , el instante t=1.7 y su valor es x=0.9039

En la tabla la variable d 2 x contiene un valor distinto de cero en el índice 16 , uno menos que el
máximo. Se utiliza find para encontrar este índice distinto de cero. Se debe escribir el código
para identificar los máximos de la función.
4 Binomio de Newton

¿ Ec. 8

Se define una función denominada pol newton que devuelve un vector que contiene los coeficientes
del polinomio resultado del desarrollo del binomio, cuando se le pasa el grado n y el valor del
término a .

Primero se calcula el número combinatorio n C m o n sobre m .

(mn )= m! (n−m)
n!
!
=
n ⋅(n−1)⋅(n−2) ⋅⋅⋅⋅(n−m+1)
1 ⋅2 ⋅⋅⋅⋅(m−1)⋅m
Ec. 9

(mn )=( 1n )( n−1


2
)⋅⋅⋅⋅(
n−m+2 n−m+1
m−1
)(
m
) Ec. 10

Se define la función combinatorio

Alternativamente, se puede utilizar la función de MATLAB nchoosek tal como se muestra en la


porción de código
¿ combinatorio (4 , 2)
ans=6
¿ nc h oosek( 4 , 2)
ans=6
Se define la función pol newton , pasándole la potencia n y el valor del coeficiente a , que
devuelve los coeficientes p del polinomio.

Al probar la función se debe obtener:


>>pol_newton(2,1)

ans=1 2 1

>>pol_newton(3,1)
ans=1 3 3 1
>>pol_newton(4,1)
ans=1 4 6 4 1
>>pol_newton(3,-1)
ans=1 -3 3 -1

También podría gustarte