Stavo confrontando alcuni dei miei codici con codici MATLAB "stock". Sono sorpreso dai risultati. Ho eseguito un codice di esempio (Sparse Matrix) n = 5000; a = diag(rand(n,1)); b = rand(n,1); disp('For a\b'); tic;a\b;toc; disp('For LU'); tic;LULU;toc; disp('For Conj Grad'); tic;conjgrad(a,b,1e-8);toc; disp('Inv(A)*B'); tic;inv(a)*b;toc; Risultati: For a\b Elapsed time is 0.052838 …
Devo valutare numericamente l'integrale di seguito: ∫∞0sinc′(xr)rE(r)−−−−√dr∫0∞sinc′(xr)rE(r)dr\int_0^\infty \mathrm{sinc}'(xr) r \sqrt{E(r)} dr dove ,x∈R+eλ,κ,ν>0. QuiKè la funzione di Bessel modificata del secondo tipo. Nel mio caso particolare hoλ=0,00313,κ=0,00825eν=0,33.E(r)=r4(λκ2+r2−−−−−−√)−ν−5/2K−ν−5/2(λκ2+r2−−−−−−√)E(r)=r4(λκ2+r2)−ν−5/2K−ν−5/2(λκ2+r2)E(r) = r^4 (\lambda\sqrt{\kappa^2+r^2})^{-\nu-5/2} K_{-\nu-5/2}(\lambda\sqrt{\kappa^2+r^2})x∈R+x∈R+x \in \mathbb{R}_+λ,κ,ν>0λ,κ,ν>0\lambda, \kappa, \nu >0KKKλ =0.00313λ=0,00,313 mila\lambda = 0.00313κ = 0,00825κ=0,00,825 mila\kappa = 0.00825ν= 0,33ν=0,33\nu = 0.33 Sto usando …
Questa domanda è stata migrata da Cross Validated perché è possibile rispondere a Exchange Stack di calcolo scientifico. Migrato 7 anni fa . Voglio imparare la programmazione MATLAB in modo da poter condurre alcune ricerche / analisi da solo e anche, in modo da poter studiare / modificare alcuni script …
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?
Vorrei simulare il comportamento di un sistema a doppio pendolo. Il sistema è un manipolatore robotico a 2 gradi di libertà che non viene attivato e, pertanto, si comporterà principalmente come un doppio pendolo influenzato dalla gravità. L'unica differenza principale con un doppio pendolo è che è composto da due …
Persone di scienza computazionale: Originariamente ho pubblicato questa domanda a Math Stack Exchange e qualcuno ha commentato che potrei ottenere risposte "molto migliori" qui: Sono alle prime armi con metodi numerici e Matlab. Sto tentando di valutare la seguente somma di due integrali tripli (ovviamente può essere scritta più semplicemente, …
Sto vivendo una certa frustrazione per il modo in cui matlab gestisce l'integrazione numerica rispetto a Scipy. Osservo le seguenti differenze nel mio codice di prova di seguito: La versione di Matlab funziona in media 24 volte più veloce del mio equivalente di Python! La versione di Matlab è in …
Ho letto che se, ad esempio, ho un doppio forciclo che corre sopra gli indici di una matrice, quindi mettere la colonna che esegue l'indice nel ciclo esterno è più efficiente. Per esempio: a=zeros(1000); for j=1:1000 for i=1:1000 a(i,j)=1; end end Qual è il modo più efficiente per codificarlo se …
Ciao e grazie per aver dato il tempo di guardare la mia domanda. Questa è una versione aggiornata della mia domanda che ho pubblicato in precedenza in physics.stackexchange.com Attualmente sto studiando uno spinatore di eccitoni 2D Bose-Einstein Condensate e sono curioso di sapere lo stato fondamentale di questo sistema. Il …
Sto usando MATLAB per risolvere un problema che comporta la risoluzione di in ogni momento, dove b cambia nel tempo. In questo momento, sto realizzando questo usando MATLAB :A x = bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}Bb\mathbf{b}mldivide x = A\b Ho la flessibilità di effettuare tutti i precomputamenti necessari, quindi mi chiedo se esiste …
Io attualmente cercando di calcolare una stima a buon mercato buon grado per una matrice . Pertanto computo un columnt che ruota la decomposizione QR usandoUNAA [Q,R,E]=qr(A) in Matlab. Stimo il rango di usandoUNAA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Funziona bene e una trama su tutte le voci diagonali …
Sto seguendo un corso di calcolo scientifico e abbiamo appena superato l'approssimazione dei minimi quadrati. La mia domanda riguarda in particolare l'approssimazione con i polinomi. Capisco che se hai n + 1 punti dati, puoi trovare un polinomio unico di grado n che descrive tutti questi punti. Ma posso anche …
Sto producendo immagini 1d e 2d usando Matlab in questo momento per il confronto della precisione con un determinato modello. Devo confrontare i miei metodi con il modello standard .wfn gaussiano e lo farò analizzando la densità di una molecola e il laplaciano. Sono interessato a produrre immagini 2D di …
Dato un insieme arbitrario di matrici complesse (numeriche) quadrate , sono interessato a calcolare la vera algebra di Lie della matrice generata da , chiamala . Cioè, vorrei una base per dove è definito in modo ricorsivo come e per .A L A L A = s p a n …
Sto implementando il documento " Trasporto di massa ottimale per la registrazione e la deformazione ", il mio obiettivo è quello di metterlo online poiché non riesco proprio a trovare alcun codice di trasporto di massa euleriano online e questo sarebbe interessante almeno per la comunità di ricerca nell'elaborazione delle …
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.