Mehrfachintegration mit Python und Sympy

Hier möchte ich über die Mehrfachintegration mit Python schreiben.

Grundform

python



from sympy import *

x = symbols('x')
y = symbols('y')

f = x**2 + y**2 + 1

integrate(f,(x, 0, 1),(y,0,1))

Apropos

python



x = symbols('x')
y = symbols('y')

Teil ist

python



x = Symbol('x')
y = Symbol('y')

Oder

python



x,y = symbols('x y')

Aber es ist okay.

Sie können auch eine gekrümmte Fläche mit z = f zeichnen, indem Sie wie folgt schreiben.

python


import numpy as np
%matplotlib inline
import matplotlib.pyplot as plt

from mpl_toolkits.mplot3d import Axes3D

x = y = np.linspace(-5,5)
X,Y = np.meshgrid(x,y)
f = X**2 + Y**2 + 1
fig = plt.figure(figsize = (10,10))
ax = fig.add_subplot(1,1,1,projection="3d")
ax.plot_surface(X, Y, f)

ダウンロード.png

Der Integrationsbereich ist dreieckig

python


from sympy import *

x = symbols('x')
y = symbols('y')

f = x**2 + x*y*2 + 1

#Integrieren Sie aus den links geschriebenen Variablen.
integrate(f,( x, 0, 2-y),(y, 0, 2))

Der Integrationsbereich ist fächerförmig

python


from sympy import *

x = symbols('x')
y = symbols('y')

f = x**2 + y**2 + 1

#Das Schreiben des Integrationsbereichs als Ungleichung funktioniert nicht.
integrate(f,(x, 0, sqrt(1-y**2)),(y,0,1))

Für den Bereich des Skalarfeldes

python



\vec{r}(u,v) = ( \cos u, \sin u, v)\\
D:0 \leq u \leq \pi,~~0 \leq v \leq 1\\
Skalarfeld in~f=\sqrt{x^2+y^2+z^2}Finden Sie den Wert für den Bereich von.

python


from sympy import *

u = symbols('u')
v = symbols('v')

r =Matrix([ cos(u), sin(u), v])

A = [0]*3
A[0] = diff(r,u)[0]
A[1] = diff(r,u)[1]
A[2] = diff(r,u)[2]
B = [0]*3
B[0] = diff(r,v)[0]
B[1] = diff(r,v)[1]
B[2] = diff(r,v)[2]

C = np.cross(A,B)
print(C)
#Daher beträgt die Länge des Außenprodukts von A und B 1.

#Python ist cos(u)**2+sin(u)**2=Sie können 1 nicht verwenden, um die Formel selbst zu transformieren.
f = sqrt(1+r[2]**2)
integrate(f,(v, 0, 1),(u,0,pi))

Recommended Posts

Mehrfachintegration mit Python und Sympy
Frage: Die Mehrfachintegration mit Python funktioniert nicht
Entwickelt und verifiziert mit mehreren Python-Versionen mit direnv
[Python] Löse Gleichungen mit Sympy
Programmieren mit Python und Tkinter
Ver- und Entschlüsselung mit Python
Python und Hardware-Verwenden von RS232C mit Python-
Python mit Pyenv und Venv
Funktioniert mit Python und R.
Lösen Sie simultane normale Differentialgleichungen mit Python und SymPy.
Kommunizieren Sie mit FX-5204PS mit Python und PyUSB
Roboter läuft mit Arduino und Python
Installieren Sie Python 2.7.9 und Python 3.4.x mit pip.
Neuronales Netzwerk mit OpenCV 3 und Python 3
AM-Modulation und Demodulation mit Python
Scraping mit Node, Ruby und Python
Scraping mit Python, Selen und Chromedriver
Spielen Sie mit Poancare-Serien und SymPy
JSON-Codierung und -Decodierung mit Python
Hadoop-Einführung und MapReduce mit Python
[GUI in Python] PyQt5-Drag & Drop-
Veröffentlichen Sie mehrere Twitter-Bilder mit Python
Lesen und Schreiben von NetCDF mit Python
Animieren Sie mehrere Standbilder mit Python
Ich habe mit PyQt5 und Python3 gespielt
Koexistenz von Python2 und 3 mit CircleCI (1.0)
[Python] Erstellen mehrerer Fenster mit Tkinter
Sugoroku-Spiel und Zusatzspiel mit Python
FM-Modulation und Demodulation mit Python
Datenpipeline-Aufbau mit Python und Luigi
Berechnen Sie das Standardgewicht und zeigen Sie es mit Python an
FM-Modulation und Demodulation mit Python Part 3
[Automatisierung] Bearbeiten Sie Maus und Tastatur mit Python
Passwortlose Authentifizierung mit RDS und IAM (Python)
Python-Installation und Paketverwaltung mit pip
Verwenden von Python und MeCab mit Azure Databricks
POST verschieden mit Python und empfange mit Flask
Bilder mit Pupil, Python und OpenCV aufnehmen
Fraktal zum Erstellen und Spielen mit Python
Ein Memo mit Python2.7 und Python3 in CentOS
CentOS 6.4, Python 2.7.3, Apache, mod_wsgi, Django
Umgang mit "Jahren und Monaten" in Python
Verarbeiten Sie mehrere Listen mit for in Python
Ich habe Numba mit Python3.5 installiert und verwendet
Tweet-Analyse mit Python, Mecab und CaboCha
Verknüpfung von Python und JavaScript mit dem Jupiter-Notizbuch
Verkehrsüberwachung mit Kibana, ElasticSearch und Python
FM-Modulation und Demodulation mit Python Part 2
Laden Sie einfach mp3 / mp4 mit Python und youtube-dl herunter!
Betreiben Sie Haushaltsgeräte mit Python und IRKit
Bereinigen Sie die Python-Umgebung mit Pythonz und virtualenv
Üben des Web-Scrapings mit Python und Selen
Einfaches Web-Scraping mit Python und Ruby
Importieren und Exportieren von GeoTiff-Bildern mit Python
Ich verwende Tox und Python 3.3 mit Travis-CI
Fröhliche GUI-Konstruktion mit Elektron und Python
Verwenden Sie Python und MeCab mit Azure-Funktionen
Verwalten Sie mehrere Python-Versionen mit Update-Alternativen (Ubuntu)
Berühren Sie AWS mit Serverless Framework und Python