如何解决程序循环次数增加后出现Inf或NAN

huangxuemei513

:@( 在matlab矩阵计算
X=A/B
前300个循环一切正常
300以后突然发散成Inf或NAN
可能是甚么原因啊？
麻烦高手指点

[ 本帖最后由 eight 于 2007-7-23 12:50 编辑 ]

hector1982

能不能把程序写多点啊！Ａ　Ｂ　是什么啊？

huangxuemei513

AX=B
A是系数矩阵，B是列矩阵
附程序如下，望各位高手不吝赐教，小女子暂且先谢过：

1. clear
   
2. yy=1e-6;zz=1e-6;tt=1e-3;yita=4.6857e-7; R=3.3333e-5; seta=4.027e+6; B=0.5;rou=2450;
   
3. N=1;
   
4. A=5e-4;
   
5. m=200;n=500;
   
6. C=zeros(2*m+1,3*n); E=zeros(2*m+1,3*n);
   
7. F=zeros(2*m+1,3*n); H=zeros(2*m+1,3*n);
   
8. D=939140;G=-935140;
   
9. X=zeros(2*m+1,3*n)+eps;
   
10. QX=zeros(2*m+1,3*n)+eps;
    
11. U=ones(2*m+1,3*n)*eps;
    
12. QU=ones(2*m+1,3*n)*eps;
    
13. 
    
14. vy=zeros(2*m+1,3*n);
    
15. vz=zeros(2*m+1,3*n);
    
16. C=-[(N*vy)/(2*yy)]-yita/yy^2; E=[(N*vy)/(yy^2)]-yita/yy^2;
    
17. F=[(N*(vz+R))/(2*zz)]+yita/zz^2;
    
18. H=-[(N*(vz+R))/(2*zz)]+yita/zz^2;
    
19. 
    
20. 
    
21. A1=zeros(2*m+1);
    
22. for j=2:n-1
    
23. i=m;
    
24. while i*yy>800*(j*zz)^2
    
25. i=i-1;
    
26. end
    
27. 
    
28.                                                 
    
29. A1(1+m-i,1+m-i)=1;
    
30. A1(1+m+i,1+m+i)=1;
    
31. for k=2+m-i:m+i
    
32. A1(k,k-1)=C(k,k-1);A1(k,k)=D;A1(k,k+1)=E(k,k+1);
    
33. end
    
34. 
    
35. AA=zeros(2*i+1);
    
36. u=1;v=1;
    
37. for h=1+m-i:1+m+i
    
38. for g=1+m-i:1+m+i
    
39. AA(u,v)=A1(h,g);
    
40. v=v+1;
    
41. end
    
42. u=u+1;
    
43. v=1;
    
44. end
    
45. 
    
46. Y=zeros(2*m+1,1);
    
47. for f=2+m-i:m+i
    
48. Y(f,1)=F(f,j-1)*X(f,j)+G*X(f,j)+H(f,j+1)*X(f,j+1)-seta*(B^2)*X(f,j)/rou;
    
49. end
    
50. Y(1+m-i,1)=2*eps/yy^2;
    
51. Y(1+m+i,1)=2*eps/yy^2;
    
52. 
    
53. Y1=zeros(2*i+1,1);
    
54. u=1;
    
55. for h=1+m-i:1+m+i
    
56. Y1(u,1)=Y(h,1);
    
57. u=u+1;
    
58. end
    
59. 
    
60. X1=zeros(2*i+1,1);
    
61. X1=(AA+eps)\Y1;
    
62. 
    
63. Xj=zeros(2*m+1,1);
    
64. u=1;
    
65. for h=1+m-i:1+m+i
    
66. Xj(h,1)=X1(u,1);
    
67. u=u+1;
    
68. end
    
69. 
    
70. X(:,j)=Xj;
    
71. end                                          
    
72. X

复制代码

xjzuo

试试pinv.

huangxuemei513

改成
X1=pinv((AA+eps)\Y1);
运行出现：

??? Index exceeds matrix dimensions
Error in   ==> Xj(h,1)=X1(u,1);

pinv试了 不好用也:@(

[ 本帖最后由 ChaChing 于 2010-4-4 20:57 编辑 ]

eight

请看看置顶贴，然后找到“常见出错问题整理”的帖子，建议先阅读一下

[ 本帖最后由 ChaChing 于 2010-4-4 20:59 编辑 ]

lxlxlxx

回复 1 # huangxuemei513 的帖子

请问，您的问题得到解决了么？我现在也遇到了同样的问题，计算过程总是应为NaN终止，比如100次的迭代，一般到10次以下就终止了，望指点

liyukui89

huangxuemei513 发表于 2007-7-21 19:29
AX=B
A是系数矩阵，B是列矩阵
附程序如下，望各位高手不吝赐教，小女子暂且先谢过：

你好，你的问题解决了吗？我在分解循环的时候也遇到矩阵的最后一行是INF，望指点，谢谢！

chybeyond

liyukui89 发表于 2014-6-19 21:44
你好，你的问题解决了吗？我在分解循环的时候也遇到矩阵的最后一行是INF，望指点，谢谢！

参考4楼方法，试试pinv