DAC 的SNR matlab help！

XD0909

那位大大帮忙解释一下以下Code的意义。
. o0 u$ r- ?( L
clear all;

N=2^12;                         % The input sequence is made by 2^12 points
) d9 s+ C8 Q8 ^0 O% {/ G! cNbit=8;                         % 8-bit DAC
Ntransient=11;                  %
Ntot=N+Ntransient;
nper=311;                       % prime integer number of sine waveforms
Fs=1;
bw=0.5;
Fin=nper*Fs/N;                  % Input signal frequency (Fin=nper*Fs/N)
3 o& T2 w6 M% M* S, m, xVref=0.5;                       % Full scale for an input sinewave
Amp_dB=-0;                      % Amplitude in dB
$ g1 s: S9 Z( [% J, q. UAmpl=10^(Amp_dB/20)*Vref;       % Input signal amplitude
7 r8 ~( \3 H& s& O
; O7 j( d! X2 j2 I+ S0 Zfinrad=Fin*2*pi;
alphadx=0.0001;                 % gradient in the resistivity
+ \" o1 G/ g! A  w$ ?den=2^Nbit-1+alphadx*2^Nbit*(2^Nbit-1)/2;

  U: y+ ]* P9 g; kfor i=1:Ntot,
: U: n! N. C" F. f8 l; d    x=Vref*sin(finrad*(i-1))+Vref;      % Unipolar range
    k=round(x*2^Nbit);
    vout=(k+alphadx*k*(k+1)/2)/den;
    y(i)=vout;
" w/ a# |6 R# q9 q% E0 lend

y=y-mean(y)                            % Output with zero average
. A; x8 z8 G/ }0 {- R; ]
%--------------------------------Graphics---------------------------------%

9 }: {0 {: x  E5 cfigure(1);
clf;
2 y9 `$ F1 ?! f4 Q" e% s9 _+ }plot(y,'r');
grid;
4 d& Y6 Z1 g1 E$ G( Z( {xlim([1600 1900]);
. }) c1 n9 V$ _7 e$ Gtitle('Sinewave with distortion');

# n7 U# [  v# K$ U5 Zw=ones(1,N);
" a7 Z* X! J2 S) w. t+ a6 R  kf=Fin/Fs;                               % Normalized signal frequency
1 S7 Y4 k1 A, x: g* H2 NfB=N*(bw/Fs);                           % Base-band frequency bins
[snr,ptot]=calcSNR(y(1:N),f,fB,w,N);
3 ]: F1 z0 F& {ptot=ptot-max(ptot);                    % Normalize total spectrum

figure(2);
% R0 z4 B( D: X! T- }clf;
$ t# D! j- E: rplot(linspace(0,Fs/2,N/2), ptot(1:N/2), 'r');
8 Z6 F: {1 k( fgrid on;
title('PSD of the Output')
5 C+ Z( R" s$ I; \! ]' L! Nxlabel('Frequency [Hz]')
# U7 G8 M0 ^2 I/ q1 k8 W4 M8 ?ylabel('PSD [dB]')
axis([0 Fs/4 -120 0]);
text(Fs/10,-20,sprintf('\\bfSNR = %3.1fdB  @  \\alpha\\DeltaX = %1.0e\n'...
8 o) v: R+ P- t- m5 F, X7 m; ]    ,snr,alphadx),'Fontsize',12);
7 r  [# w# W  x2 y; h7 P. jtext(f*2,ptot((f*N+1)*2-1)+3,sprintf('\\bf%2.1fdB',ptot((f*N+1)*2-1)),...
) e. ^; n, M) \7 Y9 m    'Fontsize',8,'HorizontalAlignment','center');