这个工作区怎么看,超出索引范围了(像图片里这样),接下来我应该怎么调试呢,真的麻烦了,还需要大家的帮助
问题相关代码,请勿粘贴截图
nx=10;
nt=20000;
tmax=10;r=1;
%计算步长
dx=r/(nx-1);
T=tmax/(nt-1);
.
.
.
.(省略了中间程序)
tshort=linspace(0,tmax,nt/100);
qxshort=zeros(nx,nt/100);
qyshort=zeros(nx,nt/100);
qzshort=zeros(nx,nt/100);
dqxshort=zeros(nx,nt/100);
dqyshort=zeros(nx,nt/100);
dqzshort=zeros(nx,nt/100);
vxshort=zeros(nx,nt/100);
vyshort=zeros(nx,nt/100);
vzshort=zeros(nx,nt/100);
dvxshort=zeros(nx,nt/100);
dvyshort=zeros(nx,nt/100);
dvzshort=zeros(nx,nt/100);
for j=1:nt/100
for i=1:nx
qxshort(i,j)=qx(i,j100); %Using true y(i,j)
qyshort(i,j)=qy(i,j100);
qzshort(i,j)=qz(i,j*100);
dqxshort(i,j)=dqx(i,j*100); %Using true dy(i,j)
dqyshort(i,j)=dqy(i,j*100);
dqzshort(i,j)=dqz(i,j*100);
vxshort(i,j)=vx(i,j*100); %Using true y(i,j)
vyshort(i,j)=vy(i,j*100);
vzshort(i,j)=vz(i,j*100);
end
end