>> zkd = rand(30,3) + i*rand(30,3);
>> R = sum(zkd * zkd');
>> size(R)
ans =
1 30
>> R
R =
Columns 1 through 2
43.7252707339352 - 14.3845494385565i 37.9012903137849 - 1.21833325457709i
Columns 3 through 4
44.0278780546839 - 5.91428978829146i 40.0850287858898 - 10.9783386248396i
The exact values are not important, but it is clear that all of the values are not the same.
If your zkd is 1000 by 5 then zkd*zkd' is going to be 1000 by 1000, and sum(zkd*zkd') is going to be 1 by 1000, in contradiction to your indication that M is 1000 by 1000.
Best Answer