Cкрытая передача информации в аудиосигналах с использованием вейвлет-преобразования

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Список литературы

 

1. Афонский А. А., Дьяконов В. П. Цифровые анализаторы спектра, сигналов и логики. Под ред. проф. В. П. Дьяконова. — М: СОЛОН-Пресс, 2009.
2. Дьяконов В. П. MATLAB 6.5 SP1/7.0/7.0 SP1 + Simulink 5/6/ Обработка сигналов и проектирование фильтров. — М: СОЛОН-Пресс, 2005.
3. Дьяконов В. П. Вейвлеты. От теории к практике. М.: СОЛОН-Пресс, 2002, 448 с..
4. Matsui K., Tanaka K., and Nakamura Y. Digital signature on a facsimile document by recursive MH coding // Symposium On Cryptography and Information Security, 1989.
5. Osborne C., van Schyndel R., Tirkel A. A Digital Watermark // IEEE Intern. Conf. on Image Processing, 1994. P. 86-90.
6. Arnold M., Kanka S. MP3 robust audio watermarking // International Watermarking Workshop. 1999.
7. Bassia P., Pitas I., Robust audio watermarking in the time domain // Department of Informatics, University of Tressaloniki.
8. Boney L., Tewfic A.H., Hamdy A.K., Digital watermarks for audio signals, Department of Electrical engineering, University of Minnesota.

 

 

 

 

 

 

 

 

 

 

 

 

 

Приложение

 

Программа, написанная на MatLab, реализующая встраивание ЦВЗ в аудиосигнал.

 

clear all

k=0.1;

N=45000;

A=unifrnd(0,1,1,N);

[C,D]=schetcd(A);

[y, Fs, nbits] = wavread('E:\track no07');

y(1000001:1000100,1);

for i=1:100;

t=(y((1000001:1000000+N)+ (i-1)*N,1));

zz(i)=nat(t,k,C,D);

f=vstrznak(t,k,C,D);

z(i)=nat(f,k,C,D);

fprintf('i=%d\n',i)

end

 

 

function X=vstrznak(x,k,C,D)

% x=dct(x);

% [x,l] = wavedec(x,5,'db5');

X=x;

for i=1:length(C)

if D(i)==1

X(i)=x(i)-k*x(i);

else

if C(i)==1

X(i)=x(i)+k*x(i);

else

X(i)=x(i);

end

end

end

% X=idct(X);

% X = waverec(X,l,'db5');

 

 

function z=nat(x,k,C,D)

%x=dct(x);

%[x,l] = wavedec(x,5,'db5');

x=x(1:length(C));

sumc=sum(abs(x(C)));

sumd=sum(abs(x(D)));

z=(sumc-sumd)/(std(x)*sum(C));

return;

 

 

function [C,D]=schetcd(A)

E=sort(A);

C=A>E(length(A)-20000);

D=A<=E(20000);