clear all
close all

R=500;     % Resistenza
C=0.001;   % Capacita' 

% Definizione delle matrici del sistema
A=-1/(R*C);
B=-A;
C=1;
D=0;

% Creazione del sistema
S=ss(A,B,C,D);

% Guadagno in continua del sistema
dcgain(S)

% Risposta a un gradino unitario
figure(1)
step(S)

% Effetto passa-basso
T=0:0.01:50;
U=sin(0.1*T')+0.05*sin(50*T'); % Sinusoide rumorosa

figure(2)
plot(T,U)
xlabel('t')
ylabel('u(t)')
title('Sinusoide rumorosa')

y0=0;
Y=lsim(S,U,T,y0); % Sinusoide filtrata

figure(3)
plot(T,Y)
xlabel('t')
ylabel('y(t)')
title('Sinusoide filtrata')