Fs = 20e3; % [Hz] Ts = 1/Fs; % [s] f0 = 1; % [Hz] V0 = 0.1; % Sweep Sine Amplitude [V] C = 2.7e-6; % Capacitance of the Load [F] Imax_rms = 0.57; % [A] f1 = 0.8*Imax_rms/(C*V0*20/sqrt(2))/2/pi; % [Hz] f1 = min(f1, Fs/4); % [Hz] t1 = 100; % [V] t = 0:Ts:t1; Vin = struct('time', t, ... 'signals', struct('values', V0*chirp(t, f0, t1, f1, 'logarithmic', -90)', ... 'dimensions', 1));