nco.zip_Code_nco MATLAB_NCO_NCO MATLAB_oscillator

clear all; close all; %cpu time %% NCO %parameter Fs = 8000; % samping frequency 8Mhz N = 8192; quaterwavesinLUTsize = N/4; % 2048, N = 8192 M = N/2; % capacity 4096 P = 16; % max 16 level level = 2^P; % 양자화 레벨 수 inc = 510; % 유효[1 4095] [0 1 8191] tic % LUT lut = quant (sin(2 * pi * (0 : (N/4) - 1) / N), 1/level); %룩업테이블을 엔트리 개수의 1/4로 줄임 addr = 0 : inc : inc * (N-1); addr = mod(addr, N); % x[N+n] = x[n] y = zeros(1,N); for i = 1 : N if addr(i) < (N* 1/4) y(i) = lut (addr(i) +1); elseif addr(i) < (N* 2/4) && addr(i) >= (N* 1/4) y(i) = lut ((N/4)-(addr(i)-(N/4))); % LUT 거꾸로 읽음 elseif addr(i) < (N* 3/4) && addr(i) >= (N* 2/4) y(i) = -lut ((addr(i)-(N* 2/4)+1)); % LUT에 - 값 else y(i) = -lut ((N/4)-(addr(i)-(N* 3/4))); % LUT 거꾸로 읽고 - end end %% DAC %sinc를 곱하자 sinc(x)=sin(pi*x)/pi*x %% channel %% ADC %% FFT F = fft(y,N); toc %% 그래프 그리기 subplot(2, 1, 1); plot(y(1, 1:1000), '.-'); %시간축 그래프 그리기 xlabel('x[n]'); ylabel('Amplitude'); grid on; hold; subplot(2,1,2); hold [P1, F1] = pwelch(y(1,:), hann(N), 0, N, Fs); plot(1,F(1:N)); plot(F1, 10*log10(P1), 'b.--'); %주파수축 그래프 그리기 axis([0 4096 -120 0]); grid xlabel('Hz') ylabel('dB')