最美情侣中文字幕电影,在线麻豆精品传媒,在线网站高清黄,久久黄色视频

歡迎光臨散文網(wǎng) 會員登陸 & 注冊

隨機秩揭方法的理論分析(一)

2023-04-18 23:17 作者:安漢少俠  | 我要投稿

本文是對?Joel A. Tropp 教授所寫的?Randomized Numerical Linear Algebra: Foundations & Algorithms 中第 11 章 the?randomized?rangefinder 的閱讀筆記.? 這一章給出了隨機秩揭方法的理論分析,其中的部分結(jié)論十分有趣且富有啟發(fā)性,故作此筆記.?

首先給出秩揭問題的定義. 設(shè)?%5Cmathbf%7BB%7D%5Cin%20%5Cmathbb%7BF%7D%5E%7Bm%5Ctimes%20n%7D?是輸入矩陣,l%5Cle%20%5Cmin%5Cleft%5C%7B%20m%2Cn%20%5Cright%5C%7D是子空間維數(shù). 找出一個列正交陣 %5Cmathbf%7BQ%7D%5Cin%20%5Cmathbb%7BF%7D%5E%7Bm%5Ctimes%20l%7D,使得 %5Cmathbf%7BQ%7D?恰好對應(yīng)于?%5Cmathbf%7BB%7D?的前?l 個左奇異向量. 對于我們求出的?%5Cmathbf%7BQ%7D,用?%5C%7C%20%5Cmathbf%7BB%7D%20-%20%20%5Cmathbf%20%7BQ%7D%5Cmathbf%7BQ%7D%5E%7B*%7D%20%5Cmathbf%7BB%7D%20%20%5C%7C%20%3D%20%5C%7C%20(%5Cmathbf%7BI%7D%20-%20%20%5Cmathbf%20%7BQ%7D%5Cmathbf%7BQ%7D%5E%7B*%7D)%5Cmathbf%7BB%7D%20%20%5C%7C?來度量結(jié)果的準(zhǔn)確度,這里所用的范數(shù)為譜范數(shù)(也就是 2-范數(shù)).? 由 Eckhart-Young 定理,不難得到該最小值應(yīng)為 %5Csigma_%7Bl%2B1%7D.?

解決這個問題的一種方法是用測試矩陣?%5Cmathbf%7B%5COmega%7D%5Cin%5Cmathbb%7BF%7D%5E%7Bn%5Ctimes%20l%7D?右乘 %5Cmathbf%7BB%7D,得到 %5Cmathbf%7BY%7D%20%3D%20%20%5Cmathbf%7BB%7D%5Cmathbf%20%7B%5COmega%7D,對?%5Cmathbf%20%7BY%7D 做 QR 分解得到?%5Cmathbf%20%7BY%7D 的正交基?%5Cmathbf%20%7BQ%7D,將?%5Cmathbf%20%7BQ%7D 作為我們要求的正交陣. 接下來我們將探究該方法在理論上的可靠性.?

首先引入一些記號. 設(shè) %5Cmathbf%7BX%7D 表示測試矩陣,即 %5Cmathbf%7BY%7D%20%3D%20%5Cmathbf%7BB%7D%5Cmathbf%7BX%7D. 設(shè)?%5Cmathbf%7BP%7D_%5Cmathbf%7BY%7D?表示到 %5Cmathbf%7BY%7D?張成的子空間上的正交投影,即?%5Cmathbf%7BP%7D_%7B%5Cmathbf%20%7BY%7D%7D%20%3D%20%5Cmathbf%20%7BQ%7D%5Cmathbf%20%7BQ%7D%5E%7B*%7D. 用?%5Cmathbf%20%7BA%7D%2F%5Cmathbf%20%7BX%7D 表示對稱(半)正定陣?%5Cmathbf%20%7BA%7D 關(guān)于?%5Cmathbf%7BX%7D?的 Schur 補,即

%5Cbegin%7Balign%7D%0A%5Cmathbf%20%7BA%7D%2F%5Cmathbf%20%7BX%7D%20%3D%20%5Cmathbf%20%7BA%7D%20-%20(%5Cmathbf%20%7BA%7D%5Cmathbf%20%7BX%7D)(%5Cmathbf%20%7BX%7D%5E%7B*%7D%5Cmathbf%20%7BA%7D%5Cmathbf%20%7BX%7D)%5E%7B%5Cdagger%7D(%5Cmathbf%20%7BA%7D%5Cmathbf%20%7BX%7D)%5E%7B*%7D.%0A%5Cend%7Balign%7D

%5Cdagger 表示廣義逆. 定義?%5Cmathbf%7BE%7D%20%3A%3D%20%5Cmathbf%7BE%7D(%5Cmathbf%7BB%7D%2C%20%5Cmathbf%7BX%7D)%20%3A%3D%20(%5Cmathbf%7BI%7D%20-%20%5Cmathbf%7BP%7D_%7B%5Cmathbf%7BY%7D%7D)%5Cmathbf%7BB%7D%20.?

據(jù)上述定義,有如下關(guān)系成立

%7C%5Cmathbf%7BE%7D%7C%5E2%20%5Ccolon%3D%20%5Cmathbf%7BE%7D%5E%7B*%7D%5Cmathbf%7BE%7D%20%3D%20(%5Cmathbf%7BB%7D%5E%7B*%7D%5Cmathbf%7BB%7D)%2F%5Cmathbf%7BX%7D.

證明??注意到?%5Cmathbf%7BP%7D_%7B%5Cmathbf%7BY%7D%7D%20%3D%20(%5Cmathbf%7BB%7D%5Cmathbf%7BX%7D)((%5Cmathbf%7BB%7D%5Cmathbf%7BX%7D)%5E%7B*%7D%20%5Cmathbf%7BB%7D%5Cmathbf%7BX%7D%20)%5E%7B%5Cdagger%7D%20(%5Cmathbf%7BB%7D%5Cmathbf%7BX%7D)%5E%7B*%7D. 將?%5Cmathbf%7BB%7D%5E%7B*%7D%5Cmathbf%7BB%7D 記為?%5Cmathbf%7BA%7D,則

%5Cmathbf%7BE%7D%5E%7B*%7D%5Cmathbf%7BE%7D%20%3D%20%5Cmathbf%7BB%7D%5E%7B*%7D(%5Cmathbf%7BI%7D%20-%20%5Cmathbf%7BP%7D_%7B%5Cmathbf%7BY%7D%7D%0A%20)%5Cmathbf%7BB%7D%20%3D%20%5Cmathbf%20%7BA%7D%20-%20(%5Cmathbf%20%7BA%7D%5Cmathbf%20%7BX%7D)(%5Cmathbf%20%7BX%7D%5E%7B*%7D%5Cmathbf%20%7BA%7D%5Cmathbf%20%7BX%7D)%5E%7B%5Cdagger%7D(%5Cmathbf%20%7BA%7D%5Cmathbf%20%7BX%7D)%5E%7B*%7D..

由定義可知結(jié)論成立.

接著,我們給出一個引理. 首先仍要定義一些記號.?設(shè)?%5Cmathcal%7BH%7D?是有限維歐式空間,M?是?%5Cmathcal%7BH%7D?的一個子空間且?%5Cmathcal%7BH%7D%20%3D%20M%20%2B%20M%5E%7B%5Cbot%7D. 則半正定陣?%5Cmathbf%7BA%7D 可表示為

%5Cmathbf%7BA%7D%20%3D%20%0A%5Cbegin%7Bbmatrix%7D%0A%5Cmathbf%7BA%7D_%7B11%7D%20%26%20%5Cmathbf%7BA%7D_%7B12%7D%20%5C%5C%0A%5Cmathbf%7BA%7D_%7B21%7D%20%26%20%5Cmathbf%7BA%7D_%7B22%7D%0A%5Cend%7Bbmatrix%7D,

其中?%5Cmathbf%7BA%7D_%7B11%7D?和?%5Cmathbf%7BA%7D_%7B22%7D 分別是?%5Cmathcal%7BM%7D 和?%5Cmathcal%7BM%7D%5E%7B%5Cbot%7D?上的線性算子,%5Cmathbf%7BA%7D_%7B21%7D 和?%5Cmathbf%7BA%7D_%7B12%7D 分別是從?%5Cmathcal%7BM%7D 到?%5Cmathcal%7BM%7D%5E%7B%5Cbot%7D 和從?%5Cmathcal%7BM%7D%5E%7B%5Cbot%7D 到?%5Cmathcal%7BM%7D 的線性算子. 定義?

%5BM%5D%5Cmathbf%7BA%7D%20%3D%20%5Cmathbf%7BA%7D_%7B11%7D%20-%20%5Cmathbf%7BA%7D_%7B12%7D%5Cmathbf%7BA%7D_%7B22%7D%5E%7B%5Cdagger%7D%5Cmathbf%7BA%7D_%7B21%7D.

用?%5Cleft%20%5Clangle%5Ccdot%2C%5Ccdot%20%20%5Cright%20%5Crangle 表示兩向量的標(biāo)準(zhǔn)內(nèi)積,用?%5Cleft%20%5Clangle%20%5Ccdot%2C%5Ccdot%20%20%5Cright%20%5Crangle_%7B%5Cmathbf%7BA%7D%7D%20 表示由?%7B%5Cmathbf%7BA%7D%7D%20 定義的準(zhǔn)內(nèi)積 (因為 %7B%5Cmathbf%7BA%7D%7D%20 不一定正定,因此只能要求是準(zhǔn)內(nèi)積). 則有

引理1? 對任意半正定陣 %5Cmathbf%7BA%7D 和子空間 %5Cmathcal%7BM%7D%5Csubset%20%5Cmathcal%7BH%7D,

%5Cleft%20%5Clangle%20(%5B%5Cmathcal%7BM%7D%5D%5Cmathbf%7BA%7D)x%20%2Cx%20%5Cright%20%5Crangle%20%3D%20%5Cinf%5Cleft%20%5C%7B%20%5Cleft%20%5C%7C%20x%20-%20y%20%20%5Cright%20%5C%7C_%7B%5Cmathbf%7BA%7D%20%7D%5C%20%5Ccolon%5C%20y%5Cin%20%5Cmathcal%7BM%7D%5E%7B%5Cbot%7D%20%20%20%5Cright%20%5C%7D%20.

下面給出一個重要推論.?

推論1?設(shè)?%20%5Cmathbf%7BB%7D%5E%7B*%7D%5Cmathbf%7BB%7D%20%5Cle%20%5Cmathbf%7BC%7D%5E%7B*%7D%5Cmathbf%7BC%7D%20,這里 %5Cle 表示?%5Cmathbf%7BC%7D%5E%7B*%7D%5Cmathbf%7BC%7D%20-%20%5Cmathbf%7BB%7D%5E%7B*%7D%5Cmathbf%7BB%7D?半正定. 則對任一測試矩陣 %5Cmathbf%7BX%7D,有

%7C%5Cmathbf%7BE%7D(%5Cmathbf%7BB%7D%2C%20%5Cmathbf%7BX%7D)%7C%5E%7B2%7D%20%3D%20(%5Cmathbf%7BB%7D%5E%7B*%7D%5Cmathbf%7BB%7D)%2F%5Cmathbf%7BX%7D%20%5Cle%20(%5Cmathbf%7BC%7D%5E%7B*%7D%5Cmathbf%7BC%7D)%2F%5Cmathbf%7BX%7D%20%3D%20%7C%5Cmathbf%7BE%7D(%5Cmathbf%7BC%7D%2C%20%5Cmathbf%7BX%7D)%7C%5E2.

證明? ?%5Cmathbf%7BS%7D%20%3D%20%5Cmathbf%7BB%7D%5E%7B*%7D%5Cmathbf%7BB%7D?及 %5Cmathbf%7BT%7D%20%3D%20%5Cmathbf%7BC%7D%5E%7B*%7D%5Cmathbf%7BC%7D.?則?%5Cmathbf%7BS%7D%20 和?%5Cmathbf%7BT%7D%20 半正定. 構(gòu)造?

%5Ctilde%7B%5Cmathbf%7BT%7D%7D%20%3D%20%20%5Cbegin%7Bbmatrix%7D%0A%5Cmathbf%7BT%7D%20%26%20%5Cmathbf%7BT%7D%5Cmathbf%7BX%7D%20%20%20%5C%5C%0A(%5Cmathbf%7BT%7D%5Cmathbf%7BX%7D)%5E%7B*%7D%20%26%20%5Cmathbf%7BX%7D%5E%7B*%7D%5Cmathbf%7BT%7D%5Cmathbf%7BX%7D%20%20%0A%5Cend%7Bbmatrix%7D 及?%5Ctilde%7B%5Cmathbf%7BS%7D%7D%20%3D%20%5Cbegin%7Bbmatrix%7D%0A%5Cmathbf%7BS%7D%20%26%20%5Cmathbf%7BS%7D%5Cmathbf%7BX%7D%20%20%20%5C%5C%0A(%5Cmathbf%7BS%7D%5Cmathbf%7BX%7D)%5E%7B*%7D%20%26%20%5Cmathbf%7BX%7D%5E%7B*%7D%5Cmathbf%7BS%7D%5Cmathbf%7BX%7D%20%20%0A%5Cend%7Bbmatrix%7D,

易驗證?%5Ctilde%7B%5Cmathbf%7BT%7D%20%7D 和?%5Ctilde%7B%5Cmathbf%7BS%7D%20%7D 都是半正定的,且?%5Ctilde%7B%5Cmathbf%7BS%7D%20%7D%5Cle%5Ctilde%7B%5Cmathbf%7BT%7D%20%7D. ?取?%5Cmathcal%7BM%7D 為算子?%5Cmathbf%7BT%7D?(或 %5Cmathbf%7BS%7D) 對應(yīng)的子空間,則對任意的??y%5Cin%20%5Cmathcal%7BM%7D%5E%7B%5Cbot%7D?及?x%5Cin%20%5Cmathcal%7BH%7D ,有

%5Cleft%20%5C%7C%20x%20-%20y%20%5Cright%20%5C%7C_%7B%5Ctilde%7B%5Cmathbf%7BS%7D%7D%7D%20%5Cle%20%5Cleft%20%5C%7C%20x%20-%20y%20%5Cright%20%5C%7C_%7B%5Ctilde%7B%5Cmathbf%7BT%7D%7D%20%7D.

由引理1 可知對任意?x%5Cin%20%5Cmathcal%7BH%7D,都有

%5Cleft%20%5Clangle%20(%5B%5Cmathcal%7BM%7D%5D%5Cmathbf%7BS%7D)x%20%2Cx%20%5Cright%20%5Crangle%20%5Cle%20%5Cleft%20%5Clangle%20(%5B%5Cmathcal%7BM%7D%5D%5Cmathbf%7BT%7D)x%20%2Cx%20%5Cright%20%5Crangle%20.

據(jù)?%5BM%5D%5Cmathbf%7BA%7D 的定義可知結(jié)論成立.?

隨機秩揭方法的理論分析(一)的評論 (共 條)

分享到微博請遵守國家法律
邵武市| 札达县| 慈利县| 开封市| 瓦房店市| 张家口市| 湘阴县| 清原| 沅江市| 石河子市| 藁城市| 永顺县| 靖边县| 澜沧| 乌兰察布市| 固阳县| 筠连县| 高密市| 涿鹿县| 化隆| 仙居县| 石渠县| 宁化县| 宁国市| 准格尔旗| 多伦县| 兴海县| 湘潭县| 武汉市| 大姚县| 黄陵县| 开鲁县| 灵宝市| 文水县| 凤山县| 林芝县| 舒城县| 平遥县| 福泉市| 同江市| 曲沃县|