%UNTITLED2 Summary of this function goes here
% Detailed explanation goes here
u=zeros(1,2511)
for t=1:1:11
u(t)=0.05*t
end
for t=12:2510
k1=0.01*u(t-1)*(1-0.1*u(t-1))+0.01*u(t-10)
k2=0.01*(u(t-1)+0.5*k1)*(1-0.1*(u(t-1)+0.5*k1))+0.01*(u(t-10)+0.5*k1)
k3=0.01*(u(t-1)+0.5*k2)*(1-0.1*(u(t-1)+0.5*k2))+0.01*(u(t-10)+0.5*k2)
k4=0.01*(u(t-1)+k3)*(1-0.1*(u(t-1)+k3))+0.01*(u(t-10)+0.5*k3)
u(t)=u(t-1)+(k1+2*k2+2*k3+k4)/6
for k=1:1:5
if t==500*k+10
u(t)=(0.5^(2*k))*u(t-1)
end
end
end
t=0:2510
plot(t,u)
end
保存名为函数名。
.m文件(即楼主所说的脚本文件)调用函数的方法与在控制台调用的方法完全一致。只是需要注意几个问题:
1、被调用的函数的文件,需与.m文件在同一路径下;
2、通常被调用函数的函数名与文件名需要一致。
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)