Calling a function inside another function with iteration

2020 2 17
1 个回答
回答已采纳
更新时间：2020 2 18
8 次查看（30 天）

0 个投票

Hello community,
I am trying to reproduce a graph, but I guess I have a problem calling a function. Could you help me?
clear all; clc; close all;
c=1;
Yi0=c;
t=linspace(0,10,101);
T=0.1;
k=4;
for ite=2:101
y(ite)=funcion(T,t(ite),Yi(ite-1),k); %(Maybe here I have an error?)
end
plot(t,y')
%%
function Yout=funcion(T,t,Yi0,k)
Yi=zeros(k,1);
Yi(1)=Yi0;
for itek=2:k
for n=1:(k-1)
Yi(itek)=(T/(itek-1))*(Yi(itek-1)-exp(t(ite))*(Yi(n)*Yi(itek-n)));
Yi(itek)=Yi(itek)+Yi(itek-1); %I want to use in my next iteration, for Yi(1)= to the sum of the previous Yi(itek), so I dont know how o put it
end
end
Yout=sum(Yi);
end
the graph should be like this:

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Adam Danz
Adam Danz 2020-2-17
编辑：Adam Danz 2020-2-17
The error message you received tells you the problem.
If you're having trouble understanding the error message, provide the full copy-pasted error message and let us know which part isn't clear.
Hi,
in your code there are several errors. Example in following line:
y(ite)=funcion(T,t(ite),Yi(ite-1),k); %(Maybe here I have an error?)
you don't initialize Yi before, thus matlab returns an error.
In following line
Yi(itek)=(T/(itek-1))*(Yi(itek-1)-exp(t(ite))*(Yi(n)*Yi(itek-n)));
when you evaluate exp(t(ite)), the variable "ite" is not defined inside function. Maybe it should be "itek" instead.
In general I advice you to choose more appropriate name variables, in this way you are less error prone.
I finish to run code after these two errors.
Do you have mathematical formulation of the curve that you want to plot? It will be easier check your code.

请先登录，再进行评论。

请先登录，再回答此问题。

 采纳的回答

Jose Sosa Lopez
Jose Sosa Lopez 2020-2-18
编辑：Jose Sosa Lopez 2020-2-18

0 个投票

%https://www.hindawi.com/journals/aaa/2014/486509/#B22
%EMHPM
clear all; clc; close all;
c=1;
yi0=c; %Initial condition
t=linspace(0,10,101); %time interval=0.1
T=0.1;
k=10;
yi=zeros(101,1); yi(1)=yi0;
for itet=2:101
yi(itet)=funcion(0.1,t(itet),yi(itet-1),k);
end
plot(t,yi','--')
hold on
%%
% Function with ode45
t=linspace(0,10,101);
yinicial = 1;
[t,y] = ode45(@(t,y) y-(exp(t))*y^2, t, yinicial);
plot (t,y)
hold off
%%
function Yout=funcion(T,t,Yi0,k)
Yi=zeros(k,1);
Yi(1)=Yi0;
for itek=2:k
for n=1:(itek-1)
Yi(itek)=(T/(itek-1))*(Yi(itek-1)-exp(t)*(Yi(n)*Yi(itek-n)));
end
end
Yout=sum(Yi);
end

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Programming 的更多信息

产品

版本

R2019b

标签

Community Treasure Hunt

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

Start Hunting!

Translated by