用matlab对动态系统绘图.m

来自「对于一个动态系统dy/dt=-y(t)+x(t-m),m为一个时间滞后常量」· M 代码 · 共 13 行

% 对于一个动态系统dy/dt=-y(t)+x(t-m),m为一个时间滞后常量，x(t)、y(t)都是关于时间
% t的函数，其中x是一个正态分布的数据集合，可随机产生，这样情况下，如何将系统情况
% 用样本点表示出来
clc;;
rand('state',0);
% Author's email: [e-mail]
m=3;
fun=inline('-y+randn*(t>m)/100','t','y','flag','m');
tspan=[0,m+1];
y0=0;
[t,y]=ode45(fun,tspan,y0,[],m);
plot(t,y);