%?Test?script?for?ADC
%?-?Josh?Carnes?20061128
%?-?Sumig?Cao?20120731
%?-?Jesaint?201211218
clear?all;
%?SET?OPTIONS?（These?need?to?be?set）
%?********************************************
DNLgo=1;??????????????????????%?Execute?INL/DNL??[yes=1/no=0]

%?SET?VALUES?FOR?FFT?（These?need?to?be?set）
%?*********************************************
N=12;??????????????????????%?Resolution?of?converter
fs_predec=90e6;???????????%?Pre-decimated?clock?frequency?[Hz]
fres=fs_predec/8192;???????%?Desired?frequency?resolution?of?FFT?[Hz]?fres=fclk/2^N=fin/M??頻率分辨率
fin=fres*2621;?????????????%?Frequency?of?interest?[Hz]?（used?to?test?if?frequency?will?be?dead?on?a?spike）
???????????????????????????%?Choose?frequency?of?resolution?times?a?prime?number選擇頻率分辨率的倍數最好是個素數
Ndelay=0;??????????????????%?Number?of?samples?to?dump?at?beginning?of?transient
Vref=2^（N-1）;??????????????%?Positive?max?voltage?swing?for?voltage?output?（Simulink?Model）
???????????????????????????%?????or?2^（N-1）?for?digital?output
fundpnts=0;???????????????%?#?of?FFT?points?on?either?side?of?fundamental?to?include
???????????????????????????%???in?fundamental?power?due?to?non-coherence（非相干性）
num_harm=5;???????????????%?#?of?harmonics?to?find?and?use?in?THD/SNR?calculation

%?GET?DATA?（These?need?to?be?set）
%?******************************
M?=?importdata（‘90-28.79.csv‘）;
indata?=?M.data（:2）;
%dataa?=?round（indata*10）/10;
dataa?=?round（indata）;
%dataa?=?indata;
data=dataa（100:8291）;
datab=dataa（100:8291）;
%datac=datab+（2^N-1）/2;%把負的提升為正的
datac=datab+2^（N-1）;

%?CALCULATE?VALUES?FOR?FFT
%?***********************************
fs=fs_predec;??????????????%?Decimated?clock?frequency
Ts=1/fs;???????????????????%?Calculate?period
Nsamples=round（fs/fres）;???%?Determine?number?of?samples?required?for?desired?FFT
???????????????????????????%???frequency?resolution
frescalc=fs/Nsamples;??????%?Calculated?freq?resolution?based?on?calculated?number?of?samples


%?DATA?&?CONFIGURATION?CHECK
%?************************************
%????Check?if?there?are?enough?samples?in?the?data?for?the?desired
%??????FFT?frequency?resolution
if?length（data）?????disp（‘ERROR:?not?enough?samples?taken?in?transient‘）;
????disp（‘??increase?transient?sampling?time?above?‘）;
????disp（（Nsamples+Ndelay）*Ts）;
????disp（‘?you?are?short?by?this?number?of?samples‘）;
????disp（Nsamples+Ndelay-length（data））;
????return;
end

%????Check?if?frequency?of?interest?is?greater?than?the?Nyquist
%??????frequency?and?folds?back?into?the?baseband
if?fs?????disp（‘WARNING?Frequency?of?interest?is?folded?back?into?base?band‘）
end

%????Check?if?the?frequency?of?interest?sits?directly?on?top?of?an?FFT?bin
%??????to?prevent?spectral?leakage/skirting
if?fin/frescalc?~=?round（fin/frescalc）
????disp（‘ERROR?frequency?of?interest?will?exhibit?skirting?spectral?leakage‘）;
????disp（‘make?the?frequency?of?interest?a?prime?multiple?of?your?frequency?resolution‘）;
????re