Clean unused files
This commit is contained in:
F5OEO
parent
264dd932c9
commit
e1d90e5a29
BIN
opera/decode
BIN
opera/decode
Binary file not shown.
BIN
opera/opera
BIN
opera/opera
Binary file not shown.
BIN
opera/opera.rfa
BIN
opera/opera.rfa
Binary file not shown.
BIN
opera/opera2
BIN
opera/opera2
Binary file not shown.
531
opera/opera2.c
531
opera/opera2.c
|
|
@ -1,531 +0,0 @@
|
|||
//******************************************************************************
|
||||
// OPERA_Coding_Test.cpp : Defines the entry point for the console application.
|
||||
//
|
||||
// Purpose : Algorithm testing by a laptop computer before implementation
|
||||
// into PIC micro processor.
|
||||
//
|
||||
// Usage: "OPERA_Coding_Test [? | d | s] ["callsign"]";
|
||||
// where : s = help, d = debug and s = sample";
|
||||
//
|
||||
// Version 1.0.4, 2015/11/29: Modified print format
|
||||
// Version 1.0.3, 2015/11/13: Delete an additional start bit for decoder
|
||||
// Version 1.0.2, 2015/11/11: Add Visual C++ directives
|
||||
// Version 1.0.1, 2015/11/10: Changed help message
|
||||
// Version 1.0.0, 2015/11/07: Initial Release
|
||||
//
|
||||
// Copyright(C) 2015 F4GCB
|
||||
// Partial copyright (C)2015 7L1RLL
|
||||
//
|
||||
// Environment : Microsoft Visual C++ 2010 Express under Windows 10.
|
||||
// Compiler version 10.0.40219.1 SP1Rel
|
||||
//
|
||||
// Acknowledgement :
|
||||
// 1)Portion of this OPERA is derived from the work of EA5HVK.
|
||||
// 2)All functions of this program are a copy of JUMA-TX500/136
|
||||
// Transmitter Controller which written by F4GCB.
|
||||
//******************************************************************************
|
||||
|
||||
//#include "stdafx.h"
|
||||
//#include <locale.h>
|
||||
//#include <tchar.h>
|
||||
#include "stdio.h"
|
||||
#include "string.h"
|
||||
#include <stdlib.h>
|
||||
#include <sys/types.h>
|
||||
#include <sys/stat.h>
|
||||
#include <fcntl.h>
|
||||
#include "stdint.h"
|
||||
#include "math.h"
|
||||
|
||||
//#define __VCpp__ TRUE
|
||||
|
||||
// Grobal Variables
|
||||
|
||||
int FileFreqTiming;
|
||||
// Test program using SNDFILE
|
||||
// see http://www.mega-nerd.com/libsndfile/api.html for API
|
||||
|
||||
void WriteTone(double Frequency,uint32_t Timing)
|
||||
{
|
||||
typedef struct {
|
||||
double Frequency;
|
||||
uint32_t WaitForThisSample;
|
||||
} samplerf_t;
|
||||
samplerf_t RfSample;
|
||||
|
||||
RfSample.Frequency=Frequency;
|
||||
RfSample.WaitForThisSample=Timing; //en 100 de nanosecond
|
||||
//printf("Freq =%f Timing=%d\n",RfSample.Frequency,RfSample.WaitForThisSample);
|
||||
if (write(FileFreqTiming,&RfSample,sizeof(samplerf_t)) != sizeof(samplerf_t)) {
|
||||
fprintf(stderr, "Unable to write sample\n");
|
||||
}
|
||||
|
||||
}
|
||||
|
||||
static const short int pseudo_sequence[51] = {
|
||||
1,1,1,0,0,0,0,1,0,1, 0,1,0,1,1,1,1,1,1,0, 0,1,1,0,1,1,0,1,0,0, 0,0,0,0,0,1,1,0,0,1,
|
||||
0,0,0,1,0,1,0,1,0,1, 1
|
||||
};
|
||||
static short int walsh_matrix[8][7] = {
|
||||
{0,0,0,0,0,0,0},{1,0,1,0,1,0,1},{0,1,1,0,0,1,1},{1,1,0,0,1,1,0},
|
||||
{0,0,0,1,1,1,1},{1,0,1,1,0,1,0},{0,1,1,1,1,0,0},{1,1,0,1,0,0,1}
|
||||
};
|
||||
static short int symbol[239];
|
||||
char call[7], call_coded[45], vector[52];
|
||||
short int vector_to_tx[51];
|
||||
short int symbol_interleaving[119], symbol_coding[119];
|
||||
short int DEBUG = 0;
|
||||
const char sampleCall[7] = "AA1AA";
|
||||
|
||||
// Declaration of functions
|
||||
|
||||
void genn_opera(void);
|
||||
void generate_call(char call[7], char call_coded[45]);
|
||||
void add_crc16(char call_coded[45], char vector[52]);
|
||||
void scramble(char vector[52], short int symbol_coding[119]);
|
||||
void Walsh_Hammered_code(short int symbol_coding[119], short int vector_to_tx[51]);
|
||||
void interleave(short int vector_to_tx[51], short int symbol_interleaving[119]);
|
||||
void ManchesterEncode(short int symbol_interleaving[119], short int symbol[239]);
|
||||
char chr_norm_opera(char bc);
|
||||
void print_help(void);
|
||||
void print_short_int(const char caption[], short int code[239], int length);
|
||||
void print_str(const char caption[250], char code[52]);
|
||||
void strcpy_w(char s1[52], char s2[52], int length);
|
||||
void strcat_w(char s1[52], char s2[52], int lenS1, int lenS2);
|
||||
void encodepitx(short int *code, int length,float Nop);
|
||||
|
||||
#ifdef __VCpp__
|
||||
//**********************************
|
||||
// main forVisual C++
|
||||
int _tmain(int argc, _TCHAR* argv[])
|
||||
//**********************************
|
||||
{
|
||||
_tsetlocale(LC_ALL, _T("")); //Change Unicode to OS-Default locale
|
||||
#else
|
||||
int main(int argc, char* argv[])
|
||||
//**********************************
|
||||
{
|
||||
#endif
|
||||
int i = 0;
|
||||
char s1 = 0x00, s2[7] = "";
|
||||
|
||||
|
||||
switch (argc)
|
||||
{
|
||||
case 1 : // Help required
|
||||
case 2 : // Help required
|
||||
{
|
||||
printf("Usage : %s CALLSIGN file.rfa \n", argv[0]);
|
||||
//print_help();
|
||||
return 0;
|
||||
}
|
||||
case 3: // 2 arguments
|
||||
{
|
||||
s1 = (char)argv[1][0];
|
||||
if ((s1 == 's') && (argv[1][1] == 0x00)) // sample callsign
|
||||
{
|
||||
printf("%s\n","A sample callsign is used.");
|
||||
DEBUG = 0; i = 0;
|
||||
while (sampleCall[i] != 0 && i < 7)
|
||||
{
|
||||
call[i] = sampleCall[i]; call[++i] = 0x00;
|
||||
}
|
||||
genn_opera();
|
||||
return 0;
|
||||
}
|
||||
else if ((s1 =='?') && (argv[1][1] == 0x00))
|
||||
{
|
||||
printf("%s\n","Help requested");
|
||||
print_help();
|
||||
return 0;
|
||||
}
|
||||
else
|
||||
{
|
||||
// range check
|
||||
if (!((argv[1][0] >= '0' && argv[1][0] <= '9') || (argv[1][0] >= 'A' && argv[1][0] <= 'Z') ||
|
||||
(argv[1][0] >= 'a' && argv[1][0] <= 'z')))
|
||||
{
|
||||
printf("%s\n","Callsign must be began with an alphan/numeric character");
|
||||
print_help();
|
||||
return 0;
|
||||
}
|
||||
DEBUG = 0; i = 0;
|
||||
while (argv[1][i] != 0 && i < 7)
|
||||
{
|
||||
call[i] = argv[1][i]; call[++i] = 0x00;
|
||||
}
|
||||
|
||||
FileFreqTiming = open( argv[2], O_WRONLY|O_CREAT, 0644);
|
||||
genn_opera();
|
||||
return 0;
|
||||
}
|
||||
}
|
||||
|
||||
default:
|
||||
{
|
||||
printf("%s\n", "Too many arguments.");
|
||||
print_help();
|
||||
break;
|
||||
}
|
||||
} // end of switch argc
|
||||
return 0;
|
||||
} // end of _tmain
|
||||
|
||||
//*******************
|
||||
void genn_opera(void)
|
||||
//*******************
|
||||
{
|
||||
printf("\nCallsign = %s\n", call);
|
||||
|
||||
generate_call(call, call_coded);
|
||||
if (DEBUG)
|
||||
print_str("call_coded =", call_coded);
|
||||
|
||||
add_crc16(call_coded, vector);
|
||||
|
||||
if (DEBUG)
|
||||
print_str("crc16 vector =", vector);
|
||||
|
||||
scramble(vector, vector_to_tx);
|
||||
if (DEBUG)
|
||||
print_short_int("vector_to_tx =", vector_to_tx, 44);
|
||||
|
||||
Walsh_Hammered_code(vector_to_tx, symbol_coding);
|
||||
if (DEBUG)
|
||||
print_short_int("symbol_coding =", symbol_coding, 119);
|
||||
|
||||
interleave(symbol_coding, symbol_interleaving);
|
||||
if (DEBUG)
|
||||
print_short_int("symbol_interleaving =", symbol_interleaving, 119);
|
||||
|
||||
ManchesterEncode(symbol_interleaving, symbol);
|
||||
print_short_int("symbol =", symbol, 239);
|
||||
encodepitx(symbol,239,1);
|
||||
|
||||
} // genn_opera
|
||||
|
||||
//****************************************************
|
||||
// Normalize characters space S..Z 0..9 in order 0..36
|
||||
char chr_norm_opera(char bc)
|
||||
//****************************************************
|
||||
{
|
||||
char cc = 0;
|
||||
|
||||
if (bc >= '0' && bc <= '9') cc = bc - '0' + 27;
|
||||
if (bc >= 'A' && bc <= 'Z') cc = bc - 'A' + 1;
|
||||
if (bc >= 'a' && bc <= 'z') cc = bc - 'a' + 1;
|
||||
if (bc == ' ') cc = 0;
|
||||
|
||||
return (cc);
|
||||
} // enf of chr_norm_opera
|
||||
|
||||
//**********************************************
|
||||
void generate_call(char *call, char *call_coded)
|
||||
//**********************************************
|
||||
{
|
||||
int i;
|
||||
unsigned long code_sum;
|
||||
|
||||
//the thired character must always be a number
|
||||
if (chr_norm_opera(call[2]) < 27)
|
||||
{
|
||||
for (i=5; i> 0; i--) call[i] = call[i-1];
|
||||
call[0]=' ';
|
||||
}
|
||||
|
||||
// the call must always have 6 characters
|
||||
for (i=strlen(call); i < 6; i++)
|
||||
call[i] = ' ';
|
||||
call[6] = 0x00;
|
||||
|
||||
if (DEBUG) printf("NormalizedCall=%s\n", call);
|
||||
code_sum = chr_norm_opera(call[0]);
|
||||
code_sum = code_sum * 36 + chr_norm_opera(call[1]) - 1;
|
||||
code_sum = code_sum * 10 + chr_norm_opera(call[2]) - 27;
|
||||
code_sum = code_sum * 27 + chr_norm_opera(call[3]);
|
||||
code_sum = code_sum * 27 + chr_norm_opera(call[4]);
|
||||
code_sum = code_sum * 27 + chr_norm_opera(call[5]);
|
||||
|
||||
if (DEBUG) printf("code_sum=%Lu\n", code_sum);
|
||||
|
||||
// merge coded callsign ino a string
|
||||
call_coded[28] = 0x00;
|
||||
call_coded[27] = (short int) ((code_sum & 0x01) + 0x30);
|
||||
for (i = 26; i >= 0; i--)
|
||||
{
|
||||
code_sum = code_sum >> 1;
|
||||
call_coded[i] = (short int)((code_sum & 0x01) + 0x30);
|
||||
}
|
||||
} // end of pack_callsign
|
||||
|
||||
//***************************************************
|
||||
void generate_crc(char *datas, char *crc, int length)
|
||||
//***************************************************
|
||||
{
|
||||
unsigned int i, j, k;
|
||||
char buffer[52]; //strlen(datas)];
|
||||
short int wcrc[16] = {0}, byte1 = 0, byte2 = 0;
|
||||
|
||||
#ifdef __VCpp__
|
||||
strcpy_s(buffer, 52, datas);// strcpy(buffer, datas);
|
||||
#else
|
||||
strcpy(buffer, datas);// strcpy(buffer, datas);
|
||||
#endif
|
||||
if (DEBUG)
|
||||
print_str("buffer_crc = ", buffer);
|
||||
|
||||
for (i = 0; i < strlen(datas); i++)
|
||||
{
|
||||
for (j = 0; j < 8; j++)
|
||||
{
|
||||
if (j > 0) buffer[i] = buffer[i] >> 1;
|
||||
byte1 = buffer[i] & 0x01;
|
||||
byte2 = byte1 ^ wcrc[0];
|
||||
wcrc[0] = byte2 ^ wcrc[1];
|
||||
for (k = 1; k < 13; k++)
|
||||
wcrc[k] = wcrc[k+1];
|
||||
wcrc[13] = byte2 ^ wcrc[14];
|
||||
wcrc[14] = wcrc[15];
|
||||
wcrc[15] = byte2;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// if msb byte crc = 0 then value at 27
|
||||
byte2 = 0;
|
||||
for (i = 0; i < 8; i++)
|
||||
#ifdef __VCpp__ // add for Visual C++ by 7L1RLL 11/07/2015
|
||||
byte2 = byte2 + (short) wcrc[i] * pow((double)2.0, (int)i);
|
||||
#else
|
||||
byte2 = byte2 + wcrc[i] * pow(2, i);
|
||||
#endif
|
||||
if (byte2 == 0) byte2 =27;
|
||||
|
||||
// if lsb byte crc = 0 then value at 43
|
||||
byte1 = 0;
|
||||
for (i = 8; i < 16; i++)
|
||||
#ifdef __VCpp__
|
||||
byte1 = byte1 + (short) (wcrc[i] * (double)pow(2.0, (int)i - 8)); // add cast by 7L1RLL
|
||||
#else
|
||||
byte1 = byte1 + (wcrc[i] * pow(2, i - 8));
|
||||
#endif
|
||||
if (byte1 == 0) byte1 = 43;
|
||||
if (DEBUG) printf("byte1 = %x, byte2 = %x\n", byte1, byte2);
|
||||
|
||||
// merge crc into a string
|
||||
for (i = 0; i < 8; i++)
|
||||
{
|
||||
if (i > 0) byte2 = byte2 >> 1;
|
||||
wcrc[7 - i] = byte2 & 0x01;
|
||||
|
||||
if ( i > 0) byte1 = byte1 >> 1;
|
||||
wcrc[15 - i] = byte1 & 0x01;
|
||||
}
|
||||
if (length > 16) length = 16;
|
||||
for (i = 16 - length; i < 16; i++)
|
||||
crc[i - (16 - length)] = wcrc[i] + 0x30;
|
||||
crc[length] = 0x00;
|
||||
} // end of genarate_crc
|
||||
|
||||
//*********************************************
|
||||
void add_crc16(char * call_coded, char *vector)
|
||||
//*********************************************
|
||||
{ // input: |28 bits|, output : |51 bits|
|
||||
char crc1[17] = "", crc2[4] = "";
|
||||
|
||||
#ifdef __VCpp__ // for wide character compiler
|
||||
char temp[52] = "";
|
||||
|
||||
_tsetlocale(LC_ALL, _T("")); //Change Unicode to OS-Default locale
|
||||
|
||||
if (DEBUG)
|
||||
print_str("call_coded in add CRC16 =", call_coded);
|
||||
strcpy_w(temp, call_coded, 28); // 28 bit
|
||||
if (DEBUG)
|
||||
print_str("temp in add CRC16=", temp);
|
||||
generate_crc(call_coded, crc1, 16);
|
||||
if (DEBUG)
|
||||
print_str("crc1 =", crc1);
|
||||
strcat_w(temp, crc1, 28, 16); // 28 + 16 = 44
|
||||
generate_crc( temp, crc2, 3);
|
||||
if (DEBUG)
|
||||
print_str("crc2 =", crc2);
|
||||
#else // PIC C compiler using ASCII
|
||||
// char crc1[17] = "", crc2[4] = "";
|
||||
generate_crc(call_coded, crc1, 16);
|
||||
if (DEBUG) printf("crc1 =%s\n", crc1);
|
||||
generate_crc(strcat(call_coded, crc1), crc2, 3);
|
||||
if (DEBUG) printf("crc2 =%s\n", crc2);
|
||||
#endif
|
||||
// |4 bits sync| + |28 bits call| + |19 bit crc|
|
||||
#ifdef __VCpp__
|
||||
strcpy_w(vector, "0000", 4); // 4
|
||||
strcat_w(vector, temp, 4, 44); // 4 + 44 = 48
|
||||
strcat_w(vector, crc2, 48, 3); // 48 + 3 = 51
|
||||
#else // not VC++ ex : PIC
|
||||
strcpy(vector, "0000"); // 4
|
||||
strcat(vector, call_coded); // 4 + 44 = 48
|
||||
strcat(vector, crc2); // 48 + 3 = 51
|
||||
#endif
|
||||
} // end of add_crc16
|
||||
|
||||
//**************************************************
|
||||
void scramble(char *vector, short int *vector_to_tx)
|
||||
//**************************************************
|
||||
{ // encoding |51 bits|
|
||||
int i=0;
|
||||
|
||||
for (i = 0; i < 51; i++)
|
||||
{
|
||||
vector_to_tx[i] = vector[i] & 0x01;
|
||||
// convert ASCII to binary
|
||||
vector_to_tx[i] = vector_to_tx[i] ^ pseudo_sequence[i];
|
||||
}
|
||||
} // end of scrambling
|
||||
|
||||
//*************************************************************************
|
||||
void Walsh_Hammered_code(short int *vector_to_tx, short int *symbol_coding)
|
||||
//*************************************************************************
|
||||
{ // order 8 walsh matrix codification : |119 bits|
|
||||
int data = 0, idx = 0, i = 0, j = 0;
|
||||
|
||||
for (i = 0; i < 51; i += 3)
|
||||
{
|
||||
data = 0;
|
||||
for (j = 0; j < 3; j++)
|
||||
data = data + (vector_to_tx[i + j] << (2 - j));
|
||||
for (j = 0; j < 7; j++)
|
||||
{
|
||||
symbol_coding[idx] = walsh_matrix[data][j];
|
||||
idx++;
|
||||
}
|
||||
}
|
||||
} //end of Walsh_Hammered_code
|
||||
|
||||
//***********************************************************************
|
||||
void interleave(short int *symbol_coding, short int *symbol_interleaving)
|
||||
//***********************************************************************
|
||||
{ // interleaving : |119 bits|
|
||||
int idx = 0, i = 0, j = 0;
|
||||
|
||||
idx = 0;
|
||||
for (i = 0; i < 7; i++)
|
||||
{
|
||||
for (j = i; j < 119; j += 7)
|
||||
{
|
||||
symbol_interleaving[idx]= symbol_coding[j];
|
||||
idx++;
|
||||
}
|
||||
}
|
||||
} // end of interleave
|
||||
|
||||
//**********************************************************************
|
||||
void ManchesterEncode(short int *symbol_interleaving, short int *symbol)
|
||||
//**********************************************************************
|
||||
{ // manchester codification : |11| + |238 bits| - |1 bit| modified by 7L1RLL 11/07/2015
|
||||
int idx = 0;
|
||||
int i = 0, j = 0;
|
||||
|
||||
symbol[0] = 1;
|
||||
for (i = 0; i < 119; i++)
|
||||
{
|
||||
if (symbol_interleaving[i] == 0)
|
||||
{
|
||||
symbol[idx + 1] = 1;
|
||||
symbol[idx + 2] = 0;
|
||||
}
|
||||
else
|
||||
{
|
||||
symbol[idx + 1] = 0;
|
||||
symbol[idx + 2] = 1;
|
||||
}
|
||||
idx += 2;
|
||||
}
|
||||
} // end of Manchester_encode
|
||||
|
||||
//***************
|
||||
void print_help()
|
||||
//***************
|
||||
{
|
||||
printf("%s\n","Usage : OPERA_Coding_Test [? | s | [d \"callsign\"]]");
|
||||
printf("%s\n"," Normal : OPERA_Coding_Test \"callsign\"");
|
||||
printf("%s\n"," Help : OPERA_Coding_Test ?");
|
||||
printf("%s\n"," Sample : OPERA_Coding_Test s");
|
||||
printf("%s\n"," Debug : OPERA_Coding_Test d \"callsign\"");
|
||||
printf("%s\n"," Callsign format shall be like \"AA1AAA\". ");
|
||||
printf("%s\n"," Third character mast be a numeric character(0..9).");
|
||||
} // end of help
|
||||
|
||||
//********************************************************************
|
||||
void print_short_int(const char *caption, short int *code, int length)
|
||||
//********************************************************************
|
||||
{ // This is service function for debugging
|
||||
int i = 0;
|
||||
|
||||
printf("%s\n", caption);
|
||||
for (i = 0; i < length; i++)
|
||||
{
|
||||
printf("%d", code[i]);
|
||||
if (((i+1) % 4) == 0) printf(" ");
|
||||
if (((i+1) % 40) == 0) printf("\n");
|
||||
}
|
||||
printf("\n");
|
||||
} // end fo print_short_int
|
||||
|
||||
|
||||
void encodepitx(short int *code, int length,float Nop)
|
||||
{
|
||||
int i = 0;
|
||||
int j=0;
|
||||
/*and each of the
|
||||
239 symbols are transmitted by keying the transmitter as CW on and off with a symbol
|
||||
rate of 0.256*n s/symbol, where n is the integer of operation mode OPn that corresponds
|
||||
with the Opera frequency recommendation: */
|
||||
//WriteTone(1*32767,1e6*(256*Nop));
|
||||
WriteTone(1*32767,1e6*(256*Nop));
|
||||
for (i = 0; i < length-1; i++)
|
||||
{
|
||||
//for(j=0;j<1000;j++)
|
||||
WriteTone(code[i]*32767,1e6*(256*Nop));
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
//********************************************************************
|
||||
void print_str(const char * caption, char * code)
|
||||
//********************************************************************
|
||||
{ // This is a service function for debugging
|
||||
int i = 0;
|
||||
|
||||
printf("%s\n", caption);
|
||||
for (i = 0; i < strlen(code); i++)
|
||||
{
|
||||
printf("%c", code[i]);
|
||||
if (((i + 1) % 4) == 0) printf(" ");
|
||||
if (((i + 1) % 40) == 0) printf("\n");
|
||||
}
|
||||
printf("\n");
|
||||
} // end fo print_str
|
||||
//*******************************************************
|
||||
void strcpy_w(char * s1, char * s2, int length)
|
||||
//*******************************************************
|
||||
{
|
||||
int i;
|
||||
|
||||
for (i = 0; i < length; i++)
|
||||
s1[i] = s2[i];
|
||||
s1[length] = 0x00;
|
||||
} // end of strcpy_w
|
||||
|
||||
//****************************************************************
|
||||
void strcat_w(char * s1, char * s2, int lenS1, int lenS2)
|
||||
//****************************************************************
|
||||
{
|
||||
int i;
|
||||
|
||||
for (i = 0; i < lenS2; i++)
|
||||
s1[i + lenS1] = s2[i];
|
||||
s1[lenS1 + lenS2]= 0x00;
|
||||
|
||||
} // end of strcat_w()
|
||||
//************** End of Program ********************************
|
||||
Binary file not shown.
121
opera/piopera.c
121
opera/piopera.c
|
|
@ -1,121 +0,0 @@
|
|||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
#include <unistd.h>
|
||||
#include <string.h>
|
||||
#include <errno.h>
|
||||
#include <stdarg.h>
|
||||
#include <stdint.h>
|
||||
#include <math.h>
|
||||
#include <time.h>
|
||||
#include <signal.h>
|
||||
#include <sys/types.h>
|
||||
#include <sys/stat.h>
|
||||
#include <fcntl.h>
|
||||
#include <sys/mman.h>
|
||||
#include <math.h>
|
||||
|
||||
#include <sndfile.h>
|
||||
|
||||
#define ln(x) (log(x)/log(2.718281828459045235f))
|
||||
#define BUFFER_LEN 1024*8
|
||||
|
||||
int FileFreqTiming;
|
||||
// Test program using SNDFILE
|
||||
// see http://www.mega-nerd.com/libsndfile/api.html for API
|
||||
|
||||
void WriteTone(double Frequency,uint32_t Timing)
|
||||
{
|
||||
typedef struct {
|
||||
double Frequency;
|
||||
uint32_t WaitForThisSample;
|
||||
} samplerf_t;
|
||||
samplerf_t RfSample;
|
||||
|
||||
RfSample.Frequency=Frequency;
|
||||
RfSample.WaitForThisSample=Timing; //en 100 de nanosecond
|
||||
//printf("Freq =%f Timing=%d\n",RfSample.Frequency,RfSample.WaitForThisSample);
|
||||
if (write(FileFreqTiming,&RfSample,sizeof(samplerf_t)) != sizeof(samplerf_t)) {
|
||||
fprintf(stderr, "Unable to write sample\n");
|
||||
}
|
||||
|
||||
}
|
||||
|
||||
int main(int argc, char **argv) {
|
||||
float data [2*BUFFER_LEN] ;
|
||||
float data_filtered[2*BUFFER_LEN] ; // we generate complex I/Q samples
|
||||
SNDFILE *infile, *outfile ;
|
||||
SF_INFO sfinfo ;
|
||||
|
||||
int readcount, nb_samples ;
|
||||
char *infilename ;
|
||||
char *outfilename ;
|
||||
int k ;
|
||||
float x ;
|
||||
|
||||
if( argc < 2 ) {
|
||||
printf("Usage : %s in.wav [out.wav]\n", argv[0]);
|
||||
return(1);
|
||||
}
|
||||
infilename = argv[1];
|
||||
if( argc == 3 ) {
|
||||
outfilename = argv[2];
|
||||
} else {
|
||||
outfilename = (char *)malloc( 128 );
|
||||
sprintf( outfilename, "%s", "out.ft");
|
||||
}
|
||||
if (! (infile = sf_open (infilename, SFM_READ, &sfinfo)))
|
||||
{ /* Open failed so print an error message. */
|
||||
printf ("Not able to open input file %s.\n", infilename);
|
||||
/* Print the error message from libsndfile. */
|
||||
puts (sf_strerror (NULL));
|
||||
return 1;
|
||||
}
|
||||
|
||||
if( sfinfo.samplerate != 48000 ) {
|
||||
printf("Input rate must be 48K.\n");
|
||||
return 1;
|
||||
}
|
||||
|
||||
FileFreqTiming = open(outfilename, O_WRONLY|O_CREAT, 0644);
|
||||
|
||||
/** **/
|
||||
printf ("Reading file : %s\n", infilename );
|
||||
printf ("Sample Rate : %d\n", sfinfo.samplerate);
|
||||
printf ("Channels : %d\n", sfinfo.channels);
|
||||
printf ("----------------------------------------\n");
|
||||
printf ("Writing file : %s\n", outfilename );
|
||||
|
||||
/* While there are.frames in the input file, read them, process
|
||||
** them and write them to the output file.
|
||||
*/
|
||||
|
||||
while ((readcount = sf_readf_float(infile, data, BUFFER_LEN)))
|
||||
{
|
||||
nb_samples = readcount / sfinfo.channels;
|
||||
for( k=0 ; k < nb_samples ; k++ ) {
|
||||
x = data[k*sfinfo.channels];
|
||||
if( sfinfo.channels == 2 ) {
|
||||
// stereo file, avg left + right
|
||||
x += data[k*sfinfo.channels+1];
|
||||
x /= 2;
|
||||
}
|
||||
//printf("%f \n",x);
|
||||
float FactAmplitude=2.0; // To be analyzed more deeply !
|
||||
/*
|
||||
double A = 87.7f; // compression parameter
|
||||
double ampf=x/32767.0;
|
||||
ampf = (fabs(ampf) < 1.0f/A) ? A*fabs(ampf)/(1.0f+ln(A)) : (1.0f+ln(A*fabs(ampf)))/(1.0f+ln(A)); //compand
|
||||
x= (int)(round(ampf * 32767.0f));
|
||||
*/
|
||||
WriteTone(x*32767*FactAmplitude,1e9/48000.0);
|
||||
|
||||
}
|
||||
|
||||
}
|
||||
|
||||
/* Close input and output files. */
|
||||
sf_close (infile) ;
|
||||
close(FileFreqTiming);
|
||||
|
||||
return 0;
|
||||
}
|
||||
Loading…
Reference in a new issue