s=tf('s');						% Variabile s di Laplace
gs=25/((s+1)^2+4);			% Funzione di trasferimento
%
figure(1); clf					% Si apre una nuova figura
Tfin=4;							% Durata temporale della simulaizone 
[y,t]=step(gs,Tfin);			% Si calcola la risposta al gradino del sistema gs
plot(t,y); hold on			% Tale risposta viene graficata
T=0.1;							% Si definisce un periodo di campionamento 
gz1=c2d(gs,T,'zoh');			% Discretizzazione con ricostr. di ordine 0
[y,t]=step(gz1,Tfin);		% Risposta al gradino del sistema discreto  
stairs(t,y,'k')				% Graficazione costante a tratti
T=0.4;							% Nuovo periodo di campionamento 
gz2=c2d(gs,T,'zoh');			% Nuova discretizzazione 
[y,t]=step(gz2,Tfin);		% Risposta al gradino
stairs(t,y,'r')				% Graficazione
xlabel('Tempo (s)')
ylabel('Funzioni gs, gz1 e gz2')
title('Campionamento ideale (T=0.1 e T=0.4) della funzione gs=25/((s+1)^2+4)')
clc; echo on 
% Funzione di trasferimento gs originaria
gs
% Funzione di trasferimento gz1 campionata: T=0.1
gz1
% Funzione di trasferimento gz2 campionata: T=0.4
gz2
%
% Si noti che negli istanti di campionamento
% il segnale originazio e quelli campionati coincidono
%
echo off
