如何用MATLAB产生wiener分布的噪声

lihaitao123

如何用MATLAB产生wiener分布的有界噪声

lihaitao123

再顶！
function [Tt,Wn]= Wiener_Process(mu,sigma,Lt)
global w
w=0.75;
dt=2*pi/w/256;inc=2*pi/w/256;
k=sigma*sqrt(dt);
p=0.5*(1+mu*sqrt(dt)/sigma);
q=1-p;sumx=0;X=[sumx];
for bt=0:inc:Lt;
    nx=inc/dt;
    ns=binornd(nx,p);nf=nx-ns;net=ns-nf;
    sumx=sumx+(net*k);X=[X sumx];
end
t=0:inc:Lt;t=[0 t];n=length(t);
[W]=zeros(n,2);
W(:,1)=t';
W(:,2)=X';
plot(W(:,1),W(:,2));
Tt=W(:,1);
Wn=W(:,2);
网上找到一个，是这样子的。但是我想把Wn放到ode45求解的状态方程中，运行显示维数不正确。急死个人