আমি খুব শোরগোলের অডিও ফাইল (.wav) এর মধ্যে লুকানো কোনও বার্তাটি বোঝার চেষ্টা করছি (আমার মনে হয় এটি অতিরিক্ত লো ড্রোন সহ সাদা শব্দ)। বার্তাটি একটি ছয় সংখ্যার নম্বর । গোলমাল সম্পর্কে আমার আর কোনও বিবরণ নেই।
আমি বেশিরভাগ উচ্চ ফ্রিকোয়েন্সি মুছে ফেলার ফলে আমার সংখ্যা শোনার সুযোগ হবে এই আশায় একটি লো-পাস ফিল্টার ব্যবহার করার চেষ্টা করেছি তবে, ভয়েসটি ভালভাবে শুনতে শুনতে আমি কম ড্রোন থেকেও মুক্তি পেতে পারিনি বলে মনে হয়। আমার প্রয়াসটি নিম্নরূপ ছিল (নিযুক্ত নিয়োগটি freq_space_low_pass_filterশেষ করা হয়েছে):
[data, SampleRate, NbitsPerSample]=wavread('noisy_msg6.wav');
y=data(:,1); % we will work only with one channel in this demo
N=length(y); %number of sample points
t=( (1:N)*1/SampleRate ).'; % time spacing is 1/SampleRate and we want column vector
Y=fft(y);
spectrum_freq=fourier_frequencies(SampleRate, N);
Freq3db=100;
[spectrum_filtered,g_vs_freq]=freq_space_low_pass_filter(Y, SampleRate, Freq3db);
y_filtered=ifft(spectrum_filtered);
y_filtered=real(y_filtered);
wavwrite(y_filtered/(0.1+max(y_filtered)), SampleRate, NbitsPerSample,
'noisy_msg6_filtered.wav');
%%%%%%%%down sampling%%%%%%%%
indexes=(abs(spectrum_freq) < 10*Freq3db);
spectrum_freq_down_sampled = spectrum_freq(indexes);
spectrum_down_sampled = spectrum_filtered(indexes);
N_down_sampled = length(spectrum_down_sampled);
spectrum_down_sampled=spectrum_down_sampled*N_down_sampled/N;
SampleRate_down_sampled=SampleRate*N_down_sampled/N;
y_down_sampled=real(ifft(spectrum_down_sampled));
t_down_sampled = ( (1:N_down_sampled)*1/SampleRate_down_sampled ).';
sound(y_down_sampled, SampleRate_down_sampled)
function [spectrum_filtered,g]=freq_space_low_pass_filter(spectrum, SampleRate, Freq3db)
%% applies low pass filter in the frequency domain
% spectrum - result of fft on time series data (column vector is expected)
% SampleRate - measured in Hz, 1/dt where dt spacing of the points in time domain
% Freq3db - desired 3db roll off point in Hz
N=length(spectrum);
function G=filter_gain(freq, Freq3db)
G=1./(1+1i*freq/Freq3db); % this corresponds to low pass RC filter
end
spectrum_freq=fourier_frequencies(SampleRate, N);
% calculate filter gain for each spectrum frequency
g=filter_gain(spectrum_freq, Freq3db);
spectrum_filtered=spectrum.*g;
end
সিগন্যাল বর্ণালী প্লট:
আপনি আমাদের সাথে কাজ করছেন এমন একটি নমুনা সরবরাহ করতে পারেন? এবং সম্ভবত আপনার ফলাফল?
—
পেনেলোপ
আমি প্লটের চিত্রগুলি আপলোড করতে পারি না (পর্যাপ্ত পয়েন্ট নয়) এবং আপনাকে ফাইলে লিঙ্ক দেওয়ার কোনও উপায় নেই
—
user1825494
আপনি যদি মন্তব্যগুলিতে লিঙ্কগুলি রাখেন তবে কেউ আনন্দের সাথে এটিকে সম্পাদনা করবে And এবং যদি আপনি কোথাও ফাইলগুলি আপলোড করার এবং লিঙ্ক সরবরাহ করার কোনও উপায় খুঁজে পান তবে সেটিও ফিরে পান।
—
পেনেলোপ
যদি গ্রাফিকটি লগ-লগের পরিকল্পনা করা হত তবে গ্রাফটি আরও অনেক তথ্যপূর্ণ হবে। যেমনটি হয়, দেখে মনে হচ্ছে আপনার প্রায় 1.5 কেএজেডজ থেকে 22.5 কেজি হার্জ পর্যন্ত প্রশস্ত ব্র্যান্ডব্যান্ড রয়েছে, প্রায় 2 কেএইচআরজে সরু ব্যান্ড টোন আছে। আমি যা অনুমান করি তা হ'ল "ভয়েস" সংকেত, যা একটি উল্লেখযোগ্য ডিসি অফসেট হিসাবে উপস্থিত বলে মনে হয়, মোটামুটি শান্ত "গর্ত" এ বসে। উপযুক্ত ব্যান্ডপাস ফিল্টার দিয়ে বিচ্ছিন্ন করা সত্যিই এতটা কঠিন হওয়া উচিত নয়।
—
ডেভ টুইট করেছেন