Question

如何把MATLAB的m文件与命令写一起

M文件中：
function y=ill(t,x)
a=1;b=0.3;
y=[a*x(1)*x(2)-b*x(1),-a*x(1)*x(2)]';

命令窗口中：
>> [t,x]=ode45('ill',[0:50],[0.02,0.98]);[t,x]
>> plot(t,x(:,1),t,x(:,2)),grid,pause

请问命令窗口中的语句是什么意思
怎么把命令窗口中的语句写到m文件中

Answer 1

>> [t,x]=ode45('ill',[0:50],[0.02,0.98]);
用ode45（四阶龙格库塔）法求解数值积分的库函数
ill是自定义函数的文件名。[0，50]积分区间，[0.02,0.98]是初始值
[t,x]
显示积分求解出的t值和x值
>> plot(t,x(:,1),t,x(:,2)),grid,pause
画出tx图，加网格，暂停
调用ode45接微分方程是不能写在同一个m文件中的。但是你可以先写函数ill保存为一个m文件，在写调用程序到另一个m文件，两个m文件在同一路径下就行了

Answer 2

用函数m文件是写不到一起的，用匿名函数吧
以下内容写道命令行或者保存在脚本m文件中运行就可以了

a=1;b=0.3;

ill=@(t,x) [a*x(1)*x(2)-b*x(1);-a*x(1)*x(2)]; %这一句代替了原来的函数m文件

[t,x]=ode45(ill,[0:50],[0.02,0.98]);[t,x]
plot(t,x(:,1),t,x(:,2)),grid,pause

Answer 3

命令窗口中能执行的语句和M文件中的是一样的啊，直接写在M文件就OK了，当然你要加上分号！
[t,x]=ode45('ill',[0:50],[0.02,0.98]);-----求解方程
[t,x]------显示数值
plot(t,x(:,1),t,x(:,2)),grid,pause------------画图、加上网格、保持
具体你可以查一下各个函数的用法，希望能帮到你

追问

直接放在一起的话
function y=ill(t,x)
a=1;b=0.3;
y=[a*x(1)*x(2)-b*x(1),-a*x(1)*x(2)]';
[t,x]=ode45('ill',[0:50],[0.02,0.98]);
plot(t,x(:,1),t,x(:,2))

这样出错：
>> ill
错误使用 ill (line 3)
输入参数的数目不足。