Алгоритм стеганографии с использованием псевдослучайной последовательности

% Инициализация
Q = 1; %амплитуда встраивания
[image, Fs] = imread('C:\Users\Лашын\Desktop\Visimag.jpg'); %считывание музыки, int16
image = image(:,:,3);
image = double(image);
a = size(image, 1); %размер столбца
b = size(image, 2); %размер строки
n = a*b; %размер контейнера
Lpsp = 100; %длина ПСП
% "Передатчик".
l = 1000; %количество отсчетов в одном участке контейнера для разбиения
image_column = reshape(image, n, 1); %контейнер в столбец
image_data_write = zeros(n, 1); %для записи файла с данными
for i = 1 : n/l %цикл по участкам контейнера
 image_data = image_column((i - 1)*l + 1 : i*l, :);%участок контейнера
 
 % Подготовка ПСП с информацией
 psp = unidrnd(2, Lpsp, 1)*2 - 3;
 data = randi([0 1], l/Lpsp, 1);
 data_psk = data*2 - 1;
 pspdata = data_psk * psp.';
 pspdata = reshape(pspdata, l, 1);
 
 % Стего
 data_to_embed = Q * pspdata; %данные для встраивания
 stego = image_data + data_to_embed; %стего
 
 image_data_write((i - 1)*l + 1 : i*l, :) = stego;
end
%figure(1)
%received_image = reshape(image_data_write,376,436);
%imshow(uint8(received_image));
for i = 1 : n/l %цикл по участкам контейнера
 received_data = image_data_write((i - 1)*l + 1 : i*l, :);%участок контейнера
 received_data = stego - image_data;
 received_data = received_data*psp.';
 received_data = sum(received_data,2);
end
received_image_data = reshape(received_data,10,100);
received_image_data = received_image_data(:,1);
Out(received_image_data > 0) = 0;
Out(received_image_data < 0) = 1;