DenoisingSignal.jl

# Signaali ja kohinan poisto, denoising
# Juhani Kaukoranta 26.1.2024
using KissSmoothing, Statistics,Plots
X = LinRange(0, pi, 1000);
Y = sin.(X.^2) .* cos.(X);
p1 = plot(X,Y,legend=:topleft,lw=2,label="originaali kohinaton signaali",ylimits=(-1.2,1.2));
TN = std(Y).*randn(length(X))./5; # lisätty kohinaa keskihajonta/5 verran
noisy_S = Y .+ TN;
# denoise poistaa kohinan N, tuottaa kohinattoman S
S, N = denoise(noisy_S); #
p2 = plot(X,noisy_S,legend=:topleft,color=:gray,lw=0.8,
label="kohinainen signaali",ylimits=(-1.2,1.2))
p3 = plot(X,S,legend=:topleft,color=:red,label = "kohinanpoiston jälkeen",ylimits=(-1.2,1.2));
plot(p1,p2,p3)