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...
Para o meu trabalho de pesquisa, desenvolvi uma pequena biblioteca C ++ destinada a facilitar a comunicação entre os códigos computacionais C ++ e o Octave / Matlab (quando este último é usado para fins de pós-processamento). Gostaria de liberar essa biblioteca sob a GPL em alguns sites de...
Estou tentando calcular momentos Zernike de ordem superior (por exemplo m=0, n=46) para alguma imagem. No entanto, estou com um problema em relação ao polinômio radial (consulte a Wikipedia ). Este é um polinômio definido no intervalo [0 1]. Veja o código MATLAB abaixo function R =...
Estou usando a caixa de ferramentas Matlab pde para resolver uma certa equação elíptica em 2D. A solução está correta, embora eu precise plotá-la ao longo de uma determinada linha, ou seja, cortar uma fatia plana da malha 3D que representa a solução. Não consigo descobrir uma maneira que envolva...
Estou escrevendo um pequeno pacote MATLAB que resolverá certos problemas numéricos de classe. Existem 3 estágios do algoritmo e o usuário tem 5 opções para cada estágio. Eu implementei todo o problema usando funções e 3 casos de switch (um para cada estágio do algoritmo). Funciona bem, mas estou...
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...
O Answer fornece um software para calcular o expoente Lyapunov condicional (CLE) para osciladores acoplados na sincronização do caos. No entanto, é difícil de seguir e não há saída gráfica do gráfico (e é em C, mais complexo). Alguém sabe como modificar a LET Toolbox, que é muito boa para o sistema...
Parece que o Matlab eigsestá me dando más aproximações dos menores vetores próprios de uma matriz. Suponho que posso usar alguns métodos mais lentos que também seriam mais precisos ... Estou procurando encontrar o segundo menor vetor próprio de uma matriz lapalciana (conhecido como vetor...
Existe alguma biblioteca / caixa de ferramentas que tenha implementação de SVD incremental no MATLAB. Eu implementei este documento , é rápido, mas não funciona bem. Eu tentei isso, mas neste erro também se propaga rapidamente (dentro de atualização de 5 a 10 pontos, o erro é...
Eu tenho uma matriz esparsa e simétrica A (nxn). O método Lanczos transforma a matriz A em matriz tridiagonal e simétrica T e os vetores Lanczos na matriz V. A partir daí, como computo k autovalores menores ou maiores e autovetores
Estou tentando modelar a condução de calor dentro de um cilindro de madeira usando métodos implícitos de diferenças finitas. A equação geral de calor que estou usando para formas cilíndricas e esféricas é: Onde p é o fator de forma, p = 1 para cilindro ep = 2 para esfera. As condições de...
Suponha que você tenha um sistema de PDEs para resolver. Pelo menos por simplicidade, vamos supor que seja independente do tempo, quase linear (linear em suas derivadas) resolvido em uma grade retangular no espaço (x, y) e com condições de contorno especificadas ao redor. Minha pergunta é mais...
O resultado em que estou interessado encontra-se na seção "Sincronização: um conceito universal em ciências não lineares" na página figura 14.3 . O fragmento peculiar também é fornecido no final deste post.33333333314,314,314.3 Portanto, basicamente existe esse acoplamento dissipativo aplicado a...
Quero começar a programação científica no Mac usando o Objective-C. Estou familiarizado com o MATLAB, que facilita o armazenamento de dados de formas de onda com valores complexos em vetores e geração de gráficos. Fiquei me perguntando se poderia haver um bom recurso que pudesse comparar como você...
Eu preciso calcular apenas um pequeno número de componentes de Fourier de baixa frequência de uma matriz bidimensional complexa. Estarei computando os mesmos componentes de Fourier repetidamente conforme a matriz de entrada muda. Claramente, no limite em que eu quero apenas um componente de...