এই প্রশ্নটি আমাকে দীর্ঘকাল ধরেও বিভ্রান্ত করেছে। @ হটপাউ 2 এর ব্যাখ্যাটি ভাল। আপনি মাতলাব ব্যবহার করে সাধারণ পরীক্ষায় আগ্রহী হতে পারেন।
https://poweidsplearningpath.blogspot.com/2019/04/dftidft.html
আপডেট তথ্য।
এই সত্যটি যাচাই করা সহজ, আমাদের কেবল সতর্কতার সাথে একটি আদর্শ (?) ব্যান্ড পাস ফিল্টারের অনুপ্রেরণামূলক প্রতিক্রিয়া বর্ণালীটি পর্যবেক্ষণ করতে হবে যা কেবল এফএফটি বিনগুলি জিরো করে ফেলেছে। আমার কেন "সাবধানে" বিজ্ঞাপনটি যুক্ত করা দরকার? আমরা যদি অনুপ্রবেশের প্রতিক্রিয়া পর্যবেক্ষণ করতে কেবল এফএফটি-র একই আকার ব্যবহার করি, তবে চিত্র 1 এ দেখানো হয়েছে বলে আমরা প্রতারণা করব । তবুও, যদি আমরা ফিল্টারটির আউটপুট পর্যবেক্ষণ করার সময় ডিএফটি এর ক্রম যুক্ত করি, অর্থাত, শূন্য প্যাডিং প্রতিক্রিয়া দেখায়, আমরা চিত্র 2 তে বর্ণিত চিত্র হিসাবে গিগস ঘটনাটি, ফ্রিকোয়েন্সি ডোমেনে রিপলগুলি খুঁজে পেতে পারি ।
বাস্তবে ফলাফলগুলি উইন্ডোংয়ের প্রভাব থেকে আসে। আপনি যদি সমস্যাটি সম্পূর্ণরূপে বুঝতে চান তবে অনুগ্রহ করে ডিএসপি (1) এর বাইবেলের অধ্যায় 7.6 এবং অধ্যায় 10.1-10.2 দেখুন। সংক্ষেপে, তিনটি মূল পয়েন্ট এখানে উল্লেখ করা হয়।
- ডিএফটি (এফএফটি) এর উইন্ডোর আকার এবং ক্রম সম্পূর্ণ স্বাধীন। তাদের একসাথে মিশ্রিত করবেন না।
- উইন্ডোর বৈশিষ্ট্য (টাইপ / আকার) ডিটিএফটি আকারে আধিপত্য করে। (উদাঃ ফ্রিকোয়েন্সি প্রতিক্রিয়াতে বৃহত্তর প্রধান লব সীসা বৃহত্তর ক্ষণস্থায়ী ব্যান্ড।)
- ডিএফটি হ'ল ফ্রিকোয়েন্সি ডোমেনে ডিটিএফটি-র নমুনা। তদুপরি, ডিএফটি-র ক্রমটি তত বেশি, ডিএফটি-র স্পেনট্রামটি হ্রাসযোগ্য।
সুতরাং, চিত্র 2-তে ডেনসার স্পেকট্রামের সাহায্যে , আমরা আদর্শ (জাল) ব্যান্ড পাস ফিল্টারটির মুখোশটি দিয়ে দেখতে পারি।
প্রতারণামূলকভাবে ফ্রেইক। প্রতিক্রিয়া।
ফ্রিক মধ্যে গীবস ঘটনা। প্রতিক্রিয়া।
(1) অ্যালান ভি। ওপেনহাইম এবং রোনাল্ড ডব্লু। শেফার। 2009. বিচ্ছিন্ন-সময় সংকেত প্রক্রিয়াকরণ (তৃতীয় সংস্করণ)। প্রেন্টিস হল প্রেস, আপার স্যাডল রিভার, এনজে, মার্কিন যুক্তরাষ্ট্র।
fps = 15;
LPF = 1;
HPF = 2;
n = -511:512;
n0 = 0;
imp = (n==n0);
NyquistF = 1/2*fps;
%% Ideal BPF
tmp_N = 512;
tmp_n = 0:1:tmp_N-1;
freq = ( n .* fps) ./ tmp_N;
F = fft(imp, tmp_N);
F_bpf = IdealBandpassFilter(F, fps, LPF, HPF);
imp_rep =[real(ifft(F_bpf))'];
% Zero padding.
imp_rep2 =[zeros(1,2048) real(ifft(F_bpf))' zeros(1,2048)];
N = 2^nextpow2(length(imp_rep));
F = fft(imp_rep,N);
freq_step = fps/N;
freq = -fps/2:freq_step:fps/2-freq_step;
freq = freq(N/2+1:end)';
figure;
plot(freq,abs(F(1:N/2)));
xlabel('freq(Hz)');
ylabel('mag');
title('Mis leading Freq Response');
N = 2^nextpow2(length(imp_rep2));
F = fft(imp_rep2,N);
freq_step = fps/N;
freq = -fps/2:freq_step:fps/2-freq_step;
freq = freq(N/2+1:end)';
figure;
plot(freq,abs(F(1:N/2)));
xlabel('freq(Hz)');
ylabel('mag');
title('Zero Padding (DFT) with more points');
%% Function
function filered_signal = IdealBandpassFilter(input_signal, fs, w1, w2)
N = length(input_signal);
n = 0:1:N-1;
freq = ( n .* fs) ./ N;
filered_signal = zeros(N, 1);
for i = 1:N
if freq(i) > w1 & freq(i) < w2
filered_signal(i) = input_signal(i);
end
end
end