FFT sur signal CW

[invite74228fda]

Bonjour,

Quelques mots de présentation avant de vous exposer ma question : je suis Fabrice, radioamateur (avec l'indicatif F4BJH), et mon club radioamateur (F8KGL) s'est lancé dans un beau projet de Radiologicielle, que nous avons appelé "le BAVAR" BAVAR = Beacon Analyser Vauréal Amitié Radio.

Il s'agit d'analyser la réception des balises RF sur les bandes radioamateurs, dans le but d'étudier les phénomènes de propagation radio (vaste sujet , pour améliorer nos perfs de concours en temps réel. Les concours sont un rassemblement des radioamateurs sur les ondes, où le but du jeu est de contacter le plus de monde possible, le plus loin.

Pour étudier les phénomènes de propagation, des radioamateurs ont installés des balises (émetteur code morse) un peu partout en europe (et dans le monde). Leur écoute permet de juger de la qualité de réception, à un instant donné, de station proche géographiquement de ces balises.

J'en viens à ma question : les balises émettent en code morse. Une succession de "1 et 0", correspondant, respectivement, à "présence" ou "absence" de porteuse.

Pour commencer, j'ai donc cherché à calculer une FFT sous Octave d'une porteuse modulée en CW. Mais le résultat donné par le graphe ne correspond pas du tout à la fréquence recherchée (en l'occurence 144,425MHz)

Avant de moduler ma porteuse RF, avec une succesion de 1 et 0, j'aurais bien aimé déjà avoir une belle porteuse à la fréquence voulue...ça me semble être le minimum. Le message est donc [1,1,1,1,1]

Le message ne contient que des "1" => présence de porteuse tout le temps => le signal entrant dans la FFT est donc un sinus pur à 144,425MHz, et je m'attends à avoir une belle raie à cette fréquence. La figure montre que non

Quelqu'un aurait-il une idée ? Où ai-je fais une erreur ?

En vous remerciant par avance pour votre aide,
73's
Fabrice, F4BJH
-----

[invite74228fda]

Et voici le code source

Merci d'avance
73's
Fabrice, F4BJH

[JPL]

Il est préférable de mettre le code dans la balise Code (# dans la barre d'outils) [Code]....[/Code]. Ainsi il sera lisible directement.

[azad]

Salut.....
Et pourquoi ne pas émettre en CW pure une séquence propre à chaque participant ? C'est à dire en tout ou rien ? Chaque participant aurait une séquence sur une dizaine de bits qu'il émettrait à intervalles réguliers, ce qui permettrait (connaissant le QTH de chacun) d'évaluer la propagation.
Ou même son indicatif, après tout ce qui éviterait d'avoir à affecter à chacun un code personnel.

[A voir en vidéo sur Futura]

[invite74228fda]

Beh, je pense que c'est déjà le cas en fait. D'après moi A1A = tout ou rien = présence ou absence de porteuse. Après, le "tout" peut être plus ou moins long, pour faire la distinction d'un trait ou d'un point...

Mais bon, de toute façon, je n'ai pas tellement le choix. Les balises sont déjà existantes, sur tout un plan de bande :

http://thf.ref-union.org/144.htm

Et la modulation des balises est déjà dans le mode A1A.

73's
Fabrice, F4BJH

[invite74228fda]

%#####################################################
%génération porteuse modulée CW A1A
%#####################################################
	
	rf_out = generate_cw_A1A(f_F5ZAM, msg_cw_k, pow_rf1_dbm);
	%rf_out = generate_cw_A1A(freq_rf1,msg_cw_k, pow_rf1_dbm) + generate_cw_A1A(freq_rf2,msg_cw_a, pow_rf2_dbm) + generate_cw_A1A(freq_rf2,msg_cw_4, pow_rf3_dbm)

%calcul de la DSP du signal de sortie
	rf_out_hat = fft(rf_out)/size(rf_out,2);
	dsp_rf_out = rf_out_hat.*conj(rf_out_hat);

function [rf_out] = generate_cw_A1A( freq_rf, msg_cw, rf_pow_dbm)

%#####################################################
%constantes
%#####################################################
	%échantillonnage
	r_sample = 10;
	f_sample = r_sample*freq_rf;
	T_sample = 1/f_sample;
	%constante CW : 1 point = 10 périodes de rf
	%dot_last = 10/freq_rf;
	dot_last = 1;
	%puissance de sortie (hypothèse : le modulateur sort 1W)
	txgain = 10^(rf_pow_dbm/10);

%#####################################################
%génération OL
%#####################################################
	t_dot = T_sample : T_sample : dot_last;
	dot_sin = sin(2*pi*freq_rf*t_dot);

%#####################################################
%modulation
%#####################################################

	msg_size = size(msg_cw,2);
	dot_sin_size = size(dot_sin,2);
	
	for (i=1:1:msg_size)
		data_out((i-1)*dot_sin_size +1 : 1 : i*dot_sin_size) = msg_cw(i).*dot_sin;
	endfor

	rf_out = txgain*data_out;

[invite74228fda]

Ahhh....je viens de comprendre : faire émettre de la CW à chaque participant du contest...

Oui, mais ce n'est pas dans le règlement...et vu qu'en plus, il y a pas loin de 100 000 participants, si vous voulez...c'est pas si simple de demander à tout le monde de changer de matériel, juste pour mon "petit plaisir personnel

mais sinon, oui...ce serait une "pure" idée

73's
Fabrice, F4BJH

[azad]

Wouha !!!!
Cent mille participants ! Pour peu que les Smètres soient à la hauteur (et cohérents entre-eux), vous allez avoir une représentation dynamique de la propagation qui va faire des envieux dans les hautes sphères. En tout cas : bravo, beau projet.