
band pass filter design on composite sinusoidal signals of 1 khz, 2khz and 3 khz
hi,
i have 3 sinusoidal signals of frequency 1Khz, 2khz and 3 khz added
up to form a composite signal. i designed one bandpass filter centered
using FDATool in matlab. i tried to give pass band range as small as
1950 to 2050 (remember i want to allow only 2khz rejecting all others)
to as large as 1.5 Khz to 2.5Khz. The designed filter is an FIR filter
with attenuation set at 80. i was able to design many
filters..changing sampling frequency and passband parameters but none
of them gave me my input 2khz signal in output (all gave some
different type of waveforms in output). i am more interested in
extracting out 2khz signal. i am not getting the output same as my
input sinusoidal 2khz signal. i feel..i am missing some basic theory
in this regard..plz do help me in this regard or guide me.
thanks
aizza


0




Reply

aizzaahmed (35)

2/8/2010 5:36:00 PM 

See related articles to this posting
On Mon, 08 Feb 2010 09:36:00 0800, aizza ahmed wrote:
> hi,
> i have 3 sinusoidal signals of frequency 1Khz, 2khz and 3 khz added
> up to form a composite signal. i designed one bandpass filter centered
> using FDATool in matlab. i tried to give pass band range as small as
> 1950 to 2050 (remember i want to allow only 2khz rejecting all others)
> to as large as 1.5 Khz to 2.5Khz. The designed filter is an FIR filter
> with attenuation set at 80. i was able to design many filters..changing
> sampling frequency and passband parameters but none of them gave me my
> input 2khz signal in output (all gave some different type of waveforms
> in output). i am more interested in extracting out 2khz signal. i am
> not getting the output same as my input sinusoidal 2khz signal. i
> feel..i am missing some basic theory in this regard..plz do help me in
> this regard or guide me.
>
> thanks
> aizza
Attenuation set at 80 what? What sort of "different type of waveforms"
were you getting? How are you implementing the filter? Have you looked
at a frequencydomain plot of the filter response, to see how well it
attenuates at 1kHz and 3kHz?
You are specifying a passband, and what sounds like an ultimate
attenuation (I assume you mean "80dB"), but you're not specifying a shape
factor  what attenuation defines your passband ripple, and at what
frequencies does your filter attain 80dB of attenuation?
You may have specified the filter incorrectly.
Matlab may have built the filter incorrectly.
You may be implementing the filter incorrectly.
You need to figure out which of the above has happened, then you can move
on to the next step.

www.wescottdesign.com


0




Reply

Tim

2/8/2010 5:59:05 PM


On Feb 8, 11:36=A0am, aizza ahmed <aizzaah...@gmail.com> wrote:
> hi,
> =A0 =A0i have 3 sinusoidal signals of frequency 1Khz, 2khz and 3 khz adde=
d
> up to form a composite signal. i designed one bandpass filter centered
> using FDATool in matlab. i tried to give pass band range as small as
> 1950 to 2050 (remember i want to allow only 2khz rejecting all others)
> to as large as 1.5 Khz to 2.5Khz. The designed filter is an FIR filter
> with attenuation set at 80. i was able to design many
> filters..changing sampling frequency and passband parameters but none
> of them gave me my input 2khz signal in output (all gave some
> different type of waveforms in output). =A0i am more interested in
> extracting out 2khz signal. i am not getting the output same as my
> input sinusoidal 2khz signal. i feel..i am missing some basic theory
> in this regard..plz do help me in this regard or guide me.
>
> thanks
> aizza
Aizza
FDATool requires more data than you are giving for a bandpass design.
You need to have 5 frequencies: stopband1, passband1, passband2,
stopband2, and the sampling frequency.
Maurice Givens


0




Reply

maury

2/8/2010 6:12:23 PM


On Feb 8, 11:12=A0pm, maury <maury...@core.com> wrote:
> On Feb 8, 11:36=A0am, aizza ahmed <aizzaah...@gmail.com> wrote:
>
>
>
>
>
> > hi,
> > =A0 =A0i have 3 sinusoidal signals of frequency 1Khz, 2khz and 3 khz ad=
ded
> > up to form a composite signal. i designed one bandpass filter centered
> > using FDATool in matlab. i tried to give pass band range as small as
> > 1950 to 2050 (remember i want to allow only 2khz rejecting all others)
> > to as large as 1.5 Khz to 2.5Khz. The designed filter is an FIR filter
> > with attenuation set at 80. i was able to design many
> > filters..changing sampling frequency and passband parameters but none
> > of them gave me my input 2khz signal in output (all gave some
> > different type of waveforms in output). =A0i am more interested in
> > extracting out 2khz signal. i am not getting the output same as my
> > input sinusoidal 2khz signal. i feel..i am missing some basic theory
> > in this regard..plz do help me in this regard or guide me.
>
> > thanks
> > aizza
>
> Aizza
> FDATool requires more data than you are giving for a bandpass design.
> You need to have 5 frequencies: stopband1, passband1, passband2,
> stopband2, and the sampling frequency.
>
> Maurice Givens
hi maurice and wescott,
thanks for response. i am giving all details. plz let me know
if you need any.
Each signal i.e. (1khz, 2khz and 3khz sinusoidal signal is sampled at
96khz and added to form a composite signal. and i need to extract 2khz
signal by passing it through bandpass filter designed below). i used
% All frequency values are in Hz.
Fs =3D 96000; % Sampling Frequency
Fstop1 =3D 1900; % First Stopband Frequency
Fpass1 =3D 1950; % First Passband Frequency
Fpass2 =3D 2050; % Second Passband Frequency
Fstop2 =3D 2100; % Second Stopband Frequency
Astop1 =3D 80; % First Stopband Attenuation (dB)
Apass =3D 1; % Passband Ripple (dB)
Astop2 =3D 80; % Second Stopband Attenuation (dB)
match =3D 'passband'; % Band to match exactly
i tried in all possible types of FIR filter based designs available in
matlab..but still couldnt get the original signal after filtering. The
signals i get are usually chirp signals and they dont show sinusoidal
shape.
thanks
aizza


0




Reply

aizza

2/9/2010 2:30:41 PM


On Feb 9, 8:30=A0am, aizza ahmed <aizzaah...@gmail.com> wrote:
> On Feb 8, 11:12=A0pm, maury <maury...@core.com> wrote:
>
>
>
>
>
> > On Feb 8, 11:36=A0am, aizza ahmed <aizzaah...@gmail.com> wrote:
>
> > > hi,
> > > =A0 =A0i have 3 sinusoidal signals of frequency 1Khz, 2khz and 3 khz =
added
> > > up to form a composite signal. i designed one bandpass filter centere=
d
> > > using FDATool in matlab. i tried to give pass band range as small as
> > > 1950 to 2050 (remember i want to allow only 2khz rejecting all others=
)
> > > to as large as 1.5 Khz to 2.5Khz. The designed filter is an FIR filte=
r
> > > with attenuation set at 80. i was able to design many
> > > filters..changing sampling frequency and passband parameters but none
> > > of them gave me my input 2khz signal in output (all gave some
> > > different type of waveforms in output). =A0i am more interested in
> > > extracting out 2khz signal. i am not getting the output same as my
> > > input sinusoidal 2khz signal. i feel..i am missing some basic theory
> > > in this regard..plz do help me in this regard or guide me.
>
> > > thanks
> > > aizza
>
> > Aizza
> > FDATool requires more data than you are giving for a bandpass design.
> > You need to have 5 frequencies: stopband1, passband1, passband2,
> > stopband2, and the sampling frequency.
>
> > Maurice Givens
>
> hi maurice and wescott,
> =A0 =A0 =A0 =A0 =A0thanks for response. i am giving all details. plz let =
me know
> if you need any.
>
> Each signal i.e. (1khz, 2khz and 3khz sinusoidal signal is sampled at
> 96khz and added to form a composite signal. and i need to extract 2khz
> signal by passing it through bandpass filter designed below). i used
>
> =A0% All frequency values are in Hz.
> =A0Fs =3D 96000; =A0% Sampling Frequency
> =A0Fstop1 =3D 1900; =A0 =A0 =A0 =A0% First Stopband Frequency
> =A0Fpass1 =3D 1950; =A0 =A0 =A0 =A0% First Passband Frequency
> =A0Fpass2 =3D 2050; =A0 =A0 =A0 =A0% Second Passband Frequency
> =A0Fstop2 =3D 2100; =A0 =A0 =A0 =A0% Second Stopband Frequency
> =A0Astop1 =3D 80; =A0 =A0 =A0 =A0 =A0% First Stopband Attenuation (dB)
> =A0Apass =A0=3D 1; =A0 =A0 =A0 =A0 =A0 % Passband Ripple (dB)
> =A0Astop2 =3D 80; =A0 =A0 =A0 =A0 =A0% Second Stopband Attenuation (dB)
> =A0match =A0=3D 'passband'; =A0% Band to match exactly
>
> i tried in all possible types of FIR filter based designs available in
> matlab..but still couldnt get the original signal after filtering. The
> signals i get are usually chirp signals and they dont show sinusoidal
> shape.
>
> thanks
> aizza Hide quoted text 
>
>  Show quoted text 
Seems to work for me. I suggest you doublecheck your procedures
Maurice Givens


0




Reply

maury

2/10/2010 9:30:25 PM


On Tue, 09 Feb 2010 06:30:41 0800, aizza ahmed wrote:
> On Feb 8, 11:12 pm, maury <maury...@core.com> wrote:
>> On Feb 8, 11:36 am, aizza ahmed <aizzaah...@gmail.com> wrote:
>>
>>
>>
>>
>>
>> > hi,
>> > i have 3 sinusoidal signals of frequency 1Khz, 2khz and 3 khz
>> > added
>> > up to form a composite signal. i designed one bandpass filter
>> > centered using FDATool in matlab. i tried to give pass band range as
>> > small as 1950 to 2050 (remember i want to allow only 2khz rejecting
>> > all others) to as large as 1.5 Khz to 2.5Khz. The designed filter is
>> > an FIR filter with attenuation set at 80. i was able to design many
>> > filters..changing sampling frequency and passband parameters but none
>> > of them gave me my input 2khz signal in output (all gave some
>> > different type of waveforms in output). i am more interested in
>> > extracting out 2khz signal. i am not getting the output same as my
>> > input sinusoidal 2khz signal. i feel..i am missing some basic theory
>> > in this regard..plz do help me in this regard or guide me.
>>
>> > thanks
>> > aizza
>>
>> Aizza
>> FDATool requires more data than you are giving for a bandpass design.
>> You need to have 5 frequencies: stopband1, passband1, passband2,
>> stopband2, and the sampling frequency.
>>
>> Maurice Givens
>
> hi maurice and wescott,
> thanks for response. i am giving all details. plz let me know
> if you need any.
>
> Each signal i.e. (1khz, 2khz and 3khz sinusoidal signal is sampled at
> 96khz and added to form a composite signal. and i need to extract 2khz
> signal by passing it through bandpass filter designed below). i used
>
> % All frequency values are in Hz.
> Fs = 96000; % Sampling Frequency
> Fstop1 = 1900; % First Stopband Frequency
> Fpass1 = 1950; % First Passband Frequency
> Fpass2 = 2050; % Second Passband Frequency
> Fstop2 = 2100; % Second Stopband Frequency
> Astop1 = 80; % First Stopband Attenuation (dB)
> Apass = 1; % Passband Ripple (dB)
> Astop2 = 80; % Second Stopband
> Attenuation (dB) match = 'passband'; % Band to match exactly
>
> i tried in all possible types of FIR filter based designs available in
> matlab..but still couldnt get the original signal after filtering. The
> signals i get are usually chirp signals and they dont show sinusoidal
> shape.
Perhaps your problem is that you are blindly applying library functions
from Matlab without really understanding what you're doing.
So you apply a filter design, and you get back this ENORMOUS vector of
numbers  right? What does the FFT of that vector look like? What you
_should_ get is something whose amplitude, after normalizing by the
vector size and sampling rate, is the desired amplitude response of your
filter. Are you getting that? And do you know why you should?

www.wescottdesign.com


0




Reply

Tim

2/11/2010 12:24:03 AM



5 Replies
644 Views
Similar Articles
[PageSpeed]
35

Similar Artilces:
just like half band filters are ther 1/3 band filters too?Hello,
Half band filters are good for interpolation and decimation by
2 as they have alternating 0's as their coefficients.
Is there something like 1/3rd band filters for doing interpolation
and decimation by 3 with inbetween 2 coeffcients being 0?
Regards
Bharat
On 2/2/2011 6:16 AM, bharat pathak wrote:
> Hello,
>
> Half band filters are good for interpolation and decimation by
> 2 as they have alternating 0's as their coefficients.
>
> Is there something like 1/3rd band filters for doing interpolation
> ... i need to design low pass , high pass, bandpass filtr with cut off frequency 6,18,20 khz. pls guide mewhen i type these functions in matlab , i get default fstop, astop, etc. how could i correlate these things with that of cut off frequency and gain of filter. pls help me.....
"nethaji anandhavalli" <mydreamprojects@yahoo.com> wrote in message <hq7fu8$rbb$1@fred.mathworks.com>...
> when i type these functions in matlab , i get default fstop, astop, etc. how could i correlate these things with that of cut off frequency and gain of filter. pls help me.....
Hi Nethaji, I assume you are talking about using a filter specification object. You do not state in your post ... Narrow pass band filter designDear all,
This is in ref. to below link:
http://www.mathworks.com/support/solutions/en/data/116LFT/?solution=116LFT
In above link, they mention 0.012 as very small. What is the threshold value required for ratio (2*fc/fs) to design a stable IIR Butterworth low pass filter? Here fc is cutoff frequency & fs is sampling frequency.
Regards,
abhishek
"Abhishek Ballaney" <denim83@yahoo.com> wrote in message <igk010$pv4$1@fred.mathworks.com>...
> Dear all,
>
> This is in ref. to below link:
> http://www.mathworks.com/support/solutions/en/dat... How can i design a bandpass filter with scipy?Hi all,
i would like to design a high pass filter with scipy.signal module.
This is the code i'm using to:
import scipy.signal as signal
import scipy
#first of all i design the lowpass fir filter. This is a 10 taps filter
with cutoff frequency =1 (as help tell me to do)
lpwindow=signal.firwin(10,1)
#with the following instruction i'm creating a band pass filter from
the low pass one
bpwindow=signal.lp2bp(lpwindow,1,0.5,0.2)
My problem is that the bandpass filter obtained with lp2bp function is
16 taps one!
How is it possible??
thanks,
Vincenzo
LabWINC wrote:
> Hi all,
>... IIR band pass filter design queryHi,
i have to design a digital Band pass filter using butter worth
4th order filter.... the requirement is that the pass band frequency i
1KHz and the sampling frequency is 48KHz. the filter should have a ver
high Q factor...no band is specified as such
my question is
1. Q factor is related to analog filter and how is it translated to
digital filter design
2. what is the best Q factor that could be got and how to measure it
3. any matlab functions to implement it like butter()
any other related material would be good...
thanks
srikk
On Aug 29, 8:03 am, "srikk123... Low pass filter design with complexsignal inputHi,
I am recently trying to use lowpass filter under Simulink. My input signal is complexvalued which is BPSK with PNsequence after AWGN channel, but I've tried "FDA tool" and lowpass filter and RF lowpass filter, the error message are all the same, as "Complex signal mismatch. Input port 1 of 'testfile/Lowpass RF Filter/RF Filter/Out' expects a signal of numeric type real. However, it is driven by a signal of numeric type complex" like that.
How can I deal with this problem? I think there should be somewhere parameter setting under the lowpass filter ... Signal processing; Fourier Transformation and Band Pass FilteringHello All;
I have a homework; in signal processing,
What i need is that i ha ve a data like;
ID miliseconds microvolts
0 1024.00 .2811
1 1022.00 .3166
2 1020.00 .0769
3 1018.00 .2151
4 1016.00 .3521
5 1014.00 .5004
6 1012.00 .3115
7 1010.00 .2304
8 1008.00 .2887
9 1006.00 .1213
......
1016 1008.00 .1873
1017 1010.00 .1467
1018 ... How to apply 1/3 octave bandpass filter of about 10Khz on a signalI would like to know if octave filtering is possible in Labview.
There is a toolkit for octave analysis. <a
href="http://sine.ni.com/apps/we/nioc.vp?cid=3122&lang=US">NI LabVIEW
Sound and Vibration Toolset</a>
Hello Pavitra,
The Sound and Vibration Toolset has many powerful tools for octave
analysis. Please refer to chapter 8 of the <a
href="http://digital.ni.com/manuals.nsf/webAdvsearch/C31B0809E0EE1934862567B5004F30F6?OpenDocument&vid=niwc&node=132100_US">Sound
and Vibration Toolset Reference Manual</a> for more information on
these... How to apply 1/3 octave bandpass filter of about 10Khz on a signal #2I am unable to retrieve the time signal after doing the 1/3 octave
analysis.
Can you post a simple example program?
...



