Dear Clement,
In your function Vb = vib + Vw*exp(t/Ta), for t=0, Vb = vib + 0.1=10.1; There are many solutions if you like to set the Vb at 10. Here you can find two possible solutions:
Option 1
clc;clear all;close all
vib = 10; Vw = 0.1;
t = [0 5 10 15 ];
Ta = [5];
for i=1:length(t)
if t(i)==0
Vb(i) = 10;
else
Vb(i) = vib + Vw*exp(t(i)/Ta);
end
end
plot(t,Vb)
Option 2
clc;clear all;close all
vib = 10; Vw = 0.1;
t = [0 5 10 15 ];
Ta = [5];
Vb = vib + Vw*exp(t/Ta);
Vb = Vb - Vw;
plot(t,Vb)
Kind regards, Christiaan
