«numerical-limitations» etiketlenmiş sorular

İki vektör arasındaki açı için klasik formülü uygularken: α=arccosv1⋅v2∥v1∥∥v2∥α=arccos⁡v1⋅v2‖v1‖‖v2‖\alpha = \arccos \frac{\mathbf{v_1} \cdot \mathbf{v_2}}{\|\mathbf{v_1}\| \|\mathbf{v_2}\|} çok küçük / akut açılar için hassasiyet kaybı olduğu ve sonucun doğru olmadığı bulunmuştur. Bu Yığın Taşması cevabında açıklandığı gibi , çözümlerden biri arktanjant kullanmaktır: α=arctan2(∥v1×v2∥,v1⋅v2)α=arctan⁡2(‖v1×v2‖,v1⋅v2)\alpha = \arctan2 \left(\|\mathbf{v_1} \times \mathbf{v_2}\|, \mathbf{v_1} \cdot \mathbf{v_2} \right) …

14 algorithms  precision  numerical-limitations 

Intel MKL'nin SVD'sini ( dgesvdSciPy aracılığıyla) kullanıyorum ve matrisim kötü koşullu / tam sıralı değilken hassasiyeti değiştirdiğimde float32ve sonuçların önemli ölçüde farklı olduğunu fark ettim float64. Ben sonuçları için duyarsız hale getirmek için eklemek gerekir regülarizasyon minimum miktarda bir rehber var mı float32> - float64değişiklik? Özel olarak, kullanıcının durumu , …

10 numerical-analysis  stability  svd  intel-mkl  numerical-limitations 

Ben (örneğin yüksek sırasını hesaplamak çalışıyorum m=0, n=46bazı görüntü için) Zernike anlar. Ancak, radyal polinom ile ilgili bir sorunla karşılaşıyorum (bkz. Wikipedia ). Bu, [0 1] aralığında tanımlanan bir polinomdur. Aşağıdaki MATLAB koduna bakın function R = radial_polynomial(m,n,RHO) R = 0; for k = 0:((n-m)/2) R = R + (-1).^k.*factorial(n-k) …

9 matlab  polynomials  numerical-limitations