版上的前辈们好，最近在做deconvolution时用到贝氏迭代法
https://i.imgur.com/4sNPcIy.jpg
但我天真的用了四个for loop 想做迭代到1000次求R，但因为a、W、R都是上百行列的矩阵，发现计算速度实在太慢了，爬了一些跟加速有关的文但没有很能理解，想请教有没有改善计算速度的方式？
如果能顺利计算的话可以站内信谈报酬(; Д `)
题外话：若可以使用快速傅立叶转换进行反卷积的话也可以站内信谈报酬><
谢谢matlab版！

矩阵向量相乘用内建语法就好啦，不用跑for-loop