Computer >> คอมพิวเตอร์ >  >> การเขียนโปรแกรม >> Python

รวมชุด Legendre และกำหนดลำดับการรวมใน Python


ในการรวมซีรี่ส์ Legendre ให้ใช้เมธอด polynomial.legendre.legint() ใน Python วิธีการส่งคืนสัมประสิทธิ์ซีรีส์ Legendre c รวม m ครั้งจาก lbnd ตามแนวแกน ในการวนซ้ำแต่ละครั้ง อนุกรมผลลัพธ์จะถูกคูณด้วย scl และค่าคงที่การรวม k ถูกเพิ่มเข้าไป ตัวประกอบมาตราส่วนมีไว้สำหรับใช้ในการเปลี่ยนแปลงเชิงเส้นของตัวแปร พารามิเตอร์ที่ 1 c คืออาร์เรย์ของชุดค่าสัมประสิทธิ์ของ Legendre หาก c มีหลายมิติ แกนที่ต่างกันจะสัมพันธ์กับตัวแปรต่างๆ โดยที่องศาในแต่ละแกนกำหนดโดยดัชนีที่เกี่ยวข้อง

พารามิเตอร์ตัวที่ 2 m คือลำดับการรวม ต้องเป็นค่าบวก (ค่าเริ่มต้น:1). พารามิเตอร์ที่ 3 kis ค่าคงที่การรวม ค่าของอินทิกรัลแรกที่ lbnd คือค่าแรกในรายการ ค่าของอินทิกรัลที่สองที่ lbnd คือค่าที่สอง ฯลฯ หาก k ==[] (ค่าเริ่มต้น) ค่าคงที่ทั้งหมดจะถูกตั้งค่าเป็นศูนย์ ถ้า m ==1 สามารถระบุสเกลาร์เดี่ยวแทนรายการได้

พารามิเตอร์ตัวที่ 4 lbnd คือขอบเขตล่างของอินทิกรัล (ค่าเริ่มต้น:0) พารามิเตอร์ที่ 5 scl คือ ascalar หลังจากการผสานรวมแต่ละครั้ง ผลลัพธ์จะถูกคูณด้วย scl ก่อนที่ค่าคงที่การรวมจะถูกเพิ่ม พารามิเตอร์ที่ 6 แกนคือแกนที่ใช้อินทิกรัล (ค่าเริ่มต้น:0).

ขั้นตอน

ขั้นแรก นำเข้าไลบรารีที่จำเป็น -

import numpy as np
from numpy.polynomial import legendre as L

สร้างอาร์เรย์ของสัมประสิทธิ์ -

c = np.array([1,2,3])

แสดงอาร์เรย์ -

print("Our Array...\n",c)

ตรวจสอบขนาด -

print("\nDimensions of our Array...\n",c.ndim)

รับประเภทข้อมูล -

print("\nDatatype of our Array object...\n",c.dtype)

รับรูปร่าง -

print("\nShape of our Array object...\n",c.shape)

ในการรวมซีรี่ส์ Legendre ให้ใช้เมธอด polynomial.legendre.legint() ใน Python -

print("\nResult...\n",L.legint(c, m = 3))

ตัวอย่าง

import numpy as np
from numpy.polynomial import legendre as L

# Create an array of coefficients
c = np.array([1,2,3])

# Display the array
print("Our Array...\n",c)

# Check the Dimensions
print("\nDimensions of our Array...\n",c.ndim)

# Get the Datatype
print("\nDatatype of our Array object...\n",c.dtype)

# Get the Shape
print("\nShape of our Array object...\n",c.shape)

# To integrate a Legendre series, use the polynomial.legendre.legint() method in Pytho
print("\nResult...\n",L.legint(c, m = 3))

ผลลัพธ์

Our Array...
   [1 2 3]

Dimensions of our Array...
1

Datatype of our Array object...
int64

Shape of our Array object...
(3,)

Result...
   [ 1.66666667e-02 -1.78571429e-02 4.76190476e-02 -1.73472348e-18
1.90476190e-02 9.52380952e-03]