1

MÉTODOS COMPUTACIONAIS
MATLAB: UMA INTRODUÇÃO

Key Fonseca de Lima

Viviana Cocco Mariani
2008
 2
SUMÁRIO

1. Introdução ao Matlab
2. Operações matemáticas elementares
3. Funções matemáticas elementares
4. Formatos de visualização numérica
5. Vetores
6. Matrizes
7. Análise de dados e funções estatísticas
8. Gráficos bidimensionais
9. Ajuste de curvas e interpolação unidimensional
10. Gráficos tridimensionais
11. Gráficos em coordenadas polares

12. Polinômios
13. Matemática Simbólica
14. Operadores Relacionais e Lógicos
15. Controles de Fluxo
16. Introdução a Programação
17. Função fprintf
 3
1. INTRODUÇÃO AO MATLAB

MATLAB ⇒ Ferramenta computacional para resolução de problemas de engenharia.

abrevia MATrix LABoratory (Laboratório de Matrizes)

O programa MATLAB foi implementado na linguagem C e oferece uma ampla biblioteca de

funções predefinidas para que a programação se torne mais simples. Essa variedade de funções faz
com que o MATLAB se torne competitivo com outras linguagens de programação como C, Fortran
etc.

1.1 Execução do Matlab

A execução do Matlab pode ser feita de duas maneiras:

I) Online ⇒ Através da janela de comandos (Command Window). (Prompt do Matlab)

Barra de Barra de
Comandos Ferramentas

II) Editor de textos ⇒ Através do Depurador de Arquivos M (Matlab Editor)

Acesso 1:
Janela de comandos ⇒ Barra de comandos ⇒ File ⇒ New ⇒ M-file.

Acesso 2:
 4

Janela de comandos ⇒ Barra de ferramentas ⇒ ícone (New M-file)

 5
2. OPERAÇÕES MATEMÁTICAS ELEMENTARES

Adição ⇒ +
Subtração ⇒ -
Multiplicação ⇒ * (asterisco)
Divisão ⇒ /
Potenciação ⇒ ^ (acento circunflexo)

2.1 Operações Elementares Utilizando a Janela de Comandos.

» 1+2
ans = ⇐ ans é a abreviação de answer (resposta)
3

» 5-8
ans =
-3

» 4*1.2 ⇐ 1.2 é equivalente a 1,2

(No Matlab utiliza-se somente ponto decimal e não vírgula)
ans =
4.8000 ⇐ Formato padrão do Matlab são 4 casas decimais

» -7/2 ⇐ Divisão pela esquerda -7÷2

ans =
-3.5000

» 5^2
ans =
25

2.2 Operações Elementares Utilizando o Depurador de Textos

⇒ MATLAB EDITOR

PROCEDIMENTO:

⇒ETAPA 1
Iniciar o Matlab Editor (ver seção 1.1)

⇒ETAPA 2

Iniciar através dos comandos

clear all ⇒ limpa a memória do MATLAB
clc ⇒ limpa a tela do computador sem limpar a memória
 6

exceto se for usar uma function então estes comandos não deverão ser usados:

⇒ETAPA 3
Salvar o arquivo:
A) Barra de comandos ⇒ File ⇒ Save As...
B) Barra de ferramentas ⇒ ícone

Salvar em: ⇒ diretório c: ⇒ Pasta: Temp

Nome do arquivo

OBSERVAÇÕES:

I) Utilizar nomes de arquivos com no máximo 8 (é possível usar mais, mas não há
necessidade) caracteres sem espaço ou acento gráfico, utilize preferencialmente nomes curtos.

II) O nome do arquivo.m não pode começar com números e utilize letras minúsculas (é possível
utilizar letras maiúsculas).

III) Salve as modificações sempre que alterar um arquivo.m, antes de executá-lo

(*) ⇒ o asterisco indica que a

alteração não foi salva.
 7

IV) Inicie um novo arquivo.m sempre que iniciar um novo cálculo ou programa, ou quando você
desejar.

⇒ETAPA 4
Voltar à Janela de Comandos (prompt do MATLAB) e digitar a pasta onde o arquivo.m
foi salvo.

Digitar cd c:\pasta onde foi salvo o arquivo.m

Se você digitar pwd na janela de comandos vai aparecer o diretório corrente

⇒ETAPA 5
Digitar o nome do arquivo.m salvo, para executá-lo.

Após digitar o nome do arquivo ⇒ <enter>

 8

Mostra o resultado do arquivo exe1.m

MATLAB espera novo comando

2.3 Criação de Variáveis Literais

A) Armazenando em Variáveis Curtas

⇒ MATLAB EDITOR ⇒ RESPOSTA

clear all a =
clc 3.5000
b=
a=3.5 ⇐ cria a variável a=3,5 5
b=5 ⇐ cria a variável b=5 c=
c=a^b ⇐ cria a variável c resultante da operação ab 525.2188

B) Armazenando em Variáveis Longas

⇒ MATLAB EDITOR ⇒ RESPOSTA

clear all var1 =
clc 7.8000

var1=7.8 var2 =
var2=2.33 2.3300

resp=var1/var2 resp =
3.3476

OBSERVAÇÕES:

I) Os nomes das variáveis não podem começar por caracteres numéricos.

II) Não se cria variáveis com letras separadas por espaço ou utilizando acento gráfico.

III) Letras maiúsculas são diferentes de letras minúsculas (o MATLAB é sensitivo).

IV) Não se deve criar variáveis com o nome de comandos já definidos pelo MATLAB.

2.4 Resolução de Expressões

As expressões são calculadas da esquerda para a direita utilizando a ordem fundamental da
matemática conforme o diagrama que segue:
 9

Potenciação

Multiplicação ou Divisão
Obedecendo a ordem de procedência

Adição e Subtração
Obedecendo a ordem de procedência

Na Janela de Comandos digite, e verifique a ordem de procedência

a=2
b=5
c=4
a+b*c^a

OBSERVAÇÕES:

I) Utiliza-se o sinal de parênteses ( ) para alterar a ordem de procedência das equações.

II) Não se utiliza colchetes ou chaves na resolução de expressões, somente parênteses.

OBSERVAÇÃO:
Para salvar a maior parte do conteúdo da Janela de Comandos basta digitar na mesma:

diary teste1.dat (onde teste1.dat é um nome de arquivo, este nome pode ser alterado)
diary off (desativa o arquivo, caso você não queira continuar salvando)
diary on (retorna a salvar no arquivo teste1.dat)

2.5 Exercícios Propostos

4
1) x = 3 + 5 − 2 +8×2
3

» clear all
» clc
» x=3+5-2*4/3+8*2 ⇒Resolução direta

x=
21.3333

OBSERVAÇÃO:
Note que ao entregar um resultado em um relatório ou trabalho os números em português são
escritos com vírgula e não com ponto, contudo para digitá-los no MATLAB você deverá usar o
ponto já que o software é na língua inglesa. Note este fato nos exercícios que seguem.
 10

3 + 2,34
2) x = 2 −
52

» clear all
» clc

» x=2-(3+2.34)/(5^2) ⇒Resolução direta

x=
1.7864

3 − 52 52,5
3) resp = +
−3 + 6 − 5 72
b
a
⇒MATLAB EDITOR
clear all
clc

a=(3-5^2)/(-3+6^-5) ⇒Resolução por partes

b=52.5/(7^2)

resp=a+b

⇒RESPOSTA
a=
7.3336

b=
1.0714

resp =

8.4051

2 2,21 2,43 + 21
4) y = +
3 − 2,3 3 + 21,44 − 3

a b

⇒MATLAB EDITOR
clear all
clc
a=(2^2.21)/(3-2.3); ⇒Resolução por partes
 11
b=(2.43+21)/(3+2^(1.44-3));

y=a+b

⇒RESPOSTA
y=
13.6264

OBSERVAÇÃO:
O ponto e vírgula (;) no final suprime a visualização da operação precedida a ele.

5) Se a = 4 e b = 3,89 encontre o valor de y da expressão abaixo:

5a − 3b 1067,44
y= −
2a − b 3( 2a − b)

⇒MATLAB EDITOR y1 y2
clear all
clc

a=4;
b=3.89;
y1=(5*a-3*b)/(2*a^(-b));
y2=1067.44/(3^(2*a-b));

y=y1-y2

⇒RESPOSTA

y=
903.7586

− 25,11 + 7,09
6) Encontre o valor de y =
1,2
3,7 ×
3

2
2,3 + 7,83  103,44 
7) Encontre o valor de z = − 
234,1 − 2 3,3  2,33 × 5,76 

3,21
8) Sendo a = 3; b = -2.1; c=3,4; resolva a expressão E = 2,4 a − b +
2c − a
 12
a + 2b 2a + b
−
− 3,32 2 4,56
9) Sendo a = 23; b = -51,2; resolva: resp =
(ab − 3a )

10) Sendo uva = 3,4 ; carro = 4; jacare = -2,11; resolva:

uva + carro jacare 5(7 + jacare) − 2
E=2 −
2,33carro uva (3,21× carro) + carro

11) Se a = 3 e b = -2,89 calcule o valor de X, sendo:

a − b + 2a 5b
X= −
5b b a − 2b

2,544 7,4 150,2

12) Calcule o valor de R onde R = +2
3 + 3,11
4,21 2,212
0,98

−5 3,72 − 3,4 2
13) Calcule o valor de Z = +
3,71 + 2 2,3 0,98

50m 50n
14) Sendo m = -4,55 e n = -5,44 calcule R = −
3 2
(mn − m n ) (2m − 3n ) 2

15) Seja a = 3 e b = 4, avalie as seguintes expressões:

5a (2a ) −5 a2 4
(a) (b) (c) (d) a 2
2b (a + b ) 2
b −a
4 4
3

16) Calcule:
−3
−2 1  1  1 3 
(a)   (b) − 2   − 
 5  4  2  5 4 

17) Calcule o valor numérico da expressão 7a2b + 4ab2 + 3a3 + (2ab - b)2 onde a = 3 e b = 2.

18) Calcule o valor do trinômio x2 – 5x + 6 para (a) x = 2 e (b) x = - 10

19) Divida x2 – 5x + 6 por x – 2 para (a) x = 4 e (b) x = 3

− 10 + 3,45 2,71 − 1,4 3

20) Calcule o valor de R = +
2,81 + 53 ,5 0,5
 13
2.6 Mensagens de ERRO

5n + m 3
Sendo n = -2,8 e m = 9,1 resolva R =
m 2 − 7,78

a) Exemplo 1
⇒ MATLAB EDITOR (salvando o arquivo como exe1.m)
clear all
clc

n=-2.8
m=9.1
Falta parêntesis no
R=(5*n+m^3)/(m^2-7.78 denominador

⇒RESPOSTA Arquivo executado

» exe1
??? (5*n+m^3)/(m^2-7.78
A barra vertical
| indica onde está
Sinal o erro
indicativo
de erro

A closing right parenthesis is missing. Mensagem de erro

Check for a missing ")" or a missing operator.

Error in ==> C:\temp\exe1.m Mostra a localização do

On line 7 ==> R=(5*n+m^3)/(m^2-7.78 arquivo onde está o erro

Indica a linha onde

está o erro

b) Exemplo 2
⇒MATLAB EDITOR
clear all
clc
Falta operador de
n=-2.8 potenciação (^)
m=9.1

R=(5*n+m^3)/(m2-7.78)

⇒RESPOSTA

n=
-2.8000

m=
 14
9.1000

??? Undefined function or variable 'm2'.

Error in ==> C:\temp\exe1.m

On line 7 ==> R=(5*n+m^3)/(m2-7.78)

c) Exemplo 3
⇒MATLAB EDITOR
clear all
clc
Falta parêntesis no
n=-2.8 denominador
m=9.1

R=(5*n+m^3)/m^2-7.78)

⇒RESPOSTA

» exe1
??? (5*n+m^3)/m^2-7.78)
|
Missing operator, comma, or semi-colon.

Error in ==> C:\temp\exe1.m

On line 7 ==> R=(5*n+m^3)/m^2-7.78)

2.7 Respostas dos Exercícios Propostos

6) y = -18.5441
7) z = -57.2782
8) E = 87.7563
9) resp = -0.0067
10) E = 0.2340
11) X = -2.8582
12) R = 122.7457
13) Z = –8.5791
14) R = 5.1370
15) (a) 1.8750 (b) 2.6245e-006 (c) 0.0514 (d) 12
16) (a) -15.6250 (b) 0.8000
17) 355
18) (a) 0 (b) 156
19) (a) 1 (b) 0
20) -0.0912
 15

3. FUNÇÕES MATEMÁTICAS ELEMENTARES

3.1 Funções Matemáticas Elementares: (Elementary math functions)

x ⇒ abs(x)
sen( x ) ⇒ sin(x)
cos( x ) ⇒ cos(x)
tg( x ) ⇒ tan(x)
arcsen(x ) ⇒ asin(x)
arccos(x ) ⇒ acos(x)
arctg( x ) ⇒ atan(x)
cos sec( x ) ⇒ csc(x)
sec( x ) ⇒ sec(x)
cot g( x ) ⇒ cot(x)
ar cos sec( x ) ⇒ acsc(x)
arc sec( x ) ⇒ asec(x)
arc cot g ( x ) ⇒ acot(x)
⇒ exp(x)
ex
log e ( x ) = ln(x ) ⇒ log(x)
log10 ( x ) ⇒ log10(x)
log 2 ( x) ⇒ log2(x)
x ⇒ sqrt(x)
a ⇒ (x^b)^(1/a)
x b = x b / a = ( x b )1 / a
π ⇒ pi
n! ⇒ factorial(n) onde n <= 21

OBSERVAÇÕES:
I) O Matlab opera com arcos trigonométricos somente em radianos, ( 180° = π rad ), ou seja, basta
multiplicar o ângulo em graus por π e dividi-lo por 180. Exemplo:

52° ⇒ 52 * pi /180

II) Pode-se acessar a lista completa de funções matemáticas elementares utilizando a janela de
ajuda (help window), através da:

a) Barra de Comandos:
 16

Exemplo:
Ir em:
Help, Help Window no lugar de MATLAB Help Topics digitar abs após digite Home
comandos
alternativos

explicação

b) Barra de Ferramentas:

Digitando no símbolo de ? você irá para a mesma janela anterior.

c) Janela de Comandos – Comando help on line

Pode-se acessar os comandos de ajuda do MATLAB utilizando o help on line na Janela de
Comandos digitando:

help + comando desejado + <enter>

 17

digitado

explicação

comandos
alternativos

3.2 Outras Funções Especiais

ceil(x) ⇒ Arredonda o número (x) para o inteiro mais próximo na direção de

mais infinito.

» ceil(2.7)
ans =
3

floor(x) ⇒ Arredonda o número (x) para o inteiro mais próximo na direção de

menos infinito.

» floor(2.4)
ans =
2

round(x) ⇒ Arredonda o número (x) para o inteiro mais próximo.

» round(2.4)
ans =
2

» round(2.7)
ans =
3

fix(x) ⇒ Arredonda o número na direção do zero.

» fix(2.4)
ans =
2

rem(x,y) ⇒ Resto de x ÷ y
 18
» rem(4,2)
ans =
0

» rem(4,3)
ans =
1

rand ⇒ Gera números aleatórios entre 0 e 1.

» rand
ans =
0.6068 !! Cuidado sua resposta não dará o mesmo número que o meu

3.3 Exercícios Propostos

7,4 + 32,8
1) Calcule x =
33,1

⇒MATLAB EDITOR
clear all
clc

x=sqrt(7.4+3^2.8)/33.1

⇒RESPOSTA
x=
0.1629

OBSERVAÇÕES (já estão na aula 1, mas vale destacar novamente):

I) Não salvar arquivos com nome de comandos, por exemplo: cos, log, sqrt...

II) Inicie um novo arquivo.m sempre que iniciar um novo cálculo ou programa.

III) Salve as modificações sempre que alterar um arquivo.m antes de executá-lo.

a

2) resp =
(
log 3,22 + 9,7 2 )
e 2, 4 + 5

b
 19
⇒MATLAB EDITOR ⇒RESPOSTA
clear all a=
clc 1.9882

a=log10(3.22+9.7^2) b=
b=exp(2.4)+5 16.0232

resp = a/b resp =

0.1241

3) Sendo a = 2,3 e b = 265°, encontre o valor de:

E1

10 2 × 3 3a
E=
cos(b)

E2

⇒MATLAB EDITOR ⇒RESPOSTA

clear all a=
clc 2.3000

b=
a=2.3 4.6251
b=265*pi/180
E1 =
E1=10^2*(3^(a))^(1/3) 232.1630
E2=cos(b)
E2 =
E=abs(E1/E2) -0.0872

E=
2.6638e+003

é equivalente a 2,6638.103

y1 y3

7 8,8
3,4 log(300 2,4 )
4) y = +
20(sen 33° − cos 44,5°) e1,5

y2
 20

⇒MATLAB EDITOR
clear all
clc

y1=sqrt((7^8.8)/3.4);
y2=20*(sin(33*pi/180)-cos(44.5*pi/180));
y3=log10(300^2.4)/exp(1.5);

y=y1/y2+y3

⇒RESPOSTA
y=
-839.6327

(1 + cos 65°)2 ln(2,8259 )

5) Calcule z = ×
(1 − cos 65°)2 3, 4

16,5 + e 2 tan(75°) 9,12 2,5

6) Avalie x = + 10 −
4
10 ln(6 ) cos 2 (20º ) + sen 2 (88°) 5,3 − 3,10,9
2

cos2(x)⇒ cos(x)^2 = (cos(x))^2 = cos(x)*cos(x)

1000 cos 25,5° − log(4 2,9 + 8,21 )

7) Calcule a =
7 + e 4,56

8) Dados a = 4,56 ; b = −98,5 ; c = 176° e d = 98,5° ; encontre o valor da expressão:

cos(c).sen(d ) log(e 3,9 )

E=
4
ab − 5b − a a / 5

9) Encontre o valor absoluto de X, sendo a = 5,32 , b = 5.10 − 3 e

3,44a + log(3,54 8,6a ) 5b 2
X= +
b 2 sec(44°) ba

2
 3a 2 − 5 − b 5 a 2 
10) Sendo a = 9, b = -4,21 e c = 0,98; calcule o valor de R = 
3 + (e ) 
 10c 
 21

[ ]
2
 cos sec(β + δ) 
11) Sendo β = 23° e δ = 107,5°, calcule o valor de α = −7,125  × sen 3 (2δ)
 cot g (2β) 
12) Encontre o inteiro mais próximo de E utilizando somente comandos do MATLAB
3
 {log (7,2)}2  e3
E= 2
 + (Utilize o help on line para log 2 ou use a aula2)
 cos(11,3°)  2
5,9 + e

13) Verifique se a função log10(x) é crescente ou decrescente.

14) Verifique se a função 2x é crescente ou decrescente.

15) Verifique se a função (1/3)x é crescente ou decrescente.

16) Resolva as expressões.

(a) -5+-2= (b) -5+8= (c) |-1/4+3||

17) Calcule
3−4 7 2
(a) (b) (c)
2− 5 3− 2 3
4 −3 2

3.4 Variáveis Especiais

a) O MATLAB trabalha com números complexos no formato ⇒ a + bi

i ou j ⇒ i = j = −1

Exemplos:
»i
ans =
0 + 1.0000i ⇒ Equivalente a −1
»j
ans =
0 + 1.0000i ⇒ Equivalente a −1

» 2 + sqrt(-45)

ans = parte real ⇒ 2.0000

2.0000 + 6.7082i
parte imaginária ⇒ 6.7082

b) Inf ⇒ infinito (divisão por zero)

 22
Exemplo:

» 1/0
Warning: Divide by zero.
ans =
Inf

c) Nan ⇒ Valor não numérico (Not a Number – Não é um número)

Exemplo:

» 0/0
Warning: Divide by zero.
ans =
NaN

3.5 Respostas dos Exercícios Propostos

5) z = 16.6888
6) x = -175.7809
7) a = 88.4543
8) E = 0.5127
9) X = 2.1853e+008
10) R = 30.1692
11) alpha = 2.4934
12) E = 567 (utilizar o comando ‘round’)
13) função crescente
14) função crescente
15) função decrescente
16) (a) 7 (b) 3 (c) 2.7500
17) (a) 9.6072 (b) 4.4142 (c) 6.1072
 23
4. FORMATOS DE VISUALIZAÇÃO NUMÉRICA
Servem para formatar a visualização dos resultados numéricos sem que haja alteração da
representação interna de um número o qual será usado sempre com o maior formato.

format ⇒ 4 dígitos após ponto decimal (default - padrão);

format short ⇒ idem a format;
format long ⇒ 16 dígitos incluindo o ponto decimal;
format short e ⇒ 5 dígitos + expoente;
format long e ⇒ 16 dígitos + expoente;
format short g ⇒ 5 dígitos;
format long g ⇒ 16 dígitos sem expoente;
format bank ⇒ formato para uso em cálculo com moedas;
format rational ⇒ formata a saída utilizando números racionais (frações).

2
1) Faça o cálculo da expressão x = com cada um dos formatos para verificar as diferenças
3
4− 2

⇒ MATLAB EDITOR ⇒ RESPOSTA

x = 11.5482
clear all x = 11.5482
clc
x=2/(4^(1/3)-2^(1/2));
x = 11.54817824703048
x = 1.1548e+001
format x = 1.154817824703048e+001
x x = 11.548
x = 11.5481782470305
format short
x
x = 11.55
x = 5393/467
format long
x

format short e
x

format long e
x

format short g
x

format long g
x

format bank
x

format rational
x
 24
Também pode-se modificar a visualização dos valores numéricos utilizando a Janela de Comandos
⇒ File ⇒ Preferences ⇒ General

OBSERVAÇÕES:

I) Quando o resultado é um número inteiro o MATLAB apresenta-o como inteiro;

II) Na Janela de Comandos em File ⇒ Preferences ⇒ Comand Window Font pode-se modificar a
formatação dos caracteres da Janela de Comandos (Veja detalhes na Figura anterior).

2) Sendo a = 3,44; b = -561.10-2 e c = 53°, encontre o valor da expressão abaixo utilizando o

formato com 16 dígitos sem expoente:

10a ⋅ tg (c)
z=
ab + a 2b − a

⇒ MATLAB EDITOR ⇒ RESPOSTA

a=
clear all 3.44
clc b=
format long g -5.61
c=
a = 3.44 0.925024503556995
b = -561e-2
c = 53*pi/180 z1 =
45.6503418637421
z1 = 10*a*tan(c)
z2 = a*b+a^(2*b-a) z2 =
-19.2983999863839

z = z1/z2 z=
-2.36549879243621
 25

4.1 Exercícios Propostos

b 3 − e 2 + 2a
1) Sendo a = 32 e b = −12 , encontre o valor de C = .
a+b

3m − n + e − m
2) Sabendo que m=-2,2 e n=-7,8, encontre o valor de E = .
n−m
20 +
2

3) Resolva o exercício 9 utilizando o formato de 16 dígitos com expoente.

4) Calcule:
3−4 7
(a) (b)
2− 5 3− 2
2 2+ a
(c) (d) para a = 9
3
4 −3 2 3 a
x x
(e) para x = 4 (f) para x = 16
x+ 2 x −3
2+ 3 x+ 8
(g) (h) para x = 10
2 16

5) Julgue as expressões abaixo com V (verdadeira) ou F (falsa), considerando a, b, c e d

números reais quaisquer.

( a + b) b
(a) =a+ ( ) (b) (a + b) 2 = a 2 + b 2 ( )
c c
c c c
(c) a2 − b2 = a − b ( ) (d) = + ( )
( a + b) a b
(e) (-2)2 = - 4 ( ) (f) (-2)2 = 4 ( )
( a + b)
(g) =1+ b ( ) (h) a 2 = a ∀ a ∈ IR. ( )
a

6) No triângulo retângulo, que segue, calcule o valor de x, onde

B = 30o sen(30o) = 0,5 cos(30o) = 0,86 tg(30o) = 0,57
C

B
A 2
 26

7) Observe o triângulo retângulo da figura a seguir e obtenha o valor de:

(a) senB C
(b) cosC
(c) tgB
(d) senC 3 cm
(e) tgC
(f) cotg(C) A 4 cm B

8) Observe o triângulo retângulo da figura a seguir e dê o valor de:

(a) senB C
(b) cosC
(c) tgB
(d) senC 4 cm
(e) tgC
(f) cotg(C)
A 5 cm B

9) Num triângulo retângulo ABC os catetos são b = 15 cm e c = 20 cm . Calcule o seno, o cosseno e

a tangente dos ângulos B e C, que são opostos respectivamente aos lados b e c.

10) Num triângulo retângulo os catetos medem 1 cm e 2 cm . Qual é o valor do seno do menor
ângulo desse triângulo?

11) Num triângulo ABC, retângulo em A , a hipotenusa é a = 25 cm e cos(B) = 0,96. Calcule o

perímetro do triângulo.

12) Determinar os catetos de um triângulo de hipotenusa 10 cm onde um dos ângulos mede 30º.

13) Calcule a distância entre os pontos dados. Use a representação em formato racional e
formato de 16 dígitos incluindo o ponto decimal.
(a) (1, 2) e (2, 3) (b) (0, 1) e (1, 3) (c) (-1, 2) e (0, 1)

4.2 Respostas dos Exercícios Propostos

1) C = 3.2000 + 2.0829i
2) X = 23.7730 - 1.9890i
3) E = 2.185326839156010e+008
4)
» format short
» (sqrt(3)-4)/(2-sqrt(5))
ans =
9.6072

» 7/(3-sqrt(2))
ans =
 27
4.4142

» 2/(4^(1/3)-2^(1/3))
ans =
6.1072

» a=9;(2+sqrt(a))/(3*sqrt(a))
ans =
0.5556

» x=4;x/(sqrt(x)+sqrt(2))
ans =
1.1716

» x=16;x/(sqrt(x)-3)
ans =
16

» (sqrt(2)+sqrt(3))/sqrt(2)
ans =
2.2247

» x=10;(sqrt(x)+sqrt(8))/sqrt(16)
ans =
1.4977

5) (a) F (b) F (c) F (d) F

(e) F (f) V (g) F (h) F

6) » x = tan(30*pi/180)*2
x=
1.1547

7)
» h=sqrt(4^2+3^2)
h=
5

» senB = 3/h
senB =
0.6000

» cosC = 3/h
cosC =
0.6000

» tanB = 3/4
tanB =
0.7500

» senC = 4/h
senC =
 28
0.8000

» tanC = 4/3
tanC =
1.3333

» cotanC=1/tanC
cotanC =
0.7500

8) » h=sqrt(4^2+5^2)
h=
6.4031

» senB=4/h
senB =
0.6247

» cosC=4/h
cosC =
0.6247

» tanB=4/5
tanB =
0.8000

» senC=5/h
senC =
0.7809

» tanC=5/4
tanC =
1.2500

» cotanC=1/tanC
cotanC =
0.8000

9) » b=15
b=
15

» c=20
c=
20

» h=sqrt(15^2+20^2)
h=
25

» senB=15/h
senB =
 29
0.6000

» cosB=20/h
cosB =
0.8000

» tanB=15/20
tanB =
0.7500

» senC=20/h
senC =
0.8000

» cosC=15/h
cosC =
0.6000

» tanC=senC/cosC
tanC =
1.3333

10) » h=sqrt(1+2^2)
h=
2.2361

» senB=1/h
senB =
0.4472

11)
h=25;
» caa=h*0.96
caa =
24

» cao=sqrt(h^2-caa^2)
cao =
7

» perimetro=caa+cao+h
perimetro =
56

12) h=10; ca=cos(30*pi/180)*h

ca =
8.6603

» co=sin(30*pi/180)*h
co =
5.0000
 30

13) format short

» d=sqrt((2-1)^2+(3-2)^2)
d=
1.4142

» d=sqrt((1-0)^2+(3-1)^2)
d=
2.2361

» d=sqrt((0+1)^2+(1-2)^2)
d=
1.4142

» format rational
» d=sqrt((2-1)^2+(3-2)^2)
d=
1393/985

» d=sqrt((1-0)^2+(3-1)^2)
d=
2889/1292

» d=sqrt((0+1)^2+(1-2)^2)
d=
1393/985

» format long
» d=sqrt((2-1)^2+(3-2)^2)
d=
1.41421356237310

» d=sqrt((1-0)^2+(3-1)^2)
d=
2.23606797749979

» d=sqrt((0+1)^2+(1-2)^2)
d=
1.41421356237310
 31
5. VETORES
Todos os cálculos considerados até este ponto envolveram números individuais chamados
escalares. As operações escalares são a base da matemática. Quando se deseja efetuar uma
operação em mais de um número de uma única vez, operações escalares repetidas são
inconvenientes e demoradas. Para resolver este problema o MATLAB utiliza-se de vetores e
matrizes.

5.1 Construção de vetores LINHA

A) Cria um vetor linha com elementos especificados entre [ ].

1) Crie o vetor x = 1 [ 2,4 2 −π ]

» x=[1 2.4 sqrt(2) -pi]

x=
1.0000 2.4000 1.4142 -3.1416

r r
2) Crie os vetores a = 5i + 3 j − k e b = −3i − k

» a=[5 3 -1]

a=
5 3 -1

» b=[-3 0 -1]

b=
-3 0 -1

OBSERVAÇÃO:
Pode-se utilizar a vírgula (,) para separar os elementos do vetor.

 3
3) x =  − 10 e 2,31 2
 4 

» x=[-10,exp(2.31),sqrt(2),3/4]

x=
-10.0000 10.0744 1.4142 0.7500

B) Cria um vetor linha x começando em início com incrementos unitários e pára em

fim.
 32
x=(início:fim) ou x=início:fim

4) Crie um vetor com início no algarismo –5, com fim no algarismo 2 e com incremento unitário.

» x=(-5:2)

x=
-5 -4 -3 -2 -1 0 1 2

5) Crie um vetor com início no algarismo 20, com fim no algarismo 50 e com incremento unitário.

» x=(20:50)

x=
Columns 1 through 12

20 21 22 23 24 25 26 27 28 29 30 31

Columns 13 through 24

32 33 34 35 36 37 38 39 40 41 42 43

Columns 25 through 31

44 45 46 47 48 49 50

C) Cria um vetor linha x começando em início com incremento especificado e pára em

fim;

x=(início:incremento:fim) ou x=início:incremento:fim

6) Crie um vetor com início no algarismo –5, com fim no algarismo 20 e com incremento igual a 5.

» x=(-5:5:20)

x=
-5 0 5 10 15 20

7) Crie um vetor com início no algarismo 3, com fim no algarismo 19 e com incremento igual a 3.

» x=(3:3:19) ⇒ Termina o vetor no algarismo mais próximo de 19 contido entre o início, 3

e o fim, 19.
x=
3 9 12 15 18
 33

8) Exemplo crie um vetor com início no algarismo 3, com fim no algarismo -15 e com incremento
igual a -1,5.

» x=(3:-1.5:-15)

x=

Columns 1 through 7

3.0000 1.5000 0 -1.5000 -3.0000 -4.5000 -6.0000

Columns 8 through 13

-7.5000 -9.0000 -10.5000 -12.0000 -13.5000 -15.0000

D) Cria um vetor linha x com n-elementos começando em início e terminando em fim.

x=linspace(início, fim , n° de elementos)

9) Crie um vetor com 15 elementos iniciando em 0 e terminando em 2π.

» x=linspace(0,2*pi,15)

x=

Columns 1 through 7

0 0.4488 0.8976 1.3464 1.7952 2.2440 2.6928

Columns 8 through 14

3.1416 3.5904 4.0392 4.4880 4.9368 5.3856 5.8344

Column 15

6.2832

OBSERVAÇÃO:
Quando não é informado o n° de elementos ao comando linspace, o vetor é construído com 100
elementos, (ver exemplo na página seguinte).

10) Crie um vetor com 100 elementos iniciando em 1 e terminando em 80.

» x=linspace(1,80)

E) Cria um vetor linha logarítmico x com n-elementos começando em 10início e

terminando em 10fim.
 34

x=logspace(início, fim, n°de elementos)

11) Crie um vetor logarítmico com 10 elementos iniciando em 100 e terminando em 102.

» x=logspace(0,2,10)

x=
Columns 1 through 7
1.0000 1.6681 2.7826 4.6416 7.7426 12.9155 21.5443

Columns 8 through 10
35.9381 59.9484 100.0000

5.2 Construção de vetores COLUNA

A) Construa o vetor coluna x dado por:

 1 
 
x = − 2,5
 11 

A.1) Com valores A.2) Separados por (;)

especificados
» x=[1; -2.5; 11] ⇒ O (;)
» x=[1 <enter> executa a mudança
-2.5 <enter> de linha
11] <enter>
x=
x= 1.0000
1.0000 -2.5000
-2.5000 11.0000
11.0000

5.3 Transposta de um vetor

(' ⇒ apóstrofo - apóstrofe)

» a=[1 -2 -3 4.5]

a=
1.0000 -2.0000 -3.0000 4.5000

» b=a' ⇒ (') aposto

b=
1.0000
-2.0000
-3.0000
 35
4.5000

» c=[1; 5; 9]

c=
1
5
9

» d=c'

d=
1 5 9

5.4 Operações com vetores

A) Operações vetor - escalar

A.1) Adição e subtração

» a=[ 1 4 -7]; ⇒ vetor

» b=5; ⇒ escalar

» x=a+b ⇒ x = [1 4 − 7] + 5
x=
6 9 -2

» y=a-b ⇒ y = [1 4 − 7] − 5
y=
-4 -1 -12

» z=-b-a ⇒ z = −5 − [1 4 − 7]
ans =
-6 -9 2

» w=3-a ⇒ w = 3 − [1 4 − 7]
w=
2 -1 10

A.2) Multiplicação e divisão

» a=[ 1 3 2.5]; ⇒ vetor

» b=4; ⇒ escalar
 36
» x=b*a ⇒ x = 4 × [1 3 2,5]
x=
4 12 10

» y=-b*a ⇒ y = −4 × [1 3 2,5]
y=
-4 -12 -10

» z=3*(-a) {
⇒ z = 3 × − [1 3 2,5] }
z=
-3.0000 -9.0000 -7.5000

[1 3 2,5]
» w=a/b ⇒ x= Divisão vetor/escalar
4
w=
0.2500 0.7500 0.6250

4
» m=b/a ⇒ m= Divisão escalar/vetor não é definida
[1 3 2,5]
??? Error using ==> /
Matrix dimensions must agree. As dimensões das matrizes
devem concordar
concordantes
Divisão escalar/vetor só é
» m=b./a definida quando utilizado o ponto
(.) após o escalar.
ans = O ponto indica que a divisão será
4.0000 1.3333 1.6000 realizada dividindo o escalar por
todos os elementos do vetor, um
de cada vez.

B) Operações vetor – vetor (elemento por elemento)

Quando dois vetores possuem mesma dimensão, as operações elementares de adição, subtração,
multiplicação, divisão e potenciação aplicam-se de elemento em elemento.

B.1) Adição e subtração

» a=[1 3 6];
» b=[2 -1.5 2];
» c=[1 2 3 4];

» x=a+b ⇒ x = [1 3 6] + [2 − 1,5 2]
x=
3.0000 1.5000 8.0000
 37

» y=a-b ⇒ y = [1 3 6] − [2 − 1,5 2]
y=
-1.0000 4.5000 4.0000

» z=2*a-b ⇒ z = 2 × [1 3 6] − [2 − 1,5 2]
z=
0 7.5000 10.0000

» m=a-3*b ⇒ m = [1 3 6] − 3 × [2 − 1,5 2]
m=
-5.0000 7.5000 0

» n=a+c ⇒ n = [1 3 6] + [1 2 3 4]

??? Error using ==> + Vetores (a) e (c) possuem

Matrix dimensions must agree. dimensões diferentes:

Vetor a ⇒dimensão 1 x 3
Vetor c ⇒dimensão 1 x 4

B.2) Multiplicacão e divisão pontuada

» a=[1 3 6];
» b=[2 -1.5 2];
» x=a.*b ⇒ x = [1 3 6] × [2 − 1,5 2]
⇒ multiplicação de cada elemento de a com seu
respectivo elemento em b
x=
2.0000 -4.5000 12.0000

» y=2*a./b ⇒ x = 2 × [1 3 6] / [2 − 1,5 2]
⇒ divisão de cada elemento de a com seu
respectivo elemento em b
y=
1 -4 6

» m=3*a.*b x = 3 × [1 3 6] × [2 − 1,5 2]
m=
6.0000 -13.5000 36.0000

B.3) Potenciação vetorial pontuada

 38

» a=[2 -4 9];
» b=[1 4 1.5];

» x=a.^b ⇒ x = [2 − 4 9][1 4 1,5]

⇒ potenciação de cada elemento de a com seu
respectivo elemento em b
x=
2 256 27

» y=2*a.^2 ⇒ x = 2 × [ 2 − 4 9] 2
⇒ mesmo que 2*a.*a
y=
8 32 162

⇒ z = [1 4 1,5] − [2 −4
2
» z=b.^2-a 9]
⇒ mesmo que b.*b-a
z=
-1.0000 20.0000 -6.7500

» m=3*a-b.^(1/2) ⇒ m = 3 × [2 −4 9] − [1 +4 1,5]1/ 2
m=
5.0000 -14.0000 25.7753

⇒ m=2
[2 4 9]
» n=2.^a
n=
4.0000 0.0625 512.0000

C) Resumo das operações vetoriais

Sendo a = [a1 a2 a3 ... an]

b = [b1 b2 b3 ... bn]
c = escalar

Adição ou subtração entre um vetor e um escalar

a ± c = [a1 ± c a2 ± c a3 ± c ... an ± c]

Multiplicação ou divisão entre um vetor e um escalar

a × c = [a1 × c a2 × c a3 × c ... an × c]
 39
a / c = [a1 / c a2 / c a3 / c ... an / c]

Divisão pontuada entre um escalar e um vetor

c./ a = [c/a1 c/a2 c/a3 ... c/an]

Adição ou subtração entre vetores

a ± b = [a1 ± b1 a2 ± b2 a3 ± b3 ... an ± bn]

Multiplicação ou divisão entre vetores

a.× b = [a1 × b1 a2 × b2 a3 × b3 ... an × bn]

a./ b = [a1 / b1 a2 / b2 a3 / b3 ... an / bn]

Potenciação envolvendo vetores

a.^ b = [a1 ^ b1 a2 ^ b2 a3 ^ b3 ... an ^ bn]

a.^ c = [a1 ^c a2^c a3^c ... an^c]

c.^ a = [c^a1 c^a2 c^a3 ... c^an]

5.5 Exercícios propostos

1) Sendo a=[1 4 6 8] e b=[-1 –2 –3 2], resolva:
x1

2a + b
x=
ab

⇒MATLAB EDITOR x2
clear all
clc

a=[1 4 6 8];
b=[-1 -2 -3 2];

x1=2*a+b;
x2=a.*b;

x=x1./x2

⇒RESPOSTA
x=
-1.0000 -0.7500 -0.5000 1.1250

2) Sendo a=[1 2 4] e b=[0,5 1,5 7], resolva:

 40

12,5 log(ab + b)
x=
5a − 1,5b

n
⇒MATLAB EDITOR
clear all
clc

a=[1 2 4];
b=[0.5 1.5 7];

m=12.5*log10(a.*b+b);
n=5*a-1.5*b;

x=m./n

⇒RESPOSTA
x=
0 1.0536 2.0317

3) Sendo x=[1 6 –2], y=[2 6 9] e z=[-2 2 1.5], resolva:

w1 w2

xyz xy − zx
w= +
2x 2 ln(2 x + 3,2 y)

w3
⇒MATLAB EDITOR ⇒RESPOSTA
clear all w=
clc -1.5301 6.7120 7.3531

x=[1 6 -2];
y=[2 6 9];
z=[-2 2 1.5];

w1=(x.*y.*z)./(2*x);
w2=abs(sqrt(x.*y-z.*x));
w3=2*log(2*x+3.2*y);
w=w1+w2./w3

4) Sendo a=[1 4 6 1] e b=[0.5 -0.5 8 1], resolva:

 41
5 a 2 b + ab 2 − 4 b
r= +
ab 3 (a ⋅ b ⋅ tg10°) 4

⇒MATLAB EDITOR
clear all
clc

a=[1 4 6 1];
b=[0.5 -0.5 8 1];

r1=5./(a.*b);
r2=(a.^2).*b+a.*(b.^2)-4.^b;
r3=((a.*b*tan(10*pi/180)).^(4)).^(1/3);

r=r1+(r2./r3)

⇒RESPOSTA
r=

1.0e+003 *

-0.0219 -0.0326 -3.7605 -0.0152

OBSERVAÇÕES:

I) A resposta acima é equivalente a:

1.0⋅103×[-0.0219 –0.0326 –3.7605 –0.0152]

II) Para uma visualização completa dos algarismos dos elementos pertencentes ao vetor r pode ser
utilizada a formatação (format long g):

r=
Columns 1 through 3

-21.8560295718585 -32.6020616677436 -3760.49503161746

Column 4

-15.2273179415593

5) Sendo a=[1 2 4 6 8], calcule o valor de R dada pela expressão abaixo.

2
2a 4,5 + 3a
R= +
aa log(10a 2 )

6) Determine o valor de Z na expressão abaixo sabendo que m = [− 1 e 2 5] e

n = [log(80) cos(12°) π] .
 42
2m + 3n mn − m
Z= ⋅
m−n 3n

7) Determine o valor K sabendo que x=[1,5 1,5 1,5 2,8] e y=[0,5 2 -4 -9].

x3 x−y
K= −
y2 3 ( 2 x − y) 2

8) Sendo c = [5 88 − 2,5] e d = [9,9 − 23,5] , encontre o valor de X.

cd − 2c − 3d 4,56
X= +
4,56 cd − 2c − 3d

9) Sendo os vetores a=[1 5 7 0]; e b=[1 4 –2,5 -1]; resolva a expressão abaixo:
2 2 (1+ ab)
ab − a e
x= +
log(5a + 10) 4 3b
10) Sendo x=[-2,5 2] e y=[3 –4,5]. Encontre o valor de A.

ln(3xy + y 2 ) x2
A= +
xy y3

11) Resolva a expressão abaixo se o vetor a inicia no algarismo 3 e termina no algarismo 14 em

incrementos de 2,2; e o vetor b inicia no algarismo 20 e termina no algarismo 13,95 em incrementos
de –1,21.
2
5a − 3b a
150,5 × −
e=
ab ab − a b
a−b

5.6 Cálculo vetorial

A) Produto escalar (interno)

r r
Sendo os vetores: a = a1i + a 2 j + a 3 k e b = b1i + b 2 j + b 3k
r r
Por definição ⇒ a • b = a1b1 + a 2 b 2 + a 3 b 3

a=[a1 a2 a3]
No MATLAB ⇒ b=[b1 b2 b3]
dot(a,b)
 43
B) Produto vetorial (externo)
r r
Sendo os vetores: a = a1i + a 2 j + a 3 k e b = b1i + b 2 j + b 3k
i j k
r r a a3 a a3 a a2
Por definição ⇒ a × b = a1 a2 a3 = i 2 −j 1 +k 1
b2 b3 b1 b 3 b1 b 2
b1 b 2 b3

a=[a1 a2 a3]
No MATLAB ⇒ b=[b1 b2 b3]
cross(a,b)

C) Norma de um vetor (módulo de um vetor)

r
Sendo o vetor: a = a1i + a 2 j + a 3 k
r
Por definição ⇒ a = a12 + a 22 + a 32

a=[a1 a2 a3]
No MATLAB ⇒
norm(a)

r r r
12) Dados os vetores a = 3i + 5 j e b = i − j + 8k , encontre a norma de x , dada pela equação
abaixo

r r
r 2( a × b )
x= r r
7 (a • b )
⇒MATLAB EDITOR
clear all
clc

a=[3 5 0];
b=[1 -1 8];

x1=2*cross(a,b);
x2=7*dot(a,b);

x=norm(x1/x2)

⇒RESPOSTA
x=
6.7612
r v
13) Sendo o vetor a = i − 4 j + 2k e b = 1,5i + 2 j − 3k , encontre a norma do vetor W.
 44
r r r r r r
r (a × b ) (a + b ) × ( b − a )
w= r r +
(a • b ) 2

⇒MATLAB EDITOR ⇒RESPOSTA

clear all w=
clc 11.7817

a=[1 -4 2];
b=[1.5 2 -3];

w1=cross(a,b)/dot(a,b);
w2=cross(a+b,b-a); w=norm(w1+w2/2)

5.7 Exercícios propostos

r r
14) Sendo x = 3i − j − 1,5k e y = 7i − j − 10k , encontre o vetor z dado pela expressão abaixo:
r r
v ( x × 3y )
z= r r
(3x • y)
r r r
15) Sendo x = 3i − 3 j − 5k e y = 7i − 2k , encontre o módulo do vetor w dada pela expressão
abaixo:
r  [2 y • x ]
r r 

w= r r r 
 (5 y) × ( x − y) 

r r
16) Encontre a medida em graus do ângulo θ formado entre os vetores u = (2,0,−3) e v = (1,1,1) .
r r
u•v
cos θ = r r ⇒ θ em radianos
u ⋅ v

r
17) Encontre a medida em graus do ângulo θ formado entre os vetores u = (1,0,1) e
r
v = (−2,10,2) .

18) Calcule a área do paralelogramo ABCD, sendo AB = (1,1,−1) e AD = (2,1,4)

B
C

area = AB × AD
A
D

19) Calcule a área do triângulo ABC, delimitada pelos vetores

 45

AC = (1,1,3) e CB = (−1,1,0) .

20) Calcule o volume do paralelepípedo delimitado pelos vetores

r r r
u = (1,0,1) , v = (0,3,3) e w = (2,1,2) . C
r r r r
volume =| ( u × v ) • w | u

A
r r D
v w
B

r r
21) Calcule o volume do tetraedro ABCD delimitado por: u = (1,1,0) , v = (0,1,1) e
r
w = (−4,0,0) .

1 r r r
volume = | ( u × v ) • w |
6

5.8 Respostas dos exercícios propostos

5) R = 9.5000 8.5541 7.6110 8.8033 10.1562
6) Z = 0.2021 -0.1520 -7.7977
7) K = 12.9571 1.3438 1.2921 1.7043
8) X = -2.7116 29.4748
9) x = 0.1155 114.2277 -3.1734 173.9700
10) A = 0.4398 0.2540
11) e = 6.6397 3.4422 1.7887 0.1570 -3.3853 433.8530

14) z = 0.2297 0.5270 0.1081

15) w = 0.3416
16) teta = 99.2143°
17) teta = 90°
18) area = 7.8740
19) area = 2.3452
20) volume = 3
21) volume = 0.6667
 46

6. MATRIZES

6.1 Construção de matrizes

A) Matriz com elementos especificados entre [ ]

1) Crie a matriz:

 2 π 4
 
x = − 1 − 1 0 
 2,5 5 e 3 
 

A.1) Forma matricial

» x=[sqrt(2) pi 4 <enter>
-1 -1 0 <enter>
2.5 5 exp(3)] <enter>

x=
1.4142 3.1416 4.0000
-1.0000 -1.0000 0
2.5000 5.0000 20.0855

A.2) Elementos separados por espaço e linhas por(;)

» x=[sqrt(2) pi 4; -1 –1 0; 2.5 5 exp(3)]

x=
1.4142 3.1416 4.0000
-1.0000 -1.0000 0
2.5000 5.0000 20.0855

A.3) Elementos separados por (,) e linhas por (;)

» x=[sqrt(2),pi,4;-1,-1,0;2.5,5,exp(3)]

x=
1.4142 3.1416 4.0000
-1.0000 -1.0000 0
2.5000 5.0000 20.0855

6.2 Dimensões de vetores e matrizes

Vetor ⇒ necessita de 1 coordenada para ser definida.
Matriz ⇒ necessita de 2 coordenadas para ser definida.
 47
A) Vetor sendo considerado como um vetor.
[
2) Sendo o vetor a = 1 0 − 5 2]

» a=[1 0 -5 2];

» a(1) ⇒vetor (a) elemento (1)

ans =
1

» a(2) ⇒vetor (a) elemento (2)

ans =
0

» a(3) ⇒vetor (a) elemento (3)

ans =
-5

» a(4) ⇒vetor (a) elemento (4)

ans =
2

B) Vetor sendo considerado como uma matriz.

[
3) Sendo o vetor a = 1 0 − 5 2]

» a=[1 0 -5 2];

» a(1,1) ⇒matriz (a) elemento (1,1)

ans =
1

» a(1,2) ⇒matriz (a) elemento (1,2)

ans =
0

» a(1,3) ⇒matriz (a) elemento (1,3)

ans =
-5

» a(1,4) ⇒matriz (a) elemento (1,4)

ans =
2
C) Matriz sendo considerada como um vetor.

As coordenadas são contadas de cima para baixo e da esquerda para a direita.

 48
 elemento 1 elemento 5 
 
elemento 2 elemento 6
a= 
elemento 3 M M 
 
elemento 4 elemento n

 1 5 0
4) Dada a matriz a =  
− 2 8 1

» a=[1 5 0; -2 8 1];

» a(1) ⇒vetor (a) elemento (1)

ans =
1
» a(2) ⇒vetor (a) elemento (2)
ans =
-2

» a(3) ⇒vetor (a) elemento (3)

ans =
5

» a(4) ⇒vetor (a) elemento (4)

ans =
8

» a(5) ⇒vetor (a) elemento (5)

ans =
0

» a(6) ⇒vetor (a) elemento (6)

ans =
1

» a(7) ⇒vetor (a) elemento (7)

??? Index exceeds matrix dimensions.
(Índice excede a dimensão da matriz)

D) Matriz sendo considerada como uma matriz.

 1 5 0
5) Dada a matriz a =  
− 2 8 1

» a=[1 5 0; -2 8 1];
 49

» a(1,1) ⇒matriz (a) elemento (1,1)

ans =
1

» a(1,2) ⇒matriz (a) elemento (1,2)

ans =
5
» a(1,3) ⇒matriz (a) elemento (1,3)
ans =
0

» a(2,1) ⇒matriz (a) elemento (2,1)

ans =
-2

» a(2,2) ⇒matriz (a) elemento (2,2)

ans =
8

» a(2,3) ⇒matriz (a) elemento (2,3)

ans =
1

» a(3,3) ⇒matriz (a) elemento (3,3)

??? Index exceeds matrix dimensions.
(Índice excede a dimensão da matriz)

6.3 Comandos de dimensão

⇒Retorna um vetor linha (s) , cujo primeiro elemento é o n° de linhas

s=size(A)
de (A) e cujo segundo elemento é o n° de colunas de(A).
⇒Retorna dois escalares (l) e (c) contendo, respectivamente, o n° de
[l,c]=size(A)
linhas e o n° de colunas de (A).
n=length(A) ⇒Retorna a maior dimensão de (A)

6) Utilize os diferentes comandos nos vetores que seguem.

» a=[1 3 5] ⇒ a = [1 3 5]
a=
1 3 5

2 3 

» b=[2 3; 4 7; 9 10] ⇒ b = 4 7

 
9 10
b=
2 3
 50
4 7
9 10

» s=size(a)
s=
1 3 ⇒ 1 linha e 3 colunas

» [l,c]=size(b)
l=
3 ⇒ 3 linhas
c=
2 ⇒ 2 colunas

» n=length(b)
n=
3 ⇒ máxima dimensão do vetor b

6.4 Operações com matrizes

A) Operações com matriz elemento por elemento

Todas as operações vetoriais são validas para matrizes desde que suas dimensões sejam iguais.

 3 3 − 4 1 − 1 8
7) Sendo a =  eb= 
4 − 2 6  2 2 3

» a=[3 3 -4; 4 -2 6]

a=
3 3 -4
4 -2 6
» b=[1 -1 8; 2 2 3]

b=
1 -1 8
2 2 3

 3 3 − 4  1 − 1 8 
» x=a+b ⇒ x=  + 
4 − 2 6  2 2 3
x=
4 2 4
6 0 9
 51
1 − 1 8
» y=b-2 ⇒ y= −2
 2 2 3
y=
-1 -3 6
0 0 1

3 3 − 4 1 − 1 8
» z=2*a-b ⇒ y = 2×   − 2 2 3
 4 − 2 6   
z=
5 7 -16

6 -6 9

 3 3 − 4
» w=-2-a ⇒ w = −2 −  
4 − 2 6 
w=
-5 -5 2
-6 0 -8

 3 3 − 4
4 − 2 6 
» m=a/2 ⇒ m=  
2
m=
1.5000 1.5000 -2.0000
2.0000 -1.0000 3.0000

2
» n=2/a ⇒ n= Divisão escalar/matriz
 3 3 − 4 não é definida
4 − 2 6 
 
??? Error using ==> /
Matrix dimensions must agree.
» n=2./a

ans =
0.6667 0.6667 -0.5000
0.5000 -1.0000 0.3333

2
 3 3 − 4
» p=a.^2 ⇒ p= 
4 − 2 6 
p=
9 9 16
16 4 36
 52
 3 3 − 4 1 − 1 8
» q=a.*b ⇒ q=  × 2 2 3
 4 − 2 6   
q=
3 -3 -32
8 -4 18

1 −1 8
 
» r=2.^b-2 ⇒ r = 2  2 2 3 −2
r=
0 -1.5000 254.0000
2.0000 2.0000 6.0000

 3 3 − 4
1 − 1 8  4 − 2 6 
» s=b.^a ⇒ s=
2 2 3
s=
1.0000 -1.0000 0.0002
16.0000 0.2500 729.0000

B) Operações matriciais

B.1) Multiplicação de matrizes (Linha por Coluna)

A x B =C
m,p p,n m,n

 1 2 2 0 2 
8) Sendo as matrizes A =   eB=  , faça A x B
 − 3 4   1 1 − 2 

A2x2 x B2x3 = C2x3

» A=[1 2; -3 4]; Multiplicação Linha x Coluna

» B=[2 0 2; 1 1 -2];

» C=A*B  1x 2 + 2 x1 1x 0 + 2x1 1x 2 + 2 x ( −2) 

A×B =  
C= − 3x 2 + 4x1 −3x 0 + 4 x1 −3x 2 + 4 x (−2)
 53
4 2 -2
-2 4 -14

B.2) Operações especiais com matrizes quadradas

1 − 6
9) Sendo b =  
3 7 

» b=[1 -6; 3 7];

» b^2 ⇒idem a (b x b)

ans =
-17 -48
24 31

» b^3 ⇒idem a (b x b x b)

ans =
-161 -234
117 73

B.3) Determinante de uma matriz quadrada ⇒ det(X)

 1 6
10) Encontre o determinante da matriz a =  
− 3 9
» a=[1 6; -3 9];
» det(a)
ans =

27

B.4) Inversa de uma matriz ⇒ inv(X) ou X^-1

 1 6
11) Encontre A-1 da matriz A =  
− 3 9
» a=[1 6; -3 9];
» inv(a)

ans =
0.3333 -0.2222
0.1111 0.0370

OBSERVAÇÃO:
Só existe inversa de uma matriz se ela não for singular, ou seja, o determinante dela não for nulo.
 54

11 2
12) Encontre o determinante da matriz b =
5,5 1

» b =[11 2;5.5 1]

b=
11.0000 2.0000
5.5000 1.0000

» inv(b)
Warning: Matrix is singular to working precision.

ans =
Inf Inf
Inf Inf

» det(b)
ans =
0

B.5) Rank de uma matriz (número de linhas não nulas) ⇒ rank(X)

13) Encontre o número de linhas não nulas do determinante de A.

1 1 1
A= 0 0 2
0 0 0
» A=[1 1 1; 0 0 2; 0 0 0];

» rank(A)

ans =
2

B.6) Transposta de uma matriz ⇒ (' aposto, apóstrofe, apóstrofo)

1 2
14) Encontre a matriz b transposta de a = 0 − 11
12 −3

» a=[1 2; 0 -11; 12 -3];

» b=a'

b=
1 0 12
2 -11 -3
 55
6.6 Resolução de sistema de equações do tipo SPD
SPD: Sistema Possível Determinado
Sendo o sistema formado por três equações com três variáveis desconhecidas:

3x +2y -z = 10
-x +3y +2z = -5
x -y -z = 1

Pode-se escrever o sistema de três equações na forma de matrizes:

 3 2 −1  x  10 
 −1 3 2   y  =  −5
    
 1 −1 −1  z   1 
onde
 3 2 − 1 x  10 
 
A = − 1 3 2  X = y B =  −5
 
 1 − 1 − 1 z   1 
 

Então, o sistema de equações pode ser escrito como:

AX = B

Para encontrar o valor das variáveis desconhecidas basta isolar o vetor X na equação acima,
portanto:

X = A −1B

No Matlab pode-se resolver de duas maneiras:

⇒MATLAB EDITOR ⇒RESPOSTA

A=[3 2 -1; -1 3 2; 1 -1 -1]; X=
B=[10; -5; 1]; -18.0000 ⇒ valor de x
X=inv(A)*B ⇒Utilizando o comando inv 15.0000 ⇒ valor de y
-34.0000 ⇒ valor de z

⇒MATLAB EDITOR ⇒RESPOSTA

A=[3 2 -1; -1 3 2; 1 -1 -1]; X=
B=[10; -5; 1]; -18.0000 ⇒ valor de x
X=A\B ⇒Utilizando a barra invertida 15.0000 ⇒ valor de y
-34.0000 ⇒ valor de z

6.7 Autovalores e Autovetores (Álgebra Linear)

Sendo a equação matricial:
 56
AX = λX
onde
A nxn é uma matriz quadrada de ordem n;
X é um vetor de n linhas;
λ é um escalar.

O valor de λ para que X seja vetor não nulo é chamado de autovalores da matriz A e X é chamado
de autovetores de A.

1 / 2 1 / 4
1) Sendo a matriz A =   , encontre seus autovalores e autovetores.
1 / 4 1 / 2

» A=[1/2 1/4;1/4 1/2];

» [avt,avl]=eig(A) ⇒Linha com o Comando

avt =
0.7071 0.7071 ⇒Colunas com os autovetores
-0.7071 0.7071

avl =
0.2500 0 ⇒Diagonal com os autovalores
0 0.7500

OBSERVAÇÃO:
Cada autovalor corresponde a um autovetor, ou seja:

 0,7071 
λ1 = 0,2500 ⇒ X1 =  
− 0,7071

0,7071
λ2 = 0,7500 ⇒ X2 =  
0,7071

6.8 Manipulação de vetores e matrizes

Vetor ⇒ 1 coordenada
Matriz ⇒ 2 coordenadas

REGRA GERAL

⇒ Operações com sinal (=) alteram ou incluem elementos em matrizes ou em vetores.

⇒ Operações sem sinal (=) somente retornam os elementos das matrizes ou dos vetores.

6.9 Matrizes especiais

A) Matrizes formadas por zeros
⇒ zeros(n° de linhas,n° de colunas)
 57
⇒ zeros(ordem)

Exemplos.

» zeros(2,4)
ans =
0 0 0 0
0 0 0 0

» zeros(3)
ans =
0 0 0
0 0 0
0 0 0

B) Matrizes formadas por 1 (um)

⇒ ones(n° de linhas,n° de colunas)
⇒ ones(ordem)
Exemplos.

» ones(4,2)

ans =
1 1
1 1
1 1
1 1

» ones(4)

ans =
1 1 1 1
1 1 1 1
1 1 1 1
1 1 1 1

C) Matriz Identidade ⇒ eye(ordem)

» eye(3)

ans =

1 0 0
0 1 0
0 0 1
 58
2) Sendo as matrizes:
1 2 3 10
a = [10 20 30] b = 4 5 6 c = 11
7 8 9 12
» a=[10 20 30];
» b=[1 2 3; 4 5 6; 7 8 9];
» c=[10;11;12];

» a(1) ⇒Retorna o elemento (1) do vetor (a)

a = [10 20 30] ans = 10

» a(1,3) ⇒Retorna o elemento (1,3) do da

a = [10 20 30] ans = 30 matriz (a)
» b(5) ⇒Retorna o elemento (5) do vetor (b)
1 2 3 ans = 5
b = 4 5 6
7 8 9
» b(3,2) ⇒Retorna o elemento (3,2) da matriz
1 2 3 ans = 8 (b)
b = 4 5 6
7 8 9

» c(1) ⇒Retorna o elemento (1) do vetor (c)

10 ans = 10
c = 11
12
» c(3,1) ⇒Retorna o elemento (3,1) da matriz
10 ans = 12 (c)
c = 11
12

» a(1)=0 a= ⇒Transforma o elemento (1) do vetor

a = [10 20 30] 0 20 30 (a) no algarismo 0 (zero)

» a(1,2)=0 ⇒Transforma o elemento (1,2) da

a = [0 20 30] a= matriz (a) no algarismo 0 (zero)
0 0 30
» b(7)=-11 ⇒Transforma o elemento (7) da do
1 2 3 b= vetor (b) no algarismo –11
b = 4 5 6
1 2 -11
4 5 6
7 8 9 7 8 9
» b(4,4)=-10 ⇒Inclui o algarismo 10 no elemento
b= (4,4) na matriz (b) e complementa os
1 2 -11 0 demais elementos não existentes com
 59
1 2 − 11 4 5 6 0 algarismo 0 (zero)
b = 4 5 6 
7 8 9 0
0 0 0 -10
7 8 9 
» c(2,2)=1 ⇒Inclui o algarismo 1 no elemento
10 c= (2,2) na matriz (c). e complementa os
c = 11
10 0 demais elementos não existentes com
11 1 algarismo 0 (zero)
12 12 0

3) Sendo a matriz A6x6:

1 2 3 4 5 6
7 8 9 10 11 12 
 
13 14 15 16 17 18 
A= 
19 20 21 22 23 24 
 25 26 27 28 29 30 
 
 31 32 33 34 35 36 

Execute as seguintes operações:

a) Crie uma matriz B formada pela segunda coluna de A;
b) Crie uma matriz C formada pela quinta linha de A;
c) Crie uma matriz D formada pela terceira e quarta coluna de A;
d) Crie uma matriz E formada pela primeira e a última linha de A;
e) Crie uma matriz linha F formada pelas duas primeiras linhas de A;

⇒MATLAB EDITOR
clear all
clc

A=[ 1 2 3 4 5 6
7 8 9 10 11 12
13 14 15 16 17 18
19 20 21 22 23 24
25 26 27 28 29 30
31 32 33 34 35 36];

B=A(:,2) ⇒B recebe todas as linhas e a 2ª coluna de A

C=A(5,:) ⇒C recebe a 5ª linha e todas as colunas de A
D=[A(:,3) A(:,4)]
E=[A(1,:); A(6,:)]
F=[A(1,:) A(2,:)]

⇒RESPOSTA
B=
2
8
14
 60
20
26
32

C=
25 26 27 28 29 30

D=
3 4
9 10
15 16
21 22
27 28
33 34

E=
1 2 3 4 5 6
31 32 33 34 35 36

F=
1 7 13 19 25 31 2 8 14 20 26 32

6.9 Exercícios propostos

1) Resolva as seguintes operações elemento por elemento envolvendo as matrizes abaixo.

 1 0 −3 
 1 2 − 4 e 2 2   2 ,3 
A=  , B =   , C = 3 5 0 ,
 − 8 π  11 0 cos(11,5°) )  2 2 ln(12)

− 2 2 3  − 3
   
D= 0 π π2  e E =  7,5  .
 0,5 − 0,5 − 3 / 4  2 

2
 2B  eB
a) X =   +
 2 + B 5B + 1

b) R = −
Alog 10A 2 ( )
2A − A 3 A

C+D CD 2
c) Y = −
C − D2 CD + 2 / 3
 61
2
3C − D 2  D 
d) E = +3  2C + 3 
log 2 (D 2 + 2)  

2) Execute, se possível, as seguintes operações matriciais utilizando as matrizes do exercício 1.

a) det(A) h) 2C-D o) D-1 v) E*D

b) det(B) i) det(2C-D) p) det(C-1) w) D*E
c) A+A j) A×B q) C×B x) D2
d) A+B k) A×D r) B×C y) det(C-2D)
e) A-3A l) A-1 s) 5B×(-C) z) det(A2-D)
f) A+C m) (2B)-1 t) A2
g) 2C+3D n) (3C-D)-1 u) B2

1 9 3
3) Sendo A = 4 1 − 9 , encontre o valor de E na expressão abaixo;
3 5 2
3 det(A) 25 det(2A −1 )
E= +
det(2A) det(A − A −1 )

0 log 10 log 100 log 1000 

 
 1 log 20 log 200 log 2000
4) Sendo D = , calcule:
2 0 2 −2 
 
3 e e2 ln(100) 

det(D 2 ) det(2D T − D −1 )
R = +
det(D −1 ) det(D / 2)
onde T simboliza a operação transposta.

1 3 5
1 5
5) Sendo as matrizes A = 3 e B = −1 3 − 6 , encontre o valor de C.
5 ln(12)
0,5 3,5 2
2 3
 det(A 2 )   det(3A) 
C=  − 
 det(B) − det(A)   det(B −1 ) 

− 11 − 9
6) Sendo X = , encontre o valor da expressão abaixo;
1 2
 62
− 3 det(X)
det(3X) + det(X 2 )
W=
det(X −1 ) − det(1,5X ) − det(X 2 )

− 3 log(73) 0
7) Sendo D = 0 e2 0,09 , encontre o valor da expressão abaixo:
5 −35 1/ 3

det(2D − D −1 ) det(5D 2 )
E= −
[
det (D ) T −1
] [ ]
det 3D −1

onde T simboliza a operação transposta.

8) Sendo a matriz abaixo, execute as operações solicitadas utilizando somente comandos do

MATLAB:

1 2 3 4
 
5 6 7 8
A= 
 9 10 11 12
 
13 14 15 16

a)Encontre a matriz B formada pelas colunas 2 e 3 da matriz A;

b)Encontre a matriz C formada pelas três primeiras linhas de A;
c)Encontre a matriz coluna D formada pelas duas primeiras linhas e colunas de A;
d)Encontre a transposta de C.

12 11 10
9) Sendo A= 1 2 3 , execute as operações abaixo:
−3 −5 1

a) Transforme o elemento A(1,2) no algarismo –5;

b) Transforme o elemento (3) no algarismo 0 (zero);
c) Inclua o algarismo -10 nos elementos A(4,2) e A(2,4);
d) Encontre a matriz B formada pela transposta de A;
e) Encontre a matriz C formada pelas duas últimas colunas de A;

10) Execute as operações abaixo utilizando os comandos do MATLAB.

 63
1 2 
 
a) Construa a matriz A = 5 − 5 ;
3 8 
b) Transforme os elementos A(1) e A(5) no algarismo 10;
c) Transforme os elementos A(2,2) e A(3,4) no algarismo 20;
d) Encontre a matriz linha formada por todas as colunas de A;
e) Encontre a matriz C formada pelas linhas ímpares de B;
f) Encontre o determinante da Matriz C se possível.

11) Execute as operações abaixo utilizando os comandos do MATLAB.

1 3 3 2
4 1 2 6
a) Construa a matriz A =
4 2 1 6
2 5 5 1

b) Transforme os elementos A(2) e A(7) no algarismo 0 (zero);

c) Transforme os elementos A(2,2), A(2,5)e A(5,3) no algarismo 8;
d) Encontre a transposta de A.
e) Encontre a matriz B formada pelas linhas ímpares e pelas colunas ímpares de A;
f) Encontre o determinante de B;
g) Encontre a matriz C formada pelas linhas pares e pelas colunas pares de A;
h) Encontre a transposta da Matriz de D que é formada pela segunda
e terceira coluna de A.

 5 4 0
− 3 2   
Sendo as matrizes A =   e B = 2 1 − 1 , resolva as expressões matriciais
 0 2,5  
 2 3 3 
abaixo:
2
 det(A) − det(B) 
12) X =   , onde I é a matriz identidade de ordem 3.
 det(I + B) 

det(A −1 ) − det(IA)
13) Y = , onde I é a matriz identidade de ordem 2 e T simboliza a operação
det(B 2 ) − det(BT )
transposta. Encontre a resposta no formato 5 dígitos com expoente, ver aula 2.

14) Resolva os sistemas de equações abaixo:

a)2x +3y = 1 b)3x +3y = 0 c) x +y +z +w = 6

-x +y = 1/2 2x -y +5z = 2,5 2x –y –z +2w = 18
y -z = -1 3x –5y +2z -w = -8
x +2y +1,5z -w = -9
 64

15) Encontre os autovalores e autovetores das matrizes abaixo:

2 0 − 1  3 5 3
− 5 2     
a)   b) 0 1 / 2 0 c) 0 4 6
 2 − 2    
1 0 4  0 0 1

6.10 Respostas dos exercícios propostos

1a) X =
1.0e+003 *
0.0160 0.0451 0.0012
1.0720 0.0010 0.0009

1b) R =
0 -1.1546
0.3565 -0.7065

1c) Y =
3.3333 -0.5000 -3.2400
1.0000 -4.6856 -0.1013
1.1286 2.3571 2.0702

1d) E =
0.1560 -0.7843 -4.2032
37.5405 1.8254 -12.4684
5.0870 5.0870 5.2838

2a) ans = 19.1416

2b) impossível (matriz B não é quadrada)

2c) ans = 2.0000 4.0000

-16.0000 6.2832

2d) impossível (A e B não possuem a mesma dimensão)

2e) ans = -2.0000 -4.0000

16.0000 -6.2832

2f) impossível (A e C não possuem a mesma dimensão)

2g) ans = -4.0000 6.0000 3.0000

25.0270 19.4248 29.6088
5.5000 2.5000 2.7198

2h) ans = 4.0000 -2.0000 -9.0000

25.0270 6.8584 -9.8696
3.5000 4.5000 5.7198
2i) ans = -107.5989
 65
2j) ans = 18.0000 7.3891 3.3741
66.5575 -59.1124 -8.2352

2k) impossível (n° de col. de A é diferente do n° de lin. de D)

2l) ans = 0.1641 -0.1045

0.4179 0.0522

2m) impossível (matriz B não é quadrada)

2n) ans = -0.2631 0.1004 -0.2640

0.5904 -0.1744 0.6537
-0.2913 0.0709 -0.2190

2o) impossível (matriz D é singular ⇒ det(D)=0

2p) ans = -0.0306

2q) impossível (n° de col. de C é diferente do n° de lin. de B)

2r) ans = 91.2914 39.7737 15.5142

12.9598 1.9598 -30.5650

2s) ans = -456.4570 -198.8685 -77.5709

-64.7992 -9.7992 152.8249

2t) ans = -15.0000 8.2832

-33.1327 -6.1304

2u) impossível (matriz B não é quadrada)

2v) ans = 27.0000

43.3012
-6.7500

2w) impossível (n° de col. de D é diferente do n° de lin. de E)

2x) ans = -5.2554 7.7158 15.6128

1.9324 13.5707 45.8059
1.3138 -1.9290 -3.9032

2y) ans = 199.5254

2z) ans = 4.5000 -3.5000 -8.2500

37.5405 -35.7350 -86.8885
4.7575 -7.0407 -21.3906

3) E = 0.3794
4) R = 7.3040e+003
5) C = -2.9557e+006
6) W = -8.5889e-004
 66
7) E = 535.4712
8) Y=0.1255
9) Y=3.6350e-001

10a) X = 10b) X = 10c) X =

-0.1000 1.2500 1.0000
0.4000 -1.2500 -0.0000
-0.2500 -2.0000
7.0000

11a) 11b)
avt = avt =

0.8944 0.4472 -0.7071 0.7071 0

-0.4472 0.8944 0 0 1.0000
0.7071 -0.7071 0
avl =
-6 0
0 -1 avl =

3.0000 0 0
0 3.0000 0
0 0 0.5000

11c)
avt =

1.0000 0.9806 0.8427

0 0.1961 -0.4815
0 0 0.2408

avl =
3 0 0
0 4 0
0 0 1
 67
7. ANÁLISE DE DADOS E FUNÇÕES ESTATÍSTICAS

⇒Realiza a soma de todos os elementos da variável

sum(x) (x).
⇒Realiza o produto de todos os elementos da
prod(x) variável (x).
⇒Encontra o máximo valor da variável (x).
max(x)

⇒Encontra o mínimo valor da variável (x).

min(x)

⇒Encontra a média dos valores da variável (x).

mean(x)

⇒ Calcula a mediana da variável (x)

median(x)

⇒ Calcula o desvio-padrão da variável (x)

std(x)

⇒ Calcula a variância, ou seja, o quadrado do

var(x) desvio-padrão da variável (x)
⇒Coloca as colunas da variável (x) na ordem
sort(x) crescente.
⇒Coloca as linhas da variável (x) na ordem
sortrows(x) crescente.
⇒Coloca as colunas da variável (x) na ordem
fliplr(sort(x)) decrescente.
⇒Coloca as linhas da variável (x) na ordem
fliplr(sortrow(x)) decrescente.

1 n
Média: x = ∑ xi
n i =1

Mediana:
valor “do meio” da distribuição, Se x1, x2, x3, ..., xn são ordenados de forma crescente,
então a mediana corresponde ao valor da observação:
n +1
se n é impar
12 n
Desvio Padrão: s = ∑ ( xi − x )
2

− 1 i =1 entre e n + 2 se n é par
nmédia n
2 2
1 n
Variância: s 2 = ∑ ( xi − x )
2

n − 1 i =1

1) Encontre a soma da PA = 1,3,5...99.

» PA = (1:2:99);
» sum(PA)

ans =
2500
 68

2) Sendo a série definida por A(m)=2m-3 onde m=1, 2, 3...10. Encontre o produto dos 10 primeiros
termos dessa série.
Indexando vetores
» m = (1:10); » m = (1:10);
» A = 2*m-3 » A(m) = 2*m-3
A=
-1 1 3 5 7 9 11 13 15 17

»prod(A)

ans =
-34459425

3) Sendo a série definida por A(m)=m2-4 onde m=0, 1, 2...10. Encontre a média, a mediana,
variância e desvio padrão dos 11 primeiros termos dessa série.

» m = (0:10); Indexando vetores

» A = m.^2-4 » m = (0:10);
» A(m+1) = m.^2-4
A=
-4 -3 0 5 12 21 32 45 60 77 96

» mean(A)

ans =
31

» median(A)

ans =

21

» var(A)

ans =

1.1858e+003

» std(A)

ans =

34.4354

4) Sendo o vetor A=[1 –3 5 –11.5 21 33.5 8 -13], pede-se:

a) encontre o máximo valor de A;
b) encontre o mínimo valor de A;
c) coloque o vetor A na ordem crescente;
 69
d) coloque o vetor A na ordem decrescente.

» A = [1 -3 5 -11.5 21 33.5 8 -13];

» max(A)
ans =
33.5000

» min(A)
ans =
-13

» sort(A)

ans =
Columns 1 through 7
-13.0000 -11.5000 -3.0000 1.0000 5.0000 8.0000 21.0000

Column 8
33.5000

» fliplr(sort(A))

ans =
Columns 1 through 7
33.5000 21.0000 8.0000 5.0000 1.0000 -3.0000 -11.5000

Column 8
-13.0000

5) Encontre o valor de F na equação.

10 2
i 12 22 32 102
F=∑ = + + + ... +
i =1 2i 2 ⋅ 1 2 ⋅ 2 2 ⋅ 3 2 ⋅ 10

» i = (1:10);
» F = (i.^2)./(2*i)

F=
Columns 1 through 7
0.5000 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000

Columns 8 through 10
4.0000 4.5000 5.0000

» F = sum(F)

F=
 70
27.5000

Ex.6)Encontre o valor de n=7! Considerando que n é um vetor.

» n = (1:7);
» prod(n)

ans =
5040

Caso contrário
» factorial(7)

ans =

5040

7.2 Indexando Vetores

Em algumas situações é necessário indexar os vetores para armazenar valores em posições pré-
estabelecidas, isto é, quando o vetor é dependente de uma variável ou posição.

Ex.7) Encontre o vetor M dado pelas equações abaixo:

2n / n 2 se n≤4
M=
n 2 / 2n se n>4
onde n=1,2,...8

⇒MATLAB EDITOR
clear all
clc

n = (1:4); !! observe que n deve ser inteiro (+)

M(n) = (2.^n)./(n.^2);

n = (5:8);
M(n) = (n.^2)./(2.^n);

⇒RESPOSTA
M=
Columns 1 through 7
2.0000 1.0000 0.8889 1.0000 0.7813 0.5625 0.3828

Column 8
0.2500

» plot(1:8,M)
» xlabel('n')
» ylabel('M(n)')
 71

1.8

1.6

1.4

1.2

M(n)
1

0.8

0.6

0.4

0.2
1 2 3 4 5 6 7 8
n

Ex.8)Resolva os exercícios 2 e 3, utilizando a técnica de indexação de vetores.

Exercício 2
Sendo a série definida por A(m)=2m-3 onde m=1, 2, 3...10. Encontre o produto dos 10 primeiros
termos dessa série.

» m = (1:10);
» A(m) = 2*m-3

A=
-1 1 3 5 7 9 11 13 15 17

»prod(A)

ans =
-34459425

Exercício 3
Sendo a série definida por A(m)=m2-4 onde m=0, 1, 2...10. Encontre a média dos 6 primeiros
termos dessa série.

» clc
» clear all
» m = (0:10);
» A(m+1) = m.^2-4

A=
-4 -3 0 5 12 21 32 45 60 77 96

» mean(A(1:6))

ans =
5.1667
 72
7.3 Administração de arquivo de dados

Exportação de arquivos

save diretório:\fname.ext file –ascii

⇒Exporta os dados da matriz file no arquivo fname.ext com o formato ascii com 8 dígitos,
armazenando-os na raiz do diretório indicado.

save diretório:\fname.ext file –ascii –double

⇒Exporta os dados da matriz file no arquivo fname.ext com o formato ascii com 16 dígitos,
armazenando-os na raiz do diretório indicado.

save diretório:\caminho\... \fname.ext file –ascii

⇒Exporta os dados da matriz file no arquivo fname.ext com o formato ascii com 8 dígitos,
armazenando-os na pasta cujo caminho é dado por diretório:\caminho\... \

save diretório:\caminho\... \fname.ext file –ascii –double

⇒Exporta os dados da matriz file no arquivo fname.ext com o formato ascii com 16 dígitos,
armazenando-os na pasta cujo caminho é dado por diretório:\caminho\... \

Observação:
Arquivos no formato ascii podem ser editados em qualquer editor comum de textos ou dados (*.txt
ou *.dat).

Observação: Quando se fala em matriz subentende-se também os vetores

Importação de Arquivos (*.txt ou *.dat)

load diretório:fname.ext

⇒Importa a(s) matriz(es) contidas no arquivo fname.txt da raiz do diretório indicado.

Para ver o que o arquivo fname.txt contém basta digitar no Matlab (após ele ter sido importado),
fname

Se você já está no diretório em que está o arquivo basta digitar

load fname.ext

load diretório:\caminho\... \fname.ext

⇒Importa a(s) matriz(es) contidas no arquivo da pasta cujo caminho é dado por
diretório:\caminho\... \
 73

7.4 Comandos de visualização

disp('texto') ⇒Imprime o texto digitado entre aspas.

disp(X) ⇒Imprime o conteúdo da variável (x)
⇒Imprime a variável (X) inteira ao lado do texto
disp(['texto',int2str(X)]) digitado, transformando-a em uma string de
caracteres.
⇒Imprime a variável (X) real ou complexa ao lado
do texto digitado, transformando-a em uma string
disp(['texto',num2str(X,N)])
de caracteres com N dígitos (N máximo=16).

⇒Quando N não é informado utiliza-se o formato

padrão do MATLAB de 4 dígitos após o ponto
disp(['texto',num2str(X)])
decimal.

1) Um aluno de engenharia mediu aleatoriamente o diâmetro externo e o interno de alguns anéis de

aço com auxílio de um paquímetro. Os valores medidos foram gravados no notepad (bloco de
notas) em um arquivo denominado anel.txt. Nos anéis estão gravados seus respectivos números de
identificação. Importe este arquivo de dados da pasta c:\temp e com auxílio do MATLAB.

⇒JANELA DE COMANDOS importa arquivo anel.txt da pasta

» load c:\temp\anel.txt c:\temp
» anel

anel = Importando Arquivos

1.0000 100.0000 50.1000

8.0000 99.9600 49.9800 1° - Importar o arquivo na janela
3.0000 99.9800 50.0200 de comandos para verificar a
10.0000 100.0200 50.0200 ordem em que se apresenta os
2.0000 100.0200 50.0200 dados.
5.0000 100.0000 50.0400
6.0000 100.0200 50.0000
7.0000 100.0400 50.0400
4.0000 99.9200 50.0800 2° - Importar os dados no MATLAB
9.0000 99.9800 50.0000 EDITOR e realizar as mudanças
necessárias.

2) Dado o exemplo anterior, faça uma tabela com os dados importados e devidamente organizados
na ordem crescente.

⇒MATLAB EDITOR
clear all
clc

load c:\temp\anel.txt
 74
%Coloca as linhas da matriz anel na ordem crescente
anel = sortrows(anel);

%Mostra a tabela
disp(' ')
disp(' n° dext[mm] dint[mm] ')
disp(' ===============================')
disp(anel)

⇒RESPOSTA
n° dext[mm] dint[mm]
===============================
1.0000 100.0000 50.1000
2.0000 100.0200 50.0200
3.0000 99.9800 50.0200
4.0000 99.9200 50.0800
5.0000 100.0000 50.0400
6.0000 100.0200 50.0000
7.0000 100.0400 50.0400
8.0000 99.9600 49.9800
9.0000 99.9800 50.0000
10.0000 100.0200 50.0200

3) Dado o exemplo anterior encontre o valor máximo, o valor mínimo e a média dos diâmetros
externo e interno.

⇒MATLAB EDITOR
clear all
clc

load c:\temp\anel.txt

%Desacopla as colunas 2 e 3 da matriz anel

dext = anel(:,2);
dint = anel(:,3);

%Encontra o valor máx., min. e a média do diâmetro externo.

dextmax = max(dext);
dextmin = min(dext);
dextmed = mean(dext);

%Encontra o valor máx., min. e a média do diâmetro interno.

dintmax = max(dint);
dintmin = min(dint);
dintmed = mean(dint);

%Imprime a saída dos dados

disp([' Diâmetro ext. máx.= ' ,num2str(dextmax,5)])
disp([' Diâmetro ext. mín.= ' ,num2str(dextmin,5)])
disp([' Diâmetro ext. médio= ',num2str(dextmed,5)])
disp([' Diâmetro int. máx.= ' ,num2str(dintmax,5)])
 75
disp([' Diâmetro int. mín.= ' ,num2str(dintmin,5)])
disp([' Diâmetro int. médio= ',num2str(dintmed,5)])

⇒RESPOSTA
Diâmetro ext. máx.= 100.04
Diâmetro ext. mín.= 99.92
Diâmetro ext. médio= 99.994
Diâmetro int. máx.= 50.1
Diâmetro int. mín.= 49.98
Diâmetro int. médio= 50.03

4) No exemplo 1, sabe-se que a medida padrão do diâmetro externo é 100,00mm e do diâmetro

interno é 50,00mm. Grave um arquivo no formato texto (txt) contendo o n° do anel, o diâmetro
externo e o diâmetro interno com suas respectivas diferenças para os valores padrão apresentados.

⇒MATLAB EDITOR
clear all
clc

load c:\temp\anel.txt

anel=sortrows(anel);

%Desacopla as colunas 1, 2 e 3 da matriz anel

n_anel = anel(:,1);
dext = anel(:,2);
dint = anel(:,3);

%Calcula a diferença entre o valor médio e o valor padrão

dif_ext = 100-dext;
dif_int = 50-dint;

%Monta a matriz a ser gravada

resp=[n_anel dext dif_ext dint dif_ext]

%Salva a matriz resp no arquivo anel2.txt com 8 dígitos

save c:\temp\anel2.txt resp -ascii
disp(' arquivo salvo ')

⇒RESPOSTA
resp =
1.0000 100.0000 0 50.1000 0
2.0000 100.0200 -0.0200 50.0200 -0.0200
3.0000 99.9800 0.0200 50.0200 0.0200
4.0000 99.9200 0.0800 50.0800 0.0800
5.0000 100.0000 0 50.0400 0
6.0000 100.0200 -0.0200 50.0000 -0.0200
7.0000 100.0400 -0.0400 50.0400 -0.0400
8.0000 99.9600 0.0400 49.9800 0.0400
9.0000 99.9800 0.0200 50.0000 0.0200
10.0000 100.0200 -0.0200 50.0200 -0.0200
 76
Arquivo anel2.txt (salvo)
1.0000000e+000 1.0000000e+002 0.0000000e+000 5.0100000e+001 0.0000000e+000
2.0000000e+000 1.0002000e+002 -2.0000000e-002 5.0020000e+001 -2.0000000e-002
3.0000000e+000 9.9980000e+001 2.0000000e-002 5.0020000e+001 2.0000000e-002
4.0000000e+000 9.9920000e+001 8.0000000e-002 5.0080000e+001 8.0000000e-002
5.0000000e+000 1.0000000e+002 0.0000000e+000 5.0040000e+001 0.0000000e+000
6.0000000e+000 1.0002000e+002 -2.0000000e-002 5.0000000e+001 -2.0000000e-002
7.0000000e+000 1.0004000e+002 -4.0000000e-002 5.0040000e+001 -4.0000000e-002
8.0000000e+000 9.9960000e+001 4.0000000e-002 4.9980000e+001 4.0000000e-002
9.0000000e+000 9.9980000e+001 2.0000000e-002 5.0000000e+001 2.0000000e-002
1.0000000e+001 1.0002000e+002 -2.0000000e-002 5.0020000e+001 -2.0000000e-002

5) Um estudante americano utilizando um sistema de medição mediu a temperatura do Deserto do

Saara num dia de verão e obteve os seguintes dados.
1 33.8 13 132.7
2 34.34 14 132.7
3 36.14 15 131.8
4 41.0 16 130.48
5 46.4 17 126.76
Salvar na pasta c:\temp
6 56.3 18 102.5
como saara1.txt
7 67.1 19 75.38
8 79.34 20 74.3
9 90.32 21 57.56
10 107.6 22 35.6
11 131.54 23 34.167
12 133.6 24 32.0

Onde as colunas 1 e 3 representam as horas em que foram feitas as leituras das temperaturas e as
colunas 2 e 4 representam as respectivas temperaturas em °F. Os resultados foram enviados para
seu amigo no Brasil na forma de um arquivo texto, saara1.txt. Pede-se:

a)Salve um novo arquivo denominado saara2.txt com as horas e com suas respectivas temperaturas
medidas em °C;
TC TF − 32
=
5 9

⇒MATLAB EDITOR
clear all
clc

%Importa o arquivo saara1.txt

load c:\temp\saara1.txt

%Transforma a 1ª e a 3ª coluna em um único vetor linha

horas=[saara1(:,1)' saara1(:,3)'];

%Transforma a 2ª e a 4ª coluna em um único vetor linha

TF=[saara1(:,2)' saara1(:,4)'];

%Transforma o vetor TF para TC

TC=(TF-32)/9*5;

%Transforma os vetores linha horas e TC em colunas

%armazenando-os na matriz M
 77
M=[horas' TC'];

%Salva a matriz M no arquivo saara2.txt com 8 dígitos

save c:/temp/saara2.txt M -ascii
disp(' arquivo salvo ')

⇒RESPOSTAS
Arquivo saara2.txt
1.0000000e+000 1.0000000e+000
2.0000000e+000 1.3000000e+000
3.0000000e+000 2.3000000e+000
4.0000000e+000 5.0000000e+000
5.0000000e+000 8.0000000e+000
6.0000000e+000 1.3500000e+001
7.0000000e+000 1.9500000e+001
8.0000000e+000 2.6300000e+001
9.0000000e+000 3.2400000e+001
1.0000000e+001 4.2000000e+001
1.1000000e+001 5.5300000e+001
1.2000000e+001 5.6444444e+001
1.3000000e+001 5.5944444e+001
1.4000000e+001 5.5944444e+001
1.5000000e+001 5.5444444e+001
1.6000000e+001 5.4711111e+001
1.7000000e+001 5.2644444e+001
1.8000000e+001 3.9166667e+001
1.9000000e+001 2.4100000e+001
2.0000000e+001 2.3500000e+001
2.1000000e+001 1.4200000e+001
2.2000000e+001 2.0000000e+000
2.3000000e+001 1.2000000e+000
2.4000000e+001 0.0000000e+000

7.5 Exercícios propostos

1) Encontre o valor de E na expressão abaixo.

2n 3 + n 2
10
E=∑
n =1 2n − 1

2) Encontre a soma, o produto e a média dos 10 primeiros termos da série formada por:
i3
S=
3i − 2i
3) Encontre a série formada pelos 12 primeiros termos de X fornecida pelas equações abaixo:

2m −1 ⇒para m impar
X= onde m=1, 2, 3,...12
m −1
−2 ⇒para m par

Obs.: Utilize a técnica de indexação de vetores.

4) No exercício 3 encontre:
a)A soma dos elementos de X;
 78
b)A média dos elementos de X;
c)Encontre o vetor Y com os elementos de X na ordem decrescente.

5) Um engenheiro fornece as medidas de um determinado furo de um componente mecânico em cm

de acordo com a tabela abaixo. Escreva este arquivo no notepad e faça sua manipulação através do
MATLAB.

1 10.120 6 10.000
2 10.120 7 10.080
3 10.240 8 10.000
4 10.120 9 9.980
5 10.100 10 10.060

Sabendo que a medida padrão do furo é 100.00 mm. Faça um programa que:

a) A partir do arquivo dado transforme as medidas de cm para mm;

b)
c) Salve um arquivo no formato xls (excel) e outro no formato txt (texto), com três colunas: a
primeira com o n° do componente, a segunda com o valor padrão em mm e a terceira com a
diferença entre o valor medido e o valor padrão.

6) Um aluno de engenharia realizou 12 medições do diâmetro e do comprimento de um cilindro de

aço. Os cilindros são numerados de 1 a 12. Infelizmente ele não tomou cuidado ao realizar as
medições e pegou os cilindros aleatoriamente tendo os resultados abaixo:

9 30.25 120.10
2 30.20 120.10
7 30.05 119.90
8 30.05 119.00
10 29.95 120.00
3 29.90 120.20
11 30.10 120.10
4 29.80 120.00
5 29.85 119.90
1 30.00 119.95
12 30.25 120.00
6 29.80 120.00

onde: a coluna 1 representa o número do cilindro, a coluna 2 representa o diâmetro e a coluna 3

representa o comprimento cilindro. (todas a dimensões mm) Sabe-se que a medida padrão do
diâmetro é 30,00mm e do comprimento 120,50mm. Pede-se:

a) Escreva este arquivo (acima) no notepad (salve com o nome cilindro1.txt) e faça sua
manipulação através do MATLAB.

b) Elabore um programa que deve salvar um arquivo com 5 colunas denominado cilindro2.txt
contendo:
coluna1 ⇒número do cilindro (na ordem crescente);
coluna2 ⇒diâmetro medido;
coluna3 ⇒diferença entre o diâmetro medido e a média dos valores medidos;
coluna4 ⇒diferença entre o comprimento medido e a média dos valores medidos.
 79

7.6 Respostas dos Exercícios Propostos

1) Resp: E=446.0666
2) Soma=6,2259; média=0,6226; produto=3,8190.106
3) X =[1 -2 4 -8 16 -32 64 -128 256 -512 1024 -2048]
4a) Soma = -1365
4b)media = -113.7500
4c) Y=[1024 256 64 16 4 1 -2 -8 -32 -128 -512 -2048]
 80

8. GRÁFICOS BIDIMENSIONAIS

8.1 Comandos básicos

⇒Plota o gráfico do vetor abscissa (x) com seus respectivos

plot(x,y) elementos no vetor ordenada (y). Os vetores (x) e (y)
devem ter a mesma dimensão.
⇒Insere o texto digitado entre apostos ('texto') na parte
title('texto')
superior de gráfico.
⇒Insere o texto digitado entre apostos ('texto') abaixo da
xlabel('texto')
abscissa do gráfico.
⇒Insere texto digitado entre apostos ('texto') verticalmente
ylabel('texto')
na ordenada do gráfico.
Grid ⇒Insere linhas de grade no gráfico.

8.2 Comandos para formatação

⇒Ajusta as escalas dos eixos (x) e (y) no

axis([xmin xmax ymin ymax])
gráfico.
⇒Plota vários tipos de linha, marcadores e
plot(x,y,'s') cores. Onde (s) é um ou dois ou três caracteres
da tabela abaixo.
⇒Ajusta espessura (v) da linha de plotagem
u=plot(x,y)
do gráfico definido pela variável (u). Onde (v)
set(u,'linewidth',v)
= 1 – 1,5 – 2,0 – 2,5 - ... ptos.

Cores Marcas Linhas

y amarelo . ponto - linha continua
m magenta o círculo : linha pontilhada
c ciano x x-mark -. Traços e pontos
r vermelho + plus -- linha tracejada
g verde * estrela
b azul s quadrado
w branco d diamante
k preto v triângulo(p/baixo)
^ triângulo(p/cima) <
triângulo(p/esquerda) >
triângulo(p/direita) p
pentagrama h
hexagrama

A construção de gráficos bidimensionais no MATLAB deve obedecer a seqüência estruturada

abaixo:
 81

CRIAR
Vetor abscissa

CRIAR
Vetor ordenada

plot(abscissa,ordenada)
Somente se há necessidade de
ajustar os eixos x e y introduzir o
comando:
axis([xmin xmax ymin ymax])
inclusão de título
title('texto')

identificar eixo abscissa

xlabel('texto')

identificar eixo ordenada

ylabel('texto')

A) CONSTRUÇÃO DE GRÁFICOS
1) Faça o gráfico da função y = x 3 − 36 x quando − 8 ≤ x ≤ 8 .

⇒MATLAB EDITOR Não esqueça de utilizar as operações pontuadas

clear all elemento por elemento quando for necessário.
clc (Ver aula 3 - pág. 30)

x=linspace(-8,8,100); % Cria o vetor x com 100 elementos.

y=x.^3-36*x;

plot(x,y)
title('Função y=x^3-36x')
xlabel('eixo x')
ylabel('eixo y')
grid
 82

B) AJUSTE DOS EIXOS COORDENADOS

2) Faça o gráfico da função y=f(θ) sendo y = sin (θ) e 0 ≤ θ ≤ 2π rad.

⇒MATLAB EDITOR
clear
clc

theta=linspace(0,2*pi);
y=sin(theta);

plot(theta,y)
title('y=sin(theta)')
xlabel('theta [rad]')
ylabel('y')
grid
 83
Caracteres usados no Matlab em legendas e títulos de figuras, isto é, usados
como um texto.
Caracteres Resultado Caracteres Resultado Caracteres Resultado
\alpha α \upsilon υ \sim ~
\beta β \phi Φ \leq ≤
\gamma γ \chi X \infty ∞
\delta δ \psi ψ \clubsuit ♣
\epsilon ε \omega ω \diamondsuit ♦
\zeta ζ \Gamma г \heartsuit ♥
\eta η \Delta ∆ \spadesuit ♠
\theta θ \Theta Θ \leftrightarrow ↔
\vartheta ϑ \Lambda Λ \leftarrw ←
\iota ι \Xi Ξ \uparrow ↑
\Kappa κ \Pi Π \rightarrow →
\lambda λ \Sigma ∑ \downarrow ↓
\mu µ \Upsilon Y \circ o
\nu ν \Phi Φ \pm ±
\xi ξ \Psi Ψ \geq ≥
\pi π \Omega Ω \propto ∝
\rho ρ \forall ∀ \partial ∂
\sigma σ \exists ∃ \bullet •
\varsigma ζ \ni ∋ \div ÷
\tau τ \cong ≅ \neq ≠
\equiv ≡ \approx ≈ \aleph ℵ
\lm ℑ \Re ℜ \wp ℘
\otimes ⊗ \oplus ⊕ \oslash ø
\cap ∩ \cup ∪ \supseteq ⊇
\supset ⊃ \subseteq ⊆ \subset ⊂
\int ∫ \in є \o
\lfloor  \lceil  \nabla ∇
\rfloor  \cdot . \dots ...
\perp ⊥ \neg ¬ \prime ´
\wedge ^ \times x \0 ∅
\rceil  \surd √ \mid |
\vee ٧ \varpi ϖ \copyright 
\langle ‹ \rangle ›

3) Faça o gráfico da função y=f(θ) sendo y = sin (θ) e 0 ≤ θ ≤ 2π rad; e ajuste seus eixos
coordenados x e y.

⇒MATLAB EDITOR
clear
clc

theta=linspace(0,2*pi);
y=sin(theta);

plot(theta,y)
axis([0 2*pi -1.5 1.5]) %altera as escalas dos eixos x e y
 84

title('y=sin( theta )')

xlabel('theta [rad]')
ylabel('y')
grid

C) MUDANÇA DE CORES, LINHAS E INCLUSÃO DE MARCADORES.

4) Faça o gráfico da função y = x 3 − 36 x quando − 8 ≤ x ≤ 8 na cor vermelha.

⇒MATLAB EDITOR
clear all
clc

x=linspace(-8,8);
y=x.^3-36*x;

plot(x,y,'r') %altera a cor da linha para vermelho

title('Função y=x^3-36x')
xlabel('eixo x')
ylabel('eixo y')
grid
 85
5) A matriz R apresenta os dados obtidos em um experimento de MRUV no Laboratório de Física:

 0 2 4 6 8 10 12 14 16 18 
R= 
10,0 − 5,5 − 13,3 − 12,0 − 2,5 15,0 41,0 75,5 118,5 170,0

A primeira linha representa o tempo t em [s] e a segunda linha representa o deslocamento X em

[cm] de um móvel. Faça o gráfico X=f(t) com linhas pontilhadas, na cor magenta e ressaltando os
pontos medidos com círculos.

⇒MATLAB EDITOR
clear all
clc

t = (0:2:18);
X = [10 –5.5 –13 –12 –2.5 15 41 75.5 118.5 170];

plot(t,X,'mo:')
title('MRUV')
xlabel('t [s]')
ylabel('X [cm]')

D) MUDANÇA DE ESPESSURA DE LINHA

6) Faça o gráfico da função y = x 3 − 36 x quando − 8 ≤ x ≤ 8 na cor vermelha com espessura

de linha 2,5 pontos.

⇒MATLAB EDITOR
clear all
clc

x=linspace(-8,8);
y=x.^3-36*x;

u=plot(x,y,'r') %altera a cor da linha para vermelho

set(u,'linewidth',2.5) %ajusta a espessura da linha para 2.5 ptos
 86
title('Função y=x^3-36x')
xlabel('eixo x')
ylabel('eixo y')
grid

8.3 Gráficos simultâneos

⇒Plota os gráficos (x1,y1), (x2,y2), (x3,y3)..., no

plot(x1,y1,x2,y2,x3,y3,...)
mesmo par de eixos coordenados.
⇒Cria uma nova janela gráfica (i). Onde i = 1, 2,
figure(i)
3....
⇒Subdivide uma janela gráfica em m-linhas e n-
subplot(m,n,p) colunas nas quais pode-se traçar gráficos nas
posições p.

Comandos auxiliares:

legend('gráf.1','gráf.2',...) ⇒Insere legenda no gráfico.

gtext('texto') ⇒Insere texto com o mouse sob o gráfico.

OBSERVAÇÃO:
Todos os comandos para formatação de linhas (ver pág.60) são válidos para gráficos simultâneos.

E) GRÁFICOS NO MESMO PAR DE EIXOS COORDENADOS

7) Faça os gráficos das funções y=sin(θ), w=cos(θ) e z=tg(θ) em um mesmo par de eixos
coordenados, onde 0 ≤ θ ≤ 360° .

⇒MATLAB EDITOR
clear all
clc
 87

theta=linspace(0,360);
y=sin(theta*pi/180);
w=cos(theta*pi/180);
z=tan(theta*pi/180);

plot(theta,y,theta,w,theta,z) %plot(x1,y1,x2,y2,x3,y3)
axis([0 360 -3 3])
legend('y=sin(theta)' ,'w=cos(theta)' , 'z=tg(theta)')

title('Funções Trigonométricas: sin - cos - tg')

xlabel('theta [graus]')
ylabel('y , w , z')

8) Repita o exercício anterior trocando o comando legend pelo comando gtext.

⇒MATLAB EDITOR
clear all
clc
teta=linspace(0,360);
y=sin(theta*pi/180);
w=cos(theta*pi/180);
z=tan(theta*pi/180);
plot(theta,y,theta,w,theta,z)
axis([0 360 -3 3])
title('Funções Trigonométricas: sin - cos - tg')
xlabel('theta [graus]')
 88
ylabel('y , w , z')
gtext('sin(theta)')
gtext('cos(theta)')
gtext('tg(theta)')

F) GRÁFICOS NA MESMA JANELA GRÁFICA

9) Faça os gráficos das funções y=sin(θ), w=cos(θ) e z=tg(θ), onde 0 ≤ θ ≤ 360° na mesma
janela gráfica e em pares de eixo coordenados diferentes.

⇒MATLAB EDITOR
clear all
clc

teta=linspace(0,360);
y=sin(theta*pi/180);
w=cos(theta*pi/180);
z=tan(theta*pi/180);
subplot(2,2,1)
plot(theta,y)
axis([0 360 –1.5 1.5])
title('y=sin(teta)')

subplot(2,2,2)
plot(theta,w)
axis([0 360 –1.5 1.5])
title('w=cos(teta)')

subplot(2,2,3)
plot(theta,z)
axis([0 360 –3 3])
title('z=tg(teta)')

subplot(2,2,4)
plot(theta,y,theta,w,theta,z)
axis([0 360 –3 3])
title('sin – cos - tg')
 89
subplot(n,m,p)
n-Colunas

Posição Posição
p=1 p=2

m-Linhas

Posição Posição
p=3 p=4

G) GRÁFICOS EM JANELAS GRÁFICAS DISTINTAS

10) Faça os gráficos das funções y=sin(θ), w=cos(θ) e z=tg(θ), onde 0 ≤ θ ≤ 360° em três
janelas gráficas distintas.

⇒MATLAB EDITOR
clear
clc
theta=linspace(0,360);
y=sin(theta*pi/180);
w=cos(theta*pi/180);
z=tan(theta*pi/180);
figure(1)
plot(theta,y)
axis([0 360 –1.5 1.5])
title('y=sin(theta)')
figure(2)
plot(theta,w)
axis([0 360 –1.5 1.5])
title('w=cos(theta)')
figure(3)
plot(theta,z)
axis([0 360 -10 10])
title('z=tan(theta)')
 90

8.4 Exercícios complementares

1) Faça o gráfico da função y = f(t), abaixo, onde 0 ≤ t ≤ 10 na cor azul espessura de linha 2
pontos.
y = −4 t 3 + t 2 + 2 t

2) Faça o gráfico da função y = f(x), abaixo, onde − 10 ≤ x ≤ 10 .

1
y=
x

3) Faça o gráfico da função P = f(x) na cor verde utilizando marcador do tipo pentagrama e com
linha pontilhada. Sabe-se que:
x5
P= e x=[-20 –16 –12... 20]
x3 +1

4) Faça o gráfico da função y=f(t) na cor preta com espessura de linha 2,0 pontos.
y = ( t 2 − 9)( t 2 − 1) e − 4 ≤ t ≤ 4

5) Faça o gráfico de y=f(z), onde:

y = (z 2 − 4)e z ⇒ −5 ≤ z ≤ 2

6) Faça o gráfico da função y=f(x) na cor magenta, onde:

y = x 4 − 4x 3 + 10 e −2≤ x ≤5

7) Um automóvel em MRUV parte de uma posição Xô = 25m com aceleração de -2,25m/s2 e

velocidade inicial de Vo = 45m/s. Faça os gráficos de X = f(t) e V = f(t) em uma mesma janela
gráfica um ao lado do outro sabendo que 0 ≤ t ≤ 30s .

8) A amplitude de oscilação de um pêndulo é dado pela equação abaixo:

− πt
[
Y = e 25 cos(2πt ) + sen(2πt ) ]
onde t é o tempo de oscilação em [s] e varia de 0 a 30s
Y é a amplitude de oscilação em [m].
Faça o gráfico de Y = f(t) na cor verde.

9) Um caminhão tem sua posição em km (s) monitorada durante 15 horas conforme a equação
s = 20t 2 − t 3 . Faça o gráfico da posição s = f(t), da velocidade v = f(t) e de sua aceleração a = f(t)
em uma mesma janela gráfica e em pares de eixos coordenados diferentes um embaixo do outro.
Lembre que v = ds/dt e a = d2s/dt2.

10) Plote no mesmo par de eixos coordenados as funções u = sec(θ), v = cossec(θ) e w =

cotg(θ), sendo 0 ≤ θ ≤ 2π rad.

11) Faça o gráfico da função X = f(t), onde − 10 ≤ t ≤ 10 .

 91

X=
[
log ( t + 12) 3 ]
( t + 15)

12) Em uma experiência em laboratório de biologia verificou-se que o crescimento de uma

determinada colônia de bactérias é dado por:

y=
2
[
1 t
(e − 3 cos(10t ) + sen(10 t ) + 2) ]
onde t é o período de horas e y é nº de bactérias em milhões.

Depois de passar 3,6 horas do início experimento é aplicado um antídoto para aniquilamento total
das bactérias. Sabe-se que o decréscimo do número de bactérias é dado pela equação abaixo:

y = −(t − 3,64 )2 + 18,9952

Quando a experiência chega em sua oitava hora verifica-se que não há mais bactérias vivas, então,
pede-se:

a) faça o gráfico que rege todo o experimento desde seu início até a oitava hora;
b) o gráfico deve ser plotado na cor vermelha com linha contínua até a aplicação do antídoto;
c) após a aplicação do antídoto o gráfico deve ser plotado na cor vermelha pontilhada.

13) Em um determinado experimento de física sobre lançamento vertical de projéteis tem-se o

alcance máximo (A) de um projétil em função da sua velocidade de disparo (Vo). A matriz R
apresenta os dados coletados, onde a 1ª linha representa a velocidade de disparo em m/s e a 2ª
coluna representa o alcance (A) em metros.

 50 75 100 125 150 175 200 

R= 
 250 550 950 1550 2220 3020 3950

Em um mesmo par de eixos coordenados, pede-se:

I) Faça o gráfico de A=f(Vo) com os dados obtidos em R, utilizando linha contínua verde
ressaltando os pontos medidos com quadrados.
II) Faça o gráfico de A=f(Vo) utilizando a equação de lançamento oblíquo de projéteis
(abaixo), sabendo que o ângulo de disparo utilizado neste experimento foi de α=45° e a
aceleração da gravidade no local é de g=9,81m/s2. Utilize linha vermelha pontilhada com
triângulos.

Vo 2
A= sin 2 (α)
g

eu − e−u
14) Plote o gráfico da função v=f(u) na cor preta espessura 2,5 pontos, sendo v = e
eu + e−u
− 5 ≤ u ≤ −5
 92

15) Plote o gráfico de y=f(x) sabendo que:

a + b 
x = (a + b) cos(θ) − b cos  θ
 b 
a + b 
y = (a + b) sen(θ) − b sen  θ
 b 

0 ≤ θ ≤ 2π

a=1 e b=1/3

16) Faça o gráfico da função y=f(t), abaixo, na cor vermelha com espessura de linha igual à 2,5ptos.

y = cos(1,5πt ) + cos(1,75πt ) ⇒ 4 ≤ t ≤ 20

17) Faça o gráfico da função h=f(x), abaixo, na cor preta e com espessura de linha 2 ptos.

1 1
h= + − 6 ⇒ −0,5 ≤ x ≤ 1,5
( x − 0,3) 2 + 0,01 ( x − 0,9) 2 + 0,04

18) Plote simultaneamente, em um mesmo par de eixos coordenados, os gráficos das funções y=f(x)
e z=f(x), dadas por:
2 e x + e− x
y= e z=
e x +e − x 2

Sendo: − 2 ≤ x ≤ 2 .

O gráfico deve de conter:

a) A curva de y dever ser plotada em vermelho pontilhado e a curva de z em preto pontilhado;
b) O gráfico deve conter legenda sabendo que y representa a secante hiperbólica (sech(x)) e z
representa o cosseno hiperbólico de x (cosh(x));

19) O lançamento oblíquo de projéteis obedece as seguintes equações:

Fomulário
 93

Onde:
hm ⇒ altura máxima atingida pelo projétil em [m]; A⇒ alcance do projétil em [m]; t⇒ tempo de
viagem do projétil em[s]; tm⇒ tempo de viagem do projétil para atingir a altura máxima; t’⇒
tempo de viagem do projétil para atingir o alcance do projétil; vx⇒componente horizontal da
velocidade do projétil; vy⇒componente vertical da velocidade do projétil; x⇒deslocamento
horizontal do projétil; y=deslocamento vertical do projétil; vo=velocidade de disparo do projétil; vox
⇒componente horizontal da velocidade de disparo; voy ⇒componente horizontal da velocidade de
disparo; α⇒ângulo de disparo do projétil em [graus]; g⇒aceleração da gravidade no local em
[m/s2].

Faça um gráfico do deslocamento (x X y), do projétil, para α=45° (alcance máximo) sabendo que
a velocidade de disparo do projétil é de Vo=200 km/h e g=9,81 m/s2.

20) Utilizando os comando do MATLAB, faça o gráfico do deslocamento do projétil para um

ângulo de disparo (α) de 15°, 30°, 45°, 60° e 75°, no mesmo par de eixos coordenados, sabendo
que:

a) Velocidade de disparo do projétil ⇒ vo=200 km/h

b) A aceleração da gravidade no local é de g=9,81 m/s2;
c) As linhas dos gráficos devem ter espessura de linha 3,0 pontos.

21) Sendo a função:

x3 − x 2 − 4x + 4
y= ⇒ −3≤ x ≤ 4
x−2

Utilizando-se dos comandos do MATLAB faça o gráfico de y=f(x) na cor preta.

22) Em um experimento no laboratório de física verificou-se que a velocidade V=f(t) de um móvel

em três intervalos diferentes é regida pelas equações abaixo:
 94

 V = 2,5t ⇒ 0 ≤ t ≤ 10s

 V = −3t + 55 ⇒ 10 ≤ t ≤ 15s
 − 10t + 400
V = ⇒ 15 ≤ t ≤ 40s
 25

onde V e t são a velocidade do móvel em [m/s] e o tempo em [s], respectivamente.

Pede-se:
d) Faça o gráfico de V=f(t) para 0 ≤ t ≤ 40s
e) Cor de Linha: preta
f) Espessura de linha: 3,0 pontos

23) Faça o gráfico de y=f(x) sabendo que:

 y = 5x ⇒ 0 ≤ x ≤ 20


 − 50 x + 3600
y = ⇒ 20 ≤ x ≤ 35
 35
 y = 50 ⇒ 35 ≤ x ≤ 55

 y = −2 x ⇒ 55 ≤ x ≤ 80

8.5 Respostas dos Exercícios Complementares

1) 2)
500 10

0 8

6
-500

4
-1000
2
-1500
0
-2000
-2
-2500
-4

-3000
-6

-3500 -8

-4000 -10
0 1 2 3 4 5 6 7 8 9 10 -10 -8 -6 -4 -2 0 2 4 6 8 10

3) 4)
 95
450 120

400 100

350
80

300
60
250

200 40

150 20

100
0

50
-20
0 -4 -3 -2 -1 0 1 2 3 4
-20 -15 -10 -5 0 5 10 15 20

5) 6)
8 140

6 120

100
4

80
2

60
0
40

-2
20

-4
0

-6 -20
-5 -4 -3 -2 -1 0 1 2 -2 -1 0 1 2 3 4 5

7) 8) Utilize t=linspace(0,30,300)
Utilize o comando (subplot)
1.5

500 60
1
400 40

300 20 0.5

200 0
0
100 -20

0 -40 -0.5
0 10 20 30 0 10 20 30

-1

-1.5
0 5 10 15 20 25 30
 96
9) 10)
5
sec(teta)
1500 4 cossec(teta)
cotg(teta)
3

1000 2
S (k m )

s x t
1
500
0

-1
0
0 5 10 15 -2
t(h)
-3
150
-4
100
-5
0 50 100 150 200 250 300 350
V(k m /h)

50
vx t
0
11)
-50
0.3

-100
0 5 10 15
t(h)
40
0.25
20
a(k m 2 /h)

0
ax t
-20
0.2
-40

-60
0 5 10 15
t(h)
-10 -8 -6 -4 -2 0 2 4 6 8 10

12) 13)
experiencia 1
20 4000

18
3500
16
3000
Nº de bactérias em milhões

14

12 2500

10
2000
8

1500
6

4 1000

2
500
0
0 1 2 3 4 5 6 7 8
Tempo de em horas 0
50 100 150 200

14) 15)Utilize o comando axis equal

após o comando plot. (consulte help axis)
 97
1.5
1.5

1 1

0.5 0.5

0 0

-0.5 -0.5

-1 -1

-1.5
-1.5
-5 -4 -3 -2 -1 0 1 2 3 4 5 -2 -1.5 -1 -0.5 0 0.5 1 1.5

16) 17)
100

2
80
1.5

1 60

0.5
40

0
20
-0.5

-1 0

-1.5
-20
-0.5 0 0.5 1 1.5
-2
4 6 8 10 12 14 16 18 20

18) 19)
2.4 250

2.2

2
200
1.8
cosh(x)
1.6
150
y z

1.4

1.2
100
1
sech(x)
0.8

0.6 50

0.4
-1.5 -1 -0.5 0 0.5 1 1.5
x
0
0 100 200 300 400 500 600 700 800 900 1000

20) 21)
 98
Lançamento Oblíquo de Projéteis
150 20
15°
30°
45°
60° 15
75°

100
10
deslocamento y

50

-5
0 -3 -2 -1 0 1 2 3 4
0 50 100 150 200 250 300 350
deslocamento x

22) 23)
25 100

90

20 80

70

15 60

50

10 40

30

5 20

10

0 0
0 5 10 15 20 25 30 35 40 0 10 20 30 40 50 60 70 80
 99
9. AJUSTE DE CURVAS e INTERPOLAÇAO UNIDIMENSIONAL

9.1 Ajuste de curvas

p=polyfit(x,y,n)

Encontra os coeficientes do polinômio interpolador p(x) de ordem n obtido através dos

vetores x e y, sendo obrigatoriamente y=f(x). O Ajuste de curvas deve seguir a seguinte estrutura:

Fornecer os dados iniciais.

1° x = [x1 x2 x3 ... xn]
y = [y1 y2 y3 ... yn]

2° Escolher o grau do polinômio.

n = grau do polinômio

3° Encontrar o polinômio interpolador.

p=polyfit(x,y,n)

4° Criar um novo vetor x com maior número de pontos.

xnovo=linspace(x1,xn,n° de pontos);

Calcular o novo vetor y com o polinômio interpolador

5° obtido em p utilizando-se dos pontos de xnovo.
ynovo=polyval(p,xnovo);

Mostra os resultados obtidos. Intertpolar o valor de y

plot(x,y,'cml',xnovo,ynovo,'cml’) para um determinado x
6° c = cor de linha
7°
x1 = valor desejado
m = marcador
y1 = polyval(p,x1)
l = estilo de linha

FIM
 100
1) Sendo considerado os dados abaixo, encontre o polinômio que melhor se ajusta a curva original y
= f(x) e faça os gráficos com os valores originais e ajustados.

X Y
0 -0,84
1 1,61
2 3,28
3 6,16
4 7,10
5 7,34
6 7,66
7 9,56
8 9,48
9 10,14
10 10,20

⇒MATLAB EDITOR
clear all
clc

%Dados originais
1° x=[0 1 2 3 4 5 6 7 8 9 10];
y=[-0.84 1.61 3.28 6.16 7.1 7.34 7.66 9.56 9.48 10.14 10.2];

%Polinômio interpolador
Grau do polinômio interpolador
2° n=2; Escolhido.
Encontra os coeficientes do
3° polinômio interpolador de grau n.
p=polyfit(x,y,n)
Cria o novo vetor x com mais
pontos para melhorar a resolução
4° xnovo=linspace(0,10,50); (precisão) do gráfico.

Calcula o novo vetor y com o

5° ynovo=polyval(p,xnovo); polinômio ajustado a curva
original.

%Mostra os resultados

plot(x,y,'o-',xnovo,ynovo,'r:') Plota o gráfico com a curva

6° title('Curva Ajustada') original (em azul) e ajustada
xlabel('x') (em vermelho).
ylabel('y=f(x)')

⇒RESPOSTA
p=
-0.1194 2.2417 -0.5133 y( x ) = −0,1194 x 2 + 2,2417 x − 0,5133
Esta equação somente é valida
no intervalo do vetor x, ou
seja, 0 ≤ x ≤ 10
 101

2) Para o exercício acima encontrar o valor de y para x=7,25 utilizando o polinômio interpolador
encontrado.

⇒MATLAB EDITOR
Acrescentar estas duas
x1 = 7.25 linhas no final do
y1 = polyval(p,x1) exemplo anterior.
7° Etapa (ver pág 132)
⇒RESPOSTA
x1 =
7.2500

y1 =
9.4647

9.2 Interpolação unidimensional

Comando:
ynovo=interp1(x,y,xnovo,'método de interpolação')

x e y ⇒ variáveis antigas (base da interpolação);

xnovo ⇒ variáveis nova a ser interpoladas;
ynovo ⇒ varíável resultante da interpolação de xnovo na base
de interpolação.

'método de interpolação' ⇒ 'linear' ou

'cubic'

3) De acordo com a tabela abaixo, encontre o valor das grandezas u e v para n1 = 1324 e n2 = 1440
utilizando interpolação linear.

N u v
1250 1,211 89,5
 102
1300 1,816 74,4
1350 2,385 69,3
1400 3,002 64,1
1450 3,650 59,0

⇒MATLAB EDITOR
clear all
clc

n=[1300 1350 1400 1450];

u=[1.816 2.385 3.002 3.650]; Constantes a serem interpoladas
v=[74.4 69.3 64.1 59.0];

u1novo=interp1(n,u,1324,'linear') Interpolação para n1 = 1324

v1novo=interp1(n,v,1324,'linear')

u2novo=interp1(n,u,1440,'linear') Interpolação para n2 = 1440

v2novo=interp1(n,v,1440,'linear')

⇒RESPOSTA
u1novo =
2.0891

v1novo =
71.9520

u2novo =
3.5204

v2novo =
60.0200

4) Em um experimento de física mediu-se o deslocamento (x) de um corpo em função do tempo (t),

e foram encontrados os seguintes valores. Faça o gráfico de x = f(t) e estime um valor para a
distância percorrida quando t = 76s e t = 174s.

t(s) x(m) t(s) x(m)

0 4,9 200 6,80
50 -4,01 250 14,26
100 -5,02 300 18,9
150 -0,8

⇒MATLAB EDITOR (ETAPA I) Primeiramente traça-se o gráfico com

os dados originais para verificar o
clear all
tipo de interpolação (aproximação)
clc necessária.
t=[0 50 100 150 200 250 300];
x=[4.9 -4.01 -5.02 -0.8 6.80 14.26 16.0];

plot(t,x,'o-')
title('x=f(t)')
xlabel('t[s] - tempo')
 103
ylabel('x[m] - distância percorrida')

Escolha do tipo de interpolação

(cubica ou linear), a partir do
⇒MATLAB EDITOR (ETAPA II)
gráfico obtido com os dados
originais na ETAPA I.
x1=interp1(t,x,76,'cubic')
x2=interp1(t,x,174,'linear') A T E N Ç Ã O

Esta operação realizada

incluindo estas duas linhas
⇒RESPOSTA após o final do programa
x1 = inicial. (exemplo anterior)
-5.3478
x2 =
2.8480

9.3 Exercícios propostos

1) Em um experimento no laboratório de física com vibração de um sistema de amortecimento
mola-amortecedor automotivo, obteve-se os seguintes dados:

T(s) A(cm) t(s) A(cm) t(s) A(cm)

0 0 16 -0,1 30 7,5
2 1,2 18 -0,6 32 8,4
4 2,3 20 -0,7 34 8,9
6 3,2 22 -0,3 36 8,6
8 3,2 23 0,6 38 8,0
10 2,8 24 1,5 40 7,0
12 1,8 26 3,4
14 0,8 28 5,5

Onde t é o tempo em [s] e A é a amplitude de vibração em [cm]. Pede-se:

I)Encontre o polinômio que melhor represente os dados experimentais no intervalo de

0 ≤ t ≤ 40s ;
II)Plote o gráfico de A=f(t) com os dados medidos e através do ajuste com o polinômio obtido em I.
 104
III)Encontre o valor da amplitude de vibração para t=32,5s;

2) Verificou-se em um teste de consumo de combustível de um determinado automóvel os seguintes

resultados:

Km rodados Consumo km rodados Consumo

combust. combust.
X1000 km/l X1000 km/l
0 7 50 11,2
5 9 55 11,1
10 10,8 60 10,8
15 11 65 10,2
20 11,2 70 8,5
25 11,2 75 8
30 10,9 80 7,1
35 11,1 85 6,8
40 11,0 90 6,2
45 11,5 95 6,0
100 5,8

Pede-se:
a) Plote simultaneamente em um mesmo gráfico :os dados originais (da tabela) e uma curva
melhor represente estes dados utilizando o ajuste de curvas;
b) Encontre o polinômio interpolador utilizado no ajuste;
c) Encontre o consumo para 27000 e 90300 km rodados.

3) Encontre o polinômio interpolador que melhor represente os dados da tabela abaixo sabendo que
Y = f(t). A partir disso, estime o valor de Y para X = 1,52.

T Y
0 6.00
0.20 4.00
0.40 3.00
0.60 2.50
0.80 2.00
1.00 1.80
1.20 1.40
1.40 1.00
1.60 0.90
1.80 0.80
2.00 0.80
2.20 0.85
2.40 0.83
2.60 0.96
2.80 1.30
3.00 2.00

4) Encontre o polinômio interpolador para os dados da tabela abaixo sabendo que x=f(z). Plote
simultaneamente em um mesmo gráfico os dados originais (da tabela) e uma curva melhor
represente estes dados utilizando o ajuste de curvas.

z x
 105
-3.00 -15.00
-2.33 -3.40
-1.67 2.05
-1.00 3.00
-0.33 1.30
0.33 -1.32
1.00 -3.00
1.67 -2.04
2.33 3.35
3.00 15.00

5) Encontre o valor de s para p1 = 122kPa e p2 = 325kPA de acordo com a tabela abaixo. Utilize os
comandos de interpolação linear.

p(kPa) S(kJ/KgK)
100 542,32
200 618,80
300 689,45
400 755,06
500 825,21

6) Utilizando os comandos de interpolação unidimensional encontre o valor de t para t1 = 0,75; t2 =

1,22; t3 = 3,65 e t4 = 4,85; sabendo que Y=f(t).

t Y
0 0.00
0.5 0.02
1.0 0.27
1.5 0.16
2.0 -2.66
2.5 -10.79
3.0 -24.06
3.5 -37.12
4.0 -39.37
4.5 -18.30
5.0 34.09

9.4 Respostas dos Exercícios Propostos

1)
p=

Columns 1 through 4

0.00000038777784 -0.00004677672376 0.00200177758275 -0.03510492679138

Columns 5 through 7

0.20180762873568 0.19557456464554 0.08291320226615

ou
3,877 ⋅ 10 − 7 x 6 − 4,677 ⋅ 10 − 5 x 5 − 2,001 ⋅ 10 − 3 x 4 − 0,0351 ⋅ 10 − 3 x 3 − 0,201x 2 + 0,195x + 0,0829
 106

Amplitude de vibração para t=32,6s ⇒ 8.7436 cm

10

-2
0 5 10 15 20 25 30 35 40

2)
p=

Columns 1 through 4

-0.0000000006348 0.00000021628002 -0.00002785051797 0.00169564150937

Columns 5 through 7

-0.05147056097900 0.75225601427762 6.80323097192733

ou
6,348 ⋅ 10 −10 x 6 + 2,1628 ⋅ 10 − 7 x 5 − 2,785 ⋅ 10 − 6 x 4 + 1,695 ⋅ 10 − 3 x 3 − 0,0514 x 2 + 0,7522 x + 6,803

12

11

10

5
0 10 20 30 40 50 60 70 80 90 100

consumo para km=27500 ⇒ 11.0286km/l

consumo para km=90300 ⇒ 6.1471lkm/l

3) p =
1.2135 -4.7252 5.2097

Y = 1,2135t 2 − 4,7252 t + 5,2097

t=1.52 e y=0.8311
 107
6

0
0 0.5 1 1.5 2 2.5 3

4)

20

15

10

-5

-10

-15

-20
-3 -2 -1 0 1 2 3

p=

1.0006 -0.0002 -4.0005 -0.0051

5)
p = 122 kPa ⇒ s = 559.1456 kJ/KgK
p = 325 kPa ⇒ s = 705.8525 kJ/KgK

6)
t1 = 0.1450 (linear)

t2 = 0.2260 (linear)

t3 = -39.3241 (cúbica)

t4 = 18.3730 (linear)
 108

10. GRÁFICOS TRIDIMENSIONAIS

10.1 Gráfico de linha

⇒ plot3(x,y,z)
1) Faça o gráfico das equações paramétricas abaixo:

x = sin ( t )
y = cos( t )
z=t
0 ≤ t ≤ 10π
⇒MATLAB EDITOR
clear all
clc

t=linspace(0,10*pi);
x=sin(t);
y=cos(t);
z=t;

plot3(x,y,z) %plot3 cria o gráfico definido por x, y e z.

title('Hélice')
xlabel('sin(t)')
ylabel('cos(t)')
zlabel('t') %insere título ao eixo z

⇒RESPOSTA
Hélice

40

30

20
t

10

0
1
0.5 1
0 0.5
0
Ordena o eixo y -0.5
-0.5 Ordena o eixo x
da direita para cos(t) -1 -1
sin(t)
da esquerda
esquerda. para direita.
Para renumerar o eixo y da
esquerda para direita utilize o
comando axis('ij') após o
comando zlabel
OBSERVAÇÃO:
 109
Todos os comandos para gráficos bidimensionais são válidos para gráficos tridimensionais.

Controle de eixos ⇒ axis([xmin xmax ymin ymax zmin zmax])

Linhas de grade ⇒ grid
Cores–Marcadores–Linhas ⇒ plot3(x,y,z,'s')
Subdivisão da janela gráfica ⇒ subplot(m,n,p)
Espessura de linhas ⇒ set(u,'linewidth',v)
Insere legenda ⇒ legend('gráf.1','gráf.2',...)
Insere texto sobre o gráfico ⇒ gtext('texto')
Criação de janelas gráficas ⇒ figure(i)
Gráficos simultaneos ⇒ plot3(x1,y1,z1,x2,y2,z2,...)

10.2 Gráficos de rede e de superfície

O MATLAB define uma superfície de rede por meio das coordenadas z dos pontos correspondentes
a uma grade retangular no plano xy. Como conseqüência o gráfico é formado através da união dos
pontos adjacentes com linhas retas.

Os gráficos de rede e superfície seguem a seguinte estrutura:

CRIAR OS VETORES BASE

x=(xinício:xincremento:xfim);
y=(yinício:yincremento:yfim);

O comando meshgrid encontra

MONTAR A GRADE RETANGULAR DE todas as combinações possíveis
BASE de x e y.
[x,y]=meshgrid(x,y);

DEFINIR A COORDENADA Z
z=f(x,y);

GRÁFICO

MALHA SUPERFÍCIE
mesh(x,y,z) surf(x,y,z)

 
sin  x 2 + y 2 
2) Plote o gráfico de z=f(x,y), onde z =
 
x 2 + y2
R
 110

Sabendo que - 7,5 ≤ x ≤ 7,5 ou - 7,5 ≤ y ≤ 7,5 .

⇒MATLAB EDITOR
clear all
clc Os pontos do intervalo devem
ser uniformemente espaçados.
x=(-7.5:0.5:7.5);
y=x;
Cria uma matriz com elementos iguais a [x,y]
onde as linhas são cópias do vetor x e as
[x,y]=meshgrid(x,y); colunas são cópias do do vetor y.

R=sqrt(x.^2+y.^2)+eps; Insere a função z=f(x,y) e adiciona

z=sin(R)./R; (eps=2.2204e-016) a R para evitar a
divisão por zero quando x=0 e y=0.
figure(1)
mesh(x,y,z) %plota o gráfico em rede de z=f(x,y)

figure(2)
surf(x,y,z) %plota o gráfico de superfície de z=f(x,y)

mesh
 111

surf

Outro exemplo, para uma malha refinada

>> x=-10:0.5:10;
>> y=x;
>> [X,Y]=meshgrid(x,y);
>> Z=1+X.^2+Y.^2;
>> mesh(X,Y,Z) %gera um gráfico de rede

No MATLAB gráficos de curvas de níveis bi e tridimensionais são gerados usando as funções

contour e contour3, respectivamente. Por exemplo,
>> contour(X,Y,Z,20) % gera 20 curvas de nível bidimensional
 112

10

-2

-4

-6

-8

-10
-10 -8 -6 -4 -2 0 2 4 6 8 10

>> contour3(X,Y,Z,50) % gera 50 curvas de nível tridimensional

200

150

100

50

0
10
5 10
0 5
0
-5
-5
-10 -10

Outra forma interessante de visualizar as informações sobre curvas de nível consiste em usar
cores para representar as alturas. A função pcolor mapeia a altura em um conjunto de cores e
apresenta a mesma informação do gráfico de curvas de nível contour, na mesma escala. Uma vez
que pcolor e contour mostram a mesma informação na mesma escala é útil superpor os dois. Assim,
a função shading será usada para mudar a aparência do gráfico.

>> pcolor(X,Y,Z) %gera o gráfico com pseudocores

 113

>> pcolor(X,Y,Z) % gera o gráfico com pseudores

>> shading interp % remove a grade de linha
>> hold on % comandos gráficos subsequentes serão adicionados ao gráfico existente
>> contour(X,Y,Z,20,'k') % gera 20 curvas de níveis na cor preta
>> hold off % volta ao padrão

O MATLAB permite especificar o ângulo a partir do qual se vê o gráfico tridimensional usando

a função view(azimute, elev), onde azimute define o ângulo de visão e elev a elevação, que descreve
o ângulo em graus a partir do qual o gráfico é observado acima do plano xy. Outra opção é estando
na janela gráfica usar o ícone clique sobre ele com o mouse e vá até a figura girando-a até a
posição desejada.

>> x=-10:10; % malha grosseira

>> y = x;
>> [X,Y]=meshgrid(x,y);
>> f=sqrt(X.^2+Y.^2)+eps;
>> Z=sin(f)./f;
>> surf(X,Y,Z)
>> view(0,0)
 114
>> view(-30,30)

Para obter o azimute e elevação atuais, por exemplo, da figura anterior, usa-se,

>> [az,el]=view
az =
-38

el =
42

A forma view([x,y,z]) define seu ponto de vista usando um vetor que contém as coordenadas
cartesianas (x, y, z) no espaço tridimensional.
>> view([2 3 4])

Existem muitas funções para tratar gráficos no MATLAB, algumas serão citadas a seguir.
A função ribbon(x,y) é semelhante a plot(x,y), a diferença é que as colunas de y são traçadas como
faixas em três dimensões.
A função contourf traça um gráfico de curvas de nível preenchido onde as áreas entre as curvas de
nível são coloridas.
A função peaks é uma função de duas variáveis usada como exemplo pelo MATLAB.

>> x=linspace(-2,2);
>> y=x.^2;
>> ribbon(x,y) % desenha o gráfico em faixas
 115
>> contourf(peaks)

45

4 40

35
3

30
2
25

1 20

15
0
2
10
1 1.4
0 1.2
5
1
-1
0.8
-2 0.6 5 10 15 20 25 30 35 40 45

A função meshc desenha o gráfico de rede e adiciona um gráfico de curvas de nível abaixo dele. A
função meshz desenha o gráfico de rede e desenha um gráfico de cortina, ou plano de referência
abaixo da rede.

>>meshc(peaks)
>> meshz(peaks)

A função waterfall é semelhante a função mesh, a diferença está nas linhas da rede que aparecem
apenas na direção x.
A função surfc desenha um gráfico de superfície e adiciona um gráfico de curvas de nível abaixo
dele. A funçãosurfl desenha um gráfico de superfície e acrescenta contrastes luminosos a partir de
uma fonte de luz. A forma geral é surfl(x,y,z,S,K) onde S é um vetor opcional em coordenadas
cartesianas, S=[Sx Sy Sz], ou em coordenadas esféricas S = [az, el] que especifica a direção da
fonte de luz, caso S não seja definido o padrão para S são 45 graus no sentido anti-horário, a partir
do ponto de vista atual. K é um vetor opcional e define a contribuição em razão da luz ambiental, da
reflexão difusa, da reflexão especular e do coeficiente de espalhamento especular, K =[ka, kb, ks,
espalhamento].

Observação: Digite uma linha por vez.

>> waterfall(peaks)
>> surfc(peaks)
>> colormap(hsv)
>> surfl(peaks),shading interp;
>> colormap(gray)
>> surfl(peaks),shading interp;
 116

10

-5

-10
60
50
40 40
30
20 20
10
0 0

As funções bar3 e bar3h são as versões tridimensionais de bar e barh e pie3 é a versão
tridimensional da função pie.

>> X=-10:10;
>> Y=X.^2;
>> bar3(X,Y)

>> pie3([2 4 3 5])

36%

14%

29% 21%

O MATLAB define um mapa de cores como uma matriz com três colunas. Cada linha da matriz
define uma cor particular, usando números no intervalo de 0 a 1. Esses números especificam o que
se costuma chamar de valores RGB, a intensidade das componentes vermelho, verde e azul de uma
cor. Alguns exemplos representativos desses valores são apresentados na tabela 6.1. Na tabela 6.2
são apresentadas algumas funções do MATLAB que geram mapas de cores predefinidos.

Tabela 10.1 Alguns valores RGB.

 117
Vermelho Verde Azul Cor
0 0 0 Preto
1 1 1 Branco
1 0 0 Vermelho
0 1 0 Verde
0 0 1 Azul
1 1 0 Amarelo
1 0 1 Magenta
0 1 1 Ciano
0,5 0,5 0,5 Cinza-médio

Tabela 10.2 Funções que geram mapas de cores predefinidos.

Função Descrição do mapa de cores
Hsv Escla com cores saturadas
Hot Preto, vermelho, amarelo, branco
Gray Escala linear de tons de cinza
None Escala de tons de cinza levemente azulados
copper Escala linear de tons acobreados
pink Tons pastéis de rosa
white Totalmente branco
flag Vermelho, branco, azul e preto alternados
jet Variante do mapa hsv
prism Prisma
cool Tons de ciano e magenta
lines Usa as mesmas cores do comando plot
colorcube Cubo colorido
summer Tons de amarelo e verde
autumn Tons de vermelho e amarelo
winter Tons de azul e verde
spring Tons de magenta e amarelo

O comando colormap(M) define a matriz M como o mapa de cores que será usado pela figura
atual. Por exemplo, colormap(cool) usa uma versão com 64 cores do mapa cool, já colormap([0 0
1]) usa a cor azul.

As funções plot e plot3 não usam mapas de cores, usam as cores da

tabela de cores do comando plot. A maior parte das demais funções
gráficas, como mesh, surf, contour, fill, pcolor e suas variantes,
usam o mapa de cores atual.

A cor pode ser usada para adicionar informação a gráficos tridimensionais se ela for empregada
para representar uma quarta dimensão. Veja o exemplo a seguir que apresenta algumas maneiras de
se usar as cores para acrescentar novas informações ou para dar ênfase a informações que já fazem
parte do gráfico.

>> x=-8:0.5:8;y=x; % cria um conjunto de dados

>> [X,Y]=meshgrid(x,y); % gera uma malha
>> R=sqrt(X.^2+Y.^2)+eps; % cria dados radiais
>> Z=sin(R)./R; % cria a função (chapéu mexicano)
>> subplot(2,2,1),surf(X,Y,Z), % varia a cor ao longo de z
>> subplot(2,2,2),surf(X,Y,Z,R), % varia a cor com o raio
 118
>> subplot(2,2,3),surf(X,Y,Z,del2(Z)), % varia a cor com a curvatura
>> [dZdx,dZdy]=gradient(Z); % calcula a inclinação da função
>> dZ=sqrt(dZdx.^2+dZdy.^2); % calcula o módulo da inclinação
>> subplot(2,2,4),surf(X,Y,Z,dZ) % varia a cor com a magnitude da inclinação

A função colorbar acrescenta à janela de figuras uma barra de cores vertical ou horizontal com
uma escala de cores, mostrando como a cor varia ao longo do eixo que é usado para defini-la. A
função colorbar(‘h’) coloca uma barra de cores horizontal abaixo do gráfico em uso, colorbar(‘v’)
coloca uma barra de cores vertical à direita de seu gráfico, já colorbar sem argumentos acrescenta
uma barra vertical de cores, se não houver barras em uso, ou atualiza a barra de cores já existente.
>> surf(peaks)
>> colorbar

10.3 Gráficos provenientes de uma matriz de dados

3) A matriz abaixo está armazenado a altura do relevo de um terreno com largura de 400 metros e
comprimento de 600 m.

Eixo X
Largura 400 metros
 119

100 99 100 99 100 99 99 99 100

100 99 99 99 100 99 100 99 99
Eixo Y 100 98 97 97 99 100 100 100 99
101 100 98 98 100 102 103 100 100
Comprimento 102 103 101 100 102 106 104 101 100
600 metros 99 102 100 100 103 108 106 101 99
97 99 100 100 102 105 103 101 100
100 102 103 102 101 101 100 99 99
100 100 101 101 100 100 100 99 99
100 100 100 100 100 99 99 99 99
100 100 100 99 99 100 99 100 99
100 100 100 99 99 100 99 100 99

⇒MATLAB EDITOR
clear all
clc

x=(0:50:400); %O eixo x varia ao longo das colunas de z

y=(0:50:600); %O eixo y varia ao longo das linhas de z

z = [100 99 100 99 100 99 99 99 100

100 99 99 99 100 99 100 99 99
99 99 98 98 100 99 100 100 100
100 98 97 97 99 100 100 100 99
101 100 98 98 100 102 103 100 100
102 103 101 100 102 106 104 101 100
99 102 100 100 103 108 106 101 99
97 99 100 100 102 105 103 101 100
100 102 103 102 101 101 100 99 99
100 100 101 101 100 100 100 99 99
100 100 100 100 100 99 99 99 99
100 100 100 99 99 100 99 100 99
100 100 100 99 99 100 99 100 99];

mesh(x,y,z)
title('Relevo');
xlabel('Largura [m]')
ylabel('Comprimento [m]')
zlabel('Altura do terreno [m]')

⇒RESPOSTA
120
 121

10.4 Interpolação Bidimensional

A interpolação bidimensional é fundamentada nas mesmas idéias básicas da interpolação
unidimensional. Entretanto, a interpolação bidimensional interpola funções de duas variáveis, ou
seja, z=f(x,y).

Comando: znovo=interp2(x,y,z,xnovo,ynovo,'método de interpolação')

x, y e z ⇒ variáveis antigas (base da interpolação);

xnovo e ynovo ⇒ variáveis novas a serem interpoladas;
znovo ⇒ varíável resultante da interpolação de
xnovo e ynovo na base de interpolação.

'método de interpolação' ⇒ A interpolação bidimensional pode ser

do tipo 'linear' ou 'cubic' ou 'nearst'

A interpolação bidimensional deve seguir a seguinte estrutura:

CRIAR OS VETORES BASE

1° x=(xinício:xincremento:xfim);
y=(yinício:yincremento:yfim);

CRIAR A MATRIZ A SER INTERPOLADA

z = [z1,1 z1,2 ... z1,n
2° z2,1 z2,2 ... z2,n
: : :
zn,1 zn,2 ... zn,n];

REFINAR OS VETORES BASE n° de elementos

3° xnovo=linspace(xinício,xfim,n° de elementos); de xnovo deve ser
ynovo=linspace(yinício,yfim,n° de elementos); proporcional a
ynovo.

MONTAR A NOVA GRADE RETANGULAR DE BASE

4° COM OS VETORES REFINADOS
[xnovo,ynovo]=meshgrid(xnovo,ynovo);

5° CÁLCULO DE znovo PROVENIENTE DA INTERPOLAÇÃO DE xnovo e ynovo

znovo=interp2(x,y,z,xnovo,ynovo,'método de interpolação');

6° GRÁFICO

MALHA SUPERFÍCIE
mesh(xnovo,ynovo,znovo) surf(xnovo,ynovo,znovo)
 122
4) Dado o exemplo anterior, (Ex.3), interpole uma nova superfície de malha mais refinada para o
relevo do terreno.

⇒MATLAB EDITOR
clear all
clc

x=0:50:400;
1° y=0:50:600;

2° z=[100 99 100 99 100 99 99 99 100

100 99 99 99 100 99 100 99 99
99 99 98 98 100 99 100 100 100
100 98 97 97 99 100 100 100 99
101 100 98 98 100 102 103 100 100
102 103 101 100 102 106 104 101 100
99 102 100 100 103 108 106 101 99
97 99 100 100 102 105 103 101 100
100 102 103 102 101 101 100 99 99
100 100 101 101 100 100 100 99 99
100 100 100 100 100 99 99 99 99
100 100 100 99 99 100 99 100 99
100 100 100 99 99 100 99 100 99];

3° xnovo=linspace(0,400,36); Eixos x e y refinados

ynovo=linspace(0,600,52);
Malha com todas as
4° [xnovo,ynovo]=meshgrid(xnovo,ynovo); combinações para xnovo e
ynovo
znovo=interp2(x,y,z,xnovo,ynovo,'cubic');
5° A interpolação
bidimensional pode ser
mesh(xnovo,ynovo,znovo) 'linear', 'cubic' ou
title('Relevo'); 'nearst'.
6° xlabel('Largura [m]')
ylabel('Comprimento [m]')
zlabel('Altura do terreno [m]')
 123

⇒RESPOSTA

5) Encontre, interpole, o valor da altura do relevo para as coordenadas (x = 122,5m,y = 433m) e (x

= 300m, y = 237m) do terreno.

⇒MATLAB EDITOR

z1=interp2(x,y,z,122.5,433,'cubic') Acrescentar estas duas linhas no

z2=interp2(x,y,z,300,237,'cubic') final do exemplo anterior.

⇒rESPOSTA

z1 =
101.7862

z2 =
103.7188

10.5 Exercícios propostos

1) Faça o gráfico de linha dada pelas equações paramétricas abaixo:
 124
x = sen 2 t

 2
 y = cos t onde − 2π ≤ t ≤ 2π
z = sen ( t ) + cos(t )


2) Faça o gráfico de linha (na cor magenta) dada pelas equações paramétricas abaixo:

x = (4 + sen (20t )) cos(t )


 y = (4 + sen (20t ))sen ( t ) onde 0 ≤ t ≤ 2π
z = cos(20 t )


3) Faça o gráfico da função z = f(x,y)

36z 2 − 9 x 2 − 4 y 2 − 360 = 0

onde − 8 ≤ x ≤ 8 e − 8 ≤ y ≤ 8

4) Faça o gráfico de linha (na cor magenta) dada pelas equações paramétricas abaixo:
r r r
r ( t ) = cos( t ) i + sen ( t ) j + sen (2 t )k ⇒ 0 ≤ t ≤ 2π

5) Faça o gráfico de linha (utilize linha pontilhada) dada pelas equações paramétricas abaixo:

u = sen (3t ) cos(t )


v = sen (3t )sen ( t ) onde 0 ≤ t ≤ 3π
w = t


6) Faça o gráfico de superfície da função z = f(x,y) definida por:

xy( x 2 − y 2 ) − 10 ≤ x ≤ 10
z= onde 
(x 2 + y 2 ) − 10 ≤ y ≤ 10

7) Faça o gráfico de superfície da função z = f(x,y) definida por:

− 10 ≤ x ≤ 10
z = e − y cos(x ) onde 
− 10 ≤ y ≤ −4
8) Faça o gráfico de superfície da função z = f(x,y) definida por:
 125
−x − 5 ≤ x ≤ 5
z= onde 
x 2 + y2 − 5 ≤ y ≤ 5

9) Faça o gráfico de superfície da função w = f(u,v) definida por:

1 − 5 ≤ u ≤ 5
w = ( u − v − u − v) onde 
2 − 5 ≤ v ≤ 5

10) Em um experimento no laboratório de física sobre condução de calor em uma placa de aço de
100cm x 90cm (comprimento x largura) obteve-se o seguinte campo de temperaturas, abaixo:

95 96 99 99 77 46 23 23 23 15
95 98 99 99 77 48 23 23 23 15
97 97 100 100 98 79 48 23 23 15
100 100 105 105 77 70 48 30 23 24
104 105 110 110 87 70 48 48 35 30
113 115 117 110 99 89 67 63 57 57
114 116 118 118 110 100 87 83 80 79
117 120 120 120 118 110 99 93 85 97
117 120 140 120 118 110 105 100 99 97
115 119 120 120 117 105 105 100 98 96
115 118 118 118 116 100 100 98 97 95

onde as linhas correspondem ao comprimento e as colunas a largura da placa.

Pede-de:
I) Faça os gráficos tridimensionais do campo de temperaturas sem e com interpolação em janelas
gráficas separadas;

II) Encontre os valores da temperatura através da interpolação bidimensional nas posições:

Comprimento = 16,5cm e largura = 31,2cm
Comprimento = 88,3cm e largura = 76,4cm

11) Faça o gráfico da função w=f(u,v):

− u 2 v2
w= + se − 3 ≤ u ≤ 3 e −4≤ v≤4
9 16
12) Plote o gráfico do (Ex.6) em linha VERDE espessura dois pontos.

13) Faça o gráfico da função Z=f(x,y), onde:

Z = sin ( x ) cos( y) + 1 e 0≤x≤π

−π≤ y≤ π

14) Faça o gráfico da função Z=f(x,y),conhecida como Função Banana de Rosenbrock, onde:

Z = 100( y − x 2 ) 2 + (1 − x ) 2 e − 1,5 ≤ x ≤ 1,5

 126
− 4 ≤ y ≤ 10

15) Faça o gráfico do exercício anterior utilizando o comando axis ij para plotar sua vista posterior.

16) Faça o gráfico do hiperbolóide de z=f(x,y), sendo:

− 10 ≤ x ≤ 10
30 + 8x 2 = −6 y 2 − 24z e
− 10 ≤ y ≤ 10

17) A matriz abaixo representa o mapeamento do relevo (altura do relevo) de uma área rural
de x= 360 metros de largura por y= 400 metros de comprimento. Sabe-se que as colunas
representam a largura e as linhas o comprimento.

10 10 10 10 10 10 10 10 10 10 10 10 10
11 11 11 11 20 20 20 12 12 12 12 12 12
11 20 20 30 40 65 30 20 20 30 30 20 12
10 20 40 20 20 20 30 40 50 75 25 15 12
12 12 12 12 0 0 15 25 25 35 40 12 12
10 20 20 20 0 0 20 20 20 30 40 15 15
10 5 5 5 10 10 10 10 10 20 30 30 20
10 5 5 5 10 10 12 12 12 12 25 35 20
15 15 15 15 20 20 30 20 20 30 25 20 20
15 15 15 15 30 50 60 30 15 15 25 25 25
15 15 10 10 10 35 30 30 15 15 25 25 25

Pede-se:
a) Plote na mesma janela utilizando o comando subplot os seguintes gráficos:
Janela 1 - Grafico 3d com os dados originais;
Janela 2 - Gráfico 3d com superfície interpolada com 40 elementos em x e 50 elementos em
y;

b) Encontre também a altura do relevo para as coordenadas: x=105m e y=325,5m.

18) Faça o gráfico de w=f(u,v), sendo:

45 − 8u 2 = −6v 2 + 15 − 24w

onde − 10 ≤ u ≤ 10 e − 10 ≤ v ≤ 10

19) Faça o gráfico da função z=f(x,y), onde:

1 −2≤ x ≤ 2
z= e
1 + x 2 +y2 −1 ≤ y ≤ 2

⇒RESPOSTAS
 127
1) 2)

2 0.8

0.6

0.4

1 0.2

-0.2
0 -0.4

-0.6

-0.8
-1 -1
5

5
-2 0
1
2
3
4

1 -2
-1
0

-3
1 -5 -5
-4

0.5
0.5
0 0

3) 4)

6 1

0.8
5.5
0.6

5 0.4

0.2
4.5
0

4 -0.2

-0.4
3.5 -0.6

-0.8
3
10 -1
1
5 10
0.5 1
0 5 0.8
0.6
0 0 0.4
-5 0.2
-5 -0.2
0
-0.5 -0.4
-10 -10 -0.6
-0.8
-1 -1

5) 6)

10

0
1
0.5 1
0 0.5
0
-0.5
-0.5
-1 -1

7) 8) Utilize o comando rotate figure

 128

9)

10)

ponto1 = 103.9722
ponto2 = 96.8261

12)
 129

0.5

-0.5

-1
5
4
0 2
0
-2
-5 -4

13) 14)

15) 16)

17)
 130

80

60

40

20

0
400
300 400
200 300
200
100
100
0 0

18) 19)
 131
EXERCÍCIOS DE REVISÃO
b 3 − e 2 + 2a + log(10b + a )
1) Sendo a = 10 e b = −1 , encontre o valor de c = .
a+b

35 m − n + e − m sen (2m) + tan(3n )

2) Sabendo que m = -10 e n = 5log(3), encontre o valor de x = + .
n−m 5m + ln(10n )
20 +
2
Solução no formato de 5 dígitos + expoente,

3) Num triângulo ABC, retângulo em A, a hipotenusa é a = 20 cm e cos(B) = 0,96. Calcule o

perímetro do triângulo.

4) Determine o valor de a sabendo que x =[1 1,2 2,5 3,5] e y=[0,4 2 5 10]. Apresente a solução com
16 dígitos sem expoente.
x3 x− y
a= 2
− + ln( x / 2 + y ) − cos( 2 / x )
y 3
( 2 x + y )2

5) Dados os vetores a = (1, -2, 3) e b = (-3,2,0).

a) Determine |-2b|, |a-b|, 3a –2b, a•b e a × b.
b) Verifique se os vetores a e b são perpendiculares.
(1+ ab )
2a + b y = ab − a
2 2
e
c) Calcule: x = , + 3b
ab log(5a + 10) 4

6) Ache o ângulo entre os vetores a = (2, 5, -4) e b = (1, -2, -3) em graus. Determine o vetor
projeção de a sobre b. Determine a norma do vetor projeção.
 r r 
 a•b 
θ = arccos r r 
 a ⋅ b 
 
r r r
r a •b
w = projb a = r r b
|| b || .|| b ||

r r
7) Sendo u = î − 2 ĵ + k̂ e v = − î + ĵ vetores, ache:
r r
a) a medida do ângulo entre os vetores u e v .
r r
b) a medida da projeção do vetor v sobre o vetor u .
r r r r
c) a área do paralelogramo limitado pelos vetores u e v . Area = || u × v ||

 y
8) Faça o gráfico da função z = f(x,y), onde z = xsen  + ysen(2 x) e 0 ≤ x ≤ 5π e 0 ≤ y ≤ 5π .
2

9) Desenhe na mesma janela gráfica utilizando o comando subplot os seguintes gráficos:

1
(a) f ( x, y ) = 2 e (b) f ( x, y ) = ln( x 2 + y 2 ) para − 5 ≤ x ≤ 5 e − 5 ≤ y ≤ 5 .
4x + y 2
 132
11. GRÁFICOS EM COORDENADAS POLARES

⇒ polar(theta,raio)

y
P(θ,r)

raio
θ

1) Desenhe o gráfico da função r=f(θ) conhecida como LEMNISCATA, dada por:

r 2 = a 2 cos(2θ)
onde a=1 e 0 ≤ θ ≤ 2π .

⇒MATLAB EDITOR
clear
clc

a=1; %raio da alça

theta=linspace(0,2*pi);
r=sqrt(a^2*cos(2*theta));

polar(theta,r)
title('LEMNISCATA')

⇒RESPOSTA
LEMNISCATA
90
1
120 60
0.8

0.6
150 30
0.4

0.2

alça
180 0

210 330

240 300
270
 133
2) Desenhe o gráfico com os dados armazenados na matriz came em coordenadas polares sabendo
que: as colunas ímpares são as coordenadas angulares em graus e as colunas pares são as
coordenadas lineares em cm.

⇒MATLAB EDITOR
clear
clc

came=[ 0 9.0 120 8.5 250 5.4

10 9.0 130 8.25 260 5.5
20 9.0 140 7.75 270 5.85
30 9.0 150 7.25 280 6.25
40 9.0 160 6.75 290 6.75
50 9.0 170 6.25 300 7.25
60 9.0 180 5.85 310 7.75
70 9.0 190 5.5 320 8.25
80 9.0 200 5.4 330 8.5
90 8.95 220 5.3 340 8.75
100 8.9 230 5.3 350 8.9
110 8.75 240 5.3 360 8.95];

%Desacopla a matriz came em dois vetores linha: theta e raio

theta = [came(:,1)' came(:,3)' came(:,5)']
raio = [came(:,2)' came(:,4)' came(:,6)']

polar(theta*pi/180,raio)
title('Came de disco')

Came de disco
90
10
120 60
8

6
150 30
4

180 0

210 330

240 300
270

3) Desenhe o gráfico r=f(θ) em coordenadas polares das equações abaixo:

a) Rosa de três folhas ⇒ r = a cos(3θ) quando a=5 e 0 ≤ θ ≤ π

b) Cardióide ⇒ r = a (1 + cos α) quando a=2 e 0 ≤ α ≤ 2π
 134
3A 3B
90 5 90 4
120 60 120 60
3

150 2.5 30 150 2 30

1

180 0 180 0

210 330 210 330

240 300 240 300

270 270

4) Desenhe o gráfico em linha vermelha (espessura 2 pontos) com os dados armazenados na matriz
dados em coordenadas polares sabendo que: as colunas ímpares são as coordenadas angulares em
graus e as colunas pares são as coordenadas lineares em cm.

dados =[ 0 10.0 60 5.8 140 5.0 230 5.0 310 5.8

5 9.9 70 5.4 150 5.0 240 5.0 320 6.5
10 9.7 80 5.2 160 5.0 250 5.0 330 7.4
15 9.3 90 5.0 170 5.0 260 5.0 340 8.4
20 8.8 100 5.0 180 5.0 270 5.0 345 8.9
30 7.9 110 5.0 190 5.0 280 5.0 350 9.3
40 7.0 120 5.0 200 5.0 290 5.2 355 9.7
50 6.4 130 5.0 220 5.0 300 5.4 360 9.9]

90 10
120 60

150 5 30

180 0

210 330

240 300
270
 135
12. Polinômios

12.1 Representação de polinômios

No MATLAB, um polinômio é representado por um vetor linha contendo seus coeficientes em
ordem decrescente.

4 3
1) p = 2 x − 12x + 25x + 116

» p=[2 -12 0 25 116]

p=
2 -12 0 25 116

5 3
2) y = −3x + x − 4 x + 12

» y=[-3 0 1 0 -4 12]

y=
-3 0 1 0 -4 12

OBSERVAÇÃO:
Complementa-se os termos faltantes com o algarismo zero.

12.2 Raízes de um polinômio

⇒ roots(x)

2 5 4 3 2
3) p = 4 x − 3x − 1 4) y = 5x − 10x + x − x + 10

» p=[4 -3 -1]; » y=[5 -10 1 -1 0 10];

» roots(p) » roots(y)
ans =
ans = 1.7258
1.0000 1.3377 Raízes
-0.2500 -0.0928 + 0.9890i complexas
-0.0928 - 0.9890i conjugadas
-0.8780

12.3 Operações básicas

a ( x ) = x 3 + 2 x 2 − 3x + 4
Sendo dois polinômios distintos:
b( x ) = x 3 + 4x 2 − 9 x + 16
A) MULTIPLICAÇÃO ⇒ conv(x1,x2)

a × b = ( x 3 + 2x 2 − 3x + 4) × ( x 3 + 4 x 2 − 9x + 16)
 136
» a=[1 2 -3 4];
» b=[1 4 -9 16];
» c=conv(a,b)

c=
1 6 -4 -10 75 -84 64

6 5 4 3 2
O resultado acima é ⇒ c( x ) = x + 6x − 4x − 10x + 75x − 84x + 64

B) ADIÇÃO E SUBTRAÇÃO

O MATLAB não possui função direta para somar polinômios. O operador padrão de soma ou
subtração podem ser utilizados se os dois vetores polinomiais a(x) e b(x) possuírem a mesma
dimensão.

a + b = ( x 3 + 2 x 2 − 3x + 4) + ( x 3 + 4x 2 − 9 x + 16)

» d=a+b
d=
2 6 -12 20 d ( x ) = 2 x 3 + 6 x 2 − 12 x + 20

a − b = ( x 3 + 2 x 2 − 3x + 4) − ( x 3 + 4 x 2 − 9 x + 16)

» e=a-b

e=
0 -2 6 -12 e( x ) = −2 x 2 + 6x − 12

OBSERVAÇÃO:
Quando dois polinômios forem de ordens diferentes, aquele que tiver menor ordem deve
ser preenchido com zeros, com o objetivo de torná-los da mesma ordem do polinômio de
ordem mais alta.

f = c + d = ( x 6 + 6 x 5 − 4 x 4 − 10 x 3 + 75x 2 − 84 x + 64) + (2 x 3 − 6 x 2 − 12x + 20)

» f=c+[0 0 0 d] c(x) possui 7 elementos e d(x) 4 elementos

f=
1 6 -4 -8 81 -96 84

6 5 4 3 2
O resultado acima é ⇒ f ( x ) = x + 6x − 4 x − 8x + 81x − 96 x + 84

C) DIVISÃO ⇒ [q,r]=deconv(x1,x2).

a ( x ) = x 3 + 2 x 2 − 3x + 4
 137
Sendo dois polinômios distintos:
b( x ) = 2 x − 1
a ÷ b = ( x 3 + 2 x 2 − 3x + 4) ÷ (2 x − 1)

» a=[1 2 -3 4];
» b=[2 -1]; a b
» [q,r]=deconv(a,b)
r q
q=
0.5000 1.2500 -0.8750 q ( x ) = 0,5x 2 + 1,25x − 0,8750

r=
3,125
0 0 0 3.1250 r(x ) =
2x − 1

OBSERVAÇÃO:
O polinômio do denominador (b), nunca pode começar com o algarismo zero.

12.4 Derivadas
A derivada dy/dx (derivada de y em relação a x), de um polinômio y=f(x) é dada pelo comando:

polyder(y)

5 4 3 dy
1) Sendo o polinômio y( x ) = 2x − 4 x + 3x − 10x + 1 , encontre o valor de :
dx
» y=[2 -4 3 0 -10 1];
» dy=polyder(y)

dy =
dy
10 -16 9 0 -10 = 10 x 4 − 16x 3 + 9x 2 − 10
dx

d2y
2) Para o mesmo polinômio do exemplo 1 calcule o valor de :
dx 2
» d2y=polyder(dy)

d2y =
dy
40 -48 18 0 = 40 x 3 − 48x 2 + 18x
dx
ou

» d2y=polyder(polyder(y))
 138

d2y =
40 -48 18 0

12.5 Cálculo de polinômios

O cálculo do valor numérico do polinômio p(x) resultante da substituição do vetor x em p(x) dado
por:

r = polyval(p,x)

onde r armazena o resultado dessa substituição.

3) Faça o gráfico da função p( x ) = x 3 + 4x 2 − 6 x − 12 , onde x varia de [-4,3];

⇒MATLAB EDITOR
clear all
clc

x=linspace(-4,3);
p=[1 4 -6 -12];
r=polyval(p,x);

plot(x,r)
title('p(x)=x^3+4x^2-6x-12')
xlabel('x')
ylabel('p(x)')

⇒MATLAB EDITOR (SOLUÇÃO SIMPLIFICADA – AULA 05)

clear all
clc
x=linspace(-4,3);
p=x.^3+4*x^2-6*x-12;
plot(x,p)
title('p(x)=x^3+4x^2-6x-12')
xlabel('x')
ylabel('p(x)')

12.6 Exercícios propostos

1) Encontre as raízes dos seguintes polinômios:

a) x 3 − 3x 2 − 7 x + 12
b) 7 x 5 + 3x 4 + 2 x − 10

2) Dados os polinômios
 139

a ( x ) = −2 x + 4 b( x ) = x 6 − 3x 4 + 3x 3 − x + 7
c( x ) = 2 x 3 − 4 x 2 + x + 9 d ( x ) = 2x − 1

execute as operações abaixo:

d ( a × c) d 2 (c )
a)a + b f)a x b j)b – c l) m)
b)b - a g)c / a d ( b) dx dx 2
c)a + d h)b / c k)
dx
d)a - d i)a + c
e)c x d

3) Faça o gráfico do polinômio abaixo com utilizando a função polyval quando − 1,5 ≤ t ≤ 2,5 .
p( t ) = t 5 − t 4 − 7 t 2 − 10

12.7 Respostas dos exercícios propostos

1.a)
ans =
4.0000
1.3028
-2.3028

1.b)
ans =
-0.9836 + 0.6647i
-0.9836 - 0.6647i
0.2914 + 0.9877i
0.2914 - 0.9877i
0.9558

2.a)
ans =
1 0 -3 3 0 -3 11 ⇒ x 6 − 3x 4 + 3x 3 − 3x + 11

2.b)
ans =
1 0 -3 3 0 1 3 ⇒ x 6 − 3x 4 + 3x 3 + x + 3

2.c)
ans =
0 3 ⇒ 0x + 3
2.d)
ans =
-4 5 ⇒ − 4x + 5

2.e)
ans =
4 -10 6 17 -9 ⇒ 4 x 4 − 10 x 3 + 6 x 2 + 17 x − 9

2.f)
 140
ans =
-4 10 -4 ⇒ − 4 x 2 + 10 x − 4

2.g)
q=
-1.0000 0 -0.5000

r=
0 0 0 11

11
⇒ − x 2 − 0,5 +
− 2x + 4

2.h)
q=
0.5000 1.0000 0.2500 -0.7500

r=
0 0 0 0 -12.2500 -2.5000 13.7500

− 12,25x 2 − 2,5x + 13,75

⇒ − 0,5x 3 + x 2 + 0,25x − 0,75 +
2x 3 − 4x 2 + x + 9
2.i)
ans =
2 -4 -1 13 ⇒ 2 x 3 − 4 x 2 − x + 13
2.j)
ans =
1 0 -3 1 4 -2 -2 ⇒ x 7 − 3x 4 + x 3 + 4 x 2 − 2 x − 2
2.k)
ans =
6 0 -12 9 0 -1 ⇒ − x 5 − 12 x 3 + 9 x 2 − 1
2.l)
ans =
-16 48 -36 -14 ⇒ − 16 x 3 + 48x 2 − 36 x − 14

2.m)
ans =
12 -8 ⇒ 12 x − 8

3)
5

-5

-10

-15

-20

-25

-30

-35

-40
-1.5 -1 -0.5 0 0.5 1 1.5 2 2.5
 141
13. MATEMÁTICA SIMBÓLICA

13.1 Representação e visualização de variáveis simbólicas

⇒Armazena em a variável simbólica a.

a=sym('a')
⇒Cria mais de uma variável simbólica, é equivalente a:
a=sym('a')
syms a b c,... b=sym('b')
c=sym('c')

⇒Faz a visualização da expressão simbólica. (É um

pretty(x) comando facultativo que ajuda a visualização de
resultados simbólicos)

2
1) Implemente a expressão simbólica z = 2x + 3y utilizando o comando sym.

⇒MATLAB EDITOR ⇒RESPOSTA

clear all z=
clc 2*x+3*y^2

x=sym('x');
y=sym('y');

z=2*x+3*y^2

2v p
2) Implemente a expressão simbólica w = + utilizando o comando syms e faça sua
u u+v
visualização utilizando o comando pretty.

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc / p \1/2
2 v/u + |-----|
syms u v p \u + v/
w=2*v/u+sqrt(p/(u+v));

pretty(w)

2x 2
3) Sendo a função simbólica f = , faça sua visualização utilizando o comando pretty.
x + 2x

⇒MATLAB EDITOR ⇒RESPOSTA

 142
clear all
clc 2
x
x=sym('x'); 2 -------------
f=(2*x^2)/(x+sqrt(2*x)); 1/2 1/2
x+2 x
pretty(f)

Observação:
Pode-se implementar uma expressão simbólica digitando-a entre apostos (' '). Entretanto, o
comando pretty não pode ser utilizado e a expressão digitada deve ser dependente de uma única
variável como por exemplo: y=f(x).

13.2 Substituição de variáveis

⇒ subs(f,antiga,nova)

A) Substituindo uma única variável.

O comando subs substitui na função simbólica f a variável nova (numérica ou não) na variável
antiga (simbólica).

2x 2
4) Sendo a função f = onde x=2,5. Encontre o valor numérico de f utilizando os
x + 2x
recursos de matemática simbólica do MATLAB.
⇒MATLAB EDITOR ⇒RESPOSTA

clear all ans =

clc 2.6393

x=sym('x');
f=(2*x^2)/(x+sqrt(2*x));

subs(f,x,2.5)

B)Substituindo mais de uma variável.

2 xy 2
5) Sendo a função f = , pede-se:
y−x

a)Faça a visualização da expressão utilizando o comando pretty;

b)Encontre o valor numérico de f utilizando os recursos de matemática simbólica do MATLAB

sabendo que x e y valem 3 e –5, respectivamente.

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc 2
 143
xy
syms x y 2 -----
f=(2*x*y^2)/(y-x); y-x

pretty(f) f1 =
6*y^2/(y-3)
f1=subs(f,x,3)
f2=subs(f1,y,-5) f2 =
-18.7500

13.3 Derivadas.

diff(f,v) ⇒Calcula a derivada da função f em relação a variável v,

df
ou seja:
dv

diff(f,v,n) ⇒Calcula a derivada de ordem n da função f em relação a

dnf
variável v, ou seja:
n
dv

13.4 Comando simplify

⇒Executa a simplificação da função f proveniente de

simplify(f)
qualquer operação simbólica.

dy 2 2
6) Encontre a da função y dada por: y = x +
dx x

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc dy =
2*x-2/x^2
x=sym('x');
y=x^2+2/x;
2
dy=diff(y,x) 2 x - ----
2
pretty(dy) x

dy x 3 − cos(x )
7) Encontre a da função y dada por: y =
dx x

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc 2 3
 144
3 x + sin(x) x - cos(x)
x=sym('x'); ------------- - -----------
y=(x^3-cos(x))/x; x 2
x
dy=diff(y,x);
pretty(dy)

8) Simplifique o resultado simbólico do exemplo acima utilizando o comando simplify.

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc 3
2 x + x sin(x) + cos(x)
x=sym('x'); ------------------------
y=(x^3-cos(x))/x; 2
x
dy=diff(y,x);

dy=simplify(dy)
pretty(dy)

d2y x +1
9) Encontre a da função y dada por: y =
dx 2 x

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc 2 x+1
- ---- + 2 -----
x=sym('x'); 2 3
y=(x+1)/x; x x

d2y=diff(y,x,2);
pretty(d2y) 2
----
d2y=simplify(d2y); 3
pretty(d2y) x

dy xa 2 + 3ax
10) Encontre a derivada cuja função y é definida por: y =
da a

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc 2
2ax+3x xa +3ax
syms x a ----------- - ------------
y=(x*a^2+3*a*x)/a; a 2
a
dy=diff(y,a);
pretty(dy) x
 145
Observe que a solução acima é equivalente a:
dy=simplify(dy);
pretty(dy)
2ax + 3x xa 2 + 3ax 2 xa 2 + 3xa − xa 2 − 3xa xa 2
− = = =x
a a 2 a 2 a 2

13.5 Comandos rats e eval

⇒Transforma a variável numérica x em um número

rats(x)
racional.
⇒Transforma a variável numérica x racional em um
eval(x)
número real.

dy
11) Encontre a ( x = 3,2) da função y definida por:
dx
x 2 + 3x + 11
y=
x −1

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc resp1 =
-1.8926
x=sym('x');
y=(x^2+2*x+11)/(x-1);
dy=diff(y,x); resp2 =
-229/121
resp1=subs(dy,x,3.2)

resp2=rats(resp1) resp3 =
-1.8926
resp3=eval(resp2)

dy
12) Encontre a ( x = 10) da função y definida por: y = log(x )
dx

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc 1125899906842624
---------------- 1/x
x=sym('x'); 2592480341699211
y=log10(x);
dy=diff(y,x); resp =
pretty(dy) 0.0434

resp=subs(dy,x,10) A derivada de y = log(x ) é equivalente a:

dy x' dy 1 1 1
= ⇒ = =
dx x log(e) dx x log(e) log(e) x

1 1125899906842624
onde = =0,4343
log(e) 2592480341699211
 146

d2y π
13) Encontre a ( x = 10, a = ) da função y definida por:
da 2 4
y = sen 2 (3ax )

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc resp =
-1800
syms a x
y=sin(3*a*x)^2;
d2y=diff(y,a,2);

resp1=subs(d2y,x,10);
resp2=subs(resp1,a,pi/4)

13.6 Integrais definidas e indefinidas

⇒Calcula a integral indefinida da função f em relação a

∫ fdv
int(f,v)
variável v, ou seja:
⇒Calcula a integral definida da função f em relação a
variável v, avaliada entre o limite inferior a e o limite
int(f,v,a,b) b
superior b, ou seja:
∫ a
fdv

∫
x
14) Encontre o valor de y = dx
x +1

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc itg =
x - log(x+1)
x=sym('x');
y=x/(x+1); Equivalente a x − ln(x + 1)

itg=int(y,x)

 1 
15) Encontre o valor de y =
∫  x + 3  dx
 x 
 147

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc 3/2 1
2/3 x - 1/2 ----
x=sym('x'); 2
y=sqrt(x)+1/x^3; x

itg=int(y,x);
pretty(itg)

16) Encontre o valor de y =

∫ cos 2 (a )da

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc
1/2 cos(a) sin(a) + 1/2 a
a=sym('a');
y=cos(a)^2;

itg=int(y,a);
pretty(itg)

3
x3
17) Encontre o valor de y = ∫ dx
x+3
−2

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc
itg =
x=sym('x'); 295/6-27*log(2)-27*log(3)
y=x^3/(x+3);

itg=int(y,x,-2,3);

295
Equivalente a: − 27 ln(2) − 27 ln(3)
6

18) Encontre o valor numérico da integral do Exemplo17

⇒MATLAB EDITOR ⇒RESPOSTA
clear all
clc
295/6 - 27 log(2) - 27 log(3)
x=sym(‘x’);
y=x^3/(x+3);
resp =
 148
itg=int(y,x,-2,3); 0.7892
pretty(itg)

resp=eval(itg)

13.7 Comando solve

A) Raízes de uma equação simbólica

19) Determine as raízes da equação do segundo grau: x 2 − 2 x − 15

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc resp =
[ -3] ⇒ raiz x1
x=sym('x'); [ 5] ⇒ raiz x2
y=x^2-2*x-15;

resp=solve(y)

20) Encontre o valor de x para que satisfaça a equação 4 cos(x ) = 0

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc resp =
1/2*pi
x=sym('x');
y=4*cos(x);

resp=solve(y)

21) Encontre o valor de x para que satisfaça a equação 4 cos(x ) + sen(x ) = 1

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc resp =
[ 1/2*pi]
x=sym('x'); [ -atan(15/8)]
y=4*cos(x)+sin(x)-1;

resp1=solve(y)

22) Encontre o valor numérico da resposta anterior.

⇒MATLAB EDITOR ⇒RESPOSTA
clear all
clc resp1 =
[ 1/2*pi]
 149
x=sym('x'); [ -atan(15/8)]
y=4*cos(x)+sin(x)-1;
resp2 =
resp1=solve(y) 1.5708
resp2=eval(resp1) -1.0808
resp3=resp2*180/pi
resp3 =
90.0000 ⇒ 90°
-61.9275 ⇒ –61.9275° ou
298.0725°

B) Resolução de um sistema de equações lineares

23) Resolva o sistema de equações abaixo:

x+y=3
4 x + y = −3

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc x=
-2
syms x y y=
5
eq1=x+y-3;
eq2=4*x+y+3; ⇒ O matlab apresenta a solução na ordem
crescente, de acordo com a ordem de
[x,y]=solve(eq1,eq2) procedência das variáveis nas equações.

24) Resolva o sistema de equações que segue:

x+y−z=4
Sistema
possível e x + 2 y + 3z = 2
determinado
2 x + 3z = −1

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc x=
1
syms x y z
y=
eq1=x+y-z-4; 2
eq2=x+2*y+3*z-2;
eq3=-3*x-3*y+3z+3; z=
-1
[x,y,z]=solve(eq1,eq2,eq3)

25) Resolva o sistema de equações que segue:

 150
x+y−z=4
Sistema
possível e 2 x + 2 y − 2z = 8
indeterminado
2 x + 3z = −1

⇒MATLAB EDITOR ⇒RESPOSTA

clear all
clc x=
-3/2*z-1/2
syms x y z
y=
eq1=x+y-z-4; 5/2*z+9/2
eq2=2*x+2*y-2*z-8;
eq3=2*x+3*z+1; z=
z
[x,y,z]=solve(eq1,eq2,eq3)

26) Resolva o sistema de equações que segue:

Sistema 2a − b = 0
impossível
2a − b = −1

⇒MATLAB EDITOR
clear all
clc

syms a b

eq1=2*a-b;
eq2=2*a-b+1;

[a,b]=solve(eq1,eq2)

⇒RESPOSTA
Warning: Explicit solution could not be found.
> In C:\MATLAB\toolbox\symbolic\solve.m at line 130
In C:\MATLAB\toolbox\symbolic\@sym\solve.m at line 49

a=

[ empty sym ]

b=

[]

13.7 Exercícios Propostos

 151
dy
1) Encontre a derivada das funções simbólicas abaixo:
dx
3x 2 − x
a) y = x 5 − 3x 2 − x + 11 b) y =
1− x
cos 2 x + 1
c) y = cos( x ) / sen( x ) d) y =
sen( x )
x −1
e) y = f) y = sec 2 ( x )
ex
cos sec( x )
g) y =
x

dy
2) Encontre a derivada das funções simbólicas abaixo:
dz
xz x 2z tan( xz 2 w )
a) y= − b) y =
z + x z −1 z

3) Dada as funções abaixo:

ax 2 + 3ax
f = g = 3 wz cos( wz) h = −a 3 x 5 + 3a 3 x 2 + ax 3
xa 2

Resolva:
df df 1
a) b) (a = ; x = 1)
dx da 2

dg dg π π
c) d) (z = ; w = )
dw dz 6 3
d 3h d 2h
e) f) (a = −1; x = 3,5)
dx 3 da 2

4) Resolvas as integrais abaixo.

( x + 1) 2 3 3 7 2 
a)
∫ dx
x
b)
∫  4 z
− z dz
8 
1
∫ ∫−3 (x
c) tan( x )dx d) 2 + x 2 + 1)dx
π

∫e
e) y cos(4 y)dy

π/3

5) Encontre o valor das variáveis (x ou θ) nas equações abaixo utilizando o comando solve:
 152
2 3 6 4 8 3 8
a) − x − x + x + =0 b) − 6 x 4 + 32x 3 − 56x 2 + 32x = 0
15 5 5 3
2
c) 3( x +1) + 3 x = 27 d) log(x + 1) + log(x − 1) = 2
log x 2
e) sen(θ 2 − θ − π / 4) = 1 / 2 f) 2 2 = 1/ 4

g) 7 − 3x + 6 = 1 / 343 h) 4 x2 +2 x 2 +3
+ 3⋅ 2 = 160
i) 3 x + 2 − j) (1 / 3)log 2 (log 5 x ) = 3
9
= 8 ⋅ 3x
3x

6) Discuta e resolva os sistema de equações abaixo, utilizando o comando solve.

4 x − y = 0
x + y + z + w = 5 3m + n − 3o = 2  x + y + 2 z = 32
  
a) 2 x − y − 2 w = −1  1 2 c) 5 x − y + 2 z = 0
b) m + n − o =
 3 3 − 3 x − 3 y − 6 z = 1
1  
3 z + 2 w =  m + 3n = 4
 2 

RESPOSTAS:
4
1a) 5x -6x-1

2
-6 x + 3 x + 1
1b) - ---------------
2
(-1 + x)

2
cos(x)
1c) -1 - -------
2
sin(x)

2
cos(x) (cos(x) - 3)
1d) - --------------------
2
-1 + cos(x)

1/2

(-3 + 2 x) exp(-x)
1e) - 1/2 ------------------
1/2
(x - 1)

sin(x)
 153
1f) 2 -------
3
cos(x)

cos(x) x + sin(x)
1g) -----------------
2 2
x (-1 + cos(x) )

2 2 2
x (2 z - 2 z + 1 + 2 x z + x )
2a) --------------------------------
2 2
(z + x) (z - 1)
2
2 2 tan(x z w)
2b) 2 (1 + tan(x z w) ) x w - -----------
2
3a) z

1/2 1/2
2 x (a x) - 3
1/2 -------------------
1/2
a x (a x)

3b)

-18.6969

3c)
(w z) (w z)
3 z log(3) cos(w z) - 3 sin(w z) z

3d)
0.7963

3e)
3 2
-60 a x + 6 a

3f)
2.9308e+003

4a)
5/2 3/2 1/2
2/5 x + 4/3 x +2x
4b)
4 3
3/16 z - 7/24 z
 154
4c)
-log(cos(x))
4d)
16.1338

4e)
2.0257

5a) [ -10]
[ 2]
[ -1]

5b) [ 0]
[ 4/3]
[ 2]
[ 2]

5c)
1.3553

5d)
10.0499

5e)
100.1873°
-42.8915°
5f) –0.5000 e 0.5000

5g) 3

5h) 1.4142135623730950488016887242097

5i) 1

5j) 2.2361

6a) SPD
x = -1/2 y = 1 z = 4 w = 1/2

6b) SPI
m = -3*n+4 n=n o = -8/3*n+10/3

6c) SI
 155
EXERCÍCIOS DE REVISÃO

1) Desenhe todos os gráficos de r = f(θ) em coordenadas polares na mesma janela gráfica.

r = cos(aθ) quando a = 2, 3 e 5, 0 ≤ θ ≤ 2π . Escolha uma cor diferente e uma linha diferente para
cada gráfico.

2) Dado os polinômios
p(x) = 2 x 4 − 3 x 2 + 12 f(x) = 5 x 5 − 2 x 4 + 3 x g ( x) = 10 x 3 + 2 x 2 + 10 x + 9

execute as operações abaixo (+, -, conv, [q,r]=deconv(a,b), polyder, polyval):

a) p + g b) p – g c) f + g d) f/g e) p x g
d ( p) d ( p × g) 2
d (g)
f) g) h) i) p/g j) p(2)
dx dx dx 2
 156
14. OPERADORES RELACIONAIS E LÓGICOS

14.1 Operadores de decisão

A finalidade dos operadores lógicos e relacionais é fornecer respostas do tipo Falso/Verdadeiro a
perguntas. A saída de todas as expressões lógicas e relacionais produz:

⇒ 1 para VERDADEIRO
⇒ 0 para FALSO

Tipo Descrição Operador

Menor que <
Menor ou igual a <=
Relacionais

Maior que >

Maior ou igual a >=
Igual a ==
Diferente de ~=
E &
Lógicos

Ou |
Não ~

» a=(1:8)
a=
1 2 3 4 5 6 7 8

» b=(2:3:20)
b=
2 5 8 11 14 17 20

» a>4 ⇒a > 4

ans =

0 0 0 0 1 1 1 1

1 2 3 4 5 6 7 8
» b>=5 ⇒b ≥ 5

ans =

0 1 1 1 1 1 1

2 5 8 11 14 17 20
 157

» (b>=5)&(b<13) ⇒ (5 ≤ b < 13)

0 1 1 1 0 0 0

2 5 8 11 14 17 20

OBSERVAÇÃO:
No exemplo acima, nunca faça:

» (5<=b<13) ⇒é diferente de (5 ≤ b < 13) visto no Cálculo

1 1 1 1 1 1 1

2 5 8 11 14 17 20

» (b<=6)|(b>=12) ⇒ (b ≤ 6) ou (b ≥ 12)

ans =

1 1 0 0 1 1 1

2 5 8 11 14 17 20

» z = (b~=8)&(b<=12) ⇒ (b ≠ 8) e (b ≤ 12) (indexando a resposta

na variável z)

z=

1 1 0 1 0 0 0

2 5 8 11 14 17 20
 158
15. CONTROLES DE FLUXO

Permitem a execução de comandos, com base em estruturas de tomada de decisões.

Loops for

Loops while
Controles de Fluxo
Estrututas if – else – end

Estruturas switch – case

15.1 Loops for

Os loops for possibilitam que uma série de comandos seja repetida por um determinado número de
vezes pré-definido. A forma geral é:

for vetor

comandos...

end

Os comandos entre as instruções for e end são executados uma vez para cada elemento do vetor.

1) Crie um vetor dado pela expressão A(m)=2m-5 onde m=1,2...10, utilizando o loop for.

⇒MATLAB EDITOR
clear all
clc

for m=1:10
A(m)=2*m-5;
end
A

⇒RESPOSTA
A=
-3 -1 1 3 5 7 9 11 13 15
 159
⇒TRADUÇÃO
clear all
clc

for m=1
A(1) = 2*1 - 5 = -3; clear all
looping clc

for m=2 for m=1:10

A(2) = 2*2 - 5 = -1; A(m)=2*m-5;
looping end
A
for m=3
A(3) = 2*3 - 5 = 1;

loopings

for m=10
A(10) = 2*10 - 5 = 15;
end
A

2) Crie a matriz M de ordem 4 dada pela expressão M(i,j)=2i-j, utilizando os loops for.

⇒MATLAB EDITOR
clear all
clc

for i=1:4 %Cria as linhas

for j=1:4 %Cria as colunas
M(i,j)=2*i-j;
end
end

⇒RESPOSTA
M=
1 0 -1 -2
3 2 1 0
5 4 3 2
7 6 5 4

⇒TRADUÇÃO
clear all
clc

for i=1
for j=1
 160
M(1,1) =2*1 – 1 = 1;
looping j

for j=2
M(1,2) =2*1 – 2 = 0;
looping j

for j=3 clear all

M(1,3) =2*1 – 3 =-1; clc
looping j
for i=1:4
for j=4 for j=1:4
M(1,4) =2*1 – 4 =-2; M(i,j)=2*i-j;
end end
end
looping i
M
for i=2
for j=1
M(2,1) =2*2 – 1 = 3;

loopings j

for j=4
M(1,4) =2*2 – 4 =-2;
end

loopings i

for i=4
for j=1
M(4,1) =2*4 – 1 = 7;

loopings j

for j=4
M(4,4) =2*4 – 4 = 4;
end

end

Ex.3) Encontre a valor de E na expressão abaixo utilizando o loop for.

10 n (n +1)
E= ∑ n!
n =0

⇒MATLAB EDITOR
 161
clear all
clc

for n=0:10
E(n+1) = n^(n+1)/factorial(n);
end

E=sum(E)

⇒RESPOSTA
E=

4.2219e+004

Se o factorial não está disponível, fazer da seguinte forma:

clear all
clc
for n=0:10
aux(n+1)=prod(1:n) % este é o factorial
E(n+1) = n^(n+1)/aux(n+1);
end

⇒TRADUÇÃO
clear all
clc

for n=0
E(0+1) = 00+1/0! = 1 clear all
looping clc

for n=1 for n=0:10

E(1+1) = 11+1/1! = 1 E(n+1) = n^(n+1)/factorial(n);
looping
end
for n=2
E(2+1) = 22+1/2! = 4 E=sum(E)
loopings

for n=10
E(10+1) = 1010+1/10! = 2.7557e+004
end

E = E(1)+E(2)+E(3)+...+E(11) = 1+1+4+...+2.7557e4 = 4.2219e+004

15.2 Contadores
Os contadores sempre são utilizados para incrementar uma determinada variável dentro de um loop
(for ou while) e geralmente possuem a seguinte estrutura:
 162

Inicia o contador
cont = 0 ⇒ se existe a necessidade de uma soma.
cont = 1 ⇒ se existe a necessidade de uma multiplicação.
cont = x ⇒ se existe a necessidade de uma subtração.
(x ≠ 0)

Inicia um loop for ou while

Incrementa o contador

cont = cont + var. à ser somada ⇒ para uma soma

cont = cont * var. à ser multiplicada ⇒ para uma multiplicação

cont = cont - var. à ser subtraída ⇒ para uma subtração

Fim do loop for ou while

4) Encontre a soma dos termos da PA definida por (1; 2,5; 4; ...; 41,5).
Utilize o loop for e um contador para realizar a soma dos termos da PA.

⇒MATLAB EDITOR ⇒RESPOSTA

clear all n=
clc
28
PA=(1:1.5:41.5); O comando length(PA)
n=length(PA) encontra a maior
dimensão do vetor PA,
ou seja encontra o n° soma =
soma=0; de termos do vetor PA.
595
for i=1:n;
soma = soma +PA(i);
end

soma

⇒TRADUÇÃO
clear all
clc
 163

PA = (1, 2.5, 4, 5.5, 7, 8.5,...41.5);

n = 28

soma=0; ⇒ inicia o contador

for i=1
soma = soma + PA(1) = 0+1 = 1; ⇒ incrementa o contador
looping para realizar a soma

for i=2 clear all

soma = soma + PA(2) = 1+2.5 = 3.5; clc
looping
PA=(1:1.5:41.5);
for i=3 n=length(PA)
soma = soma + PA(3) = 3.5+4 = 7.5;
looping soma=0;

for i=4 for i=1:n;

soma = soma + PA(4) = 7.5+5.5 = 13; soma = soma +PA(i);
loopings end

soma
for i=28
soma = soma + PA(28)= 553.5+41.5 = 595;
end

soma = 595

5) Encontre o valor de (7!). Utilize o loop for e um contador para realizar o fatorial de 7.

⇒MATLAB EDITOR
clear all
clc

n=7
prod=1;

for i=1:n;
prod = prod*i;
end

fatorial=prod

⇒RESPOSTA
n=
7

fatorial =
5040
 164

⇒TRADUÇÃO
clear all
clc

n=7

prod=1; ⇒ inicia o contador

for i=1
prod = prod * i = 1*1 = 1; ⇒ incrementa o contador
looping para realizar o produto

for i=2 clear all

prod = prod * i = 1*2 = 2; clc
looping
n=7
for i=3 prod=1;
prod = prod * i = 2*3 = 6;
looping for i=1:n;
prod = prod*i;
for i=4 end
prod = prod * i = 6*4 = 24;
fatorial=prod
loopings

for i=7
prod = prod * i = 720*7 = 5040;
end

fatorial = prod = 5040

15.3 Loops while

Os loops while permitem que um comando seja repetido um número indefinido de vezes, enquanto
uma condição lógica seja satisfeita. A forma geral é:

while expressão

comandos...

end

Os comandos entre as instruções while e end são executados enquanto todos os elementos da
expressão sejam verdadeiros.

1) Calcule o fatorial de um número qualquer utilizando o loop while.

⇒MATLAB EDITOR ⇒TRADUÇÃO
clear all clear all
clc clc

n=8
 165
n=input('Fatorial de n = '); fat=1;
fat=1;
while 8 > 1 (verdadeiro)
while n > 1 fat = fat * n = 1 * 8 = 8
fat=fat*n; n=n–1=8-1=7
n=n-1; looping
end
while 7 > 1 (verdadeiro)
fat = fat * n = 8 * 7 = 56
disp(['n! = ',num2str(fat)]) n=n–1=7-1=6
looping
⇒RESPOSTA
while 6 > 1 (verdadeiro)
Fatorial de n = 8
fat = fat * n = 56 * 6 = 336
n! = 40320 n=n–1=6-1=5

loopings

while 1 > 1 (falso)

end

n! = 40320

15.4 Estruturas if-else-end

Em diversas situações, as seqüências de comandos têm de ser executadas condicionalmente, com
base em um teste relacional. Essa lógica é implementada por meio de uma das três formas da
estrutura if-else–end.

1º Tipo if expressão

comandos...

end

Os comandos. entre as instruções if e end são executados se todos os elementos na expressão

forem verdadeiros.

1) Construa a matriz A(3,3) dada pela expressão abaixo utilizando-se das estruturas if-end:
2i + j ⇒ i ≥ j
A(i, j ) = 
2i − j ⇒ i < j

⇒MATLAB EDITOR ⇒TRADUÇÃO

clear all
clear all clc
clc for i=1
for j=1
for i=1:3 if 1 >= 1 (verdadeiro)
for j=1:3 A(1,1)=2*1+1=3
end
if i >= j if 1 < 1 (falso)
A(i,j)=2*i+j; end
 166
end
for j=2
if 1 >= 2 (falso)
if i < j end
A(i,j)=2*i-j;
end if 1 < 2 (verdadeiro)
A(1,2)=2*1-2=0
end end
end for j=3
A if 1 >= 3 (falso)
end
⇒RESPOSTA
if 1 < 3 (verdadeiro)
A= A(1,3)=2*1-3=-1
end
3 0 -1 end
5 6 1
7 8 9 for i=2
for j=1

2º Tipo
if expressão

comandos...
Com a estrutura else nunca se
else utiliza uma expressão lógica.

comandos...

end

O primeiro grupo de comandos é executado se a expressão for verdadeira; o segundo grupo de

comandos é realizado se a expressão for falsa.

2) Construa a matriz A(3,3) dada pela expressão abaixo utilizando-se das estruturas if-else-end:

2i + j ⇒ i ≥ j
A(i, j ) = 
2i − j ⇒ i < j

⇒MATLAB EDITOR ⇒TRADUÇÃO

clear all
clc for i=1
for j=1
for i=1:3
for j=1:3 if 1 >= 1 (verdadeiro)
A(1,1)=2*1+1=3
if i >= j else
A(i,j)=2*i+j; end
else
A(i,j)=2*i-j; for j=2
end
if 1 >= 2 (falso)
 167
end else
end A(1,2)=2*1-2=0
A end

⇒RESPOSTA for j=3

A= if 1 >= 2 (falso)
3 0 -1 else
5 6 1 A(1,3)=2*1-3 =-1
7 8 9 end
end

for i=2
for j=1

3) Crie um programa para verificar se três lados quaisquer formam um triângulo. Teste o programa
com: (a=3, b=4, c=5) e (a=1, b=1, c=11).
⇒MATLAB EDITOR
clear all
clc
a=input('digite o lado a =');
b=input('digite o lado b =');
c=input('digite o lado c =');
if (a+b>c)&(a+c>b)&(b+c>a)
disp('Os três lados formam um triângulo')
else
disp('Os três lados não formam um triângulo')
end

⇒RESPOSTA
digite o lado a =3 digite o lado a =1
digite o lado b =4 digite o lado b =1
digite o lado c =5 digite o lado c =11
Os três lados formam um triângulo Os três lados não formam um triângulo
3º Tipo
if expressão 1

comandos...

elseif expressão 2

comandos...

elseif expressão 3

comandos...

else
comandos...

end
 168

Somente os comandos associados à primeira expressão verdadeira são executados; as expressões

relacionais seguintes não são testadas e o resto da estrutura if-else-end é ignorada. Além disso, não
é necessário que o comando final else esteja presente.

4) Estude o sinal do discriminante (∆ = b 2 − 4ac) de uma equação do 2º grau do tipo ax2+bx+c=0

utilizando-se das estruturas if – else – end.

⇒MATLAB EDITOR (RESOLUÇÃO 1)

clear all
clc

a=input('digite o valor de a = ');

b=input('digite o valor de b = ');
c=input('digite o valor de c = ');

delta = b^2-4*a*c;

if delta >0
disp('Duas raízes reais e diferentes')

elseif delta==0
disp('Duas raízes reais e iguais')

else
disp('Duas raízes complexas')

end

⇒RESPOSTA
digite o valor de a = 2
digite o valor de b = 2 Teste para 2x2+2x-2
digite o valor de c = -2
Duas raízes reais e diferentes

digite o valor de a = 1
digite o valor de b = 4 Teste para x2+4x+4
digite o valor de c = 4
Duas raízes reais e iguais

digite o valor de a = 1
digite o valor de b = 0 Teste para x2+4
digite o valor de c = 4
Duas raízes complexas

⇒MATLAB EDITOR (RESOLUÇÃO 2)

clear all
clc

a=input('digite o valor de a = ');

b=input('digite o valor de b = ');
c=input('digite o valor de c = ');
 169

delta = b^2-4*a*c;

if delta >0
disp('Duas raízes reais e diferentes')

elseif delta==0
disp('Duas raízes reais e iguais')

elseif delta <0

disp('Duas raízes complexas')

end

15.5 Estruturas switch - case

A estrutura switch - case é utilizada quando seqüências de comandos devem ser condicionalmente
executadas, com base no uso repetido de um teste de igualdade com um argumento comum. A
estrutura tem a forma:

switch expressão

case {teste expressão 1, teste expressão 2, ...}

comandos 1...

case {...teste expressão n-1, teste expressão n}

comandos 2...

otherwise

comandos 3...

end

A expressão deve ser um escalar ou uma string de caracteres. Sendo assim, na estrutura apresentada
acima, a expressão é comparada com teste expressão 1 no primeiro comando case. Se eles são
iguais, executa-se comandos 1... e o restante dos comandos antes da declaração end são omitidos.
Se a primeira comparação não é verdadeira, a segunda é considerada e assim sucessivamente até
que seja encontrada uma igualdade. Se todas as comparações cases são falsas, o comando 3... que
segue o comando otherwise (que é opcional) é executado. Note que essa implementação da
estrutura switch-case não permite que mais de um grupo de comando seja executado.

1) Faça um programa de conversão de unidades de velocidades, de m/s para km/h, ou para ft/s, ou
para mph, sabendo que:
1m / s = 3,6km / h = 3,28084 ft / s = 2,2369mph

***Teste o programa para 45m/s ⇒mph e 45m/s⇒m/h.

 170

⇒MATLAB EDITOR
clear all
clc

mps=input('Digite o valor da velocidade em m/s = ');

unid=input('Digite a unidade final desejada = ','s');

switch unid

case{'kph','km/h','KPH'}
resp=mps*3.6;

case{'ft/s','FT/S','pps'}
resp=mps*3.28084;

case{'mph','MPH'}
resp=mps*2.2369;

otherwise
resp = 0;
end

if resp==0
disp('Unidade Desconhecida !!!')
else
disp(['Valor convertido = ',num2str(resp)])
end

⇒RESPOSTA
Digite o valor da velocidade em m/s = 45
Digite a unidade final desejada = mph
Valor convertido = 100.6605

⇒RESPOSTA
Digite o valor da velocidade em m/s = 45
Digite a unidade final desejada = m/h
Unidade Desconhecida !!!

⇒TRADUÇÃO 1 ⇒TRADUÇÃO 2
clear all clear all
clc clc
mps=45;
unid=mph; mps=45;
unid=m/h;
switch unid
switch unid
case{'kph'(F),'km/h'(F),'KPH'(F)}
case{'kph'(F),'km/h'(F),'KPH'(F)}
case{'ft/s'(F),'FT/S'(F),'pps'(F)}
 171
case{'ft/s'(F),'FT/S'(F),'pps'(F)}
case{'mph'(V),'MPH'(F)}
resp = 45*2.2369 = 100.6605; case{'mph'(F),'MPH'(F)}

otherwise otherwise
end

if resp==0 (Falso) resp = 0

else
Valor convertido = 100.6605
End end

if resp==0 (verdadeiro)
Unidade Desconhecida !!!
else
end

15.6 Comandos especiais

pause ⇒Interrompe a execução do MATLAB e espera <enter> para
reinicialização;.

pause(N) ⇒Interrompe a execução do MATLAB por N segundos.

break ⇒Interrompe a execução de qualquer controle de fluxo.

Exemplos:
1) Método da bisseção
Seja f(x) contínua em [a, b] e tal que f(a).f(b) < 0.
1) Dados iniciais:
a) intervalo inicial [a, b]
b) precisão ε
2) k = 0 (contador de iterações)
3) fa = f(a), fb = f(b)
4) Se (b - a) < ε, então escolha λ para qualquer x ∈ [a, b], fim.
a+b
5) Senão x =
2
6) fx = f(x)
7) Se fa*fx > 0, faça a = x, fa = fx,
Senão b = x, fb=fx,
8) k = k+1 e vá para o passo 4.

2) Método das cordas

Seja f(x) contínua em [a, b] e tal que f(a).f(b) < 0.
1) Dados iniciais
a) intervalo inicial [a, b]
b) precisão ε
2) k = 0 (contador de iterações)
3) fa = f(a), fb = f(b)
4) Se (b - a) < ε, então escolha λ para qualquer x ∈ [a, b], fim.
 172
a . fb − b. fa
5) x =
fb − fa
6) fx = f(x)
7) Se fa*fx > 0, faça a = x, fa = fx,
Senão b = x, fb=fx,
8) k = k+1 e vá para o passo 4.

3) Método do ponto fixo (ou da iteração linear)

Seja f(x) contínua em [a, b] e tal que f(a).f(b) < 0 e considere a equação equivalente x = F(x)
(função de iteração). Supor que a hipótese |F´(x)| < 1 para todo x ∈ [a, b] esteja satisfeita.

1) Dados iniciais
a) precisão ε
c) x (aproximação inicial da raiz)
d) Calcular Fx = F(x)
e) k = 0 (contador de iterações)
2) Se |fx| < ε então escolha λ = x, fim.
3) Senão x = F(x),
fx = f(x)
k = k+1 e volte ao passo 2

4) Método de Newton-Raphson
Seja f(x) contínua em [a, b]. Supor que as condições f(x) e f´´(x) sejam não nulas e preservem o
sinal em [a, b], f(x0)*f´´(x0) > 0, f´(x) ≠ 0 estejam satisfeitas.
1) Dados iniciais
a) precisão ε
b) x (aproximação inicial da raiz)
c) Calcular fx = f(x) e dfx = f´(x).
d) k = 0 (contador de iterações)
2) Se |fx| < ε então escolha λ = x0, fim.
fx
3) Senão x = x –
dfx
fx = f(xk)
dfx = f´(x)
k = k+1 e volte ao passo 2.

Faça um programa para cada um dos métodos anteriores (4 métodos) obtendo a raiz do exemplo que
segue, usando os comandos for, while, input etc. A seguir um programa é feito para o método da
bisseção usando while.

A capacidade calorífica (Cp) do O2 na faixa de temperatura entre 298-1500K apresenta a seguinte

equação, em função da temperatura,

0 ,4.10 5
Cp(T) = 7,16 + 10-3T -
T2

onde T está em K e Cp em cal/mol.K. Determine qual a temperatura (K) em que a capacidade

calorífica do O2 é de 8,15 cal/mol.K.
Se Cp = 8,15 a equação fica:
 173
0 ,4.10 5
f(T) = -0,99 + 10-3T - =0
T2

Solução: T = 1027,86 K

No editor de textos para o método da bisseção

clc
clear all
%% metodo da bisseção

a=input('o valor de a = ')

b=input('o valor de b = ')

fa=-0.99+1e-3*a-0.4e5/a^2;
fb=-0.99+1e-3*b-0.4e5/b^2;

k=0; %contador de iteraçoes

if(fa*fb) > 0

disp('o intervalo [a,b] não é correto, mude-o')

else
while(abs(b-a)>0.0001)
x=(a+b)/2;
fx=-0.99+1e-3*x-0.4e5/x^2;
if(fa*fx)>0
a=x;
fa=fx;
else
b=x;
fb=fx;
end
k=k+1
end
disp(['a raiz é x =',num2str(x)])
end

15.7 Exercícios propostos

1) Encontre o produto dos termos da serie definida por A(j)=2j/4 onde j=1,2,3....10 utilizando
controles de fluxo.

2) Utilizando-se de controles de fluxo monte a matriz A(i,j) de ordem 5 definida por:

i + j 2 ⇒ i < j

A(i, j ) = 2 − i ⇒ i = j

i − j ⇒ i > j
2

3) Faça um programa para determinar o fatorial de um número n qualquer utilizando o loop for para
o cálculo do fatorial, lembre-se que:
a) não existe fatorial de números negativos e números racionais;
 174
b) 0! = 1 ;
c) 1! = 1 ;
d) n ≤ 21 ;

4) No exercício anterior calcule o fatorial de n utilizando o loop while.

5) A série de Fibonacci é formada pela seguinte seqüência:

1,1,2,3,5,8,13,21,34,55,....

Mostre e calcule a soma dos 20 primeiros termos da série utilizando-se dos controles de fluxo.

6) Elabore um programa utilizando controles de fluxo para formar e calcular o valor de S nas
expressões abaixo:

1 2 3 4 5 6 10
a) S= + + + + + ... +
1 4 9 16 25 36 100

1 2 3 4 5 6 12
b) S= − + − + − ... −
1 4 9 16 25 36 124
7) Encontre a matriz literal A(8x8), sabendo que:

a ⇒ se i < 4
A(i, j) = 
b ⇒ se i ≥ 4

8) Encontre a soma de uma série de números digitados aleatoriamente utilizando controle de fluxo.

9) Encontre o maior número de uma série de números digitados aleatoriamente utilizando controle
de fluxo.

10) Separe os números ímpares dos pares de uma série de números digitados aleatoriamente
utilizando controle de fluxo.

11) Se n=1,2,3...20, encontre a soma dos 20 primeiros termos da serie abaixo. Utilize os controles
de fluxo.
Se n = 1, 3, 5... ⇒ A(n)= -2n/n2
Se n = 2, 4, 6... ⇒ A(n)= 2n/n2

12) Utilizando o loop for encontre o valor de X nas equações abaixo, sendo i=1, 2,...100

i 3 − 4i 2 − 1 (i + 1)(i − 1)
a) X = ∏ b) X = ∑
i i2 i
i

13) Faça um programa utilizando os controles de fluxo para colocar em ordem crescente três
números quaisquer.
 175
14) Faça um programa utilizando os controles de fluxo para colocar em ordem decrescente três
números quaisquer.

15) Sendo a PG dada por (1,-2,4,...32768), encontre o valor médio e a soma dos termos desta PG
utilizando controles de fluxo e contadores.

16) Faça um programa para transformação de dimensões em metros para milímetros, ou para
centímetros, ou para pés, ou para polegadas utilizando as estruturas switch – case. Sebe-se que:

1m = 1000mm = 100cm = 3,2808 ft = 39,37in

17) Dado três lados de um triângulo. Faça um programa para determinar se os três lados fornecidos
formam um triângulo. Se afirmativo, verifique se o triângulo é isósceles, escaleno ou eqüilátero e
calcule a sua respectiva área.

b
S= p(p − a )(p − b)(p − c)
a
1
p= ( a + b + c)
2
c

18) O metano apresenta a seguinte equação do calor específico em função da temperatura, na faixa
entre 298 e 1500K:

Cp = 3,381 + 18,044.10-3 T – 4,3.10-6 T2

onde T está em K e Cp em cal/mol.K. Determine a temperatura (K) para a qual a capacidade

calorífica do CH4 vale 15 cal/mol.K.

Solução:T = 794,26 K

19) Calcule o volume ocupado por 0,5 mol de CO2 a 180oC, 9 atm e R = 8,314 Pa.m3/mol.K.
(a) usando a equação dos gases ideais, PV = nRT;
(b) usando a equação de Van der Walls, expressa por:

 n 2a 
 P + 2 (V − nb ) = nRT
 V 

onde a e b são constantes específicas para cada gás, apresentadas na tabela 1:

Tabela 1 – Constantes de Van der Walls.

Gás a[Pa(m3/g-mol)2] b [m3/g-mol]
Ar 1,348x10-1 3,66x10-5
NH3 4,426x10-1 3,73x10-5
-1
CO2 3,648x10 4,28x10-5
H2 0,248x10-1 2,46x10-5
-1
CH4 2,279x10 4,28x10-5
N2 1,365x10-1 3,86x10-5
O2 1,378x10-1 3,19x10-5
-1
Vapor d’água 5,553x10 3,06x10-5
 176

É necessário mudar a temperatura para Kelvin, T(K) = T(oC) + 273,15 e a pressão para Pascal
através da relação, 1 atm = 101325 Pa.

Solução: (a) V = 0,002065679 m3 (b) V = 0,00203853 m3

20) A pressão de vapor de uma substância em função da temperatura ou a temperatura em função da

pressão de vapor pode ser dada através de várias correlações. Uma das correlações bastante simples
é a equação de Antoine, que apresenta 3 constantes empíricas (que para a água são: A = 18,3036, B
= 3816,44 e C = 4613) que relacionam uma variável (P(mmHg) ou T(K)) em função da outra (T ou
P).

B
ln( P sat ) = A −
T +C

Existem correlações mais complexas que relacionam pressão de vapor de uma substância e
temperatura. Uma destas correlações está apresentada a seguir.

B
ln(P sat ) = A + + D ln(T ) + ET F
T+C

Na primeira equação se for fornecida a pressão de saturação, pode-se obter diretamente a

temperatura. No entanto, na segunda equação, isso não é possível, sendo necessário o uso de um dos
métodos numéricos para obter zeros de funções, para encontrar a temperatura. Para a pressão em
kPa e temperatura em K, as constantes da segunda equação para a água são A = 65,9278; B = -
7,22753.103; C = 0; D = -7,17695; E = 4,031.10-6 e F = 2 e para o etano são A = 4,40101.101; B = -
2,56882.103; C = 0; D = -4,97635; E = 1,46447.10-5; F = 2 (Cutlip e Shachan, 1999). Calcule a
temperatura do etano para a qual a pressão de vapor é de 250 mmHg.
É necessário mudar a pressão para kPa através da relação, 760 mmHg = 101,325 kPa.

Solução: T = 166,0982356 K

21) Para o etano, a equação que calcula a energia livre de Gibbs em função da temperatura é
apresentada a seguir:

G = a + bT + cT2 + dT3 + eT4

Para a temperatura em K, na faixa entre 25oC e 426,85oC, a energia livre de Gibbs, em kJ/kg-mol
tem os valores das constantes a = -8,5787.104; b = 1,6858.102; c = 2,6853.10-2 e d = e = 0. Calcule a
temperatura para a qual a energia livre de Gibbs vale 1,902.103 kJ/kg-mol de etano.
É necessário mudar a temperatura para Kelvin, T(K) = T(oC) + 273,15 para obter o intervalo onde
está a raiz.

Solução: T = 483,0018 K
 177
16. INTRODUÇÃO À PROGRAMAÇÃO

Alguns destes comandos já foram utilizados anteriormente nas outras aulas. Contudo são
apresentados novamente aqui.

16.1 Comandos básicos

⇒Comando input é utilizado para entrada de valores

numéricos via teclado sempre acompanhado de (;).
x=input('texto');
Imprime o texto entre aspas e armazena o valor
digitado na variável x.
⇒Comando input é utilizado para entrada de strings
de caracteres via teclado sempre acompanhado de (;).
x=input('texto','s');
Imprime o texto entre aspas e armazena a string
digitada na variável x.
disp('texto') ⇒Imprime o texto entre aspas.
⇒Insere comentários no programa que são ignorados
%texto
pelo MATLAB, usando o símbolo de porcentagem.
⇒Imprime a variável (x) inteira ao lado do texto
disp(['texto',int2str(x)])
digitado, transfomando-a em uma string de caracteres.
⇒Imprime a variável (x) real ou complexa ao lado do
disp(['texto',num2str(x,N)]) texto digitado, transfomando-a em uma string de
caracteres com N dígitos (N máximo=17).
⇒Quando N não é informado utiliza-se o formato
disp(['texto',num2str(x)]) padrão do MATLAB de 4 dígitos após o ponto
decimal.
⇒Comando para encerrar execução de alguma
operação no MATLAB antes do fim, usado
ctrl+c
principalmente se houver algum problema de
execução.
pause(n) ⇒Pausa a execução do MATLAB por (n) segundos.
⇒Pausa a execução do MATLAB até até nova
pause
instrução via teclado seja fornecida ⇒ <enter>.
⇒Indica que a linha de comando continua na próxima
... linha.

16.2 Estrutura dos programas

Os programas realizados no MATLAB devem obedecer a seguinte estrutura:

ENTRADA DE DADOS

MANIPULAÇÃO DOS DADOS

SAÍDA DE DADOS
 178
1) Crie um programa para cálculo das raízes de uma equação do 2°grau no formato
ax 2 + bx + c = 0 , onde as raízes, (com 8 dígitos sem ponto decimal) são dadas pela equação de
Báskara:

− b ± b 2 − 4ac
x1 , x 2 =
2a

Estrutura do programa:

a
Entrada de dados ⇒ b
c

− b − b 2 − 4ac
Cálculo da 1ª raiz⇒ x1 =
2a
Manipulação ⇒
− b + b 2 − 4ac
Cálculo da 2ª raiz⇒ x 2 =
2a

Saída de dados ⇒ Raízes x1 e x2 com 8 dígitos

⇒MATLAB EDITOR (SOLUÇÃO 1)

%Programa para resolução de equações do 2° Grau
%Salvar o programa com o nome: raizes1
clear all
clc

disp('Programa para resolução de equações do 2° Grau')

disp(' no formato ax^2+bx+c=0')
disp(' ')

%--------------Entrada----------------

a = input('Entre com o valor de a= ');

b = input('Entre com o valor de b= ');
c = input('Entre com o valor de c= ');
disp(' ')

%------------Manipulação--------------

delta=b^2-4*a*c;
x1=(-b-sqrt(delta))/(2*a);
x2=(-b+sqrt(delta))/(2*a);

%---------------Saída-----------------

disp('Raízes da equação')
disp(['x1= ', num2str(x1,8)])
disp(['x2= ', num2str(x2,8)])
 179

⇒RESPOSTA
Programa para resolução de equações do 2° Grau
no formato ax^2+bx+c=0

Entre com o valor de a= 2

Entre com o valor de b= 2
Entre com o valor de c= -8

Raízes da equação
x1= -2.5615528
x2= 1.5615528

⇒MATLAB EDITOR (SOLUÇÃO 2)

%Programa para resolução de equações do 2° Grau
%Salvar o programa com o nome: raizes2
clear all
clc

disp('Programa para resolução de equações do 2° Grau')

disp(' no formato ax^2+bx+c=0')
disp(' ')

%--------------Entrada----------------

a = input('Entre com o valor de a= ');

b = input('Entre com o valor de b= ');
c = input('Entre com o valor de c= ');
disp(' ')

%------------Manipulação--------------

delta=b^2-4*a*c;
x1=(-b-sqrt(delta))/(2*a);
x2=(-b+sqrt(delta))/(2*a);

%---------------Saída-----------------

disp('Raízes da equação')
disp(['x1= ', num2str(x1,8)])
disp(['x2= ', num2str(x2,8)])

%------Rotina de Reinicialização------

opc1=5;
while (opc1~=2)
disp(' ')
disp('---------------Escolha uma opção----------------')
disp('Digite 1 para Reiniciar programa')
disp('Digite 2 para Finalizar programa')
 180
opc1 = input('Digite a opção escolhida= ');

if opc1==1
raizes2 %nome do arquivo
elseif opc==2
clc
end
end

⇒RESPOSTA
Programa para resolução de equações do 2° Grau
no formato ax^2+bx+c=0

Entre com o valor de a= 2

Entre com o valor de b= 2
Entre com o valor de c= -8

Raízes da equação
x1= -2.5615528
x2= 1.5615528

---------------Escolha uma opção---------------

Digite 1 para Reiniciar programa
Digite 2 para Finalizar programa
Digite a opção escolhida= 2

2) Crie um programa para transformação de unidades conforme o esquema abaixo:

1m ⇔ 1000mm ⇔ 39,37in (polegadas)

onde o usuário deve fornecer o valor numérico a ser convertido e sua unidade.

Estrutura do programa:

Valor numérico a ser convertido

Entrada de dados ⇒
Unidade do valor numérico

Cálculo do valor em [m]

Manipulação ⇒ Cálculo do valor em [mm]
Cálculo do valor em [in]

Valor em [m]
Saída de dados ⇒ Valor em [mm]
Valor em [in]

⇒ MATLAB EDITOR
 181
%Programa para transformação de unidades m <=> mm <=> in
%Salvar o programa com o nome: transf
clear all
clc

disp(' ')
disp(' Programa para transformação de unidades. ');
disp(' ')
disp(' Transforma m <=> mm <=> in ')
disp(' ')

%-----------------Entrada-------------------

a=input('Insira o valor a ser convertido = ');

b=input('Insira a unidade inicial= ', 's');
disp(' ')

%----------------Manipulação-----------------

switch b

case{'m', 'metros'}
m=a;
mm=a*1000;
in=a*39.37;
case{'mm', 'milimetros'}
m=a/1000;
mm=a;
in=a/1000*39.37;
case{'in', 'polegadas'}
m=a/39.37;
mm=a/39.37*1000;
in=a;
otherwise
disp('***Unidade desconhecida***')
m='???';
mm='???';
in='???';
end

%-----------------Saída----------------------

disp(['Metros [m] = ',num2str(m,17)])

disp(['Milimetros [mm] = ',num2str(mm,17)])
disp(['Polegadas [in] = ',num2str(in,17)])
 182
%-------Rotina de Reinicialização------------
opc1=5;
while (opc1~=2)
disp(' ')
disp('---------------Escolha uma opção---------------')
disp('Digite 1 para Reiniciar programa')
disp('Digite 2 para Finalizar programa')
opc1 = input('Digite a opção escolhida= ');

if opc1==1
transf %nome do arquivo mãe
elseif opc==2
clc
end
end

⇒RESPOSTA 1

Programa para transformação de unidades.

Transforma m <=> mm <=> in

Insira o valor a ser convertido = 10000

Insira a unidade inicial = mm

Metros [m] = 10
Milímetros [mm] =10000
Polegadas [in] = 393.69999999999999

---------------Escolha uma opção---------------

Digite 1 para Reiniciar programa

Digite 2 para Finalizar programa

Digite a opção escolhida=2

»

⇒RESPOSTA 2

Programa para transformação de unidades.

Transforma m <=> mm <=> in

Insira o valor a ser convertido = 10000

Insira a unidade inicial= pol

***Unidade desconhecida***
Metros [m] = ???
Milimetros [mm] = ???
Polegadas [in] = ???
 183

---------------Escolha uma opção---------------

Digite 1 para Reiniciar programa
Digite 2 para Finalizar programa
Digite a opção escolhida= 2

A SEÇÃO A SEGUIR É IMPORTANTE

16.3 Criando subrotinas - function

A) Criando uma function
Pode-se criar funções específicas (subrotinas) de seu interesse no Matlab. Tais funções são
seqüências de comandos que aceitam vários parâmetros de entrada e retornam vários parâmetros de
saída. Possuem forma geral:

function [Ps1, Ps2,...] = nome_funcao(Pe1, Pe2,...)

comandos...

Tradução:
Ps1, Ps2, ... ⇒Parâmetros de saída;

nome_função ⇒Nome da função;

Pe1, Pe2, ... ⇒Parâmetros de entrada;

comandos... ⇒Comandos que devem ser executados pela

função.

Observação:
A function deve ser salva no mesmo diretório onde está o arquivo mãe e com o mesmo
nome_funcao.

B) Acessando uma function

Para chamar uma function dentro de um programa, deve-se usar a mesma estrutura utilizada em sua
criação, ou seja:

[Ps1, Ps2,...] = nome_funcao(Pe1, Pe2,...)

1) Crie uma function para o cálculo de raízes do 2° grau e implemente-a no programa 1 – solução 2.

⇒MATLAB EDITOR
%function raiz2grau
%Retorna as raízes x1 e x2 da equação
%do segundo grau do tipo ax^2+bx+c = 0
 184

function[x1,x2]=raiz2grau(a,b,c)
Observação:
disp('Programa para resolução de equações do 2°Sempre
Grau') salvar arquivo
com o mesmo nome da
disp(' no formato ax^2+bx+c=0') function, isto é, para
disp(' ') este exemplo: raiz2grau.

delta=b^2-4*a*c;
x1=(-b-delta)/(2*a);
x2=(-b+delta)/(2*a);

2) Construa uma function utilizando os controles de fluxo para colocar em ordem decrescente 3
números quaisquer.

3) Construa uma function utilizando os controles de fluxo para colocar em ordem crescente n
números quaisquer.

4) Construa uma function para o método das cordas passando a,b e precisão. Siga o algoritmo que
segue.

Dados iniciais:
f) intervalo inicial [a, b]
g) precisão ε
h) k = 0 (contador de iterações)
i) fa = f(a), fb = f(b)
2) Se f(a)×f(b) > 0, Imprimir “A raiz não pertence a este intervalo”.
Senão
3) Enquanto (b - a) > ε
a × fb − b × fa
x=
fb − fa
fx = f(x)
Se fa×fx > 0 então
a = x, fa = fx
Senão b = x, fb=fx
k = k+1
Imprimir o valor de x após a convergência

16.4 Exercícios propostos

OBS: Quando possível use function, observe que dentro do código não pode haver clc e clear all
depois de chamar a function, caso contrário as variáveis serão limpas.

1) Elabore um programa para fornecer o peso (massa) ideal de uma pessoa. Este programa deve ter
como dados de entrada a altura e o sexo de da pessoa utilizando as seguintes fórmulas:

Para homens ⇒ P = 72,7 h − 58

Para mulheres ⇒ P = 62,1h − 44,7

Onde h é a altura da pessoa.

 185
2) Elabore um programa para verificar o desconto do imposto de renda e do recolhimento do INSS
incidente sobre o salário de uma pessoa. Sabe-se que:
Tabela 1 - IRPF
Base de cálculo (R$) Alíquota % Parcela a deduzir (R$)
Até 1.058,00 - -
De 1.058,01 até 2.115,00 15,0 158,70
Acima de 2.115,00 27,5 423,08

Tabela 2 - INSS
Alíquota para fins de recolhimento ao INSS
Salário-de-contribuição (R$)
(%)

até R$ 720,00 7,65 %

de R$ 720,01 até R$ 1.200,00 9,00 %

de R$ 1.200,01 até R$ 2.400,00 11,00 %

*** Dados de janeiro de 2004

O usuário deve fornecer a valor do salário bruto e o programa deve fornecer o valor do salário
líquido e os valores do desconto do INSS e do IRPF.

Exemplos:
1)Salário bruto: R$ 700,00
⇒INSS = 7,65% de R$700,00 = R$53,55
⇒IRPF = Salário bruto – INSS = R$646,45
IRPF ⇒ isento (tab.1)
Salário líquido: R$ 700,00 – (IRPF+INSS) = R$646,45

2)Salário bruto: R$ 1200,00

⇒INSS = 9,00% de R$1200,00 = R$108,00
⇒IRPF = Salário bruto – INSS = R$1092,00
IRPF ⇒ 15% de R$1092,00 – R$158,70 = R$5,10
Salário líquido: R$ 1200,00 – (IRPF+INSS) = R$1086,90

3) Salário bruto: R$ 2800,00

⇒INSS = 11,00% de R$2400,00 = R$264,00
⇒IRPF = Salário bruto – INSS = R$2536,00
IRPF ⇒ 27,5% de R$2536,00 – R$423,08 = R$274,32
Salário líquido: R$ 2800,00 – (IRPF+INSS) = R$2261,68

3) Elabore uma function para o cálculo do desconto de INSS.

4) Elabore uma function para o cálculo do desconto do IRPF.

5) Elabore uma function para a rotina de reinicialização.

6) Faça um programa para encontrar raízes de uma equação do terceiro grau do tipo
ax 3 + bx 2 + cx + d =0. Sabendo que as raízes são x1, x2 e x3 , dadas pelas expressões abaixo, se e
somente se a=1:

3c − b 2 9bc − 27d − 2b3

Q= R=
9 54
 186
D = Q3 + R 2 M=R+ D N=R− D

Se M e N ∈ R

Se M ≥ 0 ⇒ S = 3 M Se N ≥ 0 ⇒ T = 3 N

Se M < 0 ⇒ S = −3 M Se N < 0 ⇒ T = − 3 N

Se M e N ∉ R

S=3 M

T=3N

1
x1 = S + T − b
3
1 1 1
x2 = − ( S + T ) − b + i 3 ( S − T )
2 3 2
1 1 1
x3 = − ( S + T ) − b − i 3 ( S − T )
2 3 2

**Mais detalhes: http://mathworld.wolfram.com/CubicFormula.html

EXEMPLOS:
 x1 = 0,7100

a) x 3 + 3x 2 + 3x − 4  x 2 =-1,8550 + 1,4809i
 x =-1,8550 - 1,4809i
 3

 x1 = −12,5208

b) x + 11x − 20x − 12  x 2 =2,0000 +0,0000i
3 2
 x =-0,4792 - 0,0000i
 3
7) Elabore uma function para encontrar as raízes de uma equação do 3º grau do tipo
ax 3 + bx 2 + cx + d = 0 , do exercício 1.

8) Faça um programa para implementação da function acima. O usuário deve fornecer os valores de
a, b, c e d.

9) Elabore um programa para conversão de temperaturas de °C(graus Celsius) para °F(graus

Fahrenheit) e K (Kelvin), sabendo que:

C F − 32 K − 273
= =
5 9 5
 187

Obs.: O usuário deve fornecer o valor da temperatura inicial, a unidade da temperatura inicial e a
temperatura final desejada;

10) Sendo A, a matriz principal dos coeficientes, B a matriz das respostas, D o determinante da
matriz principal e (Dx, Dy, Dz) são os determinantes das matrizes secundárias (Ax, Ay, Az),
respectivamente. Faça um programa para resolver um sistema linear de 3 equações e 3 incógnitas
onde o usuário deve fornecer somente a matriz A e B, considerando que:

⇒Sistema é SPD se D ≠ 0;
⇒Se D=0 e Dx ou Dy ou Dz é igual a zero, o sistema é SPI;
⇒Se D=0 e Dx, Dy e Dz são diferentes de zero, o sistema é SI.

Caso o sistema SPD encontre o valor de cada incógnita.

Exemplo:
 3x + y − z = 0

SPD  2 x + 3y + z = 1
x + 2 y − 2z = −5

3 1 − 1  0
a) Entradas

⇒ A= 2 3 1 e B= 1
  
   
1 2 − 2 − 5
b) Montar as matrizes secundárias trocando-se os coeficientes de x, y e z pela matriz resposta para
encontrar as respectivas matrizes secundárias.

 0 1 − 1 3 0 − 1  3 1 0 
Ax =  1 3 1  , Ay = 2 1 1  e Az = 2 3 1 
 − 5 2 − 2 1 − 5 − 2 1 2 − 5

c) Calcular D = det(A), Dx = det(Ax), Dy = det(Ay) e Dz = det(Az)

d) Verificar se o sistema é SPD, ou SPI, ou SI. Caso seja SPD encontre os valores das incógnitas
através de:
Dx Dy Dz
x= , y= e z =
D D D

e) O sistema acima é SPD com respostas são x=1, y=-1 e z=2

***Teste também para os sistemas abaixo.

 3x + 5 y + 5z = 1  − 3x + 5 y + 5z = 1
 
SI  x − 4 y − 4z = 0 SPI  6 x − 10 y − 10z = −2
5x + 2 y + 2z = 2 − 12 x + 20 y + 20z = 4
 

11) Elabore um programa para realizar a resolução de um sistema de equações de no máximo 5

equações e 5 incógnitas conforme diretrizes do exercício 10.
 188
12) Faça um programa para encontrar as raízes de uma equação do 4º grau do tipo:
ax 4 + bx 3 + cx 2 + dx + e . Desde que a=1, as raízes são dadas por:
1 1 1
x1 = − b + R + D
4 2 2

1 1 1
x2 = − b + R − D
4 2 2

1 1 1
x3 = − b − R + E
4 2 2

1 1 1
x4 = − b − R − E
4 2 2
Onde

1 2
R= b −c+F
4

b − R 2 − 2c + (4bc − 8d − b3 )R −1
3 2 1
Se R ≠ 0 ⇒
4 4
D=
3 2
Se R = 0 ⇒ b − 2c + 2 F 2 − 4e
4

b − R 2 − 2c − (4bc − 8d − b 3 )R −1
3 2 1
Se R ≠ 0 ⇒
4 4
E=
3 2
Se R = 0 ⇒ b − 2c − 2 F 2 − 4e
4

F é a raiz real da equação de 3° grau abaixo:

m = −c
y 3 + my 2 + ny + o = 0 n = d ⋅ b − 4e
o = 4c ⋅ e − d 2 − b 2 e

**Maisdetalhes: http://mathworld.wolfram.com/QuarticEquation.html

Exemplos:
 189
 x1 = −3.6688
 x = 0,9142 + 1,6057i

x + 2x − 3x + 12x + 2  2
4 3 2
 x 3 = 0,9142 − 1,6057i
 x 4 = −0,1597

 x1 = 0
 x = −2
 2
−12x − 46x − 52x − 16x 
4 3 2
 x 3 = −1,3333
 x 4 = −0,5000

13) No programa acima utilize uma function para o cálculo de F.

14) Elabore um programa para venda de automóveis onde o usuário têm três opções de venda:
Opção 1 ⇒ Venda à vista.
a) Com 3,6% de desconto se o automóvel custar menos que R$25000,00.
b) Com 3,2% de desconto para os demais automóveis.
Opção 2 ⇒ pode vender um automóvel com 50% de entrada e o restante em 24 vezes com juros de
1,89%
Opção 3 ⇒ pode vender um automóvel com 20% de entrada e o restante em 36 vezes com juros de
2,29%

 J (1 + J ) M 
P = F 
M
 (1 + J ) − 1 

onde:
F é o montante a ser financiado;
J é o índice percentual de juros a ser aplicado no financiamento divido por 100;
M é o número de meses do financiamento;
P é a parcela mensal a ser paga.

a) Utilize como variáveis de entrada para seu programa:

PT ⇒ preço de tabela do automóvel
opcao ⇒ opção de venda (1 ou 2 ou 3)

b) Utilize como variáveis de saída para seu programa:

Considerando uma venda à vista:

PV ⇒ preço de venda a vista

Considerando uma venda à prazo:

E ⇒ o valor da entrada necessário para financiamento
F ⇒ o montante a ser financiado.
P ⇒ o valor da parcela mensal a ser paga.
PT ⇒ valor total do veículo (entrada + valor financiado)
 190
c) Teste seu programa com as opções abaixo:

1)Venda à vista com preço de tabela menor que R$25000,00

PT = R$20000,00
opcao = 1 (venda à vista)
PV = R$19280,00

2) Venda à vista com preço de tabela maior que R$25000,00

PT = R$32000,00
opcao = 2 (venda à vista)
PV = R$30976,00

3)Venda à prazo em 24 vezes

PT = R$34000,00
opcao = 1(venda à prazo em 24 vezes)
E = R$17000,00
F = R$17000,00
P = R$ 887,65
PT = R$ 38303,54

4)Venda à prazo em 36 vezes

PT = R$34000,00
opcao = 3(venda à prazo em 36 vezes)
E = R$6800,00
F = R$27200.00
P = R$1117,46
PT = R$ 47028,63

15) Elabore um programa de conversão de coordenadas polares para retangulares e vise-versa.

16) Elabore um programa para jogar o jogo da velha com o computador.

17) Dado três lados de um triângulo. Faça um programa para determinar se os três lados fornecidos
formam um triângulo. Se afirmativo, verifique se o triângulo é isósceles, escaleno ou eqüilátero e
calcule a sua respectiva área.

b
S= p(p − a )(p − b)(p − c)
a
1
p= ( a + b + c)
2
c
 191
17. FUNÇÃO FPRINTF

Esta função exibe um ou mais valores juntamente com o texto relacionado e permite ao
programador controlar a maneira como os valores são exibidos. A forma geral dessa função é,

fprintf(format, dados)

onde format é uma cadeia de caracteres descrevendo a maneira como os dados devem ser exibidos e
dados é comporto por um ou mais escalares ou matrizes para exibição.

A cadeia de caracteres format contém texto a ser apresentado, mais caracteres especiais
descrevendo o formato dos dados. Por exemplo, a função

fprintf (`O valor de pi = %f \n`, pi)

>> O valor de pi = 3.141593

fprintf(`O valor de pi = %6.2f \n`, pi)

>> O valor de pi = 3.14

Observação: A função fprintf tem uma limitação bastante significativa: ela somente exibe a porção
real de um valor complexo. Essa limitação pode levar a resultados enganosos quando os cálculos
produzem respostas complexas. Nesses casos é melhor usar a função disp.

% f são denominados caracteres de conversão

\n são caracteres de escape

Alguns são apresentados na Tabela que segue

Tabela – Caracteres de formatação para a função fprintf

Caracteres de formatação Resultados
%d Exibe valor como inteiro
%e Exibe valor em formato exponencial
%f Exibe valor em formato de ponto flutuante
%g Exibe valor em formato de ponto flutuante ou exponencial – o que
for mais curto
\n Muda de linha

17.1 Exercícios Propostos

1) O que fazem os seguintes conjuntos de declarações? Qual a saída de cada um deles?

a) raio = input(`Enter circle radius: \n`)
area = pi*raio^2;
str = [`A área é ` num2str(area)];
disp(str);

b) valor = int2str(pi);
disp([`O valor é ` valor `!`]);
 192
c) valor = 123.4567e2;
fprintf(`valor = %e \n`, valor);
fprintf(`valor = %f \n`, valor);
fprintf(`valor = %g \n`, valor);
fprintf(`valor = %2.4f \n`, valor);

2) Se uma bola estacionária é lançada da altura h0 acima da superfície da terra, com velocidade
vertical v0, a posição e a velocidade da bola como função do tempo serão dadas pelas equações

h(t) = gt2/2 + v0t + h0

v(t) = gt + v0

onde g é a aceleração da gravidade (-9,81 m/s2), h é a altura acima da superfície da terra (assumindo
ausência de atrito do ar) e v é a componente vertical da velocidade.
a) Escreva um programa que solicite ao usuário a altura inicial da bola em metros e a
velocidade da bola em m/s.
b) Desenhe a altura e a velocidade como função do tempo. Não deixe de incluir as legendas e
título dos eixos apropriados no desenho

3) Use uma function para calcular a distância entre quaisquer dois pontos (x1, y1) e (x2,y2)
especificados pelo usuário. Use-o para calcular a distância entre os pontos (2, 3) e (8,-5).

4) A força requerida para comprimir uma mola linear é dada pela equação

F = kx

onde F é a força em N e k é a constante da mola em N/m. A energia potencial armazenada na mola

comprimida é dada pela equação

E = kx2/2

onde E é a energia em J. A informação a seguir é a respeito de quatro molas,

Mola 1 Mola 2 Mola 3 Mola 4

Força (N) 20 24 22 20
k (N/m) 500 600 700 800

Determine a compressão de cada mola e a energia potencial armazenada em cada uma delas. Qual
mola tem mais energia potencial armazenada?

5) Escreva um programa para ler uma nota de aluno (numérica) e associar uma letra a ela, segundo
a tabela a seguir:

nota > 95 A
95 ≥ nota > 86 B
86 ≥ nota > 76 C
76 ≥ nota > 66 D
66 ≥ nota > 0 F
 193
6) (a) Avalie a função f(x,y) para quaisquer dois valores especificados pelo usuário para x e y, com
a função f(x,y) definida como segue. Imprima o valor obtido para a função na tela do computador.

3 x + y, x ≥ 0, y ≥ 0

 x + y , x ≥ 0, y < 0
2

f ( x, y ) =  2
 x + 2 y, x < 0, y ≥ 0
 x 2 + y 2 , x < 0, y < 0


(b) No mesmo programa avalie a função f(x,y) para -5 ≤ x ≤ 5 e -5 ≤ y ≤ 5 e faça o gráfico da função
identificando todos os elementos do Matlab já vistos a respeito de gráficos.
 194
EXERCÍCIOS DE REVISÃO

1. Primeiro resolva manualmente depois resolva no computador. Qual a diferença entre os

resultados?
Use os valores N = 2, 3, 10, 16 para o algoritmo abaixo:
X = 1/N
Para i = 1 até 100 faça
X = (N+1)X – 1
Fim para
Escreva X

2. Dado um vetor formado por 10 números positivos quaisquer. Identifique quantos números são
múltiplos de 3, 5 e 7. Separe estes elementos em conjuntos distintos assim como os elementos que
não pertencem a nenhum dos conjuntos anteriores.

3. A seguir apresenta-se a densidade (kg/m3) ) do óleo de motor (sem uso) em função da

temperatura (K). Leia estes dados a partir de um arquivo. Faça o gráfico e obtenha um polinômio
com grau adequado para ajustar-se a estes dados. Faça o gráfico novamente usando os dados
ajustados. Calcule a média e desvio padrão da viscosidade da glicerina. Obtenha o valor máximo e
mínimo do conjunto de dados.

T(K) 273 280 290 300 310 320

ρ (kg/m3) 899,1 895,3 890 884,1 877,9 871,8

4. Crie uma matriz A de ordem 5 cujos elementos da primeira linha sejam unitários, os elementos
da segunda linha sejam múltiplos de 2, os elementos da terceira linha sejam múltiplos de 3 e os
demais elementos sejam 5. Obtenha a soma, média e desvio padrão dos elementos da matriz A.
Calcule o determinante da matriz. Por que tal determinante é nulo? Explique.

5. Considere a progressão geométrica 1, 2, 4, 8, 16, 32, ... e um inteiro positivo n. Deseja-se:

a) Imprimir os n primeiros termos;
b) Calcular e imprimir a soma dos n primeiros termos da PG sem utilizar a fórmula da soma

6. A seqüência de cálculos abaixo deve ser feita manualmente (usando frações, não divida-as) e
posteriormente no computador. Existe diferença entre as duas formas de solução?
H=½
X = 2/3 – H
Y = 3/5 – H
E = (X+X+X) – H
F = (Y+Y+Y+Y+Y) – H
G = E/F

7. Elabore um programa de conversão de unidades.