Suponha que o sistema linear a seguir seja dado que é o Laplaciano ponderado conhecido como positivo definido, com um espaço nulo unidimensional abrangido por , e a variação de conversão de , ou seja, não altera o valor da função (cuja derivada é ). As únicas entradas positivas de estão na...
É mostrado (Yousef Saad, métodos iterativos para sistemas lineares esparsos , p. 260) que c o n d( A′Um ) ≈ c o n d( A )2cond(A′A)≈cond(A)2cond(A'A) \approx cond(A)^2 Isso também é verdade para ?A A′AA′AA' No caso de ser com , observo queUMAAAN× MN×MN\times MN≪ MN≪MN \ll Mc o n d( A′Um ) » c o n...
Dada uma matriz esparsa genérica com m << n (correção: m ≪ n 2 ) elementos diferentes de zero (normalmente m ∈ O ( n ) ). A é genérico no sentido de que não possui propriedades específicas (por exemplo, definição positiva), e nenhuma estrutura (por exemplo, bandagem) é assumida.A∈Rn×nA∈Rn×nA...
Eu sei que, para resolver o problema simétrico de autovalor , podemos usar a Lei de Inércia de Sylvester, que é o número de autovalores de A menor que a igual ao número de entradas negativas de D de onde a matriz diagonal D vem da fatoração LDL de um - um I = G D G T . Então, pelo método de...
VVV∥⋅∥‖⋅‖\|\cdot\|F:V→RF:V→RF : V \rightarrow \mathbb R Eu gostaria de estimar a norma de (acima e abaixo). Como é uma caixa preta, a única maneira de fazer isso é testá-lo com vetores unitários de e, com base no resultado, encontre que maximize.FFFFFFVVVv∈S1Vv∈S1Vv \in S^1...
Estou pesquisando a estrutura dos complementos de Schur e encontro um fenômeno interessante: Suponha que A seja de 5 - pt laplaciano. Se eu usar a ordenação de dissecações aninhadas e o método multifrontal para calcular a fatoração da LU e depois verificar o último bloco do complemento de schur,...
Suponha que eu tenha o sistema linear grande e esparso original: . Agora, eu não tenho pois A é muito grande para fatorar ou qualquer tipo de decomposição de , mas suponha que eu tenha a solução encontrada com uma solução iterativa.A - 1 A x 0A x0 0= b0 0Ax0=b0A\textbf{x}_0=\textbf{b}_0UMA-...
Gostaria de prever tempos de execução para operações de álgebra linear densa em uma arquitetura específica usando uma biblioteca específica. Eu gostaria de aprender um modelo que se aproxima da função Fo p: :Fop::F_{op} \;::\; tamanhos de entrada tempo de execução→→ \rightarrow para operações...
A partir da definição do número da condição, parece que é necessária uma inversão da matriz para calculá-la. Será que para uma matriz quadrada genérica (ou melhor, se a definição positiva simétrica positiva) é possível explorar alguma decomposição da matriz para calcular o número da condição em um...
Eu tenho escrito uma caixa de ferramentas do sistema de controle a partir do zero e puramente em Python3 (plug descarado:) harold. Nas minhas pesquisas anteriores, sempre tenho reclamações sobre o solucionador de Riccati care.mpor motivos técnicos / irrelevantes. Por isso, tenho escrito meu...
Como aproximar o número de condição de uma matriz grande GGG , se GGG é uma combinação de transformadas de Fourier (não uniforme ou uniforme), diferenças finitas e matrizes diagonais ?FFFRRRSSS As matrizes são muito grandes e não são armazenadas na memória e estão disponíveis apenas como...
Estou muito interessado em otimizar a solução do sistema linear para matrizes pequenas (10x10), às vezes chamadas de matrizes minúsculas . Existe uma solução pronta para isso? A matriz pode ser considerada não singular. Esse solucionador deve ser executado em excesso de 1 000 000 vezes em...
A rotina QR do LAPACK armazena Q como refletores domésticos. Ele escala o vetor de reflexão com , para que o primeiro elemento do resultado se torne , para que não precise ser armazenado. E ele armazena um vetor separado , que contém os fatores de escala necessários. Portanto, uma matriz refletora...
Muitas bibliotecas adaptáveis do FEM usam estruturas de dados de malha mais avançadas para lidar com a adição / remoção de nós, arestas, triângulos, tetraedros, etc. Por exemplo, a biblioteca p4est usa estruturas de dados octree para refinamento de malha adaptável; você normalmente não encontra...
Eu tenho um projeto em que preciso usar campos quadráticos Especificamente, números do formulário com .a+b−3−−−√a+b−3a + b \sqrt{-3}a , b ∈Qa,b∈Qa,b \in \mathbb{Q} Por exemplo, aqui estão os números primos nos números inteiros de Eisenstein : Eu não quero usar sálvia. Eu gostaria de escrever...
Estou interessado em calcular a solução de um sistema Lage de ODEs usando um método de krylov como em [1]. Esse método envolve funções relacionadas ao exponencial (as chamadas funções ). Consiste essencialmente em calcular a ação da função de matriz, construindo um subespaço de Krylov usando a...
Li essa resposta e percebi que tenho usado a diferença entre iterações sucessivas para definir um critério de parada para um método iterativo de encontrar autovalores / vetores. Quais são os bons critérios de parada para métodos iterativos que convergem para autovalores e...
Eu tenho um código de modelagem baseado em Kalman Filter que desenvolvi para um aplicativo de mapeamento ionosférico regional quase em tempo real. O código assimila dados de diferentes sensores em um mapa (descrito por um conjunto de funções básicas) usando um filtro Kalman. Estou tentando escalar...
Nos meus cursos de análise numérica, fui ensinado que a principal e principal motivação para pré-condicionar sistemas lineares de equações é aumentar a taxa de convergência de solucionadores iterativos para essa LSE. Mas, existe algum efeito na precisão da solução computada? Lembro-me de um...
Eu tenho um problema generalizado de matriz: de um método espectral em um problema de análise de estabilidade linear. Minha matriz B é diagonal e positiva semi-definida. A é não-eremita e complexo.Ax=λBxAx=λBxA x = \lambda B x Meu problema é essencialmente que, ao usar o solucionador de autovalor...