Запитання з тегом «matlab»

Я порівнював декілька своїх кодів із "запасними" кодами MATLAB. Я здивований результатами. Я запустив зразок коду (розріджена матриця) 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; Результати: For a\b Elapsed time is 0.052838 seconds. For LU Elapsed …

36 linear-algebra  performance  matlab 

Мені потрібно чисельно оцінити інтеграл нижче: ∫∞0s i n c'( x r ) r E( r )----√гr∫0∞sinc′(xr)rE(r)dr\int_0^\infty \mathrm{sinc}'(xr) r \sqrt{E(r)} dr де , і . Тут - модифікована функція Бесселя другого роду. У моєму конкретному випадку я , і .x∈R+λ,κ,ν>0Kλ=0,00313κ=0,00825ν=0,33E(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λ = …

25 matlab  quadrature  special-functions 

Це питання було переміщено з перехресної перевірки, оскільки на нього можна відповісти на обміні стека обчислювальної науки. Мігрували 7 років тому . Я хочу навчитися програмуванню MATLAB, щоб я міг самостійно провести дослідження / аналіз, а також, щоб я міг вивчити / змінити деякі сценарії MATLAB, які я знайшов в …

23 software  matlab  languages  octave 

Мені хотілося б знати, чи існує швидкий спосіб обчислити евклідову відстань двох векторів в Октаві. Здається, що для цього немає спеціальної функції, тому я повинен просто використовувати формулу sqrt?

18 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 

Я хочу імітувати поведінку системи, що нагадує подвійний маятник. Система - це маніпулятор-робот 2 ступеня свободи, який не спрацьовує, і, отже, буде вести себе, як подвійне маятник, постраждалий від сили тяжіння. Єдина основна відмінність з подвійним маятником полягає в тому, що він складається з двох жорстких тіл з властивостями маси …

15 matlab  ode  numerical-modelling 

Люди з обчислювальної науки: Я спочатку опублікував це питання на біржі Math Stack, і хтось прокоментував, що я можу отримати "набагато кращі" відповіді тут: Я новачок у числових методах та Matlab. Я намагаюся оцінити наступну суму двох потрійних інтегралів (це, очевидно, можна записати простіше, але ви все ще не можете …

15 matlab 

Я відчуваю деяке розчарування з приводу того, як matlab обробляє числову інтеграцію проти Scipy. У моєму тестовому коді нижче я спостерігаю такі відмінності: Версія Matlab працює в середньому в 24 рази швидше, ніж мій еквівалент python! Версія Matlab здатна обчислити інтеграл без попереджень, тоді як python повертається nan+nanj Що я …

13 matlab  python  quadrature  numba 

Я читав, що якщо, наприклад, у мене є подвійний forцикл, який працює над індексами матриці, то введення індексу запущеного стовпця у зовнішній цикл є більш ефективним. Наприклад: a=zeros(1000); for j=1:1000 for i=1:1000 a(i,j)=1; end end Який найефективніший спосіб кодування, якщо у мене є три або більше forциклів? Наприклад: a=zeros(100,100,100); for …

12 matlab  efficiency 

Привіт, спасибі, що надали час на розгляд мого питання. Це оновлена ​​версія мого питання, яку я опублікував раніше на сайті physics.stackexchange.com На даний момент я вивчаю 2D спінор екситону Бозе-Ейнштейна конденсату і мені цікаво основний стан цієї системи. Математичний метод потрапляння в основний стан називається методом уявного часу . Метод …

12 matlab  ode  quantum-mechanics 

Я використовую MATLAB для вирішення проблеми, яка включає вирішення на кожному кроці, де b змінюється з часом. Зараз я досягаю цього за допомогою MATLAB :Ax=bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}bb\mathbf{b}mldivide x = A\b У мене є можливість зробити стільки попередніх обчислень, скільки потрібно, тому мені цікаво, чи існує швидший і / або більш точний …

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

Я в даний час намагаюся дешево обчислити оцінку хорошого рангу для матриці . Тому я обчислюю стовпчик, що розгортається QR, використовуючиААA [Q,R,E]=qr(A) в Матлабі. Я оцінюю ранг використовуючиААA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Це добре працює, і сюжет над усіма діагональними записами R виглядає так: Якщо весь алгоритм переносять …

12 matlab  lapack 

Я беру курс на наукові обчислення, і ми просто пройшли наближення до найменших квадратів. Моє запитання стосується конкретного наближення за допомогою многочленів. Я розумію, що якщо у вас n + 1 точок даних, ви можете знайти унікальний многочлен ступеня n, який описує всі ці точки. Але я також бачу, чому …

11 matlab  least-squares 

Я зараз створюю зображення 1d та 2d, використовуючи Matlab для порівняння точності та заданої моделі. Мені потрібно порівняти свої методи зі стандартною моделлю Гаусса .wfn, і я зроблю це шляхом аналізу щільності молекули, а також лаплаціана. Мені цікаво створити 2-D зображення відмінностей у моєму наближенні до наближення .wfn (можливо, за …

11 matlab  computational-chemistry  publications 

З огляду на довільний набір (числові) квадратних комплексних матриць , Я зацікавлений в обчисленні алгебра Лі речової матриці , породжену А , назвемо його Л . Тобто, я хотів би основу для L A = s p a n R { B : B ∈ ∪ ∞ k = 1 …

11 linear-algebra  matlab  basis-set 

Я реалізую статтю " Оптимальний транспорт транспорту для реєстрації та викривлення ", моя мета - розмістити його в Інтернеті, оскільки я просто не можу знайти жодного ейлерового коду масового транспорту в Інтернеті, і це було б цікаво хоча б для дослідницької спільноти в обробці зображень. Папір можна узагальнити так: - …

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