Tôi đang cố gắng hiệu chỉnh loa siêu âm với mục đích phát ra các tín hiệu có thể dự đoán được, nhưng than ôi tôi vẫn gặp rắc rối, có lẽ là do tôi không có DSP-fu.
Một chút nền tảng
Tôi muốn có thể phát lại âm thanh gần nhất có thể với bản ghi đã hiệu chỉnh mà tôi có. Theo như tôi hiểu lý thuyết, tôi cần tìm chức năng chuyển loa và giải mã các tín hiệu tôi muốn phát ra với nó. Một cái gì đó như thế này (trong miền tần số):
X -> H -> XH
X
Tín hiệu phát ra ở đâu là H
chức năng truyền loa và XH
là X
thời gian H
. Một bộ phận ( ./
) bây giờ sẽ cho tôi H
.
Bây giờ, để phát ra tín hiệu hiệu chỉnh, nên chia cho H
:
X/H -> H -> X
Những gì đã được thực hiện
- Loa được đặt và micrô được hiệu chỉnh cách nhau 1 m trên giá ba chân.
- Ghi lại hơn 30 lần quét tuyến tính 150KHz-20KHz, dài 20ms và được ghi @ 500 KS / s.
- Các tín hiệu được căn chỉnh và trung bình với tập lệnh Matlab / Octave bên dưới, dưới tập lệnh có thể nhìn thấy tín hiệu kết quả.
files = dir('Mandag*');
rng = [1.5e6, 1.52e6];
[X, fs] = wavread(files(1).name, rng);
X = X(:,1);
for i=2:length(files)
[Y, fs] = wavread(files(i).name, rng);
sig = Y(:,1);
[x, off] = max(xcorr(X', sig'));
off = length(X) - off;
if(off < 0)
sig = [zeros(1, -off), sig(1:end+off)'];
elseif (off > 0)
sig = [sig(off:end)', zeros(1, off-1)];
end
X = X + sig';
end
X = X/length(files);
Fourier biến đổi
X
vàXH
thực hiện các tính toán được đề cập ở trên, kết quả có vẻ hợp lý. Dưới đây là một âm mưu chuẩn hóa củaH
(màu tím) vàX/H
(màu xanh lá cây).
Cốt truyện đã được cắt ngắn đến các tần số có liên quan.
Xin vui lòng cho tôi biết nếu tôi đi sai về điều này.
Câu hỏi của tôi
Sau khi tính toán X/H
tôi cần chuyển đổi nó trở lại miền thời gian, tôi cho rằng việc này sẽ đơn giản ifft(X./H)
và wavwrite
, nhưng tất cả các nỗ lực của tôi cho đến nay đều không nhận được bất kỳ câu trả lời hợp lý nào. Một vector tần số Hf
, H
và X
có thể được tìm thấy ở đây trong định dạng mat7-nhị phân.
Có lẽ tôi chỉ mệt mỏi và có một giải pháp đơn giản ở đây, nhưng hiện tại tôi không thể thấy nó. Bất kỳ trợ giúp / lời khuyên được nhiều đánh giá cao.