Algebra Linear Numerica
Esta é a página oficial do curso de Álgebra Linear Numérica (graduação) do IPRJ, contendo material potencialmente útil ao publico geral.
Informacao Basica
- Instrutor: prof. Ricardo Fabbri, Ph.D.
Horario
Atualizado 2017:
- Tercas 2:20pm-5pm Sala 211
- Quintas 1:20pm-4pm Sala 211
Avisos
Software
- Usaremos o Scilab primariamente.
- Ver funcoes svd, sva, e '\' (barra invertida)
- Python e Octave sao alternativas aceitaveis
Tarefas e Projetos
As tarefas devem ser formatadas com notacao matematica adequada, preferencialmente em Latex.
Somente serao aceitos arquivos eletronicos no formato PDF ou outro formato aberto como .odt
Quando a tarefa involver qualquer programacao, o aluno devera enviar o codigo fonte. O codigo junto com a documentacao devera estar dentro de um unico diretorio comprimido com .zip ou tar, com o nome do aluno, disciplina e data.
Enviar todo email com a tag "ALNIPRJ" no subject para garantir recebimento.
Tarefa 0 (entrega 7Set17)
- Revisar algebra linear
- Assistir aos videos do Gilbert Strang, em especial:
- The four fundamental subspaces - lecture 10 (garanta a compreensao deste video, assita aos anteriores se necessario)
- Eigenvalues and Eigenvectors - lecture 21
- SVD - lecture 29
- Assistir aos videos do Gilbert Strang, em especial:
- Digitar um resumo.
Tarefa 1 (entrega 14Set17)
- SVD + Scilab
- Plotar circulo unitario
- Escolha uma matriz A
- Transformar cada ponto por A
- SVD visual: visualizar os eixos e esticamento da transformacao efetuada por A
- Realizar SVD usando a funcao SVD propriamente dita
- Repetir com outras matrizes A, aleatorias e nao-aleatorias. Qual o poder de modelagem de matrizes na pratica que voce pode observar assim?
- Digitar o resultado, preferencialmente em Latex, e entregar por email ou impresso.
Tarefa 2 (entrega 5/Out/2017)
- Repetir Tarefa 1 para matriz 3x3 em coordenadas homogeneas - modelam rotacoes, translacoes e homografias em 2D de maneira linear em 3D
- O circulo unitario se transforma para qual objeto?
Tarefa 3: Ajuste de curvas (entrega 12/Out/2016)
Uma das principais aplicações dos métodos numéricos de álgebra linear vistos em aulas consiste na solução de sistemas por mínimos quadrados. Isto, por sua vez, é muito útil para ajustar curvas algébricas (polinomiais) a dados experimentais, por exemplo: encontrar retas, círculos, ou outras cônicas que passam próximos de dado um conjunto de pontos.
- Ajuste de reta (opcional - para aprendizado apenas)
- simule a geração de um conjunto sintético de pontos "experimentais" ao longo de uma reta e os perturbe com a função rand() ao longo de cada coordenada
- use a "barra invertida" '\' do Scilab para gerar a curva que mais se aproxima dos pontos gerados
- plote seus resultados para diferentes intensidades de ruído
- Ajuste de cônicas (obrigatorio)
- simule a geração de um conjunto sintético de pontos "experimentais" ao longo de uma curva polinomial bi-dimensional de grau 2 (círculo, parábola, elipse..) e os perturbe com a função rand() ao longo de cada coordenada
- use a "barra invertida" '\' do Scilab para gerar a curva que mais se aproxima dos dados gerados
- plote seus resultados para diferentes intensidades de ruído
Projeto de Reconhecimento de Faces
Exemplo de Album com Algumas Faces (Facebook upload)
Exemplo de Faces detectadas, a serem reconhecidas no projeto (logo apos Facebook upload)
Entrega do Trabalho
O trabalho devera ser entrege nas seguintes partes:
- Relatorio
- Data de entrega: 23Mar18 até a meia-noite
- Secoes Obrigatorias:
- Introducao
- Descricao clara do problema
- Descricao da implementacao
- Mais importante: Descricao do processo de desenvolvimento
- Experimentos, discussao, analise dos resultados e estudo detalhado dos casos de falha.
- Conclusao e ideias para trabalho futuro
- Cada membro do grupo deverá também escrever uma seção de, no mínimo, 2 páginas relatando sua contribuição e experiencia com o trabalho
- Codigo
- Data de entrega: terça dia 20Mar18 até a meia-noite
- Dividir a entrega do codigo em 2 partes:
- Codigo-fonte
- README devera incluir instrucoes detalhadas de como rodar o codigo
- Calcular a taxa de acerto do algoritmo em validacao cruzada
- Codigo-fonte
Aviso:
- Criterio de avaliacao: A nota de trabalho sera dada com base na qualidade do relatorio
- Bonus #1: O grupo que tiver melhor performance obtera nota 10.0 no trabalho.
- A performance para bonus sera medida na presenca do professor, durante processo de correcao em horario de aula no dia 22Mar18
- em taxa de reconhecimento em uma base elaborada pelo professor durante a correcao, sem o conhecimento previo do aluno.
- Observacao: Sem ter ao menos feito a P2 minimamente nao sera possivel obter qualquer bonus.
- Àqueles que fizerem P2: a prova P2 sera em grande parte baseada neste trabalho.
Projeto: Análise de Redes Sociais
Neste projeto o aluno coletara uma rede social, por exemplo a rede de amigos do facebook, ou uma rede de troca de mensagens de email, e fara uma análise dos nós mais importantes em relação a PageRank e outras características de difusão na topologia do grafo. Isso poderia ser util, por exemplo, para ordenar a sua lista de amigos por ordem de amigo mais "influente", ou para identificar maneiras de aumentar sua conexao na rede (por exemplo, numa rede de negocios).
Os passos do projeto sao:
- Obter sua rede social no formato GML
- Redes de troca de emails: http://hera.ethymos.com.br:1080/redes/python/autoRede/escolheRedes.php
- Sua rede amigos do facebook: http://snacourse.com/getnet/
- Visualizar a rede no software Gephi
- Baixe o Gephi em https://gephi.org
- Realize o tutorial basico do Gephi (do menu Help)
- Abrir o arquivo GML no Gephi
- Obter a matriz de Markov dessa rede no Scilab, considerando pesos nas arestas como no algoritmo PageRank
- Exportar a rede como matriz de adjacencias no Gephi (aprender como fazer)
- Calcular o autovetor dominante (PageRank) no Scilab
- Interpretar o que isso pode vir a dizer sobre a importância relativa dos seus amigos na rede de amizades
- Por exemplo, a pessoa tem muitos amigos que tem muitos amigos
- O que um caminhante aleatório ou quer dizer aqui? Discuta o que um processo de difusao significaria na rede de amizades. Talvez uma rede de indicacoes?
- Invente alguma outra medida nessa rede da importancia relativa dos nos
Links interessantes
- veja nossas aulas sobre grafos e matrizes
- labmacambira.sf.net/redes
- Análise de redes sociais
Conteudo aproximado
- Foco: SVD/Singular Value decomposition, auto-valores/auto-vetores, PCA/Principal Component Analysis, demais decomposicoes matriciais.
Recursos principais
- Livro-texto: nao ha livro texto. mas veja a bibliografia do curso de Analise Matricial da Pos:
- Livro classico: "Matrix Computations", Gene Golub et. al. (inventor do algoritimo mais usado de SVD).
- Outro livro: "Numerical Matrix Analysis", Ilse Ipsen
- Veja tambem os livros do Gilbert Strang.
- Notas sobre sistemas de coordenadas, matriz relativa a dadas bases, rotacoes (pdf)
Aulas: Listagem parcial
- Primeiras aulas: overview, SVD, PCA, aplicacoes, notacao, bases vs sistemas de coordenadas
- Rotacoes e Sistemas de Coordenadas
- Rotacoes e Sistemas de Coordenadas
- SVD
- Aula Informal projeto e relatorio eigenfaces/PCA
- Aula Informal PCA
- Grafos e matrizes, cadeias de Markov, laplaciano de grafos
- PageRank e outras tecnicas de Ranking, reducao de dimensionalidade, e machine learning. Ver tambem [1]
Para referencia: aulas de Analise Matricial: listagem parcial
- Matrizes, Notacao e Conceitos Elementares (~cap 1 livro)
- Sensitividade, Erros e Normas - parte 1 (cobre cap 2 livro ate normas de matrizes)
Provas
- P1: 2012 **adiada** 27/Out/2016
- P2: 22/Nov/16
- Final-Sub: 29/Nov/16
Criterio de Avaliacao
- Trabalhos: 20% da media - o ultimo trabalho (reconhecimento) valera a metade disto, ou seja, 10% da media
- O criterio final ficou (favor avisar se precisar adicionar detalhes ou corrigir no caso de erro/discrepancia):
Pessoal, conversei com o professor sobre o erro no critério de avaliação, e, como pedido por ele, estou alterando o valor da média para se passar direto.
M_p = (P1 + P2)/2 M = 0.8*M_p + 0.2*T (atualizado de 10% para 20% com acordo dos alunos), onde T é a nota dos trabalhos Se M >= 5, passou --> M (facilitando: considere T=10,0 no M `a esquerda desta desigualdade aqui) prova final - faz quem quiser, mas combinamos que teria de seria feita por quem obtiver M < 5 - facilitando: considere T=10,0 e calcule uma media M otimista, para ver se precisa fazer a final. ou seja:
- Se M_p < 6.25, precisa fazer final! *************************
M_f = 0.5*(M + P_f) = 0.5*(0.8M_p + 0.2*T + P_f) = 0.2*P1 + 0.2*P2 + 0.5*P_f + 0.1*T Se M_f >= 5, passa --> M_f Sub: repoe menor de P1, P2, P_f (apenas se alguem faltou alguma prova ou quiser melhorar nota - mas quem entregar ira substituir) M_sub = media com sub Se M_sub >= 5, passou --> M_sub
Adendo (em acordo com os alunos): a M_sub = M_f pois sera considerada a mesma prova. Quem for usar a prova como Sub ira substituir a nota independentemente do resultado.