Свертка

@banan21 · Регистрация: 01.05.2013

Студворк — интернет-сервис помощи студентам

Помогите чайнику, незнающему тонкостей(( Смысл заключается в реализации свертки текста и сигнала с построением графиков. Часть программы написал, а что нужно далее сделать для свертки? найти БПФ текста и умножить его на БПФ сигнала? Подскажите с кодом программы. Заранее спасибо.

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
close all;
 clear all;
 clc;
 
 %Обработка на стороне отправителя
 %Открытие text.txt 
 fid = fopen('C:\Users\banan\Documents\диплом\диплом\прога\тест\tex t.txt');
 stegofile = fread(fid);
 status = fclose('all');
 
 %Преобразование содержащегося текста в двоичный вид
 nfile = de2bi(stegofile,8);
 
 %Вес файла в байтах
 razmerfile=size(stegofile,1); 
 
 %Преобразование матрицы в строку и перевод в int16
 for i = 1:20
    for j = 1:8
       m = nfile (i,j);
       if (m == 1)
          v = 1 
       else 
          v = 0
       end
     x(j+8*(i-1)) = int16(v);
 end
 end
 n = size(x)
 n1 = n(2)
 
 % Открытие .wav файла и его чтение
 [y,Fs,nbits] = wavread('C:\Windows\Media\tada.wav'); % читаем вав
 sound(y,Fs);
 Y=fftshift(fft(y)); % делаем БПФ - длина зависит только от вектора у
 scale=linspace(-pi,pi,length(Y)); % шкала оси Х
 figure(1);
 plot(scale,abs(Y)); % амплитуда
 y1 = y(:,1)*2^15; %выбор левого канала 
 y2 = int16 (y1); %перевод в int16
X=fftshift(fft(nfile)); %БПФ текста
Z = conv(X, Y); %свертка

далее нужно получить опять аудио файл, только уже вместо вектора Y будет вектор Z.
Моё рассуждение правильно?

Зосима · 25.05.2013, 16:56

Дружочек, ты чуток перемудрил

Свертка во временной области будет просто:
z = conv(x, y);
Где x - отсчеты текста, y - отсчеты звука.
А если упороться с БПФ, то должно быть:
Z = X.*Y;
z = ifft(Z); (обратное фурье от произведения спектров)
Где X, Y - спектры текста и звука, без fftshift и abs!!!

Смекаешь?

@banan21 · 25.05.2013, 18:37 **[ТС]**

Понял. Значит выглядит это так?
Только вылезло
Error using .*
Matrix dimensions must agree.

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
close all;
clear all;
clc;
 
%Обработка на стороне отправителя
%Открытие text.txt 
fid = fopen('C:\Users\banan\Documents\диплом\диплом\прога\тест\text.txt');
stegofile = fread(fid);
status = fclose('all');
 
%Преобразование содержащегося текста в двоичный вид
nfile = de2bi(stegofile,8);
 
%Вес файла в байтах
razmerfile=size(stegofile,1); 
 
% Открытие .wav файла и его чтение
[y,Fs,bits] = wavread('C:\Windows\Media\tada.wav'); % читаем вав
sound(y,Fs);
X=fft(y); % делаем БПФ -  длина зависит только от вектора у
scale=linspace(-pi,pi,length(X)); % шкала оси Х
figure(1);
plot(scale,abs(X)); % амплитуда
 
 
Y=fft(nfile);
Z = X.*Y;
stegw = ifft(Z);
 
 
%запись в stegowav.wav (звуковой файл, отправляемый приемной стороне)
wavwrite(stegw,'C:\Users\banan\Documents\диплом\диплом\прога\тест\stegowav.wav'); %запись в stegowav.wav

<71296x2 complex double> размер X
<23x8 complex double> размер Y

Добавлено через 41 минуту
Подскажите в чем ошибка.

Зосима · 25.05.2013, 20:15

Как в чем? размеры массивов разные! (с уважением, Ваш Кэп)

Вот и мне теперь стало совершенно непонятно, как совместить две совершенно разные матрицы

@banan21 · 25.05.2013, 21:53 **[ТС]**

А если размерность будет
<71296x1 complex double> размер X
<23x1 complex double> размер Y
тогда возможно вычисление?

Зосима · 25.05.2013, 21:55

Да

тогда просто
z = conv(X, Y);
только это должны быть временные отсчеты, а не спектры!

@banan21 · 26.05.2013, 02:08 **[ТС]**

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
close all;
close all;
clear all;
clc;
 
%Обработка на стороне отправителя
%Открытие text.txt 
fid = fopen('C:\Users\banan\Documents\диплом\диплом\прога\тест\text.txt');
stegofile = fread(fid);
status = fclose('all');
 
 
% Открытие .wav файла и его чтение
[y,Fs,nbits] = wavread('C:\Windows\Media\tada.wav'); % читаем вав
sound(y,Fs);
X=fft(y); % делаем БПФ -  длина зависит только от вектора у
scale=linspace(-pi,pi,length(X)); % шкала оси Х
figure(1);
plot(scale,abs(X)); % амплитуда
[COLOR="Red"]y1 = y(:,1)*2^15; %выбор левого канала [/COLOR]
 
% БПФ левого канала
X=fft(y1);
 
% БПФ текста
Y=fft(stegofile);
 
%Свертка
z = conv(X,Y);
 
%Вычисление ОБПФ
stegw = ifft(z);
 
 
%Запись в stegowav.wav (звуковой файл, отправляемый приемной стороне)
wavwrite(stegw,'C:\Users\banan\Documents\диплом\диплом\прога\тест\stegowav.wav'); %запись в stegowav.wav

Как теперь сложить канал который я отделил <71296x1 complex double>
и получить похожую на исходный формат <71296x2 complex double>
чтобы вав файл получился?

Зосима · 27.05.2013, 14:09

как можно помочь человеку, который не читает твоих сообщений?!

Сообщение от Зосима

Свертка во временной области будет просто:
z = conv(x, y);
...
А если упороться с БПФ, то должно быть:
Z = X.*Y;
z = ifft(Z); (обратное фурье от произведения спектров)

Так правильно, но из-за разности длинн работать не будет:

Кликните здесь для просмотра всего текста

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
close all;
clear all;
clc;
 
%Обработка на стороне отправителя
%Открытие text.txt 
fid = fopen('C:\Users\banan\Documents\диплом\диплом\прога\тест\text.txt');
stegofile = fread(fid);
status = fclose('all');
% Открытие .wav файла и его чтение
[y,Fs,nbits] = wavread('C:\Windows\Media\tada.wav'); % читаем вав
sound(y,Fs);
X=fft(y); % делаем БПФ -  длина зависит только от вектора у
scale=linspace(-pi,pi,length(X)); % шкала оси Х
figure(1);
plot(scale,abs(X)); % амплитуда
y1 = y(:,1)*2^15; %выбор левого канала  
% БПФ левого канала
X=fft(y1);
 % БПФ текста
Y=fft(stegofile);
%Свертка
z = X.*Y;
%Вычисление ОБПФ
stegw = real(ifft(z));
 
%Запись в stegowav.wav (звуковой файл, отправляемый приемной стороне)
wavwrite(stegw,'C:\Users\banan\Documents\диплом\диплом\прога\тест\stegowav.wav'); %запись в stegowav.wav

А вот так - будет работать

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
close all;
clear all;
clc;
 
%Обработка на стороне отправителя
%Открытие text.txt 
fid = fopen('C:\Users\banan\Documents\диплом\диплом\прога\тест\text.txt');
stegofile = fread(fid);
status = fclose('all');
% Открытие .wav файла и его чтение
[y,Fs,nbits] = wavread('C:\Windows\Media\tada.wav'); % читаем вав
sound(y,Fs);
X=fft(y); % делаем БПФ -  длина зависит только от вектора у
scale=linspace(-pi,pi,length(X)); % шкала оси Х
figure(1);
plot(scale,abs(X)); % амплитуда
y1 = y(:,1)*2^15; %выбор левого канала  
 
stegw = conv(y1,stegofile);
 
%Запись в stegowav.wav (звуковой файл, отправляемый приемной стороне)
wavwrite(stegw,'C:\Users\banan\Documents\диплом\диплом\прога\тест\stegowav.wav'); %запись в stegowav.wav

Вся соль в строке 19.

@banan21 · 27.05.2013, 14:34 **[ТС]**

Спасибо за помощь!) Смысл теперь ясен, теперь пытаюсь выяснить как это реализовать на уровне спектров. Как разберусь отпишусь.

Добавлено через 15 минут
еще такой вопрос:

Matlab M
1
y1 = y(:,1)*2^15; %выбор левого канала

мы разделили сигнал на 2 канала, свернули, как теперь соеденить этот свернутый канал с каналом

Matlab M
1
y2 = y(:,2)*2^15; %выбор правого канала

размерность 1го (левого канала) получилась после свертки <71318x1 double>
а размер канала который был не тронут осталась <71296x1 double>
нам нужно получить матрицу <*****х2 double>

для сохранения wav файла

Зосима · 27.05.2013, 15:02

А ты часом не прогреваешь провода в своей аудиосистеме, для достижения более "теплого, лампового звука"?

Тебе одноканального моно файла не достаточно?

Ладно, тогда так:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
close all;
clear all;
clc;
 
%Обработка на стороне отправителя
%Открытие text.txt 
fid = fopen('C:\Users\banan\Documents\диплом\диплом\прога\тест\text.txt');
stegofile = fread(fid);
status = fclose('all');
% Открытие .wav файла и его чтение
[y,Fs,nbits] = wavread('C:\Windows\Media\tada.wav'); % читаем вав
sound(y,Fs);
X=fft(y); % делаем БПФ -  длина зависит только от вектора у
scale=linspace(-pi,pi,length(X)); % шкала оси Х
figure(1);
plot(scale,abs(X)); % амплитуда
y1 = y(:,1)*2^15; %выбор левого канала  
y2 = y(:,1)*2^15; %выбор правого канала  
 
s1 = conv(y1,stegofile);
s2 = conv(y1,stegofile);
stegw = [s1, s2]; % соединяем два столбца
 
%Запись в stegowav.wav (звуковой файл, отправляемый приемной стороне)
wavwrite(stegw,'C:\Users\banan\Documents\диплом\диплом\прога\тест\stegowav.wav'); %запись в stegowav.wav

@banan21 · 27.05.2013, 15:12 **[ТС]**

Спасибо!) даа))

Новые блоги и статьи Все статьи Все блоги /
Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .	Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .
Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .	Фото: Daniel Greenwood kumehtar 13.11.2025	Расскажи мне о Мире, бродяга kumehtar 12.11.2025 — Расскажи мне о Мире, бродяга, Ты же видел моря и метели. Как сменялись короны и стяги, Как эпохи стрелою летели. - Этот мир — это крылья и горы, Снег и пламя, любовь и тревоги, И бескрайние. . .

Зосима 5245 / 3573 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	25.05.2013, 16:56
	Дружочек, ты чуток перемудрил Свертка во временной области будет просто: z = conv(x, y); Где x - отсчеты текста, y - отсчеты звука. А если упороться с БПФ, то должно быть: Z = X.Y; z = ifft(Z); (обратное фурье от произведения спектров) Где X, Y - спектры текста и звука, без fftshift и abs!!!* Смекаешь? 1

Зосима 5245 / 3573 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	25.05.2013, 20:15
	Как в чем? размеры массивов разные! (с уважением, Ваш Кэп) Вот и мне теперь стало совершенно непонятно, как совместить две совершенно разные матрицы 0

@banan21 0 / 0 / 0 Регистрация: 01.05.2013 Сообщений: 12
	25.05.2013, 21:53 [ТС]
	А если размерность будет <71296x1 complex double> размер X <23x1 complex double> размер Y тогда возможно вычисление? 0

Зосима 5245 / 3573 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	25.05.2013, 21:55
	Да тогда просто z = conv(X, Y); только это должны быть временные отсчеты, а не спектры! 1

@banan21 0 / 0 / 0 Регистрация: 01.05.2013 Сообщений: 12
	27.05.2013, 15:12 [ТС]
	Спасибо!) даа)) 0