Costruzione e analisi di algoritmi per calcolare soluzioni discrete approssimative di problemi continui. Un esempio canonico è l'approssimazione di derivati tramite quozienti di differenza.
Ho diversi problemi di ottimizzazione globale non convessi da risolvere. Attualmente uso MATLAB's Optimization Toolbox (in particolare, fmincon()con algoritmo = 'sqp'), che è abbastanza efficace . Tuttavia, la maggior parte del mio codice è in Python e mi piacerebbe fare l'ottimizzazione anche in Python. Esiste un solutore NLP con attacchi …
Attualmente sto esaminando metodi paralleli per l'integrazione ODE. C'è molta letteratura nuova e vecchia là fuori che descrive una vasta gamma di approcci, ma non ho trovato sondaggi recenti o articoli di panoramica che descrivono l'argomento in generale. C'è il libro di Burrage [1], ma ha quasi 20 anni e …
Ho imparato a conoscere il metodo degli elementi finiti (anche un po 'su altri metodi numerici) ma non so quali siano esattamente le definizioni di questi due errori e le differenze tra loro?
Vorrei sapere se esiste un modo rapido per calcolare la distanza euclidea di due vettori in Ottava. Sembra che non ci sia una funzione speciale per questo, quindi dovrei semplicemente usare la formula con sqrt?
Un sacco di intuizione può essere acquisita dall'esperienza, mi chiedevo solo se qualcuno ha visto qualcosa di simile prima. Il diagramma mostra la condizione iniziale (verde) per l'equazione di avanzamento-diffusione, quindi la soluzione all'iterazione 200 (blu) e quindi nuovamente all'iterazione 400 (rosso). La soluzione dell'equazione avanzamento-diffusione esplode dopo alcune iterazioni. …
A seguito della mia domanda precedente, sto cercando di applicare le condizioni al contorno a questa mesh di volume finito non uniforme, Vorrei applicare una condizione al contorno di tipo Robin ai lhs del dominio ( x=xL)x=xL)x=x_L) , in modo tale che, σL=(dux+au)∣∣∣x=xLσL=(dux+au)|x=xL \sigma_L = \left( d u_x + a …
Faccio regolarmente parte dei cosiddetti "Concorsi di programmazione", in cui risolvi difficili problemi algoritmici con il tuo codice e le tue capacità di risoluzione dei problemi in un arco di tempo limitato. Per esempi referenziali di come potrebbero apparire, cerca concorsi come ad esempio Google Code Jam o ACM-ICPC. (Se …
Mi sono imbattuto in un'osservazione sconcertante sul giornale PJ van der Houwen, Lo sviluppo di metodi di Runge-Kutta per equazioni differenziali parziali, Appl. Num. Matematica. 20: 261, 1996 Alle righe 8 e seguenti a pagina 264, van der Houwen scrive: "Per i polinomi di Taylor questo implica che l'intervallo di …
Un problema comune nelle statistiche è il calcolo dell'inverso della radice quadrata di una matrice definita positiva simmetrica. Quale sarebbe il modo più efficiente per calcolare questo? Mi sono imbattuto in alcuni articoli letterari (che non ho ancora letto) e in alcuni codici R casuali qui , che riprodurrò qui …
Sto cercando di capire alcuni risultati e apprezzerei alcuni commenti generali sulla lotta ai problemi non lineari. Equazione di Fisher (una PDE non lineare di reazione-diffusione), ut= dux x+ βu ( 1 - u ) = F( u )ut=duxx+βu(1−u)=F(u) u_t = du_{xx} + \beta u (1 - u) = F(u) …
In pratica, l'autonomia di risolvere numericamente un IVP x˙(t)=f(t,x(t)) for t∈[t0,t1]X˙(t)=f(t,X(t)) per t∈[t0,t1] \dot{x}(t) = f(t, x(t)) \quad \text{ for } t \in [t_0, t_1] x(t0)=x0X(t0)=X0 x(t_0) = x_0 è spesso dominato dalla durata della valutazione del lato destro (RHS). Supponiamo quindi che tutte le altre operazioni siano istantanee (cioè …
Sto riscontrando problemi nell'implementazione numerica di una funzione. Soffre del fatto che a valori di input elevati il risultato è un numero molto grande per un numero molto piccolo. Non sono sicuro che la cancellazione catastrofica sia il termine corretto, quindi correggimi se lo è. Prova di qualcosa che non …
Sto lavorando per risolvere le equazioni di poroelasticità monodimensionali accoppiate (modello di biot), dato come: ∂−(λ+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} = 0 ∂∂t[γp+∂u∂x]−κη[∂2p∂x2]=q(x,t)∂∂t[γp+∂u∂x]−κη[∂2p∂x2]=q(x,t)\frac{\partial}{\partial t} \left[ \gamma p + \frac{\partial u}{\partial x}\right] -\frac{\kappa}{\eta}\left[\frac{\partial^2 p}{\partial x^2}\right] =q(x,t) sul dominio e con le condizioni al contorno: Ω=(0,1)Ω=(0,1)\Omega=(0,1) p=0,(λ+2μ)∂u∂x=−u0p=0,(λ+2μ)∂u∂x=−u0p=0, (\lambda …
Let e essere una funzione in queste variabili.f( → x ):[0,1]n→CX⃗ = ( 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} Esiste uno schema ricorsivo per questo integrale …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.