A very simple integration method

조회 수: 2 (최근 30일)
Yara Saavedra Ortega
Yara Saavedra Ortega 2020년 5월 16일
답변: Devineni Aslesha 2020년 5월 18일
So, I am trying to program a function which creates an array that contains a table with the value of four variables and its derivatives through some time instants. Here's what I've tried. However, it says one of my functions needs additional arguments.
function [integ] = sinteg(P,IC,t,dt)
%This function tries to employ a very simple integration method in order to
%get biomass, Carbon,nitrogen and product concentration.
% P represents a vector including the set of parameters necessary to
% calculate concentrations through the selected interval of time. CI
% represents another vector which contains a set of initial conditions
% for the process. t represents the total lenght of the interval of time
% that wishes to be evaluated. delta t corresponds to the size of the
% finite but next to infinitesimal time step.
inv=t/dt;
integ=zeros(inv,9);
integ(1,1)=0;
integ(1,2)= IC(1);
integ(1,3)= IC(2);
integ(1,4)= IC(3);
integ(1,5)= IC(4);
C=[IC(1),IC(2),IC(3),IC(4)];
aux=de_sys(C,P);
integ(1,6)=aux(1);
integ(1,7)=aux(2);
integ(1,8)=aux(3);
integ(1,9)=aux(4);
i=2;
while i<= t/dt
integ(i,1)=integ(i-1,1)+dt;
for j=1:1:5
integ(i,j)=integ(i-1,j)+integ(i-1,j+4)*dt;
end
aux=de_sys([integ(i,2),integ(i,3),integ(i,4)],P);
for j=6:1:9
integ(i,j)=aux(j-5);
end
i=i+1;
end
Here it is my other function
%This function calculates the value of each derivative for a specific point in time. It needs present conditions and a vector of parameters for the said point in time.
function dDV_dIV = de_sys(C,P)
mumax=P(1);
alpha=P(2);
Ksc=P(3);
Snm=P(4);
a1=P(5);
yxpsc=P(6);
mc=P(7);
Ksn=P(8);
Scm=P(9);
a2=P(10);
yxsn=P(11);
mn=P(12);
X=C(1);
Sc=C(2);
Sn=C(3);
dDV_dIV=zeros(1,4);
dDV_dIV(1)=(mumax*((Sc/(Ksc+Sc))*(Sn/(Ksn+Sn))*(1-(Sc/Scm)^a1)*(1-(Sn/Snm)^a2)))*X;
mu=(mumax*((Sc/(Ksc+Sc))*(Sn/(Ksn+Sn))*(1-(Sc/Scm)^a1)*(1-(Sn/Snm)^a2)));
dDV_dIV(2)=-(((1/yxsn)*mu)+mn)*X;
dDV_dIV(3)=-(((1/yxpsc)*mu)+mc)*X;
dDV_dIV(4)=(alpha*mu)*X;
end
  댓글 수: 2
Yara Saavedra Ortega
Yara Saavedra Ortega 2020년 5월 16일
Here's what the error says:
>> sinteg
Not enough input arguments.
Error in sinteg (line 11)
in=(t/dt);
Walter Roberson
Walter Roberson 2020년 5월 16일
When you run sinteg that way, what value should it use for t ? Are you expecting that it will look in the current workspace (the base workspace) to find a variable named t and to use that ?

댓글을 달려면 로그인하십시오.

채택된 답변

Devineni Aslesha
Devineni Aslesha 2020년 5월 18일
The error is due to calling the function sinteg without passing any input arguments. Call the function sinteg as shown below.
P = 3*ones(1,12);
IC = [2 1 2 1];
t = 5;
dt = 1;
integ = sinteg(P,IC],t,dt);
For more information, refer the following link.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Search Path에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by