Help With Starting Point for a Low Pass Filter
I am brand new to signal processing and having a lot of trouble understanding the documentation ot implement a low pass filter, which is what I think I need. An exmaple of what I’m working with is shown in Fig1.
The signal is oxygen amperometry changes which are necessarily slow. I’ve been attempting to use matlab find peaks function which does a decent job at finding peaks, but I cannot get it to find an appropriate width as it is considering all the smaller peaks instead of the one large one. For example, see Fig2. I’m primarily interested in peaks 1,2, and 4. The findpeaks does a good job with the widths of peaks 1 and 3, but in 4 it’s considering that a small peak. Changing factors within the findpeaks function also doesn’t appear to help as even though the smaller peaks may not appear in the output, findpeaks still considers them in the anlysis. In Fig2 I’ve been able to have findpeaks suppress the output of peaks near peak 4, but that doesn’t help its width estimation.
I’m pretty sure the answer to my issue is preprocessing, but a smoothing such as moving average doesn’t quite capture what I need and so I think I need a more advanced filter like a low pass. Can anyone suggest starting parameters in some code for me to work with ? I’m having trouble getting any implementation to work. An example dataset is in ExData.xlsx where time point 0 is the event of interest.
NOTE:-
Matlabsolutions.com provide latest MatLab Homework Help,MatLab Assignment Help for students, engineers and researchers in Multiple Branches like ECE, EEE, CSE, Mechanical, Civil with 100% output.Matlab Code for B.E, B.Tech,M.E,M.Tech, Ph.D. Scholars with 100% privacy guaranteed. Get MATLAB projects with source code for your learning and research.
I am not certain what result you want.
Here is a prototype lowpass filter design you can experiment with. It has a passband frequency of 0.08 Hz (it stops everything above that) and a reasonably short transition region. I will help you refine it to do what you want, and help with findpeaks as well.
The Filter —
[D,S] = xlsread('ExData.xlsx');
t = D(:,1);
s = D(:,2);figure
plot(t,s) % Plot Original Signal
gridTs = mean(diff(t)); % Sampling Interval
L = numel(t);
Fs = 1/Ts; % Sampling Frequency
Fn = Fs/2; % Nyquist Frequency
smc = s-mean(s); % Subtract Mean (0 Hz)
FTs = fft(smc)/L; % Fourier Transform
Fv = linspace(0,1,fix(L/2)+1)*Fn; % Frequency Vector
Iv = 1:numel(Fv); % Index Vector
SEE COMPLETE ANSWER CLICK THE LINK
