Getting NaN after division of matrix by vector

Question

Maaz Madha 2020-3-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/513154-getting-nan-after-division-of-matrix-by-vector

So for my assignment I had to form a matrix A and divide it by a column vector b to get the answer needed and plot it. Problem is that when I try to plot my answeri t gives nAn for all the values and i dont know why.

%%Initialisation 1 
alpha=0.1;
beta=0.1;
dx=(2*pi)/100;
dy=dx;
dt=0.01; %%timeskip
nu=0.1;
H=2*pi;
L=4*pi;
T_end=5;
 
 
%Mesh 
n=round((L/dx)+1);. Converts rectangle to a mesh
m=round(((H/dy)+1));
t=round((T_end/dt)+1);
x=[0:dx:L]';
y=[0:dy:H]';
v=@(x,t) 5*(1-(x-2*pi)^2/(4*pi^2)).*cos(pi*t)*sin(x);  
T=@(x,y) 20*cos(x)*sin(y);
b=zeros(n*m,1);
A=zeros(n*m);
 
 
%A matrix
for i=1:n
    for k=1:t
        an(i,k)= v((i-1)*dx,dt)*dt/(2*dy) - beta*dt/dy^2;
        as(i,k)= -v((i-1)*dx,dt)*dt/(2*dy) - beta*dt/dy^2;
    end 
end 
 
ae = -dt*alpha/dx^2;
ap = 1+(2*alpha*dt/dx^2)+(2*beta*dt/dy^2);
aw=ae;
 
 
%T1
for i=2:n-1
    for j=2:m-2
        pointer=(j-1)*n+i;
        A(pointer,pointer)=ap;
        A(pointer,pointer-n)=as(i,k);
        A(pointer,pointer-1)=aw;
        A(pointer,pointer+1)=ae;
        A(pointer,pointer+n)=an(i,k);  
      
    end 
     
end 
 
%T2
for i=1
    for j=2:m-2
        pointer=(j-1)*n+i;
        A(pointer,pointer)=ap;
        A(pointer,pointer+n)=as(i,k);
        A(pointer,pointer-1)=aw;
        A(pointer,pointer+1)=ae;
        A(pointer,pointer-m+2)=an(i,k);
       
        
    end 
end 
 
%T3
 
for i=n
    for j=2:m-2
        pointer=(j-1)*n+i;
        A(pointer,pointer)=ap;
        A(pointer,pointer+m-2)=as(i,k);
        A(pointer,pointer-1)=aw;
        A(pointer,pointer+1)=ae;
        A(pointer,pointer+n)=an(i,k);
    end 
end 
 
%T4 
 
for i=2:n-1
    for j=1
        pointer=(j-1)*n+i;
        A(pointer,pointer)=ap;
        A(pointer,pointer+n)=as(i,k);
        A(pointer,pointer+1)=(aw+ae);
        A(pointer,pointer+n)=an(i,k);
    end 
end 
%for T5
for i=2:n-1
    for j=m-1
        pointer=(j-1)*n+i;    
        A(pointer,pointer)=ap;
        A(pointer,pointer-n)=as(i,k);
        A(pointer,pointer-1)=(aw+ae);
        A(pointer,pointer+n)=an(i,k);
    end 
end 
 
%T6 
 
for i=1
    for j=1
        pointer=(j-1)*n+i;
        A(pointer,pointer)=ap;
        A(pointer,pointer+n)=as(i,k);
        A(pointer,pointer+1)=(aw+ae);
        A(pointer,pointer+2)=an(i,k);
    end 
end 
        
%T7
 
for i=n
    for j=1
        pointer=(j-1)*n+i;
        A(pointer,pointer)=ap;
        A(pointer,pointer+m-2)=as(i,k);
        A(pointer,pointer+1)=(aw+ae);
        A(pointer,pointer+n)=an(i,k);
    end 
end 
 
%T8
 
for i=1
    for j=m-1
        pointer=(j-1)*n+i;
        A(pointer,pointer)=ap;
        A(pointer,pointer+n)=as(i,k);
        A(pointer,pointer-1)=(aw+ae);
        A(pointer,pointer+2-m)=an(i,k);
    end 
end 
 
%t9
 
for i=n
    for j=m-1
        pointer=(j-1)*n+i;
        A(pointer,pointer)=ap;
        A(pointer,pointer+m-2)=as(i,k);
        A(pointer,pointer-n)=(aw+ae);
        A(pointer,pointer+n)=an(i,k);
    end
   
end 
 %matrix A produces 2301*2301
%%vector B
 
for i=1:n  
    for j=1:m
        pointer=(j-1)*n+i;
        b(pointer)=T((i-1)*dx,(j-1)*dy);
    end 
end 
 %%produces 2301*1
Temp=A/b