% Owner: Pratik Desai
% www.wright.edu/~desai.22
% Digital Controls Labs
% Use this code to calculate zeta, wn , wd from the experimental response.
clc;
clear all;
close all;
xo=input('input x0:');
xn=input('input xn:');
to=input('input t0:');
tn=input('input tn:');
n=input('input n:');
zeta=(1/(2*pi*n))*log(xo/xn);
wd=(2*pi*n)/(tn-to);
wn=wd/(sqrt(1-zeta^2));