มีทางลัดสำหรับระบบที่เป็นตัวเลขโดยประมาณของสมการเชิงอนุพันธ์สามัญเมื่อเป็นแบบอัตโนมัติหรือไม่?


10

อัลกอริทึมที่มีอยู่สำหรับการแก้ปัญหา ODE จะจัดการกับฟังก์ชันที่Y∈Rn แต่ในระบบกายภาพหลายระบบสมการเชิงอนุพันธ์จะเป็นอิสระดังนั้นdydydt=f(y,t)y∈Rn,y∈Rnโดยที่t ถูกปล่อยออกมา ด้วยข้อสมมติฐานที่ทำให้เข้าใจง่ายนี้จะเห็นการปรับปรุงอะไรบ้างในวิธีการเชิงตัวเลขที่มีอยู่ ตัวอย่างเช่นถ้าn=1ปัญหาจะเปลี่ยนเป็นt=∫dydydt=f(y)y∈Rntn=1และเราเปลี่ยนเป็นอัลกอริธึมที่แตกต่างอย่างสิ้นเชิงสำหรับการรวมอินทิกรัลหนึ่งมิติ สำหรับn>1การปรับปรุงเป็นไปได้สูงสุดคือการลดมิติของปีโดย 1 เพราะกรณีนี้ขึ้นกับเวลาที่สามารถจำลองโดยการผนวกทีจะYเปลี่ยนโดเมนของYจากRnเพื่อRn+1t=∫dyf(y)n>1ytyyRnRn+1

คำตอบ:


2

yn→yn+1=U(yn)U

∂ty=AyAU(y)=exp⁡(AΔt)y

สำหรับระบบที่ไม่ใช่เชิงเส้นนั้นไม่ใช่เรื่องง่าย แต่ขึ้นอยู่กับอัลกอริทึมการประเมินราคาแพงสามารถนำกลับมาใช้ใหม่ได้

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