%macro brown_move(ranseed,loop,record);
%do t=1 %to &loop.;
data brown_movement_&t.;
do i =1 to &record.;
if i=1 then do ;
group=&t.;
x=ranuni(&ranseed.+&t.);
y=ranuni(&ranseed.+&t.+1);
output;
end;
else ;
if i>1 then do;
group=&t.;
if ranuni(&ranseed.+&t.+2)<0.5 then x=x+ranuni(&ranseed.+&t.+2);
else x=x-ranuni(&ranseed.+&t.+2);
if ranuni(&ranseed.+&t.+3)<0.5 then y=y+ranuni(&ranseed.+&t.+3);
else y=y-ranuni(&ranseed.+&t.+3);
output;
end;
end;
drop i;
run;

%if &t. = 1 %then %do ;
data brown_movement;
set brown_movement_&t.;
run;
%end;
%else %do;
data brown_movement;
set brown_movement
brown_movement_&t.;
run;
%end;
proc sql; drop table brown_movement_&t.;quit;
%end;

%mend brown_move;

[此贴子已经被作者于2007-1-20 20:00:26编辑过]

使用方法：

在程序中添加上面的代码，运行下面的代码就可以了。这程序会产生多个数据集。如果你想把所有数据集都放在一起的话，请跟贴，我帮你写代码。
在macro中，ranseed就是一个随机种子数，随便写什么都可以，前提是数字哦。

loop是你想重复多少次试验，按你的要求就是要10000了。

record是你想在一次试验中出现多少次观测值了。按你的要求就是100了

这个程序是按照等时间观测布朗运动的结果的，数据集中有group ，x ，y3个变量；

group就是试验的次数。

x，y就是小物体运动的瞬时坐标。同一group的小物体运动是连续的，不同group的运动是独立的。这个要注意哦，如果你要画图（像物理书上的那个图）就自己用excel画画就行了。

注意：重新运行程序时要注意清理work中的原有数据集。

[此贴子已经被作者于2007-1-20 20:08:05编辑过]