用excel 更简单。

randn('state',100)           % 产生随机态
T = 1; N = 500; dt = T/N;
dW = zeros(1,N);             % 存放位置
W = zeros(1,N);              % 为了加快运算速度

dW(1) = sqrt(dt)*randn;      % 循环前的初始化
W(1) = dW(1);                % W(0) = 0 不允许，所以首先置值
for j = 2:N
   dW(j) = sqrt(dt)*randn;   % 产生序列
   W(j) = W(j-1) + dW(j);
end

plot([0:dt:T],[0,W],'r-')    % 画图。哦耶。
xlabel('t','FontSize',16)
ylabel('W(t)','FontSize',16,'Rotation',0)

参考资料：http://zhidao.baidu.com/question/158825506.html

谢谢了，虽然好像不符合题意。我上面的程序[Tt,St]=Stock_Process(s0,mu,sigma,Lt)，每做一次模拟就得到一个Tt和St，那怎么做500次使它取均值了，怎么写循环了？你设个附件吧，我把币给你