silas/imdd_silas
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Added:

Browse Source 
- Class folder 'Timing Recovery' with different methods
- Minimal example for the timing recovery in the FSO folder
- New evaluation scripts for FSO data
This commit is contained in:
magf
2026-02-02 10:15:07 +01:00
parent 3886524b0e
commit a8daf12b20
3 changed files with 10 additions and 51 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
Classes/00_signals/Signal.m
View File

@@ -482,6 +482,11 @@ classdef Signal
yticks(round(linspace(y_min, y_max, min(10, max(4, ceil(y_range/10))))));
end
grid on;
% Add legend if not already present
if isempty(get(gca, 'Legend'))
legend;
end
end
@@ -924,7 +929,7 @@ classdef Signal
options.displayname = "";
end
mode = 1;
mode = 2;
histpoints = 2048; %% verticale resolution
histpoints = floor(histpoints/2)*2+1; %% to have the eye digram centered around one point make the vertical resolution uneven
@@ -965,7 +970,9 @@ classdef Signal
hold on
plot(eye_mat(:,n),'LineStyle',':','LineWidth',0.1,'Color',col(2,:));
end
ylabel('Amplitude of Signal');
xlabel('Samples','Interpreter','latex')
ylabel('Amplitude of Signal','Interpreter','latex');
xlim([0 histpoints_horizontal])
elseif mode == 1
% generate eye diagram using histogram

1
Classes/01_transmit/Pulseformer.m
View File

@@ -99,6 +99,7 @@ function data_out = process_(obj, data_in_signal)
% For your setup (powers of 2), it will likely be integer.
racos_len = obj.pulselength; % span in symbols
h = rcosdesign(obj.alpha, racos_len, sps_filter, filtertype);
% h = gaussdesign(obj.alpha, racos_len, sps_filter);
% Matched Filter Flip (Complex Conjugate Time Reversal)
if obj.matched

49
Classes/04_DSP/Timing_Recovery.m
View File

@@ -1,49 +0,0 @@
classdef Timing_Recovery < handle
properties(Access=public)
modulation
timing_error_detector
sps
damping_factor
normalized_loop_bandwidth
detector_gain
end
methods(Access=public)
function obj = Timing_Recovery(options)
arguments(Input)
options.modulation = 'PAM/PSK/QAM';
options.timing_error_detector = 'Gardner (non-data-aided)';
options.sps = 2;
options.damping_factor = 1.0;
options.normalized_loop_bandwidth = 0.01;
options.detector_gain = 2.7;
end
fn = fieldnames(options);
for n = 1:numel(fn)
obj.(fn{n}) = options.(fn{n});
end
%obj-Initialization here%
end
function [data_out, timing_error] = process(obj, data_in)
timing_synchronization = comm.SymbolSynchronizer( ...
"Modulation", obj.modulation,...
"TimingErrorDetector", obj.timing_error_detector, ...
"SamplesPerSymbol", obj.sps, ...
"DampingFactor", obj.damping_factor, ...
"NormalizedLoopBandwidth", obj.normalized_loop_bandwidth, ...
"DetectorGain", obj.detector_gain);
data_out = data_in;
[data_out.signal, timing_error] = timing_synchronization(data_in.signal);
end
end
end