Create a notch filter for complex input data
7 次查看(过去 30 天)
显示 更早的评论
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/147981/image.png)
I have a complex sampled signal where I want to filter out the signal at 244140 Hz with an Fs = 390625.
What I tried was:
% See the pic for what the signal looks like
d = fdesign.notch('N,F0,BW,Ap',200,0.625,0.045,1);
Hd = design(d);
y = filter(Hd,myComplexData);
But
- The notch didn't show up where I thought it would
- It is not complex
Now because it is not generating a filter with complex coefficients I understand why the notch is at where it is at.
My questions is: How do I create a complex coefficient filter?
回答(2 个)
Uladzimir
2015-8-12
编辑:Uladzimir
2015-8-12
Maybe it will help you
F = 244140;
Fs = 390625;
NFFT = length(YourSignal);%or 2^n
YourSpectra = fft(YourSignal,NFFT);
df = Fs/NFFT;% or, if your signal is complex, this must be a known parameter
Ns = round(F/df)+1;
ComplexHarmonic = 2*YourSpectra(Ns);
2 个评论
Uladzimir
2015-8-13
I just showed an idea with frequency transform. If you have already complex signal, 2x isn't needed, just in a case of fourier transform. You have a signal in a frequency domain. And need a value of one bin on the definite frequency. So, you just need to know the number of this bin.
Ns = round(F/df);
ComplexValue = ComplexSignal(Ns);
Billy Kihei
2021-1-22
Multiplication in frequency domain is convolution in time domain. (this works for real or complex)
- Take x[t] -> fft -> Y[f]
- zero out the values in Y[f] you want to get rid of -> Ymanuallyfiltered[f]
- Ymanuallyfiltered[f]->ifft->x_filtered[t]
- x_filtered[t]->fft->Yfiltered[f]
Note: This may make your code take much longer to run, but it's a quick and dirty way without designing a filter. It is important that you convert it back to time domain after your manually editing. Then re-fft it to get the natural frequency of the response after filtering. Hope this helps.
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Spectral Analysis 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!