——————————————————
function practice
[t,x]=ode45(@fun1,[0,20],[100,40,6])
subplot(1,2,1)
plot(t,x(:,1),'-',t,x(:,2),'-.',t,x(:,3),':')
legend('x1(t)','x2(t)','x3(t)')
grid
subplot(1,2,2)
plot3(x(:,1),x(:,2),x(:,3))
grid
function f=fun1(t,x)
r1=1r2=0.5r3=0.6lambda1=0.1lambda2=0.02lambda3=0.06mu=0.1
f =[x(1)*(r1-lambda1*x(2))
x(2)*(-r2+lambda2*x(1)-mu*x(3))
x(3)*(-r3+lambda3*x(2))]
——————————————————
打开软件后点击File/New/Script后,会出现文件编辑窗口:在文件编辑窗口输入代码,
检查代码,确保没有错误后点击窗口中运行(Run)按钮,运行程序(如下图红色箭头所示按钮):
查看运行结果,由于利用plot函数绘制曲线,所以运行后出现如下红色箭头所指的二次曲线:
matlab软件的m文件包括两种,脚本文件和函数文件,具体的可学习教程,根据具体问题,编写出程序,保存为m文件;如是脚本文件,直接点击edit窗口的工具条“运行”,就会运行程序,如有错,会有提示出现在命令窗口,就可再修改、再运行。如是函数文件,要先定义输入参数,再调用函数名,如:[a,b]=fun(x,y),x,y就是输入参数,a,b是输出参数,就是运行的结果,运行后有错,就会有提示出现在命令窗中。这就是m文件的编写、运行、调试过程。欢迎分享,转载请注明来源:内存溢出
评论列表(0条)