版上的前輩們好，最近在做deconvolution時用到貝氏迭代法
https://i.imgur.com/4sNPcIy.jpg
但我天真的用了四個for loop 想做迭代到1000次求R，但因為a、W、R都是上百行列的矩陣，發現計算速度實在太慢了，爬了一些跟加速有關的文但沒有很能理解，想請教有沒有改善計算速度的方式？
如果能順利計算的話可以站內信談報酬(; Д `)
題外話：若可以使用快速傅立葉轉換進行反卷積的話也可以站內信談報酬><
謝謝matlab版！

矩陣向量相乘用內建語法就好啦，不用跑for-loop