Domande sugli aspetti algoritmici / computazionali dell'algebra lineare, compresa la soluzione di sistemi lineari, problemi dei minimi quadrati, problemi di autovalori e altre questioni simili.
Ho una densa matrice quadrata simmetrica reale. La dimensione è di circa 1000x1000. Devo calcolare il primo componente principale e chiedermi quale potrebbe essere l'algoritmo migliore per farlo. Sembra che MATLAB usi gli algoritmi Arnoldi / Lanczos (per eigs). Ma leggendo su di loro non sono sicuro che abbiano qualche …
A quanto ho capito, il successivo rilassamento si basa sulla scelta di un parametro 0 ≤ ω ≤ 20≤ω≤20\leq\omega\leq2 e sull'utilizzo di una combinazione lineare di un'iterazione (quasi) di Gauss-Seidel e del valore al precedente timestep ... uk + 1= ( ω ) ugSk + 1+ ( 1 - ω …
Sto cercando di risolvere un sistema lineare di equazioni che è scarso, ma privo di qualsiasi tipo di struttura a bande. Ho sentito che esiste un modo per estendere i principi di un risolutore multigrid per schemi impliciti di differenza finita a un problema lineare generale (se non sbaglio, si …
Ho un sistema lineare con matrice i cui autovalori sono distribuiti uniformemente sul cerchio unitario in questo modo: È possibile risolvere questo tipo di sistema in modo efficace con il metodo iterativo, magari con qualche precondizionatore?
Uso la funzione kernel RBF per implementare un algoritmo di machine learning (KLPP) basato sul kernel, la risultante matrice del kernel KKK K( i , j ) = exp( - ( xio- xj)2σ2m)K(i,j)=exp(−(xi−xj)2σm2)K(i,j)= \exp\left({\frac{-(x_{i}-x_{j})^2}{ \sigma_{m}^2}}\right) risulta estremamente mal condizionato. Il numero di condizione della norma L2 è di1017- 10641017−106410^{17}-10^{64} C'è …
La routine QR di LAPACK memorizza Q come riflettori Householder. Ridimensiona il vettore di riflessione vvv con 1 / v11/v11/v_1 , quindi il primo elemento del risultato diventa 111 , quindi non deve essere memorizzato. E memorizza un vettore ττ\tau separato , che contiene i fattori di scala necessari. Quindi …
Sono molto interessato all'ottimizzazione della risoluzione del sistema lineare per piccole matrici (10x10), a volte chiamate minuscole matrici. C'è una soluzione pronta per questo? La matrice può essere considerata non singolare. Questo solutore deve essere eseguito oltre 1 000 000 volte in microsecondi su una CPU Intel. Sto parlando del …
Ho scritto una toolbox del sistema di controllo da zero e puramente in Python3 (plug spudorato:) harold. Dalla mia ricerca passata, ho sempre lamentato il solutore Riccati care.mper ragioni tecniche / irrilevanti. Quindi, ho scritto il mio set di routine. Una cosa che non riesco a trovare è quella di …
Non riesco a immaginare di essere il primo a pensare al seguente problema, quindi sarò soddisfatto di un riferimento (ma una risposta completa e dettagliata è sempre apprezzata): Supponi di avere un simmetrico positivo definito . n è considerato molto grande, quindi è impossibile tenere Σ in memoria. È possibile, …
Sono interessato a calcolare la soluzione di un sistema di lage di ODE usando un metodo krylov come in [1]. Tale metodo prevede funzioni correlate all'esponenziale (le cosiddette funzioni ). Consiste essenzialmente nel calcolare l'azione della funzione matrice costruendo un sottospazio di Krylov usando l'iterazione di Arnoldi e proiettando la …
Un sacco di librerie FEM adattivi utilizzare strutture dati in rete più avanzate per gestire l'aggiunta / rimozione di nodi, bordi, triangoli, tetraedri, ecc, ad esempio, i p4est utilizza la libreria octree strutture dati per adaptive mesh refinement; non troverai spesso ocre utilizzate per i calcoli su una mesh statica. …
Ho un progetto in cui devo usare i campi quadratici In particolare i numeri del modulo con . a,b∈Qa + b - 3---√un'+B-3a + b \sqrt{-3}a , b ∈ Qun',B∈Qa,b \in \mathbb{Q} Ad esempio, ecco i numeri primi negli interi di Eisenstein : Non voglio usare la salvia. Vorrei scrivere …
Come approssimo il numero di condizione di una matrice grande GGG , se GGG è una combinazione di trasformate di Fourier FFF (non uniforme o uniforme), differenze finite RRR e matrici diagonali SSS ? Le matrici sono molto grandi e non memorizzate e sono disponibili solo come funzioni. In particolare, …
Dalla definizione del numero di condizione sembra che sia necessaria un'inversione di matrice per calcolarla, mi chiedo se per una matrice quadrata generica (o meglio se definita positiva simmetrica) è possibile sfruttare una decomposizione della matrice per calcolare il numero di condizione in un modo più veloce.
Vorrei prevedere i tempi di autonomia per operazioni di algebra lineare densa su un'architettura specifica utilizzando una libreria specifica. Vorrei imparare un modello che approssima la funzione Fop::Fop::F_{op} \;::\; dimensioni di input runtime→→ \rightarrow per operazioni come moltiplicare matrice, aggiungere elementi, risolvere triangolare, ecc. Ho il sospetto che questi runtime …
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.