Getting inf in matrix?

4 次查看（过去 30 天）
Desiree
Desiree 2020-4-19
评论： Desiree 2020-4-19
I wrote a code where I’m implementing an algorithm. The code works, but the values in my Z matrix are all inf and I don’t know why. I added a picture so you can see it. r in R^d should a sample drawn uniformly at random from the Fourier transformation of the Gaussian kernel. Gamma is just numbers uniformly at random in (0,2pi). cos(...) is at maximum 1, so I don’t understand why I get inf. Help is appreciated!
n=2000;
d=1000;
s=50;
Gamma=10;
S=randn(n,s);
U=qr(randn(s,d));
F=randn(n,d);
D=zeros(s,s);
for i=1:s
D(i,i)=1-(i-1)/d;
end
X = S*D*U + F/Gamma;
% This is the data matrix
-----------------------------------
function [G_2] = rnca(A,m)
[n,d]=size(A);
G=zeros(n,n);
for k=1:n
for j=1:n
G(k,j)=(1/2*pi)^(d/2)*exp(-norm(A(k,:)-A(j,:))^2/2);
end
end
Z=zeros(n,m);
gamma=0+(2*pi-0)*rand(1,m);
D=fft(G);
r=D(randsample(size(D,1),m),1:d);
for i=1:n
for l=1:m
Z(i,l)=sqrt(2/m)*cos(dot(r(l,:),A(i,:))+gamma(l));
end
end
G_2=Z*Z';
end
  7 个评论
显示 5更早的评论
Desiree
Desiree 2020-4-19
Now I discovered that I missed () when defining G matrix as I forgot to put () on 2*pi.
Now indeed the fft(G) reproduces lots of 0 values and Z ain't inf anymore. So I guess now it is fine?

请先登录，再进行评论。

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by