bandpass filter in multiple frequencies with a for- if- ifelse loop

1 次查看（过去 30 天）
Harry Smith
Harry Smith 2021-5-17
评论： Harry Smith 2021-5-17
Hi everyone!
I have to filter a signal in 5 frequency bands, each with its own range and calculate, for each band, an RMS value. I tried to create a loop, but RMS is overwritten every time a new signal is filtered.
the code:
for freq=1:5
if freq==1
freq_min=0,1
freq_max=0,7
fs=1000
[blow,alow]=butter(2,freq_max/(fs/2));
[bhigh,ahigh]=butter(2,freq_min/(fs/2),'high');
filt_signal=filtfilt((blow,alow,signal);
filt_signal=filtfilt(bhigh,ahigh,filt_signal);
RMS_freq1=rms(filt_signal);
elseif freq==2
freq_min=0,7
freq_max=0,9
fs=1000
[blow,alow]=butter(2,freq_max/(fs/2));
[bhigh,ahigh]=butter(2,freq_min/(fs/2),'high');
filt_signal=filtfilt((blow,alow,signal);
filt_signal=filtfilt(bhigh,ahigh,filt_signal);
RMS_freq2=rms(filt_signal);
ecc..
RMS=[RMS_freq1 RMS_freq2 RMS_freq3 RMS_freq4 RMS_freq5];
end
end
each time the script restarts with a new signal, the RMS array is overwritten.
  4 个评论
显示 2更早的评论
Sulaymon Eshkabilov
Sulaymon Eshkabilov 2021-5-17
You don't need all of these "if .. elseif ... end" conditional statements. A simple [for ... end] loop.
Harry Smith
Harry Smith 2021-5-17
thanks for your help!!
I have 5 frequency bands:
b1: 0.1-0.7 Hz
b2: 0.7-0.9 Hz
b3: 0.9-1.1 Hz
b4 1.1-1.3 Hz
b5: 1.3-1.5 Hz
and i wrote:
freq_min = [0.1; 1.3];
freq_max = [0,7; 1.5];
fs = 1e3;
ecc..
error message: Unable to perform assignment because the indices on the left side are not compatible with the size of the right side.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Matched Filter and Ambiguity Function 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by