อ่าน 8 นาที
วิธี Galerkin
ใน ทางคณิตศาสตร์ สาขา การวิเคราะห์เชิงตัวเลข วิธี การของกาเลอร์กิน เป็นกลุ่มวิธีการสำหรับการแปลงปัญหาตัวดำเนินการต่อเนื่อง เช่น สมการเชิงอนุพันธ์ ซึ่งมักอยู่ใน รูปแบบอ่อน...
วิธี Galerkin
| สมการเชิงอนุพันธ์ |
|---|
| ขอบเขต |
| การจำแนกประเภท |
| สารละลาย |
| ประชากร |
| สมการที่มีชื่อ |
ในทางคณิตศาสตร์สาขาการวิเคราะห์เชิงตัวเลขวิธีการของกาเลอร์กินเป็นกลุ่มวิธีการสำหรับการแปลงปัญหาตัวดำเนินการต่อเนื่อง เช่นสมการเชิงอนุพันธ์ซึ่งมักอยู่ในรูปแบบอ่อนให้เป็นปัญหาแบบไม่ต่อเนื่องโดยการใช้ข้อจำกัดเชิงเส้นที่กำหนดโดยชุดฟังก์ชันพื้นฐานจำนวนจำกัด วิธีการนี้ตั้งชื่อตามนักคณิตศาสตร์ชาวโซเวียตบอริส กาเลอร์กิน
บ่อยครั้งเมื่อกล่าวถึงวิธีการของ Galerkin มักจะมีการระบุชื่อวิธีการ พร้อมทั้งข้อสมมติฐานและวิธีการประมาณค่าที่ใช้โดยทั่วไปด้วย:
- วิธี Ritz–Galerkin (ตาม Walther Ritz ) โดยทั่วไปจะถือว่ารูปแบบทวิเชิงเส้นสมมาตรและเป็นบวกแน่นอน ในการกำหนดสูตรแบบอ่อนโดยที่สมการเชิงอนุพันธ์สำหรับระบบทางกายภาพสามารถกำหนดสูตรได้โดยการลดฟังก์ชันกำลังสองที่ แสดงถึง พลังงานของระบบและคำตอบโดยประมาณคือการรวมเชิงเส้นของชุดฟังก์ชันพื้นฐานที่กำหนด [ 1 ]
- วิธี Bubnov–Galerkin (ตั้งชื่อตามIvan Bubnov ) ไม่จำเป็นต้องให้รูปแบบทวิเชิงเส้นมีความสมมาตรและแทนที่การลดพลังงานด้วย ข้อจำกัด เชิงตั้งฉากที่กำหนดโดยฟังก์ชันพื้นฐานเดียวกันกับที่ใช้ในการประมาณค่าคำตอบ ในการกำหนดสมการเชิงอนุพันธ์ในรูปแบบตัว ดำเนินการ วิธี Bubnov–Galerkin สามารถมองได้ว่าเป็นการประยุกต์ใช้การฉายภาพเชิงตั้งฉากกับตัวดำเนินการ
- วิธี Petrov–Galerkin (ตาม Georgii I. Petrov [ 2 ] ) อนุญาตให้ใช้ฟังก์ชันพื้นฐานสำหรับข้อจำกัดเชิงตั้งฉาก (เรียกว่าฟังก์ชันพื้นฐานทดสอบ ) ที่แตกต่างจากฟังก์ชันพื้นฐานที่ใช้ในการประมาณคำตอบ วิธี Petrov–Galerkin สามารถมองได้ว่าเป็นส่วนขยายของวิธี Bubnov–Galerkin โดยใช้การฉายภาพที่ไม่จำเป็นต้องตั้งฉากในการกำหนดตัวดำเนินการของสมการเชิงอนุพันธ์
ตัวอย่างของวิธีการ Galerkin ได้แก่:
- วิธีGalerkin ของค่าตกค้างถ่วงน้ำหนักซึ่งเป็นวิธีที่พบบ่อยที่สุดในการคำนวณเมทริกซ์ความแข็ง ทั่วโลก ในวิธีองค์ประกอบจำกัด [ 3 ] [ 4 ]
- วิธีองค์ประกอบขอบเขตสำหรับการแก้สมการเชิงอินทิกรัล
- วิธีการของปริภูมิ ย่อยKrylov [ 5 ]
สมการเชิงเส้นในปริภูมิฮิลเบิร์ต
การกำหนดสูตรแบบอ่อนของสมการเชิงเส้น
ขอแนะนำวิธีการของ Galerkin โดยใช้ปัญหาเชิงนามธรรมที่กำหนดในรูปแบบอ่อนบนปริภูมิฮิลเบิร์ต กล่าวคือ
- จงหาค่าที่ทำให้สำหรับทุกค่า
ในที่นี้เป็นรูปแบบทวิเชิงเส้น (ข้อกำหนดที่แน่นอนเกี่ยวกับจะระบุในภายหลัง) และเป็นฟังก์ชันเชิงเส้นที่มีขอบเขตบน
การลดมิติแบบกาเลอร์กิน
เลือกซับสเปซที่มีมิติnแล้วแก้ปัญหาที่ฉายออกมา:
- จงหาค่าที่ทำให้สำหรับทุกค่า
เราเรียกสมการนี้ว่าสมการกาเลอร์กินสังเกตว่าสมการยังคงไม่เปลี่ยนแปลง มีเพียงปริภูมิเท่านั้นที่เปลี่ยนแปลง การลดปัญหาให้เหลือปริภูมิย่อยเวกเตอร์มิติจำกัดทำให้เราสามารถคำนวณค่าได้โดยใช้การรวมเชิงเส้นจำกัดของเวกเตอร์ฐานในปริภูมิย่อยนั้น
ความตั้งฉากแบบกาเลอร์กิน
คุณสมบัติสำคัญของวิธีการของ Galerkin คือ ข้อผิดพลาดนั้นตั้งฉากกับปริภูมิย่อยที่เลือกไว้ เนื่องจาก เราจึงสามารถใช้เป็นเวกเตอร์ทดสอบในสมการดั้งเดิมได้ เมื่อลบทั้งสอง เราจะได้ความสัมพันธ์เชิงตั้งฉากของ Galerkin สำหรับข้อผิดพลาดซึ่งก็คือข้อผิดพลาดระหว่างคำตอบของปัญหาดั้งเดิมและคำตอบของสมการ Galerkin
สมการของ Galerkin ในรูปแบบเมทริกซ์
เนื่องจากเป้าหมายของวิธีของกาเลอร์กินคือการสร้างระบบสมการเชิงเส้นเราจึงสร้างรูปแบบเมทริกซ์ของวิธีดังกล่าว ซึ่งสามารถนำมาใช้ในการคำนวณหาคำตอบโดยใช้อัลกอริทึมได้
ให้เป็นฐานสำหรับจากนั้นก็เพียงพอที่จะใช้ฐานเหล่านี้ทีละฐานเพื่อทดสอบสมการกาเลอร์กิน กล่าวคือ หาที่ทำให้
เราขยายโดยใช้ฐานนี้แล้วแทนค่าลงในสมการข้างต้น เพื่อให้ได้
สมการก่อนหน้านี้เป็นระบบสมการเชิงเส้นโดยที่
สมมาตรของเมทริกซ์
เนื่องจากนิยามของสมาชิกในเมทริกซ์ เมทริกซ์ของสมการกาเลอร์กินจะเป็นเมทริกซ์สมมาตรก็ต่อเมื่อรูปแบบทวิเชิงเส้นเป็นเมทริกซ์สมมาตร เท่านั้น
การวิเคราะห์วิธีการของ Galerkin
ในที่นี้ เราจะจำกัดตัวเองอยู่เฉพาะรูปแบบทวิเชิงเส้นสมมาตรเท่านั้น นั่นคือ
แม้ว่านี่จะไม่ใช่ข้อจำกัดที่แท้จริงของวิธีการกาเลอร์กิน แต่การประยุกต์ใช้ทฤษฎีมาตรฐานจะง่ายขึ้นมาก นอกจากนี้ อาจจำเป็นต้องใช้ วิธีการเปตรอฟ-กาเลอร์กินในกรณีที่ไม่สมมาตร
การวิเคราะห์วิธีการเหล่านี้ดำเนินการเป็นสองขั้นตอน ขั้นแรก เราจะแสดงให้เห็นว่าสมการกาเลอร์กินเป็นปัญหาที่มีคำตอบแน่นอนในความหมายของฮาดามาร์ดและดังนั้นจึงมีคำตอบเดียว ในขั้นตอนที่สอง เราจะศึกษาคุณภาพของการประมาณค่าของคำตอบกาเลอร์กิน
การวิเคราะห์ส่วนใหญ่จะอาศัยคุณสมบัติสองประการของรูปแบบทวิเชิงเส้นได้แก่
- ขอบเขต: สำหรับทุกการยึดถือ
- สำหรับค่าคงที่บางค่า
- ความรี: สำหรับทุกการจับยึด
- สำหรับค่าคงที่บางค่า
ตามทฤษฎีบท Lax-Milgram (ดูการกำหนดรูปแบบอ่อน ) เงื่อนไขทั้งสองนี้บ่งชี้ถึงความถูกต้องของปัญหาดั้งเดิมในการกำหนดรูปแบบอ่อน บรรทัดฐานทั้งหมดในส่วนต่อไปนี้จะเป็นบรรทัดฐานที่ความไม่เท่าเทียมกันข้างต้นเป็นจริง (บรรทัดฐานเหล่านี้มักเรียกว่าบรรทัดฐานพลังงาน)
ความถูกต้องของสมการกาเลอร์กิน
เนื่องจากความมีขอบเขตและความเป็นวงรีของรูปแบบทวิเชิงเส้นนั้นใช้ได้กับดังนั้น ความถูกต้องของปัญหา Galerkin จึงสืบทอดมาจากความถูกต้องของปัญหาดั้งเดิม
การประมาณกึ่งที่ดีที่สุด (บทแทรกของ Céa)
ข้อผิดพลาดระหว่างผลลัพธ์เดิมกับผลลัพธ์ของวิธี Galerkin ยอมรับการประมาณค่า
นั่นหมายความว่า เมื่อพิจารณาถึงค่าคงที่แล้ววิธีแก้ปัญหาแบบกาเลอร์กิน จะใกล้เคียงกับวิธีแก้ปัญหาดั้งเดิมมากที่สุดเท่ากับเวกเตอร์อื่นๆ ใน โดยเฉพาะอย่างยิ่ง การศึกษาการประมาณค่าโดยใช้ปริภูมิ จะเพียงพอแล้วโดยไม่จำเป็นต้องคำนึงถึงสมการที่กำลังแก้เลย
การพิสูจน์
เนื่องจากบทพิสูจน์นั้นง่ายมากและเป็นหลักการพื้นฐานเบื้องหลังวิธีการของ Galerkin ทั้งหมด เราจึงรวมไว้ในที่นี้: โดยอาศัยความเป็นวงรีและความมีขอบเขตของรูปแบบทวิเชิงเส้น (อสมการ) และความตั้งฉากของ Galerkin (เครื่องหมายเท่ากับตรงกลาง) เราจะได้ว่า สำหรับค่าใดๆ:
การหารด้วยและการหาค่าต่ำสุดเหนือค่าที่เป็นไปได้ทั้งหมด จะได้บทพิสูจน์ย่อยนี้
คุณสมบัติการประมาณค่าที่ดีที่สุดของ Galerkin ในบรรทัดฐานพลังงาน
เพื่อความง่ายในการนำเสนอในส่วนด้านบน เราได้สมมติว่ารูปแบบทวิเชิงเส้นนั้นสมมาตรและเป็นบวกแน่นอน ซึ่งหมายความว่ามันเป็นผลคูณสเกลาร์และนิพจน์นั้นเป็นบรรทัดฐานเวกเตอร์ที่ถูกต้อง เรียกว่าบรรทัดฐานพลังงานภายใต้สมมติฐานเหล่านี้ เราสามารถพิสูจน์คุณสมบัติการประมาณค่าที่ดีที่สุดของกาเลอร์กินในบรรทัดฐานพลังงานได้อย่างง่ายดาย
โดยใช้คุณสมบัติการตั้งฉากแบบ Galerkin และอสมการ Cauchy–Schwarzสำหรับบรรทัดฐานพลังงาน เราจะได้
การหารด้วยและการหาค่าต่ำสุดเหนือค่าที่เป็นไปได้ทั้งหมดพิสูจน์ได้ว่า การประมาณแบบกาเลอร์กินเป็นการประมาณที่ดีที่สุดในบรรทัดฐานพลังงานภายในปริภูมิย่อย กล่าวคือเป็นเพียงการฉายภาพเชิงตั้งฉากโดยสัมพันธ์กับผลคูณสเกลาร์ของคำตอบไปยังปริภูมิย่อย
วิธี Galerkin สำหรับโครงสร้างแบบขั้นบันได
I. Elishakof , M. Amato, A. Marzani, PA Arvan และ JN Reddy [ 6 ] [ 7 ] [ 8 ] [ 9 ] ศึกษาการประยุกต์ใช้วิธี Galerkin กับโครงสร้างแบบขั้นบันได พวกเขาแสดงให้เห็นว่าฟังก์ชันทั่วไป ได้แก่ ฟังก์ชันขั้นบันไดหน่วย ฟังก์ชันเดลต้าของ Dirac และฟังก์ชันคู่ จำเป็นสำหรับการได้ผลลัพธ์ที่แม่นยำ
ประวัติศาสตร์
โดยทั่วไปแล้ววิธีการนี้ได้รับการยกย่องให้แก่Boris Galerkin [ 10 ] [ 11 ] วิธีการนี้ได้รับการอธิบายให้ผู้อ่านชาวตะวันตกเข้าใจโดย Hencky [ 12 ]และ Duncan [ 13 ] [ 14 ]เป็นต้น การลู่เข้าของวิธีการนี้ได้รับการศึกษาโดย Mikhlin [ 15 ]และ Leipholz [ 16 ] [ 17 ] [ 18 ] [ 19 ]ความสอดคล้องกับวิธีการ Fourier ได้รับการแสดงให้เห็นโดยElishakoff et al. [ 20 ] [ 21 ] [ 22 ]ความเท่าเทียมกันของวิธีการนี้กับวิธีการของ Ritz สำหรับปัญหาอนุรักษ์ได้รับการแสดงให้เห็นโดย Singer [ 23 ] Gander และ Wanner [ 24 ]แสดงให้เห็นว่าวิธีการของ Ritz และ Galerkin นำไปสู่วิธีการไฟไนต์เอเลเมนต์สมัยใหม่ได้อย่างไร Repin ได้กล่าวถึงการพัฒนาของวิธีการนี้มาเป็นเวลาหนึ่งร้อยปีแล้ว[ 25 ] Elishakoff, Kaplunov และ Kaplunov [ 26 ]แสดงให้เห็นว่าวิธีการของ Galerkin ไม่ได้ถูกพัฒนาโดย Ritz ซึ่งขัดแย้งกับคำกล่าวของ Timoshenko
ดูเพิ่มเติม
ลิงก์ภายนอก
- "วิธีของกาเลอร์กิน" , สารานุกรมคณิตศาสตร์ , EMS Press , 2001 [1994]
- วิธีการของ Galerkin จาก MathWorld
สรุปเนื้อหา
ข้อมูลสำคัญจากบทความ
ข้อมูลสำคัญเกี่ยวกับ วิธี Galerkin
ใน ทางคณิตศาสตร์ สาขา การวิเคราะห์เชิงตัวเลข วิธี การของกาเลอร์กิน เป็นกลุ่มวิธีการสำหรับการแปลงปัญหาตัวดำเนินการต่อเนื่อง เช่น สมการเชิงอนุพันธ์ ซึ่งมักอยู่ใน รูปแบบอ่อน...
การกำหนดสูตรแบบอ่อนของสมการเชิงเส้น
ขอแนะนำวิธีการของ Galerkin โดยใช้ปัญหาเชิงนามธรรมที่กำหนดใน รูปแบบอ่อน บน ปริภูมิฮิลเบิร์ต กล่าวคือ วี {\displaystyle V}
การลดมิติแบบกาเลอร์กิน
เลือกซับสเปซที่มีมิติ n แล้วแก้ปัญหาที่ฉายออกมา: วี n ⊂ วี {\displaystyle V_{n}\เซตย่อย V}
ความตั้งฉากแบบกาเลอร์กิน
คุณสมบัติสำคัญของวิธีการของ Galerkin คือ ข้อผิดพลาดนั้นตั้งฉากกับปริภูมิย่อยที่เลือกไว้ เนื่องจาก เราจึงสามารถใช้เป็นเวกเตอร์ทดสอบในสมการดั้งเดิมได้ เมื่อลบทั้งสอง เราจะได้ความสัมพันธ์เชิงตั้งฉากของ Galerkin...