Нужно подправить программу по соответствующей теме описанной выше. Ниже текст программы.
close all
[s, fs] = audioread('D:\Desktop\Zil1\1.wav');
soundsc(s, fs);
figure(1);
plot(s);
grid on
n=1000;
t=1000/(2*fs);
i=0:t:length(s);
z=diff(s);
soundsc(z, fs);
figure(2);
plot((fs/2)*(1:length(z))/length(z),z);
title('ДИКМ')
xlabel('x');
ylabel('y');
grid on
p=cumsum(z);
soundsc(p, fs);
figure(3);
plot(p);
title('Восстановленный сигнал ')
xlabel('y');
ylabel('а');
grid on