samplerate=150e6/(26*16); bits=16; settleoffset=8; bipolar=1; adc=importdata('myFile1.txt'); twos=ones(4096,1); if (bipolar==1) for i=settleoffset:length(adc) if (adc(i)<32768) twos(i-settleoffset+1)=adc(i)+32768; elseif(adc(i)>32768) twos(i-settleoffset+1)=adc(i)-32768; end end else twos=adc; end adcdata=twos/65535; M = length(adcdata); w = blackman(M); yw=w.*adcdata; x=20*log10((1/2^0.5)./(abs(fft(yw))).^2); %plot(0:samplerate/M:(samplerate/2-samplerate/(M)),x(1:M/2),'g') figure(1) snr(adcdata,samplerate,6) %snr %figure(2) %thd(adcdata,samplerate) %sinad(adcdata,samplerate)