การลู่เข้าที่ไม่ใช่แบบโมโนโทนิกในปัญหาจุดคงที่


13

พื้นหลัง

ฉันกำลังแก้ไขตัวแปรของสมการOrnstein-Zernikeจากทฤษฎีของเหลว abstractly ปัญหาสามารถแสดงเป็นการแก้ปัญหาจุดคง( R ) = C ( R )ที่เป็นผู้ดำเนิน Integro-เกี่ยวกับพีชคณิตและ( R )เป็นฟังก์ชั่นการแก้ปัญหา (ฟังก์ชั่นความสัมพันธ์ OZ โดยตรง) ฉันกำลังแก้ไขโดยการทำซ้ำของ Picard ซึ่งฉันได้เตรียมโซลูชันทดลองใช้เบื้องต้นc 0 ( r )และสร้างโซลูชันทดลองใช้ใหม่โดยโครงการ c j + 1 = α (Ac(r)=c(r)Ac(r)c0(r) ที่ αเป็นพารามิเตอร์ที่ปรับค่าได้ซึ่งควบคุมการผสมของ cและ A c ที่ใช้ในโซลูชันทดลองใช้ถัดไป สำหรับการสนทนานี้สมมติว่าค่าของ αนั้นไม่สำคัญ ฉันทำซ้ำจนกว่าซ้ำลู่ไปภายในความอดทนต้องการ ε : Δ J + 1 d R | c j + 1 ( r ) - c

cj+1=α(Acj)+(1α)cj ,
αcAcαϵ ในตัวแปรของปัญหา Aขึ้นอยู่กับพารามิเตอร์ λและคำถามของฉันเกี่ยวกับการลู่ของ A c = cขึ้นอยู่กับพารามิเตอร์นี้
Δj+1dr|cj+1(r)cj(r)|<ϵ .
AλAc=c

สำหรับค่าที่หลากหลายสำหรับชุดการวนซ้ำด้านบนจะมาบรรจบกันอย่างรวดเร็วแบบทวีคูณ อย่างไรก็ตามในขณะที่ฉันลดระดับλในที่สุดฉันก็มาถึงระบอบการปกครองที่คอนเวอร์เจนซ์นั้นไม่ใช่แบบโมโนโทนิกดังภาพด้านล่าง λλการโจมตีของการลู่เข้าที่ไม่ใช่โมโนโทน

คำถามสำคัญ

ในการแก้ปัญหาแบบวนซ้ำเพื่อแก้ไขปัญหาจุดคงที่การลู่เข้าที่ไม่ใช่โมโนโพนิคมีความสำคัญเป็นพิเศษหรือไม่? มันเป็นสัญญาณว่ารูปแบบการวนซ้ำของฉันใกล้จะขาดเสถียรภาพหรือไม่? สิ่งสำคัญที่สุดคือการบรรจบที่ไม่ใช่โมโนโทนทำให้ฉันสงสัยว่าการแก้ปัญหาแบบ "แปรสภาพ" ไม่ใช่วิธีแก้ปัญหาจุดคงที่ที่ดีหรือไม่?

คำตอบ:


1

xx=f(x)xfx(x)αα<1x

  1. λ

  2. หากวิธีการแก้ปัญหาของคุณมาบรรจบกันภายในความคลาดเคลื่อนที่ยอมรับได้อย่างถูกต้องซึ่งรวมถึงตัวเลขขนาดเล็กก็มีเช่นกัน


คุณช่วยชี้แจงจุดที่สองของคุณได้หรือไม่?
Endulum

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