File: resamp2_crcf_interp_recreate_test.c

package info (click to toggle)
liquid-dsp 1.3.1-1
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 7,376 kB
  • sloc: ansic: 93,650; sh: 3,208; makefile: 1,171; python: 95; asm: 11
file content (93 lines) | stat: -rw-r--r-- 2,952 bytes parent folder | download | duplicates (5)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
//
// resamp2_crcf_interp_recreate_test.c
//
// Halfband interpolator (complex), recreate() example
//

#include <stdio.h>
#include <complex.h>
#include <math.h>

#include "liquid.h"

#define OUTPUT_FILENAME "resamp2_crcf_interp_recreate_test.m"

int main() {
    unsigned int m0=5;
    unsigned int m1=4;
    unsigned int h0_len = 4*m0+1; // initial filter length
    unsigned int h1_len = 4*m1+1; // filter length (after recreate)
    float fc=0.10f;
    unsigned int N=32;
    float As=60.0f;

    resamp2_cccf f = resamp2_cccf_create(h0_len,0.0f,As);

    resamp2_cccf_print(f);

    FILE*fid = fopen(OUTPUT_FILENAME,"w");
    fprintf(fid,"%% %s: auto-generated file\n",OUTPUT_FILENAME);
    fprintf(fid,"clear all;\nclose all;\n\n");
    fprintf(fid,"h0_len = %u;\n", h0_len);
    fprintf(fid,"h1_len = %u;\n", h1_len);
    fprintf(fid,"N      = %u;\n", N);

    unsigned int i;
    float theta=0.0f, dtheta=2*M_PI*fc;
    unsigned int ix=0, iy=0;
    float complex x, y[2];
    for (i=0; i<N/2; i++) {
        x = cexpf(_Complex_I*theta);
        theta += dtheta;

        resamp2_cccf_interp_execute(f, x, y);

        ix++;
        iy++;
        fprintf(fid,"x(%3u) = %8.4f + j*%8.4f;\n", ix+1,   crealf(x),    cimagf(x));
        fprintf(fid,"y(%3u) = %8.4f + j*%8.4f;\n", 2*iy+1, crealf(y[0]), cimagf(y[0]));
        fprintf(fid,"y(%3u) = %8.4f + j*%8.4f;\n", 2*iy+2, crealf(y[1]), cimagf(y[1]));

        printf("y(%3u) = %8.4f + j*%8.4f;\n", 2*iy+1, crealf(y[0]), cimagf(y[0]));
        printf("y(%3u) = %8.4f + j*%8.4f;\n", 2*iy+2, crealf(y[1]), cimagf(y[1]));
    }

    // re-create interpolator
    printf("recreating filter...\n");
    f = resamp2_cccf_recreate(f,h1_len,0.0f,As);

    // TODO: push additional values?


    for (i=0; i<N/2; i++) {
        x = cexpf(_Complex_I*theta);
        theta += dtheta;

        resamp2_cccf_interp_execute(f, x, y);

        ix++;
        iy++;
        fprintf(fid,"x(%3u) = %8.4f + j*%8.4f;\n", ix+1,   crealf(x),    cimagf(x));
        fprintf(fid,"y(%3u) = %8.4f + j*%8.4f;\n", 2*iy+1, crealf(y[0]), cimagf(y[0]));
        fprintf(fid,"y(%3u) = %8.4f + j*%8.4f;\n", 2*iy+2, crealf(y[1]), cimagf(y[1]));

        printf("y(%3u) = %8.4f + j*%8.4f;\n", 2*iy+1, crealf(y[0]), cimagf(y[0]));
        printf("y(%3u) = %8.4f + j*%8.4f;\n", 2*iy+2, crealf(y[1]), cimagf(y[1]));
    }


    fprintf(fid,"nfft=512;\n");
    fprintf(fid,"X=20*log10(abs(fftshift(fft(x.*hamming(length(x))',nfft))));\n");
    fprintf(fid,"Y=20*log10(abs(fftshift(fft(y.*hamming(length(y))',nfft))));\n");
    fprintf(fid,"f=[0:(nfft-1)]/nfft-0.5;\n");
    fprintf(fid,"figure; plot(f/2,X,'Color',[0.5 0.5 0.5],f,Y,'LineWidth',2);\n");
    fprintf(fid,"grid on;\nxlabel('normalized frequency');\nylabel('PSD [dB]');\n");
    fprintf(fid,"legend('original','interpolated',1);");

    fclose(fid);
    printf("results written to %s\n",OUTPUT_FILENAME);

    resamp2_cccf_destroy(f);
    printf("done.\n");
    return 0;
}