为什么我算的结果（P2）是Nan，分母也不是0呀

1. clc
   
2. clear
   
3. %%
   
4. 
   
5. %%
   
6. m=100;
   
7. n=100;
   
8. H=ones(m,n);
   
9. A=ones(m,n);
   
10. B=ones(m,n);
    
11. C=ones(m,n);
    
12. D=ones(m,n);
    
13. E=ones(m,n);
    
14. F=ones(m,n);
    
15. % Psi=5*1e-2;
    
16. % lamda=2*z/l;
    
17. d_lamda=2/n;
    
18. phi=2*pi;%任意转角15
    
19. epslion=linspace(-1,1,m);
    
20. theta=linspace(pi/2,3*pi/2,n);
    
21. d_phi=phi/m;
    
22. l=0.18;%m
    
23. d=l/0.8;%m
    
24. H=1+epslion.*cos(phi-theta)';
    
25. surf(theta,epslion,H);
    
26. %% 初始条件
    
27. % p0=2*ou*miu/Psi^2;
    
28. P=zeros(m,n);
    
29. for i=1:m
    
30.     for j=1:n
    
31.   epslion1(i,j)=epslion(i);
    
32.     end
    
33. end
    
34. for i=1:m
    
35.   theta1(i,:)=theta(1,:);
    
36. end
    
37. 
    
38. 
    
39. %% 边界条件
    
40. P(1,:)=ones(1,n)*0;%最上面第一行全部是1
    
41. P(m,:)=ones(1,n)*0;%最下面第一行全部是1
    
42. for i=1:m
    
43.     P(:,1)=0;%最左边第一列为1,
    
44.     P(:,n)=0;%最右边最后一列为1，
    
45. end
    
46. %% 压力分布计算
    
47. P2=ones(m,n);
    
48. for i=2:m-1
    
49.     for j=2:n-1
    
50. A(i,j)=((H(i+1,j)-H(i,j))/2)^3;
    
51. B(i,j)=((H(i,j)-H(i-1,j))/2)^3;
    
52. C(i,j)=(d/l*d_phi/d_lamda)^2*((H(i,j+1)-H(i,j))/2)^3;
    
53. D(i,j)=(d/l*d_phi/d_lamda)^2*((H(i,j)-H(i,j-1))/2)^3;
    
54. E(i,j)=A(i,j)+B(i,j)+C(i,j)+D(i,j);
    
55. F(i,j)=3*d_phi*(H(i+1,j)+H(i,j))/2-(H(i,j)-H(i-1,j))/2;  
    
56. 
    
57.     end
    
58. end
    
59. 
    
60. 
    
61. 
    

复制代码

上面的没问题，最后迭代求解出来的结果就是nan了，P2

1. <blockquote>maxt=1;

复制代码

迭代格式P2(i,j)=(A(i,j)*P(i+1,j)+B(i,j)*P(i-1,j)+C(i,j)*P(i,j+1)+D(i,j)*P(i,j-1)-F(i,j))/E(i,j)    的分母E也不是0，我感觉是不是数太小了导致的，希望指点迷津

跟你用啥版本没多大关系，按while循环开始前 k = 0; 计，当 k =  179, i =   51, j =   43 时，B(i,j) 是0，P(i-1,j) 是inf，这俩的乘积没法算，导致出现NaN，至于为啥，这得你自己回去检查。反正这时候P2矩阵里面的数的绝对值都大到夸张，多个都是1e+250以上，正常的迭代一般不会出现这么夸张的数值，我都怀疑它发散。

我用的版本是2020b

跟你用啥版本没多大关系，按while循环开始前 k = 0; 计，当 k =  179, i =   51, j =   43 时，B(i,j) 是0，P(i-1,j) 是inf，这俩的乘积没法算，导致出现NaN，至于为啥，这得你自己回去检查。反正这时候P2矩阵里面的数的绝对值都大到夸张，多个都是1e+250以上，正常的迭代一般不会出现这么夸张的数值，我都怀疑它发散。

TouAkira 发表于 2022-5-18 12:42
跟你用啥版本没多大关系，按while循环开始前 k = 0; 计，当 k =  179, i =   51, j =   43 时，B(i,j) 是0，P(i-1,j) 是inf，这俩的乘积没法算，导

多谢指点，我去检查一下

TouAkira 发表于 2022-5-18 12:42
跟你用啥版本没多大关系，按while循环开始前 k = 0; 计，当 k =  179, i =   51, j =   43 时，B(i,j) 是0 ...

谢谢老哥，我找到问题了，现在能出结果了（抱拳）