Một tiêu chí dừng tốt khi sử dụng phương pháp lặp để tìm giá trị bản địa là gì?


8

Tôi đọc câu trả lời này và nhận ra rằng tôi đã sử dụng sự khác biệt giữa các lần lặp thành công để xác định một tiêu chí dừng cho một phương pháp lặp để tìm giá trị / vectơ.

Các tiêu chí dừng tốt cho các phương pháp lặp hội tụ đến giá trị riêng và hàm riêng là gì?

Câu trả lời:


8

Cuốn sách Phương pháp số của Youssef Saad cho các vấn đề Eigenvalue lớn, ấn bản thứ 2 sử dụng định mức của vectơ dư để xác định tiêu chí hội tụ. Ông định nghĩa vectơ còn lại như sau trên trang 59:

Cho một ma trận , một giả định eigenvalue ~ λC và eigenvector giả định ~ uC n kết hợp với ~ λ , vector dư r gắn liền với cặp ( ~ λ , ~ u )MộtCn×nλ~Cbạn~Cnλ~r(λ~,bạn~)

r= =Mộtbạn~-λ~bạn~.

Nhiều kết quả lỗi trong cuốn sách của Saad được nêu theo định mức của vectơ dư (nói chung là 2 chỉ tiêu) và anh ta sử dụng định mức của vectơ còn lại như một thước đo để hội tụ bất cứ khi nào anh ta trình bày kết quả bằng số. Dựa trên bằng chứng đó, tiêu chí dừng lại sẽ là

r<ε.

SLEPc (dựa trên PETSc ), dường như sử dụng một tiêu chuẩn hội tụ tương tự trong họ thực hành các bài tập (họ sử dụng trong bài tập 1 và 2 thay vì).r/λ~bạn~

Tuy nhiên, LAPACK không nhất thiết phải sử dụng số liệu đó để hội tụ (ví dụ, trong ghi chú làm việc LAPACK (LAWN) # 15 , sử dụng phương pháp Jacobi để tính toán hàm riêng và giá trị riêng của ma trận xác định dương đối xứng). Các LUẬT khá dày đặc (tha thứ cho trò chơi chữ) và kỹ thuật, nhưng nếu bạn muốn xem những triển khai chất lượng cao nào được sử dụng để hội tụ, có lẽ chúng đáng để đọc chi tiết.


Theo những gì tôi hiểu, đối với các hệ tuyến tính , nếu A bị điều hòa, phần dư có thể nhỏ hơn nhiều so với sai số thực tế của vectơ giải pháp x . Điều này có đúng với các vấn đề eigenvalue không? Nếu vậy, chúng ta có thể làm gì để đảm bảo rằng uλ trong trường hợp của chúng ta được hội tụ hợp lý khi phần dư nhỏ? Mộtx= =bMộtxbạnλ
nukeguy
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.