%
% f.m -- Reed-Solomon-Visualisierung mit FFT
%
% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule

N = 64;
b = 32;
l = N + b;

signal = zeros(l,1);
signal(1:N,1) = round(10 * rand(N,1));
signal

plot(abs(signal));
xlim([1, l]);
title("Signal");
pause()

codiert = fft(signal)

plot(abs(codiert));
xlim([1, l]);
title("Codiert");
pause()

fehler = zeros(l,1);
fehler(21,1) = 2;
fehler(75,1) = 1;
fehler(7,1) = 2;

plot(fehler);
xlim([1, l]);
title("Fehler");
pause()

empfangen = codiert + fehler;

plot(abs(empfangen));
xlim([1, l]);
title("Empfangen");
pause()

decodiert = ifft(empfangen)
plot(abs(decodiert));
xlim([1, l]);
title("Decodiert");
pause()

syndrom = decodiert;
syndrom(1:N,1) = zeros(N,1)
plot(abs(syndrom));
xlim([1, l]);
title("Syndrom");

pause()

locator = abs(fft(syndrom))

plot(locator);
xlim([1, l]);
title("Locator");
pause()


writematrix([transpose(counter), abs(signal)], 'signal.txt')
writematrix([transpose(counter), abs(codiert)], 'codiert.txt')
writematrix([transpose(counter), fehler], 'fehler.txt')
writematrix([transpose(counter), abs(empfangen)], 'empfangen.txt')
writematrix([transpose(counter), abs(decodiert)], 'decodiert.txt')
writematrix([transpose(counter), abs(syndrom)], 'syndrom.txt')
writematrix([transpose(counter), locator], 'locator.txt')