模糊函数归一化多普勒频移频率

最近在复现一篇算法优化序列设计与雷达模糊函数相结合的文章
文章中模糊函数的定义如下

其中Rn为相关函数矩阵，Zk是大小为[N M]的矩阵，即M条长为N的序列，通过定义一个矩阵J相乘后得到的相关性矩阵Rn

文章中将多普勒频移归一化为[-0.2 0.2]rad ，均匀离散分为L份，即

请问图一中模糊函数表达式中的θ在matlab中我应该怎么处理，以及这个模糊函数我是应该Rn一项一项乘频移exp（jkθ）求和吗？还是怎么处理呢？

L赋个实际的数作为长度，theta就＝linespace(-0.2 0.2 L)，取L个delta theta的意思；
模糊函数有没有发现非常像DFT的定义式，你可以用for循环分别计算k层R(n)的k次方和指数因子的乘积，并累加求和，或直接用matlab的求和函数（一时半会想不起来叫啥了，好像是series吧，查查mathworks官方文档，我手机打字也不太方便）编写定义式，可以找找matlab实现DFT的自定义函数例程，照葫芦画瓢改一下。

补充一下优化的目标函数