从两行对比看:
Hfy(j+1)=Hby(j+1)-dt(Ebx(j+1)-Ebx(j))(1/u0)(1/dz)
如果像你所说,只有Hfy(26)=NaN,问题可能出在Hby、Ebx,你应该把这两项也输出来检查一下。
也可以调试程序,h_loop:doj=1,s 里面的1,改成25,在后面一行加上断点,逐行调试,在调试过程中查看每一个变量的值,就能找到谁是Nan的来源了
很显然,你建立的接迭代格式不收敛,发散了。这要从你自己的算法上找原因,你可以找分母上的数,是否很接近于0,之类的。比如我在我的matlab平台下监控了Sij的分母上的数随着迭代次数的增加变化情况。迭代进行500多次后就已经出现了不合理的值,你可以试试下面的代码。
clear;
m = 32;
n = 128;
N0 =5;
X0 = rand(m,N0);
C0 = cov(transpose(X0));
[V0,D0] = eig(C0);
XX = rand(m,n);
ppp=zeros(1);
iii=1
for K = 1:n
X=XX(:,K);
lmda = 1/K;
alfa=V0X;
C0=(1-lmda)C0+lmdaXtranspose(X);
P=zeros(32);
for i=1:m
p = (alfa(i)^2-D0(i,i))/K;
P(i,i) =p;
end;
Q=zeros(32);
for i=1:m-1
for j=i+1:m
if j==i
Sij=0;
else
ppp(iii)=(D0(j,j)-D0(i,i)+P(j,j)-P(i,i));
iii=iii+1;
Sij = (alfa(i)alfa(j))/(D0(j,j)-D0(i,i)+P(j,j)-P(i,i))lmda;
end
Q(i,j)=Sij;
Q(j,i)=-Sij;
end
end
V0 = V0(eye(32)+Q);
D0 = D0+P;
T1=zeros(m);
for i=1:m
V1=V0(:,i);
T1(i,i)=-norm(V1);
end
V0=V0T1;
T2=zeros(m);
for i=1:m
V=V0(:,i);
T2(i,i)=norm(V)^2;
end
D0=D0T2;
end
figure
plot(1:iii-1,ppp)
figure
plot(1:800,ppp(1:800))
matlab默认循环次数是500次,matlab本身擅长矩阵运算,不擅长循环。而且由于计算机浮点系统运算特点,所有收敛的序列相加肯定会有一个特定的数,不会出现inf、nan。而根据题目中“循环次数增加后”说明楼主自己增加了matlab的循环次数,而如果matlab循环超过设置的次数,会显示循环超过限制。因而分析,原因很可能是:你加的序列不收敛。建议检查程序代码,是否有哪里出错;如果没错,建议该用C++试试,这个擅长循环。
以上就是关于循环过程中出现输出NaN怎么办全部的内容,包括:循环过程中出现输出NaN怎么办、Matlab编辑程序里出来NaN,这是怎么回事怎么修改、matlab运行中如何解决程序循环次数增加后出现Inf或NAN等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)