mathlab - сглаживание сигнала

[tranvay]

Привет всем!
пожалуйста помогите решить задачу!!
необходимо сгладить сигнал, который имеет случайные импульсные помехи, короткие замыкания и обрывы.
Сигнал, который все это создает вот:

F=1000; %массив отсчетов датчика
M=14.7; %математическое ожидание
Dsh=3; %дисперсия шума
x=randn(F,1)+M; %randn(m,n) - формирует массив размера m х n, элементами которого являются случайные величины,
%распределенные по нормальному закону с математическим ожиданием 0 и
%среднеквадратическим отклонением 1.


figure (2);
plot(x);

% Имитация импульсной помехи

for m = 1:2
Pom=M*5; %амплитуда помехи
t =150; %индекс элемента массива, где будет помеха

if (round(rand(1)*10)>5) %положение помехи относительно математического ожидания
x(t) = Pom
else x(t) = 2*M-Pom;
end
end


% Имитация короткого замыкания


KZ=M*3; %амплитуда короткого замыкания
tz =500; %индекс элемента массива, где будет короткое замыкание

for m= 1:50
x(tz) = KZ
tz=tz+1;
end


% Имитация обрыва

Obr=0; %амплитуда обрыва
tz =750; %индекс элемента массива, где будет обрыв
for m= 1:50
x(tz) = Obr
tz=tz+1;
end

figure (3);
plot(x);
D1 = cov(x);


plot(x)
grid on
text(1, 75, strcat('Dispersiya:', num2str(D1)))

Помогите пожалуйста со сглаживанием!