t/WirelessNetworkingTechnologies
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

task2 implemented

Browse Source
This commit is contained in:
ti_mo
2026-06-08 14:35:03 +02:00
parent eff95ac6e2
commit 1383556946
6 changed files with 1689 additions and 89 deletions
Download Patch File Download Diff File Expand all files Collapse all files
lab_4/lab_4.m
+13 -89
View File
@@ -2,7 +2,8 @@ carrier_frequency = 2.4e9;
wifi_be_subcarrier_spacing = 78.125e3;
leo_altitude = 600e3;
bandwidth = 20;
mcs_values = [2];
mcs_values = [0 2 4 8 10 12 13];
mcs_t1 = 2;
apep = 1000;
tx_antenna_count = 1;
rx_antenna_count = 1;
@@ -11,7 +12,6 @@ snrs = (0:1:40);
packets_per_snr = 1024;
maxNumErrors = 50;
% elevation_angle = 50;
elevationAngles = [ ...
linspace(10, 90, ceil(packets_per_snr / 2)), ...
linspace(90, 10, floor(packets_per_snr / 2)) ...
@@ -26,94 +26,18 @@ seed = 666;
ntn_tdl_delay_profiles = ["A" "B" "C" "D"];
[cfgEHT, simParameters, channel, chInfo, maxChDelay, satelliteDopplerShift, chanBW] = ...
create_baseline_configuration(ntn_tdl_delay_profiles(3), bandwidth, false, "None", ...
carrier_frequency, 0, ...
0, leo_altitude, random_stream, ...
seed, tx_antenna_count, ...
rx_antenna_count, apep);
use_t1_as_t4 = false;
task1;
task2;
meanSlantRange = mean(slantRanges);
meanSatelliteDopplerShift = mean(dopplerShifts);
meanNormalizedCFO = meanSatelliteDopplerShift / wifi_be_subcarrier_spacing;
disp("Mean Elevation Angle: " + meanElevationAngle + "°")
disp("Mean Slant Range: " + meanSlantRange + "m")
disp("Mean Satellite Doppler Shift: " + meanSatelliteDopplerShift + "Hz")
disp("Normalized CFO: " + meanNormalizedCFO)
errorCounts = zeros(numel(mcs_values), numel(snrs));
packetCounts = zeros(numel(mcs_values), numel(snrs));
errorRates = zeros(numel(mcs_values), numel(snrs));
cfrMagnitudesBySnr = cell(numel(mcs_values), numel(snrs));
for imcs = 1 : numel(mcs_values)
cfgEHT.User{1}.MCS = mcs_values(imcs);
errorCountsForMcs = zeros(1, numel(snrs));
packetCountsForMcs = zeros(1, numel(snrs));
errorRatesForMcs = zeros(1, numel(snrs));
cfrMagnitudesForMcs = cell(1, numel(snrs));
parfor snr_idx = 1 : length(snrs)
local_channel = clone(channel);
[errorCount, packetCount, cfrMagnitudes] = simulateTransmission(cfgEHT, local_channel, maxChDelay, snrs(snr_idx), packets_per_snr, maxNumErrors, simParameters, chanBW, snr_idx, elevationAngles);
errorCountsForMcs(snr_idx) = errorCount;
packetCountsForMcs(snr_idx) = packetCount;
errorRatesForMcs(snr_idx) = errorCount / packetCount;
cfrMagnitudesForMcs{snr_idx} = cfrMagnitudes;
end
errorCounts(imcs, :) = errorCountsForMcs;
packetCounts(imcs, :) = packetCountsForMcs;
errorRates(imcs, :) = errorRatesForMcs;
cfrMagnitudesBySnr(imcs, :) = cfrMagnitudesForMcs;
end
markers = 'ox*sd^v><ph+ox*sd^v><ph+';
figure;
for imcs = 1 : numel(mcs_values)
semilogy(snrs, errorRates(imcs, :).',['-' markers(imcs)]);
hold on;
end
grid on;
xlabel('SNR (dB)');
ylabel('PER');
% task 1
dataStr = arrayfun(@(x)sprintf('MCS %d', x), mcs_values, 'UniformOutput', false);
% task 4
% dataStr = arrayfun(@(x)'MCS Sweep Values: [0 2 4 8 10 12 13]', 'UniformOutput', false);
legend(dataStr, 'Location', 'NorthEastOutside');
title(['PER (EHT MU), ' num2str(cfgEHT.ChannelBandwidth) ', NTN-TDL-C, ' num2str(tx_antenna_count) '-by-' num2str(rx_antenna_count)]);
cfrMagnitudes = [];
for snr_idx = 1 : numel(snrs)
cfrMagnitudes = [cfrMagnitudes cfrMagnitudesBySnr{1, snr_idx}];
end
%{ for testing purposes
if ~isempty(cfrMagnitudes)
figure;
packetPlots = plot(cfrMagnitudes, 'Color', [0.75 0.75 0.75]);
hold on;
meanCfrPlot = plot(mean(cfrMagnitudes, 2), 'b', 'LineWidth', 2);
grid on;
xlabel("Subcarrier index");
ylabel("|H[k]|");
legend([packetPlots(1) meanCfrPlot], ["Packets" "Mean CFR"], 'Location', 'best');
title("CFR magnitude for collected packets");
end
%}
% extract a representative packet for cfr plot
representativeSnr = 25;
representativeSnrIdx = find(snrs == representativeSnr);
representativeCfrs = cfrMagnitudesBySnr{1, representativeSnrIdx};
representativeCfr = representativeCfrs(:, 1);
apep1000Idx = find(apep_values == 1000);
figure;
plot(representativeCfr);
semilogy(snrs, errorRates(1,:), "-o");
hold on;
semilogy(snrs, results(apep1000Idx,:), "-x");
grid on;
xlabel("Subcarrier index");
ylabel("|H[k]|");
title("CFR magnitude for representative packet at 25 dB SNR");
xlabel("SNR (dB)");
ylabel("PER");
legend("Task 1 APEP 1000", "Task 2 APEP 1000");
title("Task 1 vs Task 2 consistency check");