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

2 次查看（过去 30 天）

Harry Smith 2021-5-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/832638-bandpass-filter-in-multiple-frequencies-with-a-for-if-ifelse-loop

评论： 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

each time the script restarts with a new signal, the RMS array is overwritten.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Sulaymon Eshkabilov 2021-5-17

You don't need all of these "if .. elseif ... end" conditional statements. A simple [for ... end] loop.

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 个）

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

类别

MATLAB Graphics 2-D and 3-D Plots Line Plots Two y-axis

在 Help Center 和 File Exchange 中查找有关 Two y-axis 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by