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


10

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

คำตอบ:


2

ynyn+1=U(yn)U

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

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

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