jens
04e1b6fcbf
git-svn-id: http://moon:8086/svn/vhdl/trunk@215 cc03376c-175c-47c8-b038-4cd826a8556b
44 lines
1.2 KiB
Matlab
44 lines
1.2 KiB
Matlab
% Read data
|
|
function tb_results()
|
|
fa = 100.0;
|
|
R = TEXTREAD('ratio.txt');
|
|
N = 3;
|
|
M = 1;
|
|
fa_dec = fa/R
|
|
N_pfir = 63;
|
|
|
|
I_out_dec = TEXTREAD('i_dec.txt')';
|
|
Q_out_dec = TEXTREAD('q_dec.txt')';
|
|
I_out_fir = TEXTREAD('i_fir.txt')';
|
|
Q_out_fir = TEXTREAD('q_fir.txt')';
|
|
|
|
N_dec = length(I_out_dec)-2*N_pfir;
|
|
N_fft_dec = N_dec
|
|
N_fft2_dec = N_fft_dec/2;
|
|
|
|
fft_y2_fir = 2/N_fft_dec*abs(fft(I_out_fir(2*N_pfir+1:N_fft_dec+2*N_pfir).*hann(N_fft_dec)'));
|
|
|
|
pfir_b = sinc(0.2*(-(N_pfir-1)/2:1:(N_pfir)/2));
|
|
pfir_b = pfir_b/sum(pfir_b).*blackman(N_pfir)';
|
|
I_out_dec = I_out_dec;
|
|
Q_out_dec = Q_out_dec;
|
|
fft_y2_dec = 2/N_fft_dec*abs(fft(I_out_dec(2*N_pfir+1:N_fft_dec+2*N_pfir).*hann(N_fft_dec)'));
|
|
|
|
wavwrite(0.45*[I_out_fir'-mean(I_out_fir) Q_out_fir'-mean(Q_out_fir)], 8000, 16, 'iq_cfir.wav');
|
|
wavwrite(0.45*[I_out_dec'-mean(I_out_dec) Q_out_dec'-mean(Q_out_dec)], 8000, 16, 'iq.wav');
|
|
|
|
% Output
|
|
close all;
|
|
|
|
plot(1:length(I_out_dec), I_out_dec, 1:length(Q_out_dec), Q_out_dec); grid;
|
|
|
|
figure;
|
|
cic_plot(M, N, R);
|
|
|
|
figure;
|
|
plot(fa/(R*N_fft_dec)*(0:N_fft2_dec),20*log10(fft_y2_fir(1:N_fft2_dec+1)), '-', fa/(R*N_fft_dec)*(0:N_fft2_dec),20*log10(fft_y2_dec(1:N_fft2_dec+1)), '-');
|
|
legend('CFIR', 'no CFIR');
|
|
xlabel('MHz');
|
|
grid;
|
|
|