Universidade Estadual Paulista

FCT – Faculdade de Ciências e Tecnologia

Departamento de Cartografia

Tutorial do GNU Octave / GNU Octave Tutorial

Autores:
Mauricio Galo
Paulo de Oliveira Camargo

Presidente Prudente - SP, 2023

1
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Tutorial do GNU Octave / GNU Octave Tutorial

Tutorial do GNU Octave / GNU Octave Tutorial

1. INTRODUÇÃO
2. INICIANDO O GNU Octave
3. OPERAÇÕES COM MATRIZES E VETORES
4. GRÁFICOS
5. ARQUIVOS SCRIPT
6. CONTROLE DE FLUXO
7. ABERTURA E VISUALIZAÇÃO DE IMAGENS
8. IMPORTAÇÃO / EXPORTAÇÃO DE DADOS

Autores:
9. CRIAÇÃO E USO DE FUNÇÕES
Mauricio Galo
2023 Paulo de Oliveira Camargo
Departamento de Cartografia / FCT - UNESP

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 2

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Este material foi preparado a partir das seguintes referências (em ordem alfabética): 1. INTRODUÇÃO

CORAL, A. M.; SANTOS, M. P.; BASTOS, T. D. A.; BORBA, M. Curso de Matlab. Universidade Federal de
Santa Catarina, Dep. de Eng. de Produção, Programa Especial de Treinamento – PET. Florianópolis – SC, 1999. 1.1. O que é GNU Octave
29p.
EATON, J. W.; BATEMAN, D.; HAUBERG, S. GNU Octave – Edition 3 for Octave version 3.2.3, July, 2007.
672p.
EATON, J. W.; BATEMAN, D.; HAUBERG, S.; WEHBRING, R. GNU Octave – Edition 4 for Octave version O GNU Octave é um aplicativo que foi originalmente desenvolvido com o propósito
4.0.0 – Free Your Numbers, March, 2015. 966p.
didático, mais especificamente para o projeto de reatores químicos e surgiu a partir da
HANSELMAN, D.; LITTLEFIELD, B. Matlab - Versão do Estudante: guia do usuário (Tradução). São Paulo:
Makron Books, 1997. 305p. intenção de criar um aplicativo no qual a programação fosse mais rápida do que nas
McANDREW, A. Introduction to Digital Image Processing with MATLAB®. Thomson Course Technology,
2004. 509p. ISBN: 0-534-40011-6
demais linguagens.
PAGAMISSE, A.; SOUZA, L. H. G. Introdução ao Software Octave. Semana de Cursos de Matemática,
Estatística e Computação, 25-29 de agosto de 2003, FCT/UNESP, Presidente Prudente, 2003. 54p.
Segundo os autores, os alunos poderiam dedicar mais tempo na análise e solução dos
SIGMON, K. MATLAB Primer – Third Edition. Department of Mathematics, University of Florida, Gainesville,
1993. 35p. Disponível em [Link] Acesso em Março/2021. problemas de química, do que especificamente na implementação.
ZERI, L. M. M. Apostila de Matlab. Instituto Nacional de Pesquisas Espaciais – INPE, 2001. 19p.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 3 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 4
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

O seu desenvolvimento começou por volta de 1988 e um de seus idealizadores foi Características básicas:
John W. Eaton, além de James B. Rawlings da Universidade de Wisconsin-Madison e
John G. Ekerdt da Universidade do Texas. Atualmente o desenvolvimento é feito por - Domínio público;
vários profissionais, de diferentes países, sendo a sua distribuição feita de acordo com a - Possui vários comandos que são semelhantes ao MATLAB® e SciLab®;
licença GLP (GNU General Public License). - É um software multiplataforma uma vez são disponíveis versões para diferentes
sistemas operacionais: Linux, Unix, Windows, Mac, etc.
Onde encontrar (páginas oficiais):
Últimas versões (estáveis):
[Link]
[Link]
6.4.0 (Lançada em 30/Outubro/2021)
7.3.0 (Lançada em 03/Novembro/2022)
8.1.0 (Lançada em 03/Março/2023)
8.3.0 (Lançada em 09/Agosto/2023)

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 5 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 6

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

1.2. Ambiente de trabalho e comandos básicos do sistema

Desde o inicio do desenvolvimento do GNU Octave, em 1988, foram lançadas

inúmeras versões, para as diferentes plataformas1. Na maior parte das versões a única
possibilidade de trabalho era via uma interface não gráfica. Há poucos anos a opção de
utilização de uma interface gráfica foi possível, tornando mais prática e amigável sua
atualização. Na figura abaixo são mostradas duas interfaces (versões 3.2.4 e 4.0.0, resp.).

Interface não gráfica – GNU Octave versão 3.2.4.

1
Consulte a página “[Link] para ver algumas destas versões.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 7 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 8
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Interface gráfica do GNU Octave versão 6.4.0, lançada em outubro de 2021. Interface gráfica do GNU Octave versão 4.2.2, lançada em março de 2018.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 9 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 10

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Interface gráfica do GNU Octave versão 4.4.0, lançada em abril de 2018. Interface gráfica do GNU Octave versão 5.1.0, lançada em março de 2019.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 11 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 12
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Interface gráfica do GNU Octave versão 5.2.0, lançada em janeiro de 2020. Interface gráfica do GNU Octave versão 6.2.0, lançada em fevereiro de 2021.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 13 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 14

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Diretório de
trabalho
Janela de
comandos Editor

Interface gráfica do GNU Octave versão 8.3.0, lançada em agosto de 2023. Interface gráfica do GNU Octave versão 6.2.0.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 15 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 16
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Além dessas versões, que podem ser instaladas e executadas no modo off-line, também Alguns comandos do GNU Octave são básicos e semelhantes a outros usados em
é possível executar o GNU Octave de modo online (ver link abaixo). sistemas operacionais como DOS, Unix, Linux, dentre outros. A tabela seguinte mostra
alguns destes comandos, que podem se acionados a partir do prompt do GNU Octave.

cd Troca o diretório de trabalho atual

dir Lista o conteúdo do diretório atual
pwd Mostra o path do diretório atual de trabalho
delete Exclui arquivo
Exemplo de uso: delete [Link]
type Mostra o conteúdo de arquivos texto
Exemplo de uso: type [Link]
what Lista arquivos específicos do GNU Octave no diretório
corrente
help Ajuda do GNU Octave
Exemplo de uso: help comando
quit ou exit Fechar o aplicativo e sair do sistema
Interface gráfica do GNU Octave executado no modo online.
Fonte: [Link] Alguns outros comandos: mkdir, rmdir, etc.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 17 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 18

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

1.3. Pacotes / Bibliotecas adicionais do GNU Octave

Lista de pacotes/bibliotecas disponíveis para o usuário ao instalar o aplicativo
GNU Octave versão 8.3.0.
O GNU Octave é um software em constante desenvolvimento e após a instalação do
audio instrument-control queueing
aplicativo, outros pacotes adicionais podem ser instalados. Para ter uma lista dos pacotes communications interval signal
disponíveis para o GNU Octave sugere-se a consulta à seguinte página control io sockets
data-smoothing linear-algebra sparsersb
[Link] database lssa specfun
dataframe ltfat splines
dicom mapping statistics
Antes de instalar qualquer pacote adicional, sugere-se ativar o comando “pkg list” no financial miscellaneous stk
fits nan strings
prompt do GNU Octave para ver todos os pacotes que são pré-instalados junto com a fuzzy-logic-toolkit netcdf struct
ga nurbs symbolic
instalação do aplicativo. general ocs tisean
generate_html odepkg tsa
geometry optim video
gsl optiminterp windows
image quaternion zeromq

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 19 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 20
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Para a instalação de algum pacote adicional, de nome package_file_name.[Link], 2. INICIANDO O GNU Octave
basta ativar o seguinte comando no prompt do GNU Octave:
2.1. Variáveis

O GNU Octave trabalha essencialmente com um tipo de variável: matriz, que pode
pkg install package_file_name.[Link]
conter números (complexos ou não) e textos. Em alguns casos, um tratamento especial é
dado a uma matriz 1 x 1 (escalar) ou a matrizes 1 x n ou n x 1 (vetores).

2.1.1. Entrando com valores

No GNU Octave não é necessário declarar as variáveis e os respectivos tipos (inteiro,

char, double, etc,) para iniciá-las, como é feito em outras linguagens de programação
(C/C++, dentre outras). Ao atribuir valores numéricos (ou alfanuméricos) a uma variável,
o programa aloca a memória automaticamente.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 21 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 22

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

A maneira mais fácil de entrar com uma pequena quantidade de valores é digitando ou
diretamente os dados: » A=[1, 2, 3; 4, 5, 6; 7, 8, 9]
- envolva os elementos com colchetes, [ ]; (Para matrizes e vetor)
- separe cada elemento com espaços ou vírgulas; Resultado: A =
- use ponto-e-vírgula (;) para indicar fim da linha. 1 2 3
4 5 6
Por exemplo, para entrar com a matriz abaixo na memória do computador, e guardá-la
7 8 9
na variável A:

1 2 3
A  4 5 6 OBS: Para que o computador realize a operação e não mostre a saída, basta terminar a

7 8 9 expressão com ponto-e-vírgula (;). Isto é muito útil para evitar que o computador mostre
resultados de cálculos intermediários desnecessários, e para acelerar as operações.
Basta digitar:
» A=[1 2 3; 4 5 6; 7 8 9]

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 23 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 24
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

2.1.2. Variáveis reservadas 2.2. Expressões e comandos básicos

Existem algumas variáveis que são intrínsecas ao GNU Octave, como por exemplo: O GNU Octave, assim como o MATLAB, são linguagens de alto nível e interpretadas.
Nestas linguagens o sistema avalia as expressões digitadas, que são geralmente na forma:
help Ajuda do GNU Octave
ans Nome de variável padrão usado para resultados variável = expressão
eps Precisão do número real (ponto flutuante)
nan Not a Number (indeterminação: 0/0) Os números reais são escritos em notação decimal e para criar números complexos
realmax Maior número real positivo utilizável basta escrever i (ou j) depois da parte imaginária. Alguns exemplos de números permitidos
realmin Menor real positivo utilizável são mostrados abaixo:
inf Infinito: 1/0
computer Tipo de computador 1/3 -99 .0001
pi 3,14159265358979 9.63973 1.602E-20 6.025E23
i, j Unidade imaginária (i2=-1)
3 + 2i -3.1459i 3E5j
version Versão do GNU Octave

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 25 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 26

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Todo cálculo é realizado com todas as casas decimais (eps), embora os números O formato + é uma maneira compacta de mostrar matrizes de grandes dimensões. Os
mostrados ocultem algumas dessas casas. Para mudar o formato de saída dos números, símbolos "+", "-", e "espaço em branco" são mostrados, respectivamente para elementos
pode-se usar o comando format. positivos, elementos negativos e zeros.
O formato “default”, chamado de formato short, mostra aproximadamente 5 dígitos
significativos ou usam notação científica. Exemplo: Podem-se construir expressões com os operadores aritméticos usuais:
+ adição /e\ Divisão
>> x = [4/3 1.2345e-6]
- subtração ^ Potenciação
é mostrada, para cada formato usado, da seguinte maneira: * multiplicação ’ matriz transposta

format short 1.3333e+000 1.2345e-006

format long 1.333333333333333e+000 1.234500000000000e-006 O GNU Octave possui um conjunto de funções matemáticas elementares, com seno
format hex 3ff5555555555555 3eb4b6231abfd271 (sin), tangente (tan), logaritmo (log10), etc. Por exemplo, para calcular o seno de 30 e
format rat 4/3 9/7290401 (Aproximação racional) guardar na variável x pode-se escrever:
format bank 1.33 0.00
format + ++ » x=sin(30)
x=-0.98803 (para format short)

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 27 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 28
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Funções Elementares (colocadas em ordem alfabética) Funções Elementares (continuação)

mean(x) Média de x
abs(x) Valor absoluto ou módulo de um número
min(x) Menor elemento em x
acos(x) Arco co-seno
sin(x) Função seno
angle(x) Calcula o ângulo de fase (em radianos) para números complexos
sqrt(x) Raiz quadrada
asin(x) Arco seno
std(x) Desvio padrão
atan(x) Arco tangente
sum(x) Soma dos elementos de x
cos(x) Cosseno
tan(x) Tangente
cross(a,b) Produto vetorial dos vetores a e b
exp(x) Exponencial (ex)
inv(x) Matriz inversa da matriz x Algumas dessas operações podem ser aplicadas a vetores ou matrizes. Para a função
log(x) Logaritmo natural (loge x) std(x), por exemplo, se o argumento x for um vetor o resultado será o desvio padrão desse
log10(x) Logaritmo na base 10
max(x) Maior elemento em x vetor. Se x for uma matriz, será calculado o desvio padrão para cada coluna dessa matriz.
O mesmo ocorre com as funções sum, max e mean, entre outras.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 29 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 30

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Usando os comandos vistos, determinar a média, desvio padrão, a soma dos 2.2.1. Comandos
elementos, o valor máximo e mínimo do vetor V abaixo:
Ao sair do GNU Octave (através do comando quit ou exit) todas as variáveis do

» V=[1 2 3 4 5 6 7 8 9] workspace são perdidas, a menos que sejam salvas usando o comando save.

Como usar o comando save, bem como outros relacionados:

save nome_de_arquivo nome_da(s)_variável(is)
Salva somente as variáveis especificadas.
load nome_de_arquivo
Carrega as informações salvas.
clear nome_da_variável
Apaga a variável especificada do workspace.
O clear (sem variável) apaga todas as variáveis.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 31 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 32
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Um dos comandos mais úteis no GNU Octave, bem como em diversos aplicativos, é o
2.2.2. Exercícios
help, que fornece ajuda on-line sobre qualquer outro comando. Por exemplo, para obter
ajuda sobre o comando who:
1) Armazene no workspace os seguintes valores:
» help who
a = 3
WHO – Lista as variáveis no workspace atual. b = -2.123
WHOS – Lista as variáveis e informações adicionais sobre elas. c = 4*3
d = [1 2 3 4]
Outro comando útil: lookfor. Este comando procura e lista todas as funções do
GNU Octave que contém a palavra-chave especificada. e = [1; 2; 3; 4]

» lookfor mean
Utilize os comandos who e whos e observe as diferenças.
Resultado: mean, meansq e mean2.
Ative o comando clear, e digite novamente o comando whos.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 33 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 34

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

2) Armazene no workspace os seguintes valores: Comandos: save, quit (ou exit), load, clear, whos.
a = 3.452
b = -25.123 4) Atribua as seguintes expressões às variáveis:
c = 4*pi a) 4.12 a – pi/c para x
b) log(a + 40.1353) para y
3) Verifique o resultado das seguintes das seguintes operações: c) log(a) para z
a) (a + b) c c) tan(c) – asin(b)
b) (c – b) a / b d) sin(cos(tan(c))) 5) Salve as variáveis x, y, z em um arquivo chamado [Link].

Efetuar os cálculos (formato padrão), depois mudar para o formato long 6) Saia do GNU Octave, entre novamente e carregue as variáveis salvas anteriormente.

7) Apague a variável z e verifique se ela foi realmente apagada.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 35 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 36
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Outros exemplos com incrementos diferentes de um:

3. OPERAÇÕES COM MATRIZES E VETORES

> y = 0:pi/4:pi
O GNU Octave permite a manipulação de linhas, colunas, elementos individuais e
partes de matrizes.
que resulta em:

3.1. Geração de vetores y =

0.0000 0.7854 1.5708 2.3562 3.1416
A geração de uma seqüência de números pode ser feita utilizando:

Incrementos negativos também são possíveis.

> x = 1:5 ou > x = 1:1:5

> z = 6 : -l : l
gera um vetor linha contendo os números de 1 a 5 com incremento unitário. z =
6 5 4 3 2 1

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 37 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 38

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Pode-se também gerar vetores usando a função linspace. 3.2. Elementos das matrizes
Sintaxe: linspace (inicio, fim, número de elementos)
Um elemento específico da matriz pode ser acessado especificando a linha e a coluna
do elemento desejado, fazendo A(linha, coluna). Por exemplo, dada a matriz A:
Exemplo de uso:

A =
> k = linspace (0, l, 6)
l 2 3
4 5 6
k= 7 8 9
0 0.2000 0.4000 0.6000 0.8000 1.0000
Qual o resultado da seguinte operação?

>> A(3,3) = A(1,3) + A(3,l)

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 39 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 40
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Resulta em: “Extração” de submatrizes de uma dada matriz

A = A =
l 2 3 92 99 11 18
4 5 6 98 80 17 14
7 8 10 14 81 88 20

>> B = A(2:3,4) Armazena em B os elementos das linhas 2 e 3 da coluna 4.

B =
14
20

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 41 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 42

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

>> C = A(2:3,2:4) Geração de algumas matrizes especiais:

C =
80 17 14 Exemplo 1: B uma matriz 10x10 unitária.
81 88 20 >> B = ones (10)

C é uma submatriz 2x4, formada pelas linhas 2 e 3 e colunas 2 a 4 da matriz A. Exemplo 2: C uma matriz de dimensão 8x8 com elementos nulos.
>> C = zeros (8)
>> D = triu(A) << Extrai a triangular superior.
D =
92 99 11 18 Matriz identidade: eye
00 88 17 14 Matriz diagonal: diag << Extrai a diagonal de uma dada matriz.
00 00 88 20 Matriz de números randômicos: rand

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 43 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 44
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

» m+100
3.3. Operações com matrizes ans =
108 101 106
Operações matemáticas simples (adição, subtração, divisão e multiplicação)
103 105 107
envolvendo matrizes são semelhantes às operações com escalares.
104 109 102
» m=[8 1 6; 3 5 7; 4 9 2];
Nas operações entre matrizes devem ser respeitadas as regras usuais da matemática
» 3*m quanto ao número de linhas e colunas que duas matrizes devem ter para serem somadas,
multiplicadas, etc. No entanto existem operações especiais. Sendo A=[a1 a2 ... an] e B=[b1
ans = b2 ... bn] duas matrizes, então:
24 3 18
9 15 21 • A./B = [a1/b1 a2/b2 ... an/bn];

12 27 6 • A.*B = [a1·b1 a2·b2 ... an·bn];

• A.^B = [a1^b1 a2^b2 ... an^bn];

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 45 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 46

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

3.4. Exercícios
Alguns comandos adicionais para a operação com matrizes:
1) Dadas as seguintes matrizes e vetores: » A = [1 2 3; 4 10 6; 7 8 19]
» B = [4 5 6; 1 2 3; 8 7 6]
inv Matriz inversa de uma dada matriz
Sintaxe: inv(matriz) » C = [4 5 6]
det Calcula o determinante (D) de uma matriz bem como o » L = [4; 5; 6]
recíproco do número de condição (RCOND)
Sintaxe: [D, RCOND]=det(matriz)
2) Calcule:
trace Traço de uma matriz
Sintaxe: trace(matriz) a) D = A+B g) Determinar a inversa de A e salvar em D
cond Número de condição de uma matriz. b) E = A*B h) Produto D vezes A
Sintaxe: cond(matriz)
rank Rank de uma matriz c) F = A-B i) Autovalor e Autovetor de A
Sintaxe: rank(matriz) d) G = BT-A j) Calcule o traço e o determinante de A
T
e) H = A A k) Resolver o sistema AX = L
T
f) I = A C

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 47 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 48
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

4. GRÁFICOS
Caso o usuário queira saber qual a interface em uso basta digitar graphics_toolkit no
prompt do GNU Octave. Para modificar a interface de saída o comando
A construção de gráficos no GNU Octave é bem simples e a maior parte dos
graphics_toolkit (opção) pode ser usado.
comandos é muito parecida com os usados pelo MATLAB, sendo possível a geração de
As imagens abaixo mostram o mesmo gráfico gerado com as opções “fltk” e
gráficos bidimensionais ou tridimensionais.
“gnuplot”, respectivamente.

Existe no GNU Octave um grande conjunto de comandos que permite a construção de

gráficos.

Na versão atual a saída gráfica é, por default, criada usando o OpenGL e Qt. Caso as
bibliotecas do OpenGL não estejam disponíveis o GnuPlot é utilizado.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 49 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 50

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

4.1. Comandos gráficos básicos

Na tabela seguinte são mostrados alguns comandos básicos que permitem a geração
de gráficos bidimensionais no GNU Octave. Se y é um vetor, plot(y) produz um gráfico linear com o índice dos elementos de y na
abscissa e os elementos de y na ordenada.
plot Visualização de gráficos lineares no plano
bar Gráfico de barras Por exemplo, para plotar os números [0.0, 0.48, 0.84, 1.0, 0.91, 0.6, 0.14], entre com o
stem Sequência discreta
stairs Plotar em degraus vetor y e execute o comando plot, como mostrado:
errorbar Plotar barra de erros
hist Plotar histograma >> y = [0.0 0.48 0.84 1.0 0.91 0.6 0.14];
comet Plotar com trajetória de cometa
>> plot(y)

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 51 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 52
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Resultado:
Se t e y são vetores com dimensões iguais, o comando plot(t,y) produz um gráfico
bidimensional dos elementos de t versos os elementos de y, por exemplo

>> t = 0:0.05:4*pi;
>> y = sin(t);
>> plot(t,y)

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 53 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 54

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Resultado: O GNU Octave pode também plotar múltiplas linhas em apenas um gráfico.

Exemplo de uso do comando plot com múltiplos argumentos.

Exemplo:

>> plot(t,sin(t),t,cos(t),t,0.4*sin(t+pi),t,0.6*cos(t+pi))

Obs : Observe que as informações são fornecidas aos pares, sendo um desles o que será
representado na abscissa e o seguinte na ordenada, e assim sucessivamente.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 55 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 56
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Resultado: 4.2. Estilos de Linha e Símbolo

1
Os tipos de linhas, símbolos e cores usadas para mostrar os gráficos podem ser
0.8
controlados. Seguem alguns exemplos, seguidos dos respectivos resultados (próxima
0.6
página):
0.4 Símbolo: *
0.2
Cor: Preto
>> x = 0:0.05:1;
0
>> subplot(1,3,1); Símbolo: -
-0.2 Cor: Azul
-0.4
>> plot(x,x.^2,'*k')
-0.6 >> subplot(1,3,2);
-0.8 >> plot(x,log(x+0.1),'-b') Símbolo: +
-1 >> subplot(1,3,3); Cor: Vermelha
0 2 4 6 8 10 12 14

>> plot(x,-x.^3,'+r')

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 57 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 58

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Resultado: Outros tipos de linhas, pontos e cores também podem ser usados:

Tipo de Linha Tipo de Ponto Cores

- __________   y amarelo

  m magenta

c cian
 
Mais opções: '--', ':', '-.'
r vermelho
+ +++++++++
g verde
X xxxxxxx
b azul
^ ^^^^^^^
w branco

Mais opções: 's', 'd', 'p', 'h'. k preto

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 59 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 60
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

4.3. Exercícios Resultado:

>> t = 0:0.5:4*pi;
>> y = sin(t);

Criar os seguintes gráficos:

Gráfico simples (plot)

Gráfico de barras (bar)
Sequência discreta (stem)
Escada (stairs)

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 61 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 62

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Dado o vetor: y= [ 1 2 3 7 3 8 8 4 10 1 3 7 8 9 10 6] construa o gráfico de barras bem Resultados: Gráfico de barras e histograma de frequências, respectivamente.
como o histograma (hist) de frequências.

>> bar(y);
>> grid('on');

>> hist(y,50);
>> grid('on');

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 63 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 64
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Verificar o resultado do uso dos seguintes comandos: Resultado final:

>> t=1:pi/10:6*pi;
>> y=(1-t).*sin(t);
>> comet(y)

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 65 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 66

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

4.4. Geração de Gráficos Tridimensionais e Isolinhas

O comando mesh(x,y,z) cria uma perspectiva tridimensional plotando os elementos da
Estes são alguns comandos para gerar gráficos tridimensionais e contornos. matriz z em relação ao plano definindo pelas matrizes x e y. Exemplo de uso:

plot3 Mostrar gráfico no espaço 3D >> t = -2:0.2:2;

contour Permite visualizar curvas de nível no plano
contour3 Permite visualizar curvas de nível no espaço 3D >> [x, y] = meshgrid(t, t);
meshgrid Permite a criação de estruturas (2D e 3D) para representação 3D >> z = x.* exp(-x.^2 - y.^2);
mesh Permite plotar malhas 3D >> mesh(x,y,z)
meshc Combinação mesh/contour
surf Plotar superfície 3D
surfc Combinação surf/contour
slice Plot visualização volumétrica
cylinder Gerar cilindro
sphere Gerar esfera

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 67 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 68
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Resultado:

Verifique as seguintes opções:

>> surf(x, y, z) Alguns outros paletes: autumn, bone, cool, copper, cubehelix, flag, gray, hot, hsv, jet,
>> colormap(gray); lines, ocean, pink, prism, rainbow, spring, summer, white, winter, …
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 69 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 70

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Outros exemplos de superfícies (ainda usando a mesma função) no qual se considera que a Resultados:
superfície é iluminada por uma fonte de luz, dando o efeito de "sombreamento". Seguem
algumas das alternativas:

>> colormap(bone)
>> surfl(x, y, z), shading faceted

>> colormap(bone)
>> surfl(x, y, z), shading flat
(a)... ..(b) (c)
>> colormap(bone) (a) Sombreamento interpolado, com as linhas de transição visíveis (shading faceted)
>> surfl(x, y, z), shading interp (b) Sombreamento interpolado mostrando as faces planas (shading flat)
(c) Sombreamento interpolado (shading interp)

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 71 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 72
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

4.5. Anotações no Gráfico

Curvas de nível
O comando contour(z,10) mostra a projeção da superfície definida (pela função z), no O GNU Octave possui alguns comandos que permitem adicionar informações em um
plano cartesiano xy, com 10 iso-linhas. A figura seguinte mostra o resultado para a gráfico, como por exemplo:
superfície anterior.
title Título do gráfico
xlabel Título do eixo X
ylabel Título do eixo Y
zlabel Título do eixo Z
text Inserir anotação no gráfico
gtext Inserir anotação com o “mouse”, de modo interativo
grid Linhas de grade

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 73 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 74

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Exemplos:

>> fplot(’sin’, [-pi pi]);

>> title ’Gráfico da função f(x)=seno(x), -pi<x<pi’;
>> xlabel ’x’;
>> ylabel ’f(x)’;
>> grid (’on’);
>> text(1,.2,’Curva sin(x)’);
>> gtext ’Sin(x)’;
Exemplo do comando fplot para múltiplas funções

fplot('[tan(x),sin(x),cos(x)]',2*pi*[-1 1 -1 1])

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 75 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 76
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Resultado:
Uma outra possibilidade para a visualização de múltiplas funções:

>> figure(2)
>> t=linspace(-2*pi,2*pi,50);
>> y1=tan(t);
>> y2=sin(t);
>> y3=cos(t);
>>
>> plot(t,y1,"r");
>> hold on
>> plot(t,y2,"g");
>> plot(t,y3,"b");
>> legend('y1','y2','y3');

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 77 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 78

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

5. ARQUIVOS SCRIPT (.m)

Os comandos do GNU Octave são normalmente digitados na área de trabalho, onde

cada linha de comando é introduzida e imediatamente processada. O GNU Octave é
também capaz de executar sequências de comandos armazenadas em arquivos ASCII com Diretório de
extensão m. Estes arquivos podem ser criados por alguns aplicativos como Notepad++, trabalho
Notepad, Textpad, Sublime Text, etc. Janela de
Uma outra possibilidade para criar e editar estes arquivos é por meio do editor do comandos Editor
próprio ambiente de trabalho do GNU Octave. Como pode-se notar na próxima figura,
tem-se na parte inferior da tela a aba Editor. É sugerido que antes de criar o arquivo, com
extensão “m”, que seja escolhido e definido o diretório de trabalho, como indicado na
próxima figura.
GUI do GNU Octave versão 6.2.0, onde são destacados alguns elementos.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 79 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 80
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Exemplo de arquivo script

Na sequência são apresentados alguns comandos e declarações especiais para serem
#Mostrar a funcao y=ax^2 + bx + c no intervalo -5<x<5
usados nos arquivos script, por exemplo: clear;
aux='s';
while aux=='s',
Comando Função clc
% ou # Inserir um comentário no arquivo “.m” disp(' ');
clear Apagar todos os dados da memória disp(' Dada a equação y=ax^2 + bx + c ...'); disp(' ');
input Usado para a entrada de dados a partir da janela a=input('Digite o valor de a = ');
b=input('Digite o valor de b = ');
de comando c=input('Digite o valor de c = ');
pause Provoca uma pausa na execução do arquivo até x=-5:0.1:5;
que qualquer tecla seja digitada y=a*x.^2+b*x+c;
clc Limpar a janela de comando figure(1)
figure(n) Abrir e mostrar a janela gráfica de número n plot(x,y)
grid on
disp(‘...’) Mostrar no ambiente de trabalho o texto clc
colocado entre ‘...’ disp(' ');
close Fechar todas as janelas gráficas aux=input(' Deseja ver outra curva ? (s/n) ==> ','s');
endwhile
disp('Programa finalizado.');
close;

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 81 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 82

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

6. CONTROLE DE FLUXO É muito comum construções em que conjuntos de laços for são usados,
principalmente em operações envolvendo vetores e matrizes:
Os comandos que controlam o fluxo especificam a ordem em que as operações são
realizadas. No GNU Octave estes comandos são semelhantes aos usados na linguagem C for i=1:10;
bem como em outras linguagens. for j=1:10;
A(i,j)=i+j;
6.1. Laço for B(i,j)=i-j;
O laço for é o controlador de fluxo mais simples e usado na programação endfor
GNU Octave. Exemplo de uso do comando for: endfor
A
>> for i = 1:n;
>> for j = 1:n; B
>> ... C=A+B
>> endfor
>> endfor

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 83 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 84
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

6.2. Comando while

No exemplo abaixo o laço while é executado se a condição testada for verdadeira.

Operadores relacionais: Operadores lógicos
< Menor do que & e
a = 1; b = 15; - Inicialmente são atribuídos valores para a e b.
> Maior do que | ou
while a<b - A condição a<b é testada.
<= Menor ou igual a ~ não
a = a+1 - Se ela for verdadeira o corpo do laço será executado
>= Maior ou igual a
b = b-1 e o procedimento é repetido.
== Igual
endwhile - Quando o teste se tornar falso o laço terminará, e a
~= Diferente
disp(’fim do loop’) execução continuará no comando que segue o laço,
após o endwhile.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 85 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 86

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

6.3. Declarações if, else

7. ABERTURA E VISUALIZAÇÃO DE IMAGENS
A seguir, é apresentado um exemplo do uso da declaração if no GNU Octave.
São vários os comandos destinados à visualização e processamento de imagens. Boa
for i = 1:5; parte destes comandos é incorporada em bibliotecas (ou pacotes) específicas como a
for j = 1:5;
if i == j denominada “image”, disponível na página [Link]
A(i,j) = 2; O propósito central desta seção é apresentar apenas alguns comandos básicos, em
else
if abs(i-j) == 1 função das diversas possibilidades de processamento. Dada a diversidade de comandos
A(i,j) = 1; um, ou mais cursos poderiam ser ministrados com este propósito.
else
A(i,j) = 0; Pode-se notar que boa parte dos comandos são similares ao da biblioteca MATLAB.
endif Para detalhes adicionais sobres os diversos comandos sugere-se McANDREW (2004) e
endif
endfor EATON et al (2015).
endfor
A

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 87 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 88
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Abertura de imagem / Visualização / Conversão RBG para Tom de Cinza

A tabela seguinte mostra alguns comandos básicos, que permitem a operação com
imagens no GNU Octave: Crie um arquivo (.m), digite, salve e execute os comandos abaixo:

warning('off');
imread Faz a leitura de imagens pkg load image;
image Permite a visualização de uma matriz como uma imagem nome='[Link]';
hist Calcula o histograma de um conjunto de dados DATA=imread(nome);
rgb2gray Converte imagem colorida (RGB) para tons de cinza
figure(1);
axis Controla a aparência e a escala dos eixos image(DATA);
disp Visualização de um vetor no modo texto xlabel([ 'Imagem: ' nome ]);
figure Abre uma nova janela gráfica disp([ 'Imagem entrada: ' nome ]);
colormap Permite definir um mapa de cores
Exemplo: colormap(gray(256)) figure(2);
imfinfo Função que retorna uma estrutura que contem diversas DATAg=rgb2gray(DATA);
informações sobre determinada imagem image(DATAg);
Exemplo de uso: imfinfo(“arquivo_imagem”) colormap(gray(256));
xlabel([ 'Imagem (tom de cinza) ' ]);
disp([ 'Imagem (tom de cinza) : ok' ]);
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 89 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 90

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Abertura de imagem / Visualização / Conversão RBG para Tom de Cinza (Cont.) Abertura de imagem / Visualização / Conversão RBG para Tom de Cinza (Cont.)

No mesmo arquivo criado antes, escreva um código que permita criar a imagem No mesmo arquivo criado antes, escreve um código que permita criar a imagem
negativa em tom de cinza. Mostre esta imagem: negativa em tom de cinza. Mostre esta imagem:

... ...

figure(3); figure(3);
DATAn=255.-DATAg; DATAn=255.-DATAg;
image(DATAn); image(DATAn);
colormap(gray(256)); colormap(gray(256));
xlabel([ 'Imagem (negativa)' ]); xlabel([ 'Imagem (negativa)' ]);
disp([ 'Imagem (negativa) : ok' ]); disp([ 'Imagem (negativa) : ok' ]);

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 91 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 92
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Leitura / Cálculo do Histograma de Frequência / Visualização do Histograma

Resultado do processamento anterior, onde é mostrada a imagem original (colorida), a

warning('off');
convertida para tons de cinza e a negativa. pkg load image;
nome='[Link]';
DATA=imread(nome);
DATAg=rgb2gray(DATA);
50 50 50

100 100 100

figure(1);
150 150 150

200 200 200

image(DATAg);
250 250 250
colormap(gray(256));
300 300 300
xlabel([ 'Imagem: ' nome ]);
350 350 350
disp([ 'Imagem (original): ' nome ]);
400 400 400
50 100 150 200 250 300 350 400 50 100 150 200 250 300 350 400 50 100 150 200 250 300 350 400
Imagem: [Link] Imagem (tom de cinza) Imagem (negativa) [freq,tom]=hist(DATAg(:),0:255);
figure(2);
plot(tom,freq);
title([ 'Histograma' ]);

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 93 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 94

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Resultado do processamento anterior, onde é mostrada a imagem lida e respectivo Alguns comandos adicionais relacionado ao processamento de imagens:
histograma.
imshow Permite a visualização de uma imagem
Histograma
3500 imagesc Permite a visualização de uma imagem após a aplicação de uma escala,
de modo que todo o mapa de cores seja utilizado
50 3000
imwrite Permite salvar uma imagem em arquivo (imwrite(imagem, “[Link]”))
100
2500 Formatos aceitos: jpg, tif, gif, pgm, bmp, png, etc
150 filter2(B,X) Realiza a correlação do dado em X usando o filtro 2D armazenado em B
2000
200
conv2 (A,B) Realiza a convolução da matriz A sobre a matriz B
1500 uint8 Converte os elementos para inteiros sem sinal (8 bits)
250
int8 Converte os elementos para inteiros de 8 bits
1000
300

350 500

400 0
50 100 150 200 250 300 350 400 0 50 100 150 200 250 300

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 95 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 96
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Leitura / Aplicação de um filtro a uma imagem / Visualização Leitura / Aplicação de um filtro a uma imagem / Visualização (Continuação)

pkg load image; ...

nome='[Link]';
DATA=imread(nome); f55=(1/25).*[ 1 1 1 1 1; 1 1 1 1 1;
DATA=rgb2gray(DATA); 1 1 1 1 1; 1 1 1 1 1;
1 1 1 1 1 ];
figure(1); PARCIALb=conv2(DATA,f55,'same');
image(DATA); figure(3);
colormap(gray(256)); imagesc(PARCIALb);
title([ 'Imagem: ' nome ]); colormap(gray(256));
title([ 'Imagem: ' nome ' [Após fitro média 5x5]' ]);
f33=(1/9).*[ 1 1 1;
1 1 1; f2da=[ -1 0 1; -2 0 2; -1 0 1 ];
1 1 1 ]; PARCIAL=filter2(f2da,DATA);
PARCIALa=filter2(f33,DATA,'same'); figure(4);
figure(2); imshow(uint8(PARCIAL));
image(PARCIALa); title([ 'Imagem: ' nome ' [Apos filtro Sobel (borda
colormap(gray(256)) vertical]' ]);
title([ 'Imagem: ' nome ' [Após fitro media 3x3]' ]);
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 97 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 98

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Imagem: [Link] Imagem: [Link] [Após fitro média 3x3]

Abertura de imagem / Armazenamento em outros formatos
50 50

100 100 No mesmo arquivo script anteriormente escrito, incorpore no código o

150 150

200 200
armazenamento das imagens após aplicação dos filtros de suavização 3x3 e 5x5, em
250 250
arquivos com diferentes formatos:
300 300

350 350

400 400
50 100 150 200 250 300 350 400 50 100 150 200 250 300 350 400

Imagem: [Link] [Após fitro média 5x5] Imagem: [Link] [Após filtro Sobel (borda vertical]
...
50

100 imwrite(uint8(PARCIALa),"[Link]");
150
imwrite(uint8(PARCIALb),"[Link]");
200
disp(['Imagens [Link] e [Link] salvas.']);
250

300

350

400
50 100 150 200 250 300 350 400

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 99 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 100
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Abertura de imagem / Armazenamento em outros formatos 8. IMPORTAÇÃO / EXPORTAÇÃO DE DADOS

No mesmo arquivo script anteriormente escrito, incorpore no código o

Neste tópico são apresentados alguns comandos adicionais relacionados à
armazenamento das imagens após aplicação dos filtros de suavização 3x3 e 5x5, em
manipulação de arquivos e ao uso de funções criadas pelo usuário.
arquivos com diferentes formatos:

8.1. Importação e Exportação de Dados

...

imwrite(uint8(PARCIALa),"[Link]"); Os dados (matrizes, vetores, escalares, cadeia de caracteres, etc) disponíveis no

imwrite(uint8(PARCIALb),"[Link]"); Ambiente de Trabalho do GNU Octave, como foi visto, podem ser armazenados em
disp(['Imagens [Link] e [Link] salvas.']);
arquivos no modo texto, binário (e outros formatos), utilizando o comando save.
Existem diversas maneiras de utilizar este comando, como pode-se ver nos exemplos
mostrados na sequência:

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 101 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 102

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Por exemplo suponha que um programa qualquer criou arquivo “[Link]” que
save [Link] X Y Z –binary salva as variáveis X, Y e Z no arquivo binário contém a seguinte matriz.
“[Link]”
save [Link] X Y Z -ascii salva as variáveis X., Y e Z no arquivo texto
“[Link]” com 8 casas decimais 11
save [Link] X Y Z -ascii -double salva as matrizes X., Y e Z no arquivo texto 31
“[Link]” com 16 casas decimais 33
13
11
Os dados obtidos por outros programas podem ser importados pelo GNU Octave,
desde que estes dados sejam gravados em formato apropriado. Se os dados são Ao executar o comando:
armazenados no formato ASCII, e no caso de matrizes, com colunas separadas por >> load [Link]
espaços e cada linha da matriz em uma linha do texto, o comando load pode ser usado.

o GNU Octave importa a matriz, que é armazenada com o nome polig1.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 103 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 104
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

O GNU Octave pode também importar (através do comando load) os dados que foram
anteriormente exportados por ele. Por exemplo, para importar as variáveis X, Y e Z,
» polig1 anteriormente exportadas usando o comando save, pode-se fazer:
polig1 =
1 1 load arq1 save arq1 X Y Z
load [Link] save [Link] X Y Z -ascii
3 1
load [Link] save [Link] X Y Z -ascii -double
3 3
1 3 Deve-se ressaltar que o comando save, quando usado para exportar os dados do

1 1 GNU Octave em formato texto, exporta apenas um bloco contendo todas as variáveis.
Quando os dados são importados através do comando load, apenas uma variável com
nome do arquivo é criada.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 105 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 106

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Exemplo:

>> x = rand(3,3) >> save [Link] X Y -ascii

x = >> clear
0.2190 0.6793 0.5194 >> whos
0.0470 0.9347 0.8310 >> load [Link]
0.6789 0.3835 0.0346 >> arq2
arq2 =
>> y = rand(3,3)
0.2190 0.6793 0.5194
y =
0.0470 0.9347 0.8310
0.0535 0.0077 0.4175
0.6789 0.3835 0.0346
0.5297 0.3835 0.6868 0.0535 0.0077 0.4175
0.5297 0.3834 0.6868
0.6711 0.0668 0.5890
0.6711 0.0668 0.5890

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 107 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 108
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

9. CRIAÇÃO DE FUNÇÕES
Na sequência são mostrados dois exemplos de funções. Uma delas foi criada com o
9.1. Criação de funções objetivo de, dado um ângulo (alfa) em radianos, montar a matriz de rotação M:

O uso de funções é muito útil tanto na execução de algumas tarefas repetitivas

 cos() sin() 
quanto no aproveitamento do código em diferentes aplicações. Pode-se, deste modo, fazer M ( )   
a criação bibliotecas de funções destinadas a solução de alguns problemas específicos.
 sin() cos() .

A estrutura geral de uma função do GNU Octave é a seguinte: A segunda função tem o propósito de, dada uma matriz de rotação M e um ponto de
coordenas (x,y), aplicar a matriz de rotação de modo que se tenha:

function [ retorno ] = teste (variáveis)

...  x ' x 
 y '  M (  )  y 
endfunction    

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 109 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 110

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Função “f_matriz_rotacao” armazenada no arquivo f_matriz_rotacao.m Função “f_aplica_rotacao” armazenada no arquivo f_aplica_rotacao.m

function Mrot=f_matriz_rotacao(alfa) function [xrot,yrot]=f_aplica_rotacao(M,x,y)

# F_MATRIZ_ROTACAO Calcula a matriz de rotação no plano # F_APLICA_ROTACAO Faz a rotação de um ponto (x,y) usando
# Entrada: Ângulo em graus # a matriz de rotação M
# Retorno: Matriz de rotação 2x2 # Retorno: Coordenadas rotacionadas do
# # ponto (x,y)
# Autores: Mauricio Galo e Paulo Camargo #
# UNESP, Presidente Prudente, SP # Autores: Mauricio Galo e Paulo Camargo
# UNESP, Presidente Prudente, SP
alfa=alfa*pi/180;
Mrot=[ cos(alfa) sin(alfa); xrot=x*M(1,1) + y*M(1,2);
-sin(alfa) cos(alfa)]; yrot=x*M(2,1) + y*M(2,2);

Endfunction endfunction

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 111 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 112
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Leitura dos dados / Visualização dos polígonos

Como exemplo de aplicação das funções criadas para o GNU Octave, inicialmente
escreva um script que permite ler nas variáveis [Link] e [Link], bem como permita
a visualização destes vetores, usando o comando plot: load '[Link]'
load '[Link]'
Conteúdo do Conteúdo do [lin1,col1]=size(polig1);
arquivo [Link] arquivo [Link] [lin2,col2]=size(polig2);
0 0 2 2
1 0 3 3 x1=polig1(1:lin1,1);
1 1 2 3 y1=polig1(1:lin1,2);
0 1 2 2 x2=polig2(1:lin2,1);
0 0 y2=polig2(1:lin2,2);

figure(1);
plot(x1,y1,x2,y2);

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 113 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 114

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Cálculo da Matriz de Rotação / Rotação / Visualização dos polígonos após rotação Visualização dos polígonos antes e após rotação

No mesmo script anterior, calcule e matriz de rotação, aplique-a a todos os pontos 3

dos dois polígonos e faça a visualização.
2.5
... 2
ROT = f_matriz_rotacao(30);
1.5
[x1r,y1r]=f_aplica_rotacao(ROT,x1,y1);
[x2r,y2r]=f_aplica_rotacao(ROT,x2,y2);
1
figure(2);
plot(x1r,y1r,x2r,y2r); 0.5

figure(3); 0
plot(x1,y1,x2,y2,x1r,y1r,x2r,y2r);
-0.5
0 1 2 3 4 5

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 115 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 116
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

PAGAMISSE, A.; SOUZA, L. H. G. Introdução ao Software Octave. Semana de Cursos de

REFERÊNCIAS Matemática, Estatística e Computação, 25-29 de agosto de 2003, FCT/UNESP, Presidente Prudente,
2003. 54p.
CORAL, A. M.; SANTOS, M. P.; BASTOS, T. D. A.; BORBA, M. Curso de Matlab. Universidade SIGMON, K. Matlab Primer – Third Edition. Department of Mathematics, University of Florida,
Federal de Santa Catarina, Dep. de Eng. de Produção, Programa Especial de Treinamento – PET. Gainesville, 1993. 35p. Disponível em [Link] Acesso em
Florianópolis – SC, 1999. 29p. Março/2021.
EATON, J. W.; BATEMAN, D.; HAUBERG, S. GNU Octave – Edition 3 for Octave version 3.2.3, SIQUEIRA, A. F. Octave – Seus primeiros passos na programação científica. São Paulo: Casa do
July, 2007. 672p. Código, 2015. 202p.
EATON, J. W.; BATEMAN, D.; HAUBERG, S.; WEHBRING, R. GNU Octave – Edition 4 for ZERI, L. M. M. Apostila de Matlab. Instituto Nacional de Pesquisas Espaciais – INPE, 2001. 19p.
Octave version 4.0.0 – Free Your Numbers, March, 2015. 966p.
HANSELMAN, D.; LITTLEFIELD, B. Matlab - Versão do Estudante: guia do usuário (Tradução).
São Paulo: Makron Books, 1997. 305p.
McANDREW, A. Introduction to Digital Image Processing with MATLAB®. Thomson Course
Technology, 2004. 509p. ISBN: 0-534-40011-6
QUARTERONI, A.; SALERI, F.; GERVASIO, P. Scientific Computing with MATLAB and Octave.
Fourth Edition. Text in Computational Science and Engineering. BARTH, T. J.; GRIEBEL, M.;
KEYES, D. E.; NIEMINEN, R. M.; ROOSE, D.; SCHLICK. T. (Eds.). Lausanne: Springer, 2014.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 117 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 118

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo I
ANEXOS %{
Exemplo de uso de alguns comandos do GNU Octave
Comandos (principais): subplot, bar, hist.
Nesta seção serão apresentados alguns exemplos adicionais de comandos e
respectivos resultados, que podem ser úteis em algumas aplicações na visualização de Autores: M. Galo e Paulo de O. Camargo
dados em algumas áreas das ciências exatas e engenharias, usando o aplicativo Unesp, Presidente Prudente, SP
GNU Octave. %}
clear
y = [1 2 3 7 3 8 8 4 10 1 3 7 8 9 10 6];
figure;
subplot(2,1,1);
bar(y,'g');
xlabel('Grafico de barras');
grid
subplot(2,1,2);
hist(y,50,'b')
xlabel('Histograma de frequencia');
grid
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 119 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 120
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo II

%{
Exemplo de uso de alguns comandos do GNU Octave
Comandos (principais): bar, legend.

Autores: M. Galo e Paulo de O. Camargo

Unesp, Presidente Prudente, SP
%}

clear
figure;
periodo = 2012:1:2016;
dados = [ 19.1 23.5 23.7 24.7 24.4;
23.5 21.5 23.5 23.7 23.3;
18.5 20.8 23.0 21.3 21.3;
17.0 19.8 19.7 21.9 27.5;
14.0 17.1 12.5 15.2 0;
Resultado gráfico do script anterior. 14.3 14.3 14.5 13.9 0;
(continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 121 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 122

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

12.9 7.0 14.1 15.1 0;

18.8 9.0 15.5 16.2 0;
15.1 12.0 18.4 17.4 0;
18.0 19.2 20.4 19.3 0;
22.8 21.5 22.1 22.2 0;
21.1 25.0 22.9 22.0 0 ];
B=bar(dados);
legend(B,'2012','2013','2014','2015','2016');
legend('Boxoff');
meses=({'Jan','Fev','Mar','Abr','Mai','Jun','Jul','Ago','Set
','Out','Nov','Dez'});
set(gca,'XTickLabel',meses);
ylabel('Temperatura minima mensal [graus C]');
Resultado gráfico do script anterior.

Fonte dos dados: CIIAGRO – Centro Integrado de Informações Agrometeorológicas (site:

Obs.: Fonte dos dados: [Link]. Os dados referentes ao mês de abril/2016 [Link])
correspondem à coleta de dados dos primeiros 15 dias de abril, uma vez que os dados
foram obtidos neste site em 16/Abril/2016.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 123 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 124
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo III

%{
Exemplo de uso de alguns comandos do GNU Octave
Comandos (principais): pie.

Autores: M. Galo e Paulo de O. Camargo

Unesp, Presidente Prudente, SP
%}

PAIS = {'CHINA 49%' 'JAPAO 12%' 'ALEMANHA 6,7%' 'COREIA DO

SUL 3,8%' 'FRANCA 2,6%' 'HOLANDA 2,2%' 'ITALIA 2,6%' 'OUTROS
21,1%'};
DATA = [ 49 12 6.7 3.8 2.6 2.2 2.6 21.1 ];
figure;
pie(DATA,PAIS);
S=['PRINCIPAIS COMPRADORES DE MINERIO DE FERRO BRASILEIRO EM
2010. Fonte: IBRAM'];
title(S,'fontsize',11);
Gráfico resultante do script anterior.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 125 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 126

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo IV

%{
Exemplo de uso de alguns comandos do GNU Octave
Comandos (principais): plot, legend, uso de letras gregas.

Autores: M. Galo e Paulo de O. Camargo

Unesp, Presidente Prudente, SP
%}

alpha = 2;
beta = 2.4;
Delta = pi;
t=0:0.1:2*pi;
plot(t,alpha*sin(t),t,beta*sin(t/2),t,Delta*sin(t/2.5));
xlabel('Angulo t [rad], \alpha = 2, \beta = 2.4, \Delta =
\pi');
legend({'\alpha*sin(t)','\beta*sin(t/2)','\Delta*sin(t/2.5)'
},'location','southwest');
grid Gráfico resultante do script anterior.
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 127 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 128
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo V S = [ 'Angulo t [rad], Fator de escala: \lambda='

num2str(lambda) ];
%{ xlabel(S);
Exemplo de uso de alguns comandos do GNU Octave
Comandos (principais): num2str, linewidth. Gráfico resultante
do script anterior.
Autores: M. Galo e Paulo de O. Camargo
Unesp, Presidente Prudente, SP
%}

lambda=input('Entre com o fator de escala: ');

figure(1)
t=0:0.2:6*pi;
y=(1-t).*sin(t);
y1=lambda*(1-t).*sin(t);
plot(t,y,'linewidth',2,t,y1,'linewidth',2);
ylabel('f(t)');
legend('y(t)','\lambda.y(t)','location','northwest');
legend('boxoff');
(continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 129 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 130

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo VI [lin,col]=size(xy);
x=xy(1:lin,1);
%{ y=xy(1:lin,2);
Exemplo de uso de alguns comandos do GNU Octave figure
Comandos (principais): size, delaunay, triplot, axis.
subplot(1,2,1)
Autores: M. Galo e Paulo de O. Camargo plot(x,y,'sr');
Unesp, Presidente Prudente, SP title('Pontos da amostra.');
%} axis('square');
grid;
close all;
clear; subplot(1,2,2)
clc; plot(x,y,'sr');
tri=delaunay(x',y');
xy = [ hold on;
6.9 4.9 1.1 1.2 7.9 5.7 2.4 4.6 3.1 3.4 5.5 9.1 7.1 8.5 2.3; triplot(tri,x,y);
7.1 5.7 7.7 2.6 1.2 0.8 0.6 3.0 4.7 7.3 8.7 6.1 4.6 2.4 title('Pontos da amostra e Triangulacao de Delaunay.');
5.5]; axis('square');
xy=xy'; grid;
(continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 131 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 132
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo VII

%{
Exemplo de uso de alguns comandos do GNU Octave
Comandos (principais): subplot, bar, hist.

Autores: M. Galo e Paulo de O. Camargo

Unesp, Presidente Prudente, SP
%}

clear
load [Link]
[lin,col]=size(terreno);
x=terreno(1:lin,1);
y=terreno(1:lin,2);
z=terreno(1:lin,3);
xn = x - min(x);
Gráficos resultantes do script anterior. yn = y - min(y);
zn = z;
(continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 133 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 134

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

figure; contour(zf,50,'-b');
plot(xn,yn,'.r');
xlabel('X(m)'); xlabel('X(m)');
ylabel('Y(m)'); ylabel('Y(m)');

figure; figure;
plot3(xn,yn,zn,'.r'); surfc(xf, yf, zf);
xlabel('X(m)'); title('Superficie topografica');
ylabel('Y(m)'); xlabel('X(m)');
zlabel('Z(m)'); ylabel('Y(m)');
zlabel('Z(m)');
limX = min(xn):1.:max(xn);
limY = min(yn):1.:max(yn); Na sequência são mostrados o arquivo [Link], usado no script anterior, bem
[xf,yf] = meshgrid(limX, limY); como os quatro gráficos resultantes da execução deste script.
zf = griddata(xn, yn, zn, xf, yf);

figure;
plot(xn,yn,'.r');
hold on
(continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 135 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 136
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Conteúdo do arquivo [Link].

10605.777 14253.405 847.374 10654.248 14300.966 832.379 10645.834 14366.175 850.873
10599.733 14226.811 849.188 10677.961 14273.010 821.125 10569.665 14349.903 855.487
10615.293 14198.717 843.856 10695.797 14296.792 830.147 10566.565 14352.056 855.569
10634.561 14206.868 835.203 10687.212 14311.491 835.304 10551.639 14366.817 853.152
10630.045 14231.201 834.698 10702.598 14296.371 825.103 10598.039 14338.327 855.792
10626.659 14263.897 834.841 10727.848 14302.539 825.620 10586.884 14319.117 855.908
10646.331 14232.702 827.488 10722.865 14257.941 816.945 10573.056 14334.192 856.122
10650.098 14268.751 824.985 10757.893 14332.396 825.310 10555.576 14312.625 853.111
10665.250 14256.447 818.154 10724.310 14329.433 829.970 10558.450 14312.015 853.120
10664.420 14255.853 818.118 10720.952 14330.090 832.536 10535.919 14352.970 850.769
10664.845 14255.081 818.081 10702.665 14337.463 837.446 10527.406 14327.464 849.723
10665.840 14255.707 818.115 10721.728 14311.994 828.226 10555.024 14277.472 853.647
10673.666 14244.882 817.800 10752.754 14350.186 825.874 10552.173 14277.784 853.858
10673.593 14246.394 817.840 10719.296 14350.783 829.237 10557.990 14241.363 854.902
10675.106 14247.072 817.834 10724.365 14348.249 834.469 10560.678 14241.972 854.754
10671.091 14244.626 818.551 10693.065 14355.301 843.320 10568.406 14262.295 854.681
10678.698 14242.686 814.706 10734.297 14354.692 830.744 10576.697 14305.307 855.266
10688.459 14268.374 821.455 10695.266 14380.014 843.335 10523.849 14268.045 851.384
10691.294 14259.813 816.981 10617.750 14403.847 851.262 10587.444 14202.788 849.155
10695.599 14275.018 823.062 10621.073 14404.658 851.418 10586.427 14199.660 848.979
10697.511 14260.744 815.977 10609.384 14441.500 854.433 10579.011 14232.853 854.422
10687.189 14244.902 812.722 10587.327 14446.788 851.514 10547.881 14213.814 850.228
10695.143 14239.535 805.741 10598.051 14374.378 853.819 10595.544 14177.230 844.544 Gráficos 1 e 2 resultantes do script anterior.
10675.229 14232.082 816.265 10599.717 14372.089 853.823 10598.582 14179.058 844.836
10687.963 14217.071 809.833 10618.061 14389.550 852.020 10544.620 14192.465 846.335
10595.478 14287.307 851.135 10615.075 14390.768 851.944 10589.274 14162.826 842.976
10614.218 14318.957 850.357 10587.427 14423.912 848.785 10587.583 14166.114 842.864
10642.839 14346.278 847.213 10572.092 14397.026 852.387 10695.918 14217.893 806.531
10678.690 14326.350 838.391 10652.161 14396.553 853.625
# Fim do arquivo
(continua) (continua)

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 137 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 138

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo VIII

%{
Exemplo de uso de alguns comandos do GNU Octave
Comandos (principais): meshgrid, surf, interp2.

Autores: M. Galo e Paulo de O. Camargo

Unesp, Presidente Prudente, SP
%}

clear
Lmin = -2*pi;
Lmax = 2*pi;
[x,y] = meshgrid(Lmin:Lmax);
z = (1/3)*x.*cos(y/2);
Gráficos 3 e 4 resultantes do script anterior.
figure;
surf(x,y,z);

(continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 139 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 140
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

[xi,yi] = meshgrid(Lmin:0.5:Lmax);
zi = interp2(x, y, z, xi, yi);
figure;
surf(xi,yi,zi);

zi1 = interp2(x, y, z, xi, yi,'nearest');

figure;
surf(xi,yi,zi1);
a) b)
zi2 = interp2(x, y, z, xi, yi,'cubic');
figure;
surf(xi,yi,zi2);
Superfícies 3 a 5 resultantes do script
zi3 = interp2(x, y, z, xi, yi,'spline');
figure; anterior, onde os interpoladores “nearest”
surf(xi,yi,zi3); (a), “cubic” (b) e “spline” (c) foram
% utilizados.
% Fim do script
% c)

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 141 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 142

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo IX
% Média e desvio da amostra 1
%{ am1_med = 0;
Visualização da curva normal am1_des = 1;
Comando (principal): normpdf
% Média e desvio da amostra 2
Autores: M. Galo e Paulo de O. Camargo am2_med = -1.3;
Unesp, Presidente Prudente, SP am2_des = 2.5;
%}
% Posição correspondente ao desvio-padrao dado
clear all; L1 = [ am1_med-am1_des; am1_med+am1_des ];
clc L2 = [ am2_med-am2_des; am2_med+am2_des ];
pkg load statistics
figure(1);
disp(' ') x = -9:0.1:9;
disp(' ############################################ ');
disp(' Visualizacao da curva normal no GNU Octave '); % Visualização da curva 1
disp(' ############################################ ');
disp(' ') (continua)
(continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 143 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 144
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

plot(x,normpdf(x,am1_med,am1_des),'-b','linewidth',2);
hold on
stem(L1,normpdf(L1,am1_med,am1_des),':b','linewidth',2);
hold on

% Visualização da curva 2

plot(x,normpdf(x,am2_med,am2_des),'-g','linewidth',2);
hold on
stem(L2,normpdf(L2,am2_med,am2_des),':g','linewidth',2);
grid
title('Distribuicao Normal');
hold off

%
% Fim do script "Curvas_normais.m"
%
Resultado do script mostrando as duas curvas normais.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 145 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 146

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo X % Distribuição Normal Padronizada (z)

dist_z = norminv(1-NS);
%{
Exemplo de uso de comandos do GNU Octave, % Distribuição t
que permitem o cálculo de distribuições dist_t = tinv(1-NS, GL);
estatísticas (Normal, t e Qui-quadrado)
Comandos: norminv, tinv, chi2pdf. % Distribuição Qui-Quadrado
dist_QQ = chi2inv(1-NS, GL);
Autores: M. Galo e Paulo de O. Camargo
Unesp, Presidente Prudente, SP disp(' ')
%} disp(' ############################################ ');
disp(' Exemplo de uso de algumas* ');
clear all; clc; disp(' funcoes estatisticas no GNU Octave ');
pkg load statistics disp(' ############################################ ');
disp(' ')
% NS - Nivel de Significância (%) S=sprintf(' Nivel de significancia [%%]: %5.3f',NS);
% GL - Graus de Liberdade disp(S)
NS = 0.1; S=sprintf(' Graus de liberdade : %d',GL);
GL = 25; disp(S)
(continua) (continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 147 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 148
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

S=sprintf(' z teorico : %8.3f',dist_z);

disp(S) Resultado do script anterior.
S=sprintf(' t teorico : %8.3f',dist_t);
disp(S)
S=sprintf(' QQ teorico : %8.3f',dist_QQ);
disp(S)
disp(' ')
disp(' * Compare os valores determinados com os');
disp(' obtidos por uma tabela estatistica. ');
disp(' ')

%
% Fim do script 'Distribuicoes_exemplo.m'
%

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 149 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 150

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo XI for i=1:linhas1

x1(i)=poligonoAt1(i,1);
%{ y1(i)=poligonoAt1(i,2);
Exemplo de uso de comandos do GNU Octave endfor
Comando (principal): fill
[linhas2,colunas]=size(poligonoAt2);
Autores: M. Galo e Paulo de O. Camargo for i=1:linhas2
Unesp, Presidente Prudente, SP x2(i)=poligonoAt2(i,1);
%} y2(i)=poligonoAt2(i,2);
endfor
clear all;
close all; figure(1)
clc; plot(x1,y1,'r',x2,y2,'g');
axis('square');
load [Link] title('Polígonos - Sem preenchimento');
load [Link]
figure;
[linhas1,colunas]=size(poligonoAt1); h=fill(x1,y1,'r',x2,y2,'g');

(continua) (continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 151 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 152
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

axis('square'); Conteúdo dos arquivos [Link] e [Link].

title('Polígonos - Com preenchimento');

%
% Fim do script 'Preenchimento.m'
%

Resultados do script preenchimento.m.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 153 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 154

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo XII for i=1:linhas1

x1(i)=poligonoAt1(i,1);
%{ y1(i)=poligonoAt1(i,2);
Exemplo de uso de comandos do GNU Octave endfor
Comandos (principais): fill, facealpha
[linhas2,colunas]=size(poligonoAt2);
Autores: M. Galo e Paulo de O. Camargo for i=1:linhas2
Unesp, Presidente Prudente, SP x2(i)=poligonoAt2(i,1);
%} y2(i)=poligonoAt2(i,2);
endfor
clear all;
close all; figure(1);
clc; h=fill(x1,y1,'r',x2,y2,'g');
set(h,'facealpha',0.9);
load [Link] axis('square');
load [Link] title('Polígonos (Opacidade=0.9)');

[linhas1,colunas]=size(poligonoAt1); figure;

(continua) (continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 155 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 156
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

h=fill(x1,y1,'r',x2,y2,'g'); Resultados do script transparencia.m para diferentes valores do parâmetro

set(h,'facealpha',0.5); facealpha.
axis('square');
title('Polígonos (Opacidade=0.5)');
facealpha = 0 (Transparência máxima / Opacidade mínima)
figure; facealpha = 1 (Transparência mínima / Opacidade máxima)
h=fill(x1,y1,'r',x2,y2,'g');
set(h,'facealpha',0.1);
axis('square');
title('Polígonos (Opacidade=0.1)');

%
% Fim do script 'transparencia.m'
%

Resultados para facealpha=0.9, facealpha= 0.5 e facealpha = 0.1, respectivamente.

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 157 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 158

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo XIII
% Determinação da dimensão do arquivo
%{ [pontos,colunas]=size(residuos_calibracao);
Exemplo de uso de comandos do GNU Octave
Comandos: quiver % Colunas 1 e 2 – Coordenadas x e y
% Colunas 3 e 4 – Distorção em x e y, respectivamente
Autores: M. Galo e Paulo de O. Camargo for i=1:pontos
Unesp, Presidente Prudente, SP x(i) = residuos_calibracao(i,1);
%} y(i) = residuos_calibracao(i,2);
dx(i) = residuos_calibracao(i,3);
clear all; dy(i) = residuos_calibracao(i,4);
close all; end
clc;
figure(1);
% Leitura dos dados em arquivo plot(x,y,'+b');
load residuos_calibracao.dat title('Distribuição dos pontos');
axis([ -17 17 -13 13 ]);
% Fator de escala a ser aplicado nos vetores dos resíduos axis('equal');
esc = 1;
(continua) (continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 159 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 160
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

figure;
quiver(x,y,dx,dy,esc,'b'); Trecho do arquivo residuos_calibracao.dat.
title('Vetores de distorção'); % Coordenadas (x,y) e distorção* ...
axis([ -17 17 -13 13 ]); % #1 #2 - Coordenadas x e y (mm) ...
% #3 #4 - Distorções em x e y (mm) 3.150 8.400 -0.019 -0.045
axis('equal'); % 6.300 8.400 -0.049 -0.062
% *Distorção Radial Simétrica 9.450 8.400 -0.104 -0.089
12.600 8.400 -0.194 -0.127
% -15.750 -10.500 0.313 0.212 15.750 8.400 -0.331 -0.175
-12.600 -10.500 0.185 0.158 -15.750 10.500 0.320 -0.224
% Fim do script 'Vetores_no_Plano.m' -9.450 -10.500 0.101 0.117 -12.600 10.500 0.191 -0.168
-6.300 -10.500 0.049 0.088 -9.450 10.500 0.105 -0.125
% -3.150 -10.500 0.018 0.071 -6.300 10.500 0.052 -0.095
0.000 -10.500 -0.003 0.067 -3.150 10.500 0.019 -0.078
3.150 -10.500 -0.026 0.075 0.000 10.500 -0.003 -0.074
Na sequência é mostrado trecho do arquivo de dados usado no script 6.300 -10.500 -0.063 0.095 3.150 10.500 -0.028 -0.082
9.450 -10.500 -0.123 0.128 6.300 10.500 -0.066 -0.103
Vetores_no_Plano.m. 12.600 -10.500 -0.220 0.174 9.450 10.500 -0.128 -0.137
15.750 -10.500 -0.364 0.231 12.600 10.500 -0.226 -0.184
-15.750 -8.400 0.277 0.149 15.750 10.500 -0.371 -0.244
-12.600 -8.400 0.157 0.107
-9.450 -8.400 0.080 0.074 % Fim do arquivo
-6.300 -8.400 0.036 0.050
...
...

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 161 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 162

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

Anexo XIV
As imagens, na sequência, mostram a distribuição dos pontos (a) e os vetores das
distorções para cada pontos (b), a partir dos dados lidos no arquivo %{
residuos_calibracao.dat, resultantes do script Vetores_no_Plano.m. Raios de Curvatura

Autores: M. Galo e Paulo de O. Camargo

Unesp, Presidente Prudente, SP
%}

clear all; close all; clc;

Parametros_do_Elipsoide;

dFonte_Numeros = 15;
dFonte_Texto = 18;
fTexto = 'arial';

(a) (b) % Determinacao das constantes do elipsoide

Distribuição dos pontos em uma imagem (a) e vetores das distorções em cada pontos (b). a = a/1000.; % Conversão para [km]
b = a*(1. - f);
(continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 163 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 164
 TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

%legend('boxoff');
E2 = f*(2. - f); title('Raios de curvatura - Elipsoide', 'fontsize',
dFonte_Texto, 'fontname', fTexto);
ind=0; xlabel('Latitude [graus]','fontsize', dFonte_Texto,
for lat=0:1:90 'fontname', fTexto);
ind=ind+1; ylabel('Raios de curvatura [km]', 'fontsize', dFonte_Texto,
t(ind)=lat; 'fontname', fTexto);
faux = sqrt(1. - E2*(sin(lat*pi/180)**2)); grid;
fN(ind) = a / faux;
fM(ind) = a*(1. - E2) / ( faux**3 ); %
fMN(ind) = a*((1. - E2)**0.5)/(faux**2); % Fim do script 'Raios_de_Curvatura.m'
endfor %

figure(1); Na sequência é mostrado o arquivo “Parametros_do_Elipsoide.m”, que

h=plot(t,fN,t,fM,t,fMN); contém os parâmetros do elipsoide utilizado. Pode-se observar que o conteúdo desse
set(h,"linewidth",1); arquivo é lido pelo script a partir da seguinte linha de comandos:
set(gca, 'fontsize', dFonte_Numeros);
set(gca, 'fontname', fTexto, 'xtick',0:15:90); Parametros_do_Elipsoide;
legend('N','M','Raio Medio');
legend('location','southeast');
(continua)
Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 165 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 166

TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia TUTORIAL do GNU Octave / GNU Octave TUTORIAL UNESP / Departamento de Cartografia

O resultado do script anterior é mostrado na figura abaixo.

Conteúdo do arquivo Parametros_do_Elipsoide.m.

%{
Arquivo: Parametros_do_Elisoide.m
%}

0;

% Elipsoide SAD 69
% a: semi-eixo [m]
% f: achatamento

a = 6378160;
f = 1./298.25;

% Fim do arquivo

Gráfico mostrando os raios de curvatura do elipsoide de revolução.

### Fim do arquivo ###

Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 167 Mauricio Galo e Paulo de Oliveira Camargo, Presidente Prudente – SP, 2023 168