Question

我希望将这两个文件组合成一个m.file，而不是使用函数来定义＆＃34; monod＆＃34;传递参数，我想使用传递参数的匿名函数。所以要实现两个文件的组合。

%(file 1)
function dcdt = monod(t,c,k,ks,y,b)
dcdt = zeros(2,1);
dcdt(1) = -k*c(2)*c(1)/(ks+c(1));
dcdt(2) = y*k*c(2)*c(1)/(ks+c(1))-b*c(2);
% (file 2) ODE45
k = 3.7;ks=30;y=0.03;b=0.01;
options = odeset('Reltol',1.e-10,'AbsTol',1.e-10);
[t,c] = ode45(@monod, [0,200],[200,1],options,k,ks,y,b);
% plot 
C = c(:,1);
Xa = c(:,2);
figure(1);
grid on;
subplot(2,1,1);
plot(t,C);
title('Substrate aqueous concentration vs time(ODE45)');
xlabel('time');ylabel('Substrate aqueous concentration C')
subplot(2,1,2);
plot(t,Xa);
title('Active-cell concentration vs time');
xlabel('time');ylabel('Active-cell concentration Xa(ODE45)');

Answer 1

您只需在文件中定义嵌套函数即可。完全可以接受的语法。但是，您需要使File＃2成为实际函数，因为您无法使用脚本文件定义嵌套函数。要做到这一点，只需使它成为一个不接受输入并且不返回任何内容的函数。：

function run_ode %// Change here

    %// Include monod function here - watch the end keyword
    function dcdt = monod(t,c,k,ks,y,b)
        dcdt = zeros(2,1);
        dcdt(1) = -k*c(2)*c(1)/(ks+c(1));
        dcdt(2) = y*k*c(2)*c(1)/(ks+c(1))-b*c(2);
    end %<----

%// Begin File #2
k = 3.7;ks=30;y=0.03;b=0.01;
options = odeset('Reltol',1.e-10,'AbsTol',1.e-10);
[t,c] = ode45(@monod, [0,200],[200,1],options,k,ks,y,b);
% plot 
C = c(:,1);
Xa = c(:,2);
figure(1);
grid on;
subplot(2,1,1);
plot(t,C);
title('Substrate aqueous concentration vs time(ODE45)');
xlabel('time');ylabel('Substrate aqueous concentration C')
subplot(2,1,2);
plot(t,Xa);
title('Active-cell concentration vs time');
xlabel('time');ylabel('Active-cell concentration Xa(ODE45)');

end %// Take note of this end too as we now have nested functions

将上述代码复制并粘贴到名为run_ode.m的文件中，然后在MATLAB命令提示符中输入run_ode并按 ENTER 。

>> run_ode

你应该得到你想要的结果。

或者，如果您想使用问题标题中引用的匿名函数，则可以改为：

%// Change here
monod = @(t,c,k,ks,y,b) [-k*c(2)*c(1)/(ks+c(1)); y*k*c(2)*c(1)/(ks+c(1))-b*c(2)];

k = 3.7;ks=30;y=0.03;b=0.01;
options = odeset('Reltol',1.e-10,'AbsTol',1.e-10);
[t,c] = ode45(monod, [0,200],[200,1],options,k,ks,y,b); %// Change here too
% plot 
C = c(:,1);
Xa = c(:,2);
figure(1);
grid on;
subplot(2,1,1);
plot(t,C);
title('Substrate aqueous concentration vs time(ODE45)');
xlabel('time');ylabel('Substrate aqueous concentration C')
subplot(2,1,2);
plot(t,Xa);
title('Active-cell concentration vs time');
xlabel('time');ylabel('Active-cell concentration Xa(ODE45)');

monod现在是一个匿名函数，它接收6个输入，并输出一个两元素列向量，以便在ode45中使用。请注意，ode45现已更改，以便删除@。 monod现在已经是匿名函数的句柄，因此不需要使用@。

在解决ODE时，不是在matlab中创建两个m.files，而是如何使用匿名函数传递参数？

1 个答案: