Mein MATLAB Forum

function pt1_demo
IC=[1];
T=[0 5];
R=1;
C=1;
t_u=linspace(0,5,100);
u=sin(2*pi*t_u);
%u(1:100)=2;  % konstant
% u(1:20)=0;  % Sprung
% u(21:100)=1;
[t,v] = ode45(@(t,v) pt1(t,v,t_u,u),T,IC);
plot(t,v);
figure;
plot(t_u,u);

function vp = pt1(t,v,t_u,u)
u=interp1(t_u,u,t);
vp=u-v/(R*C);
end
end