Perguntas com a marcação «numerical-analysis»

Eu tenho vários problemas desafiadores de otimização global não convexa para resolver. Atualmente, uso o Optimization Toolbox do MATLAB (especificamente, fmincon()com o algoritmo = 'sqp'), o que é bastante eficaz . No entanto, a maior parte do meu código está em Python, e eu adoraria fazer a...

Atualmente, estou procurando métodos paralelos para integração com ODE. Há muita literatura nova e antiga por aí descrevendo uma ampla variedade de abordagens, mas não encontrei nenhuma pesquisa recente ou artigo de visão geral descrevendo o tópico em geral. Existe o livro de Burrage [1], mas tem...

Qual é o estado da arte na aproximação de integrais altamente oscilatórias em uma dimensão e dimensões superiores à precisão

Aprendi sobre o método dos elementos finitos (também um pouco sobre outros métodos numéricos), mas não sei exatamente qual é a definição desses dois erros e as diferenças entre

Seguindo a minha pergunta anterior , estou tentando aplicar condições de contorno a essa malha de volume finito não uniforme, Gostaria de aplicar uma condição de limite do tipo Robin aos lhs do domínio ( x = xeu)x=xeu)x=x_L) , de modo que, σeu= ( dvocêx+ a u ) ∣∣∣x =...

Pode-se obter uma grande quantidade de insights sobre a experiência, apenas queria saber se alguém já viu algo semelhante a isso antes. O gráfico mostra a condição inicial (verde) para a equação de advecção-difusão, depois a solução na iteração 200 (azul) e depois novamente na iteração 400...

Gostaria de saber se existe uma maneira rápida de calcular a distância euclidiana de dois vetores na oitava. Parece que não há função especial para isso, então devo usar apenas a fórmula com

Estou tentando entender alguns resultados e gostaria de receber alguns comentários gerais sobre como resolver problemas não-lineares. Equação de Fisher (um PDE de difusão da reação não linear), vocêt= dvocêx x+ βu ( 1 - u ) = F( U )vocêt=dvocêxx+βvocê(1 1-você)=F(você) u_t = du_{xx} + \beta u...

Me deparei com uma observação intrigante no jornal PJ van der Houwen, O desenvolvimento de métodos de Runge-Kutta para equações diferenciais parciais, Appl. Num. Matemática. 20: 261, 1996 Nas linhas 8ff na página 264, van der Houwen escreve: "Para os polinômios de Taylor, isso implica que o...

Eu concordo regularmente nos chamados "Concursos de programação", em que você resolve problemas algorítmicos difíceis com seu próprio código e habilidades para resolver problemas durante um período de tempo limitado. Para obter exemplos de referência de como isso pode ser, pesquise competições...

Um problema comum em estatística é calcular a raiz quadrada inversa de uma matriz definida positiva simétrica. Qual seria a maneira mais eficiente de calcular isso? Encontrei alguma literatura (que ainda não li) e algum código R incidental aqui , que reproduzirei aqui por conveniência # function...

Na prática, o tempo de execução de resolver numericamente um IVP é frequentemente dominada pela duração da avaliação do lado direito (RHS) . Vamos, portanto, assumir que todas as outras operações são instantâneas (ou seja, sem custo computacional). Se o tempo de execução geral para resolver o IVP...

Estou tendo problemas para implementar uma função numericamente. Ele sofre com o fato de que, com valores de entrada grandes, o resultado é um número muito grande vezes um número muito pequeno. Não tenho certeza se o cancelamento catastrófico é o termo correto, por favor, corrija-me se for....

Estou trabalhando na resolução das equações de poroelasticidade unidimensional acopladas (modelo de biot), dadas como: −(λ+2μ)∂2u∂x2+∂p∂x=0−(λ+2μ)∂2u∂x2+∂p∂x=0-(\lambda+ 2\mu) \frac{\partial^2 u}{\partial x^2} + \frac{\partial p}{\partial x} =

Vamos e ser uma função nessas variáveis.x⃗ =(x1,x2,…,xn)∈[0,1]nx→=(x1,x2,…,xn)∈[0,1]n\vec{x} = (x_1, x_2, \dots, x_n) \in [0,1]^nf(x⃗ ):[0,1]n→Cf(x→):[0,1]n→Cf(\vec{x}): [0,1]^n \to \mathbb{C} Existe um esquema recursivo para essa integral iterada? ∫[0,1]n∏dxif(x⃗ )∫[0,1]n∏dxif(x→)\int_{[0,1]^n}...

Quando o FEM discretiza e resolve um problema de difusão da reação, por exemplo, com 0 < ε « 1 (perturbação singular), a solução do problema discreta irá tipicamente exibem camadas oscilatórios próximo da fronteira. Com Ω = ( 0 , 1 ) , ε = 10 - 5 e elementos finitos lineares, a solução u h se...

Estou tentando inverter numericamente a seguinte transformação integral: F( y) = ∫∞0 0yexp[ - 12( y2+ x2) ] I0 0( x y) f( X )d xF(y)=∫0∞yexp⁡[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Portanto, para um dado F(...

Começando com a equação de advecção em forma de conservação. ut=(a(x)u)xut=(a(x)u)x u_t = (a(x)u)_x onde é uma velocidade que depende do espaço, e é uma concentração de uma espécie que é conservada.ua(x)a(x)a(x)uuu Discretizar o fluxo (onde o fluxo , é definido nas bordas das células entre os...

No livro de Nocedal & Wright sobre otimização numérica, há uma declaração na seção 2.2 (página 27): "De um modo geral, é mais fácil preservar a invariância de escala para algoritmos de pesquisa de linha do que para algoritmos de região confiável". Na mesma seção, eles falam sobre ter novas...

Como posso substituir o método de Euler pela 4ª ordem de Runge-Kutta para determinar o movimento de queda livre em magnitude gravimétrica não constante (por exemplo, queda livre de 10.000 km acima do solo)? Até agora escrevi uma integração simples pelo método Euler: while() { v +=...