Вопросы с тегом «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 …

36 linear-algebra  performance  matlab 

Мне нужно численно оценить интеграл ниже: ∫∞0с я н с'(xr)rE(r)−−−−√dr∫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})Икс∈R+x∈R+x \in \mathbb{R}_+λ,κ,ν> 0λ,κ,ν>0\lambda, \kappa, \nu >0КKKλ = 0,00313λ=0.00313\lambda = 0.00313κ = 0,00825κ=0.00825\kappa …

25 matlab  quadrature  special-functions 

Этот вопрос был перенесен из Cross Validated, потому что на него можно получить ответ на бирже вычислительных стеков. Мигрировал 7 лет назад . Я хочу изучать программирование на MATLAB, чтобы я мог проводить некоторые исследования / анализ самостоятельно, а также, чтобы я мог изучать / модифицировать некоторые скрипты MATLAB, которые …

23 software  matlab  languages  octave 

Я хотел бы знать, есть ли быстрый способ вычисления евклидова расстояния двух векторов в октаве. Кажется, что для этого нет специальной функции, так что я должен просто использовать формулу с 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 

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

15 matlab  ode  numerical-modelling 

Вычислительные люди науки: Первоначально я разместил этот вопрос на Math Stack Exchange, и кто-то сказал, что я мог бы получить «гораздо лучшие» ответы здесь: Я новичок в численных методах и 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 

Привет и спасибо, что нашли время взглянуть на мой вопрос. Это обновленная версия моего вопроса, который я ранее разместил на сайте физика.stackexchange.com. В настоящее время я изучаю двумерный экситонный спинор Бозе-Эйнштейна, и мне интересно узнать основное состояние этой системы. Математический метод перехода в основное состояние называется методом мнимого времени . …

12 matlab  ode  quantum-mechanics 

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

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

Я в настоящее время пытаюсь дешево вычислить оценку хорошего ранга для матрицы . Поэтому я вычисляю разложение QR-кода с помощью Columnt, используяAAA [Q,R,E]=qr(A) в Matlab. Я оцениваю ранг используяAAA 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, и я собираюсь сделать это путем анализа плотности молекулы, а также лапласиана. Я заинтересован в создании двумерных изображений различий в моем приближении к …

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