Re: [分析] 雙向遞迴數列的推導問題
這篇回應關於 pseudo inverse 的想法。
還是用這個舉例
: y_n - y_{n-1} = x_n
: y_n = y_{-∞} + (u*x)_n
: 或
: y_n = y_{∞} + [(u-1)*x]_n
: 在 y_{-∞} 存在的情況下,u 可以是脈波響應。
1 0 0 0 0
-1 1 0 0 0
0 -1 1 0 0
0 0 -1 1 0
0 0 0 -1 1
以上矩陣是差分算子的矩陣表示,然後我砍掉(truncate)了一些東西。
其偽反矩陣:
1 0 0 0 0
1 1 0 0 0
1 1 1 0 0
1 1 1 1 0
1 1 1 1 1
這個矩陣與 u* 的作用相似。
: 而在 y_{∞} 存在的情況下則換成可以拿 u-1 當作脈波響應。
-1 1 0 0
0 -1 1 0
0 0 -1 1
0 0 0 -1
以上矩陣也是差分算子的矩陣表示,然後我又砍了一些東西。
其偽反矩陣:
-1 -1 -1 -1
0 -1 -1 -1
0 0 -1 -1
0 0 0 -1
這次矩陣跟 (u-1)* 長得比較像了。
: 兩個脈波響應只差了 1(此為齊次解),都同樣能有作用,
: 即使 u*x 與 (u-1)*x 算出來不一樣,但只要存在,他們就的確都是特解。
: 是相應於不同邊界條件的特解。
: 而此時會受到初始條件影響的東西,是 y_{-∞} 和 y_{∞}。
: 也就是說,不同初始條件與特解(那個捲積)的計算過程無涉,
: 但是會改變齊次項。
砍矩陣的時候下刀的方式不同,其實就是在考慮不同的邊界條件。
第一個差分矩陣左上角是 1,代表在擴大矩陣的時候,y_n 最前面的那項要一直能算。
第二個差分矩陣右下角是 -1,代表最後項要一直存在。
這與我們的邊界條件概念上是相符的。
其實還有其他下刀的方式,出來的偽反矩陣我還沒想到解釋,就放生了。
至於一般的線性算子要怎麼算偽反矩陣我沒去研究,最近沒空。
但是數列到數列、函數到函數的這種,一定要考慮邊界行為。
跟提問「iD 是 Hermitian 算子嗎?」一樣要考慮邊界。
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 163.13.112.58 (臺灣)
※ 文章網址: https://www.ptt.cc/bbs/Math/M.1663081163.A.9A2.html
※ 編輯: Vulpix (163.13.112.58 臺灣), 09/13/2022 23:32:14
→
09/14 00:52,
3年前
, 1F
09/14 00:52, 1F
→
09/14 00:52,
3年前
, 2F
09/14 00:52, 2F
→
09/14 00:55,
3年前
, 3F
09/14 00:55, 3F
→
09/14 00:55,
3年前
, 4F
09/14 00:55, 4F
→
09/14 00:55,
3年前
, 5F
09/14 00:55, 5F
推
09/14 01:10,
3年前
, 6F
09/14 01:10, 6F
→
09/14 01:10,
3年前
, 7F
09/14 01:10, 7F
推
09/14 01:16,
3年前
, 8F
09/14 01:16, 8F
→
09/14 01:16,
3年前
, 9F
09/14 01:16, 9F
→
09/14 01:16,
3年前
, 10F
09/14 01:16, 10F
→
09/14 01:16,
3年前
, 11F
09/14 01:16, 11F
→
09/14 01:16,
3年前
, 12F
09/14 01:16, 12F
→
09/14 01:16,
3年前
, 13F
09/14 01:16, 13F
→
09/14 01:16,
3年前
, 14F
09/14 01:16, 14F
→
09/14 01:16,
3年前
, 15F
09/14 01:16, 15F
→
09/14 05:59,
3年前
, 16F
09/14 05:59, 16F
→
09/14 05:59,
3年前
, 17F
09/14 05:59, 17F
推
09/14 11:04,
3年前
, 18F
09/14 11:04, 18F
※ 編輯: Vulpix (163.13.18.218 臺灣), 09/14/2022 13:20:44
→
09/16 01:21,
3年前
, 19F
09/16 01:21, 19F
推
09/16 10:36,
3年前
, 20F
09/16 10:36, 20F
其他看起來比較合理的下刀方式:
-1 1 0 0 0
0 -1 1 0 0
0 0 -1 1 0
0 0 0 -1 1
偽反矩陣:
-0.8 -0.6 -0.4 -0.2
0.2 -0.6 -0.4 -0.2
0.2 0.4 -0.4 -0.2
0.2 0.4 0.6 -0.2
0.2 0.4 0.6 0.8
這個應該下不同的條件就可以讓他看起來收斂到適當的 u+c 上吧。
→
09/16 18:35,
3年前
, 21F
09/16 18:35, 21F
※ 編輯: Vulpix (163.13.112.58 臺灣), 10/05/2022 02:45:37
討論串 (同標題文章)
完整討論串 (本文為第 6 之 6 篇):
分析
14
72