循环过程中出现输出NaN怎么办

循环过程中出现输出NaN怎么办,第1张

从两行对比看:

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等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!

欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/zz/10079311.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-05-05
下一篇 2023-05-05

发表评论

登录后才能评论

评论列表(0条)

保存