data shuju;
length ztl$8 ybl$2 zs$2 ybbl$12;
obs=_n_;
call symput("vars",compress(attrn(open("work.shuju",'i'),'nvars')));
input ztl$  ybl$ zs$ ybbl$ @@;
cards;
20/20/20 2  3  0.1/0.9/0.1
10/10/10 2  3  0.1/0.9/0.1
15/15/15 3  4  0.2/0.8/0.3
12/12/12 4  5  0.3/0.5/0.2
;
run;
proc transpose data=shuju out=shuju1;
var _all_;
run;
data _null_;
set shuju1;
nobs=_n_;
call symput("name"||left(nobs),_name_);
run;
%macro a(row);
data _null_;
set shuju;
%do i=1 %to %eval(&vars.-1);
call symput("n"||compress(&i.)||left(obs),&&name&i.);
%end;
%put &&n1&row. &&n2&row. &&n3&row. &&n4&row.;
run;
%mend a;
%a(1);
%a(2);

n11 n21  n31 n41
n12 n22  n32 n42
n13 n23 n33 n43
n14 n24 n34 n44
这个是宏变量矩阵

其实我的水平有限，没怎么看懂你的程序，不过还是有提示作用，我用自己能理解的方式把程序编出来了。