%================================================================
%	nco.m is a program that generate N bits rom sin output
%	Copyright :blutea 2004_10_20
%   vionsion:  1.0
%================================================================
function nco_out=ncosinwave(n,m)

tsnum=4*2^n;
%---------產生一個周期和1/4周期的正弦信號---------
t=[0:tsnum-1]*2*pi/tsnum;
nco_out=sin(t);
nco_out=nco_out*(2^m-1);
%-----------沒有取整的輸出--------------------------
nco_out_prim=nco_out;
%------------取整的輸出------------------------------
nco_out=round(nco_out);
%---------一個周期的正弦信號-------------------------
nco_out_whole=nco_out;
%----------1/4周期的正弦信號-------------------------
nco_out=nco_out(1:tsnum/4);
%----------1/4周期組成的整周期正弦信號----------------
nco_out_quta=[nco_out,fliplr(nco_out),-nco_out,-fliplr(nco_out)];

%---------用FFT分析動態范圍---------
plotin=[nco_out_quta];
figure(101);
subplot(2,1,1),plot(plotin);grid on
title('nco time analyse');
xlabel('time ');
fft_num=tsnum;
nco_fft =fft(plotin,fft_num);
nco_fft =nco_fft(1:fft_num/2);
nco_fft = abs(nco_fft);
nco_fft=20*log10(nco_fft+1e-4);%由于nco_out的對稱性FFT后隔一位產生一個0
subplot(2,1,2),plot(nco_fft);grid on
title('Frequency content of nco');
xlabel('frequency');
ylabel('Magtitude(dB)');