Funções Predefinidas

Ver Todos os Exercícios

Teoria

E depois da tempestade....

Iaê galerinha.... Aula passada foi difícil né?

Sei como é, quando vi aquilo pela primeira vez também achei, mas vamos para outro papinho?

Juro que essa conversa aqui vai ser mais light e beeeeeeeeeeeeem mais prática.

Bem pessoal, hoje finalizaremos a nossa introdução ao Matlab e aqui iremos estudar algumas funções predefinidas que servem para nos auxiliar em diversas áreas e aposto que algumas delas você irá amar, eeeeeeeeee vamos lá.

Bem, já foi falado na primeira aula que o Matlab funciona como uma espécie de calculadora, então se você olhar na sua calculadora existe um bocado de função com algumas matemáticas e se a sua for mais sofisticada pode ter até diferenciação e integração, mas... o “Matlab tem isso? ”

Para a NOOOOOOOOOOOOSA ALEGRIA tem sim, e são sobre essas funções que iremos conversar hoje.

Funções Matemáticas

Assim como toda calculadora, o Matlab também possui funções matemáticas e elas são diversas, aqui iremos listar algumas das mais utilizadas na engenharia para que possamos nos situar.

Por exemplo, vamos supor que você quer retirar a raiz quadrada de 4, basta escrever o seguinte código:

A = 4 ;

B = s q r t ( A )

E sabe o que vai sair? Isso ae B = 2 .

Outro exemplo é se quisermos calcular a raiz de -4.

OPAA, olha só, o Matlab também processa números complexos, e a forma de escrevê-los é do jeito que está ai, primeiramente a parte real é uma constante seguindo as mesmas lógicas e depois a parte imaginária se multiplica por i.

Função de integração.

No Matlab existem dois tipos de integração a numérica e a simbólica, a integração numérica é uma espécie de algorítmo que serve para realizar a integração de pontos experimentais e obter uma expressão aproximada, isso será visto em cursos futuros.

Aqui no momento o nosso foco é a integração simbólica, que é aquela integração que você aprendeu no cálculo diferencial e essas funções de integração simbólica só operam em funções simbólicas, para nossa sorte só existe um único comando que é o comando i n t .

Para entendermos a sintaxe vamos fazer uns exemplos...

Primeiramente devemos definir as variáveis simbólicas:

S y m s   a   b   x

Depois vamos inserir a função simbólica:

F = a / s q r t ( x ^ 2 + b ^ 2 )

E agora iremos integrar.

I = i n t ( F , x )

Então...

Então, de uma forma geral, a sintaxe da função i n t é :

I n t ( f u n ç ã o , v a r i á v e l   d e   i n t e g r a ç ã o )

Se a integral for definida ainda podemos generalizar mais

I n t f u n ç ã o , v a r i á v e l   d e   i n t e g r a ç ã o , l i m i t e   i n f e r i o r   d e   i n t e g r a ç ã o ,   l i m i t e   s u p e r i o r   d e   i n t e g r a ç ã o

Funções de derivação.

Assim como na integração, existe a diferenciação numérica e simbólica, mas o foco aqui será a diferenciação simbólica belezinha?

Então vamos lá:

Na derivação as mais importantes para a engenharia é a função d i f f e a função g r a d i e n t e igualmente a integração ambas operam sobre variáveis simbólicas.

Galera, a função diff ela tira a derivada da função, para deixar mais claro vamos fazer um exemplo:

Da mesma forma iremos escrever primeiramente as variáveis simbólicas e as funções:

S y m s   a   b   c   x

F = a * x ^ 2 + b * x + c

Para derivar precisamos fazer D = d i f f F , x e se você fez tudo certinho deve aparecer algo assim:

E assim podemos obter qualquer derivada cavernosa, e isso funciona também como derivada parcial, olha só:

S y m s   a   b   c   x   y

F = a * y * x 2 + b * y 2 * x + c + y

D = d i f f ( F , x )

Isso é possível porque ele considera tudo símbolo excerto a variável que está na função diferenciação ou integração.

Outra função muito utilizada é a famosa função gradiente, que é aquele mesmo vetorzinho que você vê em cálculo diferencial, aqui existe uma função destinada a seu cálculo que é a função g r a d i e n t .

Vamos resolver a seguinte função

S y m s   x   y

F = y * x 2 + y 2 * x + y

E simplesmente vamos utilizar D = g r a d i e n t ( F , x , y ) .

Fora essas duas existem mais três que são utilizadas em cursos avançados de cálculo, mas é pouco vista em programação destinada a engenharia, elas são a divergência e o laplaciano, o laplaciano iremos ver aqui, porém, para o cálculo do divergente precisamos entrar com um vetor simbólico e a criação de vetores será um tópico das próximas aulas.

Então, faremos um exemplo com laplaciano, ok?

S y m s   x   y

F = y * x 2 + y 2 * x + y

D 2 = l a p l a c i a n ( F )

Logo, para a derivada, a sintaxe geral é:

D = d i f f ( f u n ç ã o   a   d e r i v a r ,   V a r i á v e i s   a   d e r i v a r )

Se quer tirar o gradiente, divergente ou laplaciano já que deriva todas as variáveis só precisa do

D = ( g r a d i e n t ,   d i v e r g e n c e   o u   l a p l a c i a n ) ( F u n ç ã o   a   d e r i v a r )  

Opa galera, existem outras funcionalidades no Matlab, mas acredito que essas aqui são as mais importantes, como já foi falado essa parte é bem prática e precisa realmente treinar, devido a isso selecionamos alguns exercícios com todo carinho do mundo para vocês treinar, então nos vemos lá.

Função de integração.

Funções de derivação.

Exercícios Resolvidos

Exercício Resolvido #1

Calcule as funções abaixo utilizando o Matlab;

  1. Exponencial de 5
  2. Raiz quadrada de 4
  3. Raiz cúbica de 8
  4. Cosseno de π 4 (obs: π   no matlab é escrito como pi)

Passo 1

E aí galerinha? Vamos resolver mais uma questão de programação?

Essa questão é bem prática e, basicamente, é brincar com a tabela de funções. Essa tabelinha bonita aqui!

Passo 2

Pra vocês verem exatamente como fica no Matlab, já vou colocar os cálculos na janela do command window.

  1. A exponencial de x fica exp(x). Então a exponencial de 5 fica exp(5)
  2. Passo 3

  3. Pra achar a raiz quadrada de x, fazemos sqrt(x). Então a raiz quadrada de 4 é sqrt(4).
  4. Passo 4

  5. Não tem função pra raíz cúbica. Então vamos elevar ao expoente ⅓.
  6. Passo 5

  7. O cosseno é bem simples ne, é bem como a gente escreve na matemática mesmo. E o nosso π vai ser do jeito que a gente fala (pi).

Resposta

  1. 148,4132
  2. 2
  3. 2
  4. 0,7071

Exercício Resolvido #2

Em termodinâmica podemos calcular o trabalho da expansão de um gás utilizando a expressão:

w = ∫ p * d V

Utilize a equação de Estado de Wan der Walls para calcular o trabalho de expansão de um gás de Wan der Walls

Obs: Equação de wan der walls     P = R * T V - b - a V 2

Passo 1

E aí, pessoal? Tudo bom?

Nessa questão, vamos precisar fazer uma integração simbólica. O primeiro passo é a inserção da função como uma função simbólica. Para isso primeiro precisa introduzir as variáveis simbólicas.

Olhando o problema temos que as variáveis são a, b, V, R e T, sendo assim:

S y m s   a   b   V   R   T

E a função é dada por:

P = ( R * T ) / ( V - b )   - a / V ^ 2

Se tudo foi feito direitinho devemos ter algo do tipo:

Passo 2

O passo 2 é integrar a função e, uma vez que a integral é indefinida, apenas precisamos utilizar a forma padrão da função i n t e a variavel de integração é o V.

I = i n t P , V

No Matlab fica assim:

SIIIIIIIIM uma coisa importante o Matlab não conta a constante de integração +C então se estiver fazendo uma prova de cálculo ou algo que precise dele não esquece de somar ela tá?

Daí fica assim:

w = a V + R   T log ⁡ V - b + C

Resposta

w = a V + R   T log ⁡ V - b + C

Exercício Resolvido #3

Em eletrodinâmica podemos calcular a força entre duas partículas a partir do potencial elétrico utilizando a seguinte expressão:

F = - d V d R

Utilize a expressão do potencial elétrico para calcular a força entre duas partículas carregadas

Obs: Equação Força é dada por   V = k * Q 1 Q 2 R

Passo 1

Fala aí galerinha! O primeiro passo, para fazer uma diferenciação simbólica, é a inserção da função como uma função simbólica. E, para isso, primeiro precisamos introduzir as variáveis simbólicas.

Mas como sabemos que a questão é de variável simbólica? Ou até mesmo quais variáveis são simbólicas?

Bem pessoal, existem diversas formas de deduzir isso, mas a forma mais comum é você perceber que o que a questão está te pedindo é uma expressão.

Em geral, o Matlab gera como saída algum número, então quando a questão te pede uma letra ou uma expressão pode confiar que o problema tem alguma base nas variáveis simbólicas, belezinha?

E para completar precisamos definir quais variáveis são simbólicas né?

As variáveis simbólicas são todas aquelas variáveis pela qual a sua função simbólica depende, ou seja, nessa questão corresponde às variáveis Q 1 , Q 2 , K , R .

Passo 2

Utilizando as informações do passo 1 precisamos implementar as variáveis simbólicas no Matlab, que corresponde a:

S y m s   Q 1   Q 2   K   R

Onde a função syms é a função que utilizamos para indicar que as variáveis que vem depois dela é tudo variável simbólica, ok?

consequentemente a função simbólica é dada por:

V = k * Q 1 * Q 2 / R

Se tudo foi feito direitinho devemos ter algo do tipo:

Passo 3

O próximo passo é derivar a função e, para isso, apenas precisamos utilizar a forma padrão da função d i f f que vimos na aula. A variável de derivação é o R.

F = - d i f f V , R

Se você digitou tudo certinho no seu Matlab deve mostrar algo mais ou menos assim:

Resposta

Exercício Resolvido #4

Em eletrodinâmica podemos calcular a força entre duas partículas a partir do potencial elétrico utilizando as seguintes expressões:

F =   g r a d V

F x = d V d x

F y = d V d y

Utilize a expressão do potencial elétrico para calcular a força entre duas partículas carregadas

Obs: Equação do potencial é dado por   V = k * Q 1 Q 2 R   e lembrar que   R = x 2 + y 2

Passo 1

Fala aí, pessoal! O primeiro passo pra fazer uma diferenciação simbólica é a inserção da função como uma função simbólica e, para isso, primeiro precisa introduzir as variáveis simbólicas.

Olhando o problema, temos que as variáveis são Q1, Q2, K e o R, sendo assim:

S y m s   Q 1   Q 2   K   R

E a função é dada por:

V = k * Q 1 * Q 2 / s q r t ( x 2 + y 2 )

Se tudo foi feito direitinho devemos ter algo do tipo:

Passo 2

O passo 2 é obter o gradiente da função, que pode ser obtido apenas utilizando a função gradient, logo:

I = g r a d i e n t V

Socorro, que resposta enorme é essa? Parece que temos 5 derivadas, sendo que nós temos apenas duas componentes de força....

Relaxaaaaaa você não fez nada de errado. É que a função gradient deriva a função em relação a todas as variaveis adicionadas. Então toma cuidado com isso ok?

Mas lembra as demais expressões: F x = d V d x

F y = d V d y

Então so iremos pegar as derivadas em relação a x e a y que corresponde as duas ultimas linhas do nosso resultado.

Resposta

Ei, a resposta está no passo a passo :)

Exercício Resolvido #5

Quando estudamos Mecânica dos materiais vimos de forma aprofundada como calcular a distribuição de força por unidade de área em sistemas deformáveis, uma dessas forças é chamada de tensão de cisalhamento e pode ser calculada para um solido linear pela expressão:

F =   ∫ ∫ ( ∇ 2 u * d x ) d y = ∫ ∫ ( L a p l a c i a n ( u ) * d x ) d y

Utilize a seguinte expressão para o cálculo dessa força de cisalhamento:

u = x * y 2 + x 2 * y + x y

Passo 1

O primeiro passo como foi visto na teoria é a inserção da função como uma função simbólica. Para isso, primeiro precisamos introduzir as variáveis simbólicas.

Olhando o problema temos que as variáveis são x e y sendo assim:

S y m s   x   y

E a função é dada por:

u = x * y ^ 2 + x ^ 2 * y + x y

Se tudo foi feito direitinho devemos ter algo do tipo:

Passo 2

O segundo passo é tirar o laplaciano da expressão anterior

I = l a p l a c i a n u

Se tudo foi feito certinho, vai aparecer algo assim na sua tela:

Passo 3

Agora que temos as expressões basta integrar em relação a uma variável e depois a outra, logo:

I 1 = i n t ( I , x )

I 2 = i n t ( I 1 , y )

E finalmente ira aparecer a resposta final na sua tela:

E essa é a função que vai descrever a distribuição de força ao longo de um solido linear deformável.

Resposta

F = x   x + 2 y 2 4

Exercício Resolvido #6

A distância d   de um ponto ( x 0 , y 0 ) a reta A x + B y + C = 0 é dada por:

d = | A x 0 + B y 0 + C | A 2 + B 2

Determine a distância do ponto ( 2 , - 3 ) a reta 3 x + 5 y - 6 = 0 . Primeiramente, declare as variáveis A , B , C , x 0   e   y 0 , então calcule d .

Passo 1

E aí galerinha? Vamos resolver mais uma questão de programação?

Essa é uma questão de cálculo simples, onde primeiramente precisamos implementar a entrada de dados, que são as coordenadas do ponto e os coeficientes da reta. E então, a partir da fórmula dada pelo enunciado, iremos calcular a distância da reta até o ponto, usando funções matemáticas predefinidas existentes no Matlab.

É um tipo de questão bastante útil para o cálculo vetorial, então essa é uma ferramenta bastante útil para te ajudar a resolvê-las, ok?

Então precisamos seguir os seguintes passos para criar o programinha:

1.Entrada de dados – Inserir as variáveis A , B , C , x 0   e   y 0 .

2.Cálculo do numerador

3.Cálculo do denominador

4.Cálculo da distância d

Passo 2

Agora que sabemos o passo a passo precisamos transformá-los em uma linguagem de programação, beleza?

Vamos inserir o primeiro passo que é a entrada de dados, que basicamente definir os valores das 5 variáveis, sendo assim:

A = 3

B = 5

C = - 6

x 0 = 2

y 0 = - 3

Agora que entramos com a entrada de dados precisamos criar o numerador e denominador.

Passo 3

Vamos começar com o numerador, que é um módulo, então para isso precisamos calcular o valor absoluto da expressão e essa função no matlab é a função a b s , logo:

n u m = a b s A * x 0 + B * y 0 + C

Agora vamos implementar o denominador, onde vemos a raiz dos coeficientes da equação da reta, logo para o cálculo do denominador precisamos utilizar a função s q r t para calcular a raiz quadrada, logo:

d e m = s q r t ( A ^ 2 + B ^ 2 )

Por fim, a distância entre o ponto e a reta é

d = n u m / d e m

Passo 4

Se você digitou tudo certinho no Matlab vai aparecer as seguintes informações:

Resposta

Se tudo foi inscrito certinho o resultado é d = 2.5725

Exercício Resolvido #7

Declare duas variáveis: α = 5 π 9 e β = π 7 e com isso calcule ambos os lados da seguinte igualdade:

c o s c o s   α   - c o s c o s   β   = 2 * s i n s i n   α + β 2   * s i n ⁡ ( α - β 2 )

Após calcular cada lado da equação confirme a existência da igualdade desta identidade trigonométrica.

Passo 1

E aí galerinha? Vamos resolver mais uma questão de programação?

Essa é uma questão de cálculo simples, onde primeiramente precisamos implementar a entrada de dados e a partir de funções predefinidas existentes no Matlab iremos calcular essa igualdade trigonométrica.

Então para isso precisamos seguir os seguintes passos para criar o programinha:

1.Entrada de dados – Inserir as variáveis α , β .

2.Cálculo do termo a esquerda da igualdade

3.Cálculo do termo a direita da igualdade

Passo 2

Agora que sabemos o passo a passo precisamos transformá-los em uma linguagem de programação, belezinha?

Vamos inserir o primeiro passo que é a entrada de dados, que basicamente é programar os dados disponibilizados, sendo assim:

a l p h a = 5 * p i / 9

b e t a = p i / 7

Agora que entramos com a entrada de dados precisamos calcular o lado esquerdo e direito da identidade acima, então para isso precisamos calcular cada etapa.

Sendo assim, primeiramente vamos começar com o lado esquerdo e para isso precisamos calcular o valor das diferenças do cosseno, sendo assim:

F = c o s c o s   a l p h a   - c o s ⁡ ( b e t a )

Agora vamos implementar o lado direito da equação, onde vamos calcular os produtos de seno. Lembrando que o seno no Matlab é chamado de s i n , sendo assim:

G = s i n s i n   a l p h a + b e t a 2   * s i n s i n   a l p h a - b e t a 2   * 2

Passo 3

Se você digitou tudo certinho no Matlab vai aparecer as seguintes informações:

Resposta

Logo analisando o Resultado vemos que essa igualdade não existe para esse valor de a l p h a e b e t a , logo essa igualdade não é válida.

Exercício Resolvido #8

Em cálculo de reatores batelada uma das grandezas mais importantes para o dimensionamento é o tempo de reação.

Um reator Batelada é modelo de reator cuja a reação varia no tempo, em geral, é bastante utilizado em sistemas de produção de medicamentos e processos de produção de etanol, a equação que dita qual é o tempo do processo é dada por:

t = ∫ 0 t f d C a R a

Onde C a é a concentração do reagente limitante e R a é uma função que dita como a concentração do componente varia dentro do reator.

Através de medidas experimentais um grupo de alunos encontrou a seguinte expressão para o R a .

R a = C a 3 * C b * k

Onde k é uma constante de proporcionalidade e a concentração do componente b é dada por:

C b = C o n s t a n t e - C a

Utilizando essas informações encontre uma expressão para o cálculo do tempo do processo.

Passo 1

E aí galerinha? Vamos resolver mais uma questão de programação?

Essa é uma questão que é bastante vista em cursos de reatores químicos, em geral o tempo do processo é a chave para dimensionar e projetar o reator, belezinha?

Aqui você vai ver como a gente calcula esse tempo a partir de uma equação estimada através de dados experimentais, sendo assim precisamos de duas coisas aqui.

Primeiramente é as variáveis simbólicas e segundamente é a função i n t .

Então para isso precisamos seguir os seguintes passos para criar o programinha:

1.Entrada de dados – Inserir as variáveis simbólicas e as funções disponibilizadas no problema.

2.Cálculo da integral através da função i n t do Matlab.

Passo 2

Agora que sabemos o passo a passo precisamos transforma-los em uma linguagem de programação, belezinha?

Vamos inserir o primeiro passo que é a entrada de dados, que basicamente é programar os dados disponibilizados, sendo assim:

s y m s   C a   C b   k

C b = c o n s t a n t e - C

r a = C a ^ 3 * C b ^ 2

Agora que entramos com a entrada de dados precisamos calcular a integral para assim encontrar uma expressão para o tempo, portanto:

t = i n t 1 r a , C a

Passo 3

Se você digitou tudo certinho no Matlab vai aparecer as seguintes informações:

Resposta

Logo, a sua resposta deve ser algo parecido com isso aqui: