Sinyal işlemede çok fazla kullanılan uygulamalardan biri örnekleme frekansının artırılması işlemidir.  Bu işlemi gerçekleştirilmesi için farklı yöntemler geliştirilmiştir. Bu yazımızda FFT kullanarak yukarı örnekleme işlemine ait uygulamayı MATLAB kodları ile gerçekleştireceğiz. Aşağıda 20 Hz'de örnkelenmiş 1 Hz'lik sünüs işareti gösterilmiştir.



Yukarı örnekleme işlemlerinde öncelikle sinyalin FFT alma işlemleri gerçekleştirilir. Daha sonra FFT simetrik olduğundan dolayı FFT sinyalinin ilk kısmı ile son kısmının arasına yukarı örnekleme yapılacak adım sayısının bir eksiği kadar FFT' boyunun katı kadar sıfırlar eklenir. Örneğin 3 ile yukarı örneklenecekse 32'lik FFT için 64 adet sıfır eklenmesi gerekmektedir. Daha sonra bu sinyalin ters FFT'si alınır ve gerçel kısımları yeni sinyalimiz olarak alınır. Aşağıda yukarıda verilen giriş sinyalinin 3 ile yukarı örneklenmesini sağlayan MATLAB kodu ve çıkışta elde edilen sinyal verilmiştir. 20 Hz ile örneklenen siynaln yeni örnekleme frekansı 60 Hz olmuştur.

Kod: Matlab M
  1.  
  2.  
  3. f_s = 20;
  4. f_I = 1;
  5.  
  6. n_t = 1 / f_s : 1 / f_s : 1;
  7.  
  8. s_Input = sin(2 * pi * f_I * n_t);
  9. figure, plot(n_t, s_Input);
  10. xlabel('Zaman(s)')
  11. title('20 Hzde örneklenmiş sinüs')
  12.  
  13. l_Input = length(s_Input);
  14. l_FFT = 2^(ceil(log2(l_Input)));
  15.  
  16. fft_I = fft(s_Input, l_FFT);
  17. up_L = 3;
  18.  
  19. Zero_Buff = zeros(1, (up_L - 1) * l_FFT);
  20. fft_I_new = up_L * [fft_I(1 : l_FFT / 2)...
  21.    Zero_Buff fft_I(l_FFT / 2 + 1 : l_FFT)];
  22. i_Input = real(ifft(fft_I_new, up_L * l_FFT));
  23. u_Input = i_Input(1 : up_L * l_Input);
  24. n_t = 1 / (up_L * f_s) : 1 / (up_L * f_s) : 1;
  25. figure, plot(n_t, u_Input)
  26. xlabel('Zaman(s)')
  27. title('3 ile yukarı örneklenmiş sinüs')
  28.  



1 Hz'lik sinüs işaretimize 0.1 genlik seviyesine sahip bir sinüs işareti eklenerek giriş işareti aşağıdaki gibi olmaktadır.



Aşağıda yukarıda verilen giriş sinyalinin 3 ile yukarı örneklenmesini sağlayan MATLAB kodu ve çıkışta elde edilen sinyal verilmiştir. 20 Hz ile örneklenen siynaln yeni örnekleme frekansı 60 Hz olmuştur.

Kod: Matlab M
  1.  
  2.  
  3. f_s = 20;
  4. f_I1 = 1;
  5. f_I2 = 5;
  6.  
  7. n_t = 1 / f_s : 1 / f_s : 1;
  8.  
  9. s_Input = sin(2 * pi * f_I1 * n_t) + 0.1 * sin(2 * pi * f_I2 * n_t);
  10. figure, plot(n_t, s_Input);
  11. xlabel('Zaman(s)')
  12. title('20 Hzde örneklenmiş 1 Hz ve 5 Hz sinüs toplamı')
  13.  
  14. l_Input = length(s_Input);
  15. l_FFT = 2^(ceil(log2(l_Input)));
  16.  
  17. fft_I = fft(s_Input, l_FFT);
  18. up_L = 3;
  19.  
  20. Zero_Buff = zeros(1, (up_L - 1) * l_FFT);
  21. fft_I_new = up_L * [fft_I(1 : l_FFT / 2)...
  22.    Zero_Buff fft_I(l_FFT / 2 + 1 : l_FFT)];
  23. i_Input = real(ifft(fft_I_new, up_L * l_FFT));
  24. u_Input = i_Input(1 : up_L * l_Input);
  25. n_t = 1 / (up_L * f_s) : 1 / (up_L * f_s) : 1;
  26. figure, plot(n_t, u_Input)
  27. xlabel('Zaman(s)')
  28. title('3 ile yukarı örneklenmiş 1 Hz ve 5 Hz sinüs toplamı')
  29.  


« Son Düzenleme: 26 Aralık 2016, 14:26:22 Gönderen: alicavuslu »