การปรับขนาดตัวแปรจำเป็นสำหรับการแก้ปัญหา PDE หรือไม่?


15

ในการจำลองเซมิคอนดักเตอร์มันเป็นเรื่องธรรมดาที่สมการจะถูกปรับอัตราส่วนเพื่อให้พวกเขามีค่าปกติ ยกตัวอย่างเช่นในกรณีที่ความหนาแน่นของอิเล็กตรอนในเซมิคอนดักเตอร์อาจแตกต่างกันไปมากกว่า 18 ลำดับความสำคัญและสนามไฟฟ้าสามารถเปลี่ยนรูปร่างได้โดยมีขนาดของคำสั่งมากกว่า 6 (หรือมากกว่า)

อย่างไรก็ตามเอกสารไม่เคยให้เหตุผลในการทำเช่นนี้จริงๆ โดยส่วนตัวฉันมีความสุขในการจัดการกับสมการในหน่วยจริงมีข้อได้เปรียบเชิงตัวเลขที่จะทำเช่นนี้เป็นไปไม่ได้หรือไม่? ฉันคิดว่าด้วยความแม่นยำสองเท่าจะมีตัวเลขเพียงพอที่จะรับมือกับความผันผวนเหล่านี้


คำตอบทั้งสองนั้นมีประโยชน์มากขอบคุณมาก!


1
"สามารถเปลี่ยนแปลงคำสั่งได้มากกว่า 18 คำสั่ง" - และหากคุณพิจารณาว่ามีจำนวนตัวเลขกี่หลักในความแม่นยำสองเท่าคุณจะเห็นว่า "ด้วยความแม่นยำสองเท่าจะมีตัวเลขเพียงพอที่จะรับมือกับความผันผวนเหล่านี้" เป็นจริง ...
JM

1
และปัญหาที่แท้จริงเริ่มต้นเมื่อคุณป้อนตัวเลขเหล่านี้เป็นอัลกอริธึมเชิงตัวเลข: นำสี่เหลี่ยมและทันใดนั้นคุณมี 36 คำสั่งที่แตกต่างกันของขนาด ...
Christian Clason

คำตอบ:


9

การแก้ PDE (เชิงเส้น) ประกอบด้วยการลดทอนสมการเพื่อให้เกิดระบบเชิงเส้นซึ่งจะถูกแก้ไขโดยตัวแก้เชิงเส้นซึ่งการลู่เข้า (อัตรา) ขึ้นอยู่กับจำนวนเงื่อนไขของเมทริกซ์ การปรับขนาดตัวแปรมักจะลดจำนวนเงื่อนไขนี้ซึ่งเป็นการปรับปรุงการบรรจบกัน (สิ่งนี้โดยทั่วไปจะใช้จำนวนผู้ที่มีเงื่อนไขเป็นเส้นทแยงมุมดูความแม่นยำและความเสถียรของอัลกอริธึมเชิงตัวเลขของ Nicholas Higham )

การแก้ PDE ที่ไม่ใช่เชิงเส้นนอกจากนี้ยังต้องใช้วิธีการในการแก้สมการไม่เชิงเส้นเช่นวิธีการของนิวตันที่การปรับสเกลสามารถส่งผลต่อการลู่เข้า

เนื่องจากการทำให้เป็นมาตรฐานทุกอย่างมักจะใช้ความพยายามน้อยมากจึงเป็นความคิดที่ดี


ฉันแน่ใจว่า @ArnoldNeumaier มีมากกว่าที่จะพูดในหัวข้อนี้
Christian Clason

จำนวนสภาพของการฝึกอบรมฉันใช้ (ตัวแปร unscaled) เป็น~ 1.25 สิ่งนี้ดูสมเหตุสมผลหรือไม่? สิ่งนี้คำนวณโดยใช้วิธี 2-norm ( docs.scipy.org/doc/numpy/reference/generated/ ...... )
boyfarrell

κ2=1

1
@ boyfarrell: ฉันทำงานเป็นประจำกับหมายเลขเงื่อนไขที่มีขนาดใหญ่ถึง 10 ^ 7 พร้อมผลลัพธ์ที่ยอมรับได้ อย่างไรก็ตามฉันไม่ยอมรับหมายเลขเงื่อนไขที่สูงกว่า 10 ^ 9 มาก
jvriesem

9

εΔu+u=0 on Ω,u=1 on Ω.

ที่กล่าวว่าไม่มีการปรับขนาดของตัวแปรหรือโดเมนที่ขจัดปัญหานี้

uα

α2Δu=fα on αΩ
αu1
Δu=f on Ω.
uα(x):=u1(x/α)αα

4
และพารามิเตอร์ที่เหลือจะต้องมีความสำคัญสำหรับการกำหนดพฤติกรรมเชิงคุณภาพของการแก้ปัญหา - นี่คือเหตุผลที่จำนวน Reynolds มีความสำคัญในการเปลี่ยนแปลงของของไหล กระบวนการนี้เรียกว่าNondimensionalization
Christian Clason

แน่นอนการหาค่าพารามิเตอร์ที่เท่าเทียมกันนั้นเป็นปัญหาหลักในการหากลุ่มสมมาตรของ PDE ซึ่งเป็นปัญหาที่ยากโดยทั่วไป
lurscher

2

การจัดการกับจำนวนจุดลอยตัวอาจเป็นกลอุบายในการลบตัวเลขขนาดเล็กออกจากจำนวนที่มากขึ้นรวมถึงแง่มุมอื่น ๆ อีกมากมาย ฉันอยากจะแนะนำให้อ่าน John D. Cooks โพสต์บนบล็อกเช่น

กายวิภาคของจำนวนจุดลอยตัว

เช่นเดียวกับ Oracle

สิ่งที่นักวิทยาศาสตร์คอมพิวเตอร์ทุกคนควรรู้เกี่ยวกับเลขคณิตจุดลอยตัว

อัลกอริธึมเชิงตัวเลขบางอย่างสำหรับการย่อขนาดหรือการขยายให้ใหญ่สุดต้องการการทำให้เป็นมาตรฐานสำหรับความเสถียรเชิงตัวเลข

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.