y(:,:,i)=sin(a(i)+(2*pi/N)*(k-1)*(p-1)/10);
s(:,k,:)=y;
w(p,:,:)=s;
等价于
y(i)=sin(a(i)+(2*pi/N)*(k-1)*(p-1)/10);%y是1×4的向量
s(:,k)=y;%s是4×(N+N/2)的矩阵
w(p,:,:)=s;%w是〔((N/2)+1)*10〕×4×(N+N/2)的矩阵;