clc;
clear?all;
N=50;
M=6;
K=5;
J=100;
X=zeros（NJ）;
T=[];
P=randperm（N）;
for?j=1:1:J
????X（:j）=zeros（N1）;
????X（P（（1:K））j）=randn（K1）;
????T（:（（j-1）*N+1）:（j*N））=randn（MN）;
????y（:j）=T（:（（j-1）*N+1）:（j*N））*X（:j）;
end
for?i=1:N
?????sum=0;
????for?j=1:1:J
????????sum=s