Làm thế nào để sửa đổi thứ hạng thấp ảnh hưởng đến sự hội tụ phương pháp Krylov?


14

Giả sử tôi có hệ thống tuyến tính , hội tụ nhanh chóng bằng phương pháp Krylov phù hợp (như CG hoặc GMRES) cho tất cả . Nếu là ma trận có thứ hạng thấp , liệu phương thức Krylov tương tự trên hệ thống cũng hội tụ nhanh chóng (lý tưởng với số lần lặp thêm phụ thuộc gần như chỉ phụ thuộc vào )?Mộtx= =bB r ( A + B ) x = b rbBr(Một+B)x= =br

Một ví dụ về một hệ thống như vậy sẽ là độ co giãn và uốn cong của màng tiền điều kiện cộng với các điều khoản áp suất không khí không phân đoạn với cấu trúc sản phẩm bên ngoài dày đặc.

Lưu ý rằng vấn đề là như nhau có hoặc không có chuẩn bị trước, vì là một cấp bậc r sửa đổi của P Một Q .P(Một+B)Q= =PMộtQ+PBQrPMộtQ

Câu trả lời:


7

Nếu không gian con Krylov của bạn dựa trên quyền hạn của , việc hội tụ sẽ bị trì hoãn bởi một số lần lặp nhiều nhất là thứ hạng của hiệu chỉnh. Nếu nó dựa trên quyền hạn của A T A thì nhiều nhất là gấp đôi con số này.MộtMộtTMột

Tôi chưa thấy một tuyên bố như vậy trong tài liệu. Nhưng để thấy tính hợp lệ trong trường hợp đầu tiên, đủ để chỉ ra rằng không gian Krylov thứ của ma trận A + U S V T trong đó U , V có các cột r được chứa trong không gian tương ứng mà không cần hiệu chỉnh thứ hạng thấp nhưng với một tương ứng chỉ số cao hơn k + r . Điều này là đơn giản để xác minh.kMột+BạnSVTBạn,Vrk+r


Bạn có thể giải thích những gì bạn có nghĩa là "dựa trên sức mạnh của "? Người giải Krylov chỉ được cung cấp thông tin về A + B chứ không phải A trực tiếp. MộtMột+BMột
Geoffrey Irving

Đừng bận tâm: có lẽ bạn có nghĩa là sức mạnh của ma trận trong câu hỏi, vì vậy trong trường hợp này. Một+B
Geoffrey Irving

Đúng. Phương pháp này có một ma trận như một tham số, và ma trận này thường được biểu hiện bằng . Một
Arnold Neumaier

Có lẽ sự quan tâm hơn nữa bạn có thể viết lại phương trình của bạn (hoặc giải pháp) với một số yêu cầu về để x = ( E + Σ k = 1 ( A - 1 B ) k ) Một - 1 b mà có thể giúp đỡ nếu B là lũy linh hay A - 1 B định mức nhỏ. Người ta cũng nhận ra sự phụ thuộc vào giải pháp của vấn đề. Bx= =(E+Σk= =1(Một-1B)k)Một-1bBMột-1Bundisturbed
Bastian Ebeling
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.