how to do a loop and call a function
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
function [g1]=func2(L,D,P,r,c,a,t)
g1=(((3*c)/(P*r))*((((D*(t-a))/L))-(L/6)));
end
function [g2]=func3(L,D,P,r,c,Y,t)
g2=1-((c/P)*exp(((-3*D)/(r*L))*(t-Y)));
end
hi guys. can someone help me? i have this to funtion and i have to set all the variable L,D,P,r,c,a and Y as a input and i have to do a for loop for t from 0 to 100 using this eqution and have to set the solution when a<t<=Y the it will use the first function and when t>Y it will use the second function and when t<a the fuction will be g=0 .but i don't know how to do it. so please help me
2 Kommentare
Antworten (2)
David Hill
am 9 Feb. 2020
c=;
P=;
r=;
D=;
a=;
L=;
Y=;
t1=0:a;
t2=a+1:Y;
t3=Y+1:100;
g=[zeros(1,length(t1)),(((3*c)/(P*r))*((((D*(t1-a))/L))-(L/6))),1-((c/P)*exp(((-3*D)/(r*L))*(t3-Y)))];
t=0:100;
0 Kommentare
Stephen23
am 9 Feb. 2020
Bearbeitet: Stephen23
am 9 Feb. 2020
As described in the question, you will have to do something like this:
L = ???;
D = ???;
P = ???;
r = ???;
c = ???;
a = ???;
Y = ???;
t = 0:???:100;
%
N = numel(t);
g = zeros(1,N); % this already fulfills the t<a requirement.
for k = 1:N
if t(k)>Y
g(k) = func3(??????, Y, t(k));
elseif t(k)>a
g(k) = func2(??????, a, t(k));
end
end
I left a few trivial things for you to complete (you need to replace all of the ???).
Both functions will need to be saved somewhere where they are visible to this code.
6 Kommentare
Stephen23
am 10 Feb. 2020
"then what if i want to plot r from 0 to 1.42 and fix the t"
If you only want to vary the r value then do the same like I showed you for L.
If you want to vary multiple different values then use ndgrid to generate all permutations.
Siehe auch
Kategorien
Mehr zu Loops and Conditional Statements finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!