ตรวจสอบ Heuristic ของความมั่นคงเชิงตัวเลข


12

สมมติว่าฉันมีฟังก์ชั่นที่มีคุณค่าจริง ๆของตัวแปรบางตัวx iที่ฉันต้องการประเมินตัวเลข โดยทั่วไปสูตรสำหรับfสามารถมีผลิตภัณฑ์, ปันส่วน, ฟังก์ชั่น trancendental เป็นต้นและจะใช้เวลานานในการตรวจสอบความเสถียรเชิงตัวเลขของการวิเคราะห์ หรืออย่างน้อยก็ต้องใช้เวลาในการฝึกฝน สมมติว่าฉันไม่ได้เทียบเท่าสั้นกว่าที่มีความมั่นคง guaruanteed มีขั้นตอนวิธีการในการวิเคราะห์เสถียรภาพเชิงตัวเลขของf หรือไม่f(x1,,xN)xiff. ฉันคิดว่าการเปรียบเทียบกับผลการตกตะกอนโดยพลการที่ได้รับโดยใช้ระบบพีชคณิตของคอมพิวเตอร์ สมมติว่าฟังก์ชั่นนี้จะใช้งานใน C โดยใช้ฟังก์ชั่น stdlib และความแม่นยำเดียวหรือสองครั้ง ฉันควรเปรียบเทียบปริมาณใดกับปริมาณที่ประมาณคุณภาพของการประมาณค่าที่ จำกัด ล่วงหน้า? ฉันจะกำหนดค่าวิกฤตของตัวแปรได้อย่างไร ฉันจะเลือกคอมไพเลอร์และการปรับแต่งคอมไพเลอร์เพื่อให้ผู้อื่นสามารถสร้างผลลัพธ์ได้อย่างง่ายดายได้อย่างไร? ... ฉันรู้ว่าการตั้งค่าปัญหาน่าจะเป็นการทั่วไปเพื่อให้คำตอบที่ดี แต่ฉันก็ยังคิดว่านี่เป็นปัญหาที่พบบ่อยในวิทยาการคอมพิวเตอร์และสงสัยว่ามีการอ้างอิงที่เสนอมาตรฐานเพื่อทำการวิเคราะห์ดังกล่าว

คำตอบ:


7

สิ่งที่คุณกำลังมองหาคือสิ่งที่เรียกว่า "การวิเคราะห์ข้อผิดพลาดอัตโนมัติ" และเป็นเรื่องของบทที่ 26 ของหนังสือ Higham "ความแม่นยำและความเสถียรของอัลกอริธึมเชิงตัวเลข" 2nd ed. สำนักพิมพ์ SIAM

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

ฉันขอแนะนำให้คุณได้รับสำเนาของหนังสือเล่มนี้อ่านบทเกริ่นนำและบทที่ 26 และการอ้างอิงในนั้น



0

|f(x+ε)f(x)|C|ε|
Cfx,ϵ

และสิ่งที่สามารถทำได้ถ้าฟังก์ชั่นจะแตกต่างกันอย่างดุเดือดเหนือโดเมนหรือถ้าไม่มีอนุพันธ์ที่เป็นไปได้? มีเทคนิคอื่น ๆ หรือเราจะจบด้วยวิธีมอนติคาร์โล?
André

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