radio95/fm95
0
1
Fork 0
mirror of https://github.com/radio95-rnt/fm95.git synced 2026-02-27 03:23:54 +01:00
Code Issues Activity

test

Browse Source
This commit is contained in:
KubaPro010
2025-01-01 02:06:26 +01:00
parent 0bf0cb9674
commit ed01827bcb
2 changed files with 19 additions and 30 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
lib/hilbert.c
View File

@@ -1,8 +1,24 @@
#include "hilbert.h"
static float generate_coefficient(int i, int nzeros) {
if (i == nzeros/2) return 0.0f;
float n = i - nzeros/2;
// Basic Hilbert transform coefficient
float h = (i == nzeros/2) ? 0.0f : (2.0f/(PI * n));
// Apply Blackman window for better stopband attenuation
float w = 0.42f - 0.5f * cosf(2.0f * PI * i / nzeros) +
0.08f * cosf(4.0f * PI * i / nzeros);
return h * w;
}
void init_hilbert(HilbertTransformer *hilbert) {
hilbert->delay = calloc(D_SIZE, sizeof(float));
hilbert->dptr = 0;
for(int i = 0; i < NZEROS; i++) {
hilbert->coeffs[i] = generate_coefficient(i, NZEROS);
}
}
void apply_hilbert(HilbertTransformer *hilbert, float sample, float *output_0deg, float *output_90deg) {
@@ -11,7 +27,7 @@ void apply_hilbert(HilbertTransformer *hilbert, float sample, float *output_0deg
hilbert->delay[hilbert->dptr] = sample;
hilb = 0.0f;
for(int i = 0; i < NZEROS/2; i++) {
hilb += (xcoeffs[i] * hilbert->delay[(hilbert->dptr - i*2) & (D_SIZE - 1)]);
hilb += (hilbert->coeffs[i] * hilbert->delay[(hilbert->dptr - i*2) & (D_SIZE - 1)]);
}
*output_0deg = hilbert->delay[(hilbert->dptr - 99) & (D_SIZE - 1)];