a1(:,1)=ones(33,1); %a1的第一列全1，共33个 a1(:,2:8)=A(:,1:7); % a1的第2~第8列等于矩阵A的第1~第7列b1=B(:,1); % b1等于矩阵B的第1列，即是个列向量y11=a % y11等于a
