Pertanyaan yang diberi tag «matlab»

Saya membandingkan beberapa kode saya dengan "stok" kode MATLAB. Saya terkejut dengan hasilnya. Saya menjalankan kode sampel (Matriks Jarang) 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; Hasil: For a\b Elapsed time is 0.052838 seconds. For LU …

36 linear-algebra  performance  matlab 

Saya perlu mengevaluasi secara integral angka-angka di bawah ini: ∫∞0sinc′(xr)rE(r)−−−−√dr∫0∞sinc′(xr)rE(r)dr\int_0^\infty \mathrm{sinc}'(xr) r \sqrt{E(r)} dr di mana , dan . Di sini adalah fungsi Bessel yang dimodifikasi dari jenis kedua. Dalam kasus khusus saya, saya memiliki , dan .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.00313\lambda = 0.00313κ=0.00825κ=0.00825\kappa = …

25 matlab  quadrature  special-functions 

Pertanyaan ini dimigrasikan dari Cross Validated karena dapat dijawab di Computational Science Stack Exchange. Bermigrasi 7 tahun yang lalu . Saya ingin belajar pemrograman MATLAB sehingga saya dapat melakukan beberapa penelitian / analisis saya sendiri dan juga, sehingga saya dapat mempelajari / memodifikasi beberapa skrip MATLAB yang telah saya temukan …

23 software  matlab  languages  octave 

Saya ingin tahu apakah ada cara cepat untuk menghitung jarak Euclidean dari dua vektor dalam Oktaf. Sepertinya tidak ada fungsi khusus untuk itu, jadi haruskah saya menggunakan rumus saja sqrt?

16 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

Saya ingin mensimulasikan perilaku sistem seperti pendulum ganda. Sistem ini adalah manipulator robot 2 derajat kebebasan yang tidak digerakkan dan oleh karena itu, akan berperilaku seperti pendulum ganda yang dipengaruhi oleh gravitasi. Satu-satunya perbedaan utama dengan pendulum ganda adalah bahwa ia terdiri dari dua benda kaku dengan massa dan sifat …

15 matlab  ode  numerical-modelling 

Orang-orang Ilmu Komputasi: Saya awalnya memposting pertanyaan ini di Math Stack Exchange dan seseorang berkomentar bahwa saya mungkin mendapatkan jawaban "jauh lebih baik" di sini: Saya seorang pemula dalam metode numerik dan Matlab. Saya mencoba untuk mengevaluasi jumlah berikut dari dua integral tiga (ini jelas dapat ditulis lebih sederhana, tetapi …

15 matlab 

Saya mengalami beberapa frustrasi atas cara matlab menangani integrasi numerik vs Scipy. Saya mengamati perbedaan berikut dalam kode pengujian saya di bawah ini: Versi Matlab berjalan rata-rata 24 kali lebih cepat dari yang setara dengan python saya! Versi Matlab mampu menghitung integral tanpa peringatan, sementara python kembali nan+nanj Apa yang …

13 matlab  python  quadrature  numba 

Saya telah membaca bahwa jika, misalnya, saya memiliki forloop ganda yang berjalan di atas indeks matriks, maka menempatkan indeks berjalan kolom di loop luar lebih efisien. Sebagai contoh: a=zeros(1000); for j=1:1000 for i=1:1000 a(i,j)=1; end end Apa cara paling efisien untuk mengkodekannya jika saya memiliki tiga forloop atau lebih ? …

12 matlab  efficiency 

Hai dan terima kasih telah meluangkan waktu untuk melihat pertanyaan saya. Ini adalah versi terbaru dari pertanyaan saya yang saya posting sebelumnya di physics.stackexchange.com Saat ini saya sedang mempelajari spin exciton 2D Bose-Einstein Condensate dan saya ingin tahu tentang keadaan dasar sistem ini. Metode matematika untuk mencapai keadaan dasar disebut …

12 matlab  ode  quantum-mechanics 

Saya menggunakan MATLAB untuk memecahkan masalah yang melibatkan penyelesaian di setiap timestep, di mana berubah seiring waktu. Saat ini, saya menyelesaikan ini menggunakan MATLAB :bAx=bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}bb\mathbf{b}mldivide x = A\b Saya memiliki fleksibilitas untuk membuat sebanyak mungkin perhitungan yang diperlukan, jadi saya bertanya-tanya apakah ada metode yang lebih cepat dan / …

12 linear-algebra  matlab  linear-solver  matrix  linear-system 

Saat ini saya mencoba untuk murah menghitung perkiraan yang baik peringkat untuk matriks . Oleh karena itu saya menghitung penggunaan QR decompostion columnt pivotingSEBUAHAA [Q,R,E]=qr(A) di Matlab. Saya memperkirakan peringkat menggunakanSEBUAHAA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Ini berfungsi dengan baik dan plot pada semua entri diagonal R terlihat seperti: …

12 matlab  lapack 

Saya mengambil kursus perhitungan ilmiah, dan kami baru saja mendekati pendekatan kuadrat. Pertanyaan saya secara khusus tentang perkiraan menggunakan polinomial. Saya mengerti bahwa jika Anda memiliki titik data n +1, Anda dapat menemukan polinomial tingkat n yang unik yang menggambarkan semua titik ini. Tapi saya juga bisa melihat mengapa ini …

11 matlab  least-squares 

Saya menghasilkan gambar 1d dan 2d menggunakan Matlab sekarang untuk perbandingan akurasi terhadap model yang diberikan. Saya perlu membandingkan metode saya dengan model Gaussian .wfn standar dan saya akan melakukannya dengan menganalisis kepadatan molekul serta laplacian. Saya tertarik untuk menghasilkan gambar 2-D dari perbedaan perkiraan saya dengan perkiraan .wfn (mungkin …

11 matlab  computational-chemistry  publications 

Mengingat set sewenang-wenang (numerik) persegi kompleks matriks , Saya tertarik dalam menghitung matriks nyata aljabar Lie yang dihasilkan oleh A , sebut saja L A . Artinya, saya ingin dasar untuk L A = s p a n R { B : B ∈ ∪ ∞ k = 1 C …

11 linear-algebra  matlab  basis-set 

Saya mengimplementasikan makalah " Transportasi Massal Optimal untuk Registrasi dan Warping ", tujuan saya adalah untuk membuatnya online karena saya tidak dapat menemukan kode transportasi massal euler online dan ini akan menarik setidaknya untuk komunitas riset dalam pemrosesan gambar. Makalah ini dapat diringkas sebagai berikut: - temukan peta awal uuu …

11 pde  finite-difference  matlab  fluid-dynamics  image-processing