Science computationnelle

Questions-réponses pour les scientifiques utilisant des ordinateurs pour résoudre des problèmes scientifiques

3
Question d'approximation des moindres carrés
Je prends un cours sur le calcul scientifique, et nous venons de passer en revue l'approximation des moindres carrés. Ma question concerne spécifiquement l'approximation à l'aide de polynômes. Je comprends que si vous avez n + 1 points de données, vous pouvez trouver un polynôme unique de degré n qui …
1
Méthodes numériques pour inverser les transformées intégrales?
J'essaie d'inverser numériquement la transformation intégrale suivante: F(y)=∫∞0yexp[−12(y2+x2)]I0(xy)f(x)dxF(y)=∫0∞yexp⁡[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Donc, pour un j'ai besoin d'approximer où:F(y)F(y)F(y)f(x)f(x)f(x) f(x)f(x)f(x) et sont réels et positifsF(y)F(y)F(y) (ce sont des distributions de probabilité continues) x,yx,yx,y sont réels et positifs (ce sont des grandeurs) J'ai une méthode très compliquée et brutale …
2
Comment puis-je tracer la surface d'un tracé 4D?
J'essaie de tracer la fonction d'onde d'une particule dans une boîte 3D. Cela me demande de tracer 4 variables: les axes x, y, z et la fonction de densité de probabilité. La fonction de densité de probabilité est: abs((np.sin((p*np.pi*X)/a))*(np.sin((q*np.pi*Y)/b))*(np.sin((r*np.pi*Z)/c)))**2 J'utilise np.arange()pour les X, Y et Z. J'ai lu que pour …
3
Étant donné un système linéaire tridiagonal SPD, pouvons-nous précalculer de sorte que trois indices quelconques puissent être liés en temps O (1)?
Considérons un système linéaire tridiagonal défini symétrique positif où et . Étant donné trois indices , si nous supposons uniquement des lignes d'équation strictement entre et tenir, nous pouvons éliminer les variables intermédiaires pour obtenir une équation de la forme où . Cette équation relie la valeur de à indépendamment …
3
Meilleures méthodologies pour gérer un maillage en calcul par éléments finis parallèles?
Je développe actuellement une méthode de décomposition de domaine pour la solution du problème de diffusion. Fondamentalement, je résous un système de BVP Helmholtz de manière itérative. Je discrétise les équations en utilisant la méthode des éléments finis sur des maillages triangulaires ou tétraédriques. Je développe le code vers ma …
3
Maillage de données de surface 3D en python
J'ai un ensemble de données de points en 3 dimensions pour lesquels je voudrais construire un maillage, en utilisant python. Tous les logiciels que j'ai vus exigent que vous fournissiez les bords. Existe-t-il un programme en python qui prend en entrée un ensemble de points en 3D et sort un …
3
Tester si deux matrices 12x12 ont le même déterminant
On me donne une matrice symétrique, inversible, définie positive et dense. Je dois tester si où J est la matrice des uns.12×1212×1212 \times 12QQQJdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Je fais actuellement cela avec la bibliothèque de tatou mais cela s'avère trop lent. Le fait est que je dois le …
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.