фильтр Баттерворта

@crupps · Регистрация: 12.01.2014

Студворк — интернет-сервис помощи студентам

Добрый вечер!
Подскажите пожалуйста, с программой столкнулась в первые.
у Меня есть текстовый файл в котором лежит сигнал ЭКГ, который необходимо очистить от шума в 50 Гц и от 1ГЦ. Для этого я выбираю полосовой фильтр Баттерворта третьего порядка. Как написать фильтр я поняла:

Matlab M
1
2
3
4
5
6
7
8
9
Rp=2;
Rs=45;
Ws=[7 314];
Wp=[32 252];
[n,Wn]=buttord(Wp,Ws,Rp,Rs,'s'];
w0=Wn;
[b,a]=butter(n,w0,'s');
freqs(b,a);
end.

а вот как вставить сигнал из текстового файла не знаю(( подскажите пожалуйста,как это сделать и правильно ли написан у меня фильтр.
заранее ОГРОМНОЕ спасибо!

@R2D2 · 14.01.2014, 15:42

Сообщение от crupps

что бы убрать частоты 0.05 гц и 50гц

Так все таки 0,05 И 50 или от 0.05 до 50?

@crupps · 14.01.2014, 19:49 **[ТС]**

Сообщение от R2D2

Так все таки 0,05 И 50 или от 0.05 до 50?

оставить от 0.05 до 50 Гц. Все остальное вырезать

@R2D2 · 14.01.2014, 20:17

Вот, пожалуйста. Только фильтр я взял не Баттерворта. Сойдет?

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
clc; close all; clear;
y = dlmread('04.txt');
y(end) = [];
fs = 260;
L = length(y);
nfft = 2^nextpow2(L);
t = (0:L-1)/fs;
fc = [0.05 50];
Rp = 1;
Rs = 90;
[b1, a1] = ellip(5, Rp, Rs, 2*fc/fs, 'bandpass');
Y = filter(b1, a1, y);
%% Plots
plot(t, y, t, Y); grid on

Но что то мне подсказывает, что нужно не это...

@crupps · 15.01.2014, 09:10 **[ТС]**

такое и надо было, просто к изначальному сигналу нужно добавить 0.01 герцовую синусойду( что бы было видно, что мы ее вырезаем потом, а сейчас этого просто не заметно, так как ее нет)
у меня теперь встает вопрос:
как к сигналу y добавить синусоиду с частотой 0.01 Гц?(в итоге должно получится график ЭКГ наложенный на синусоиду(волнистый))

я писала:
f=0.01;
y1=sin(f*t);
а потом пыталась сложить сигналы:
Y=y+y1;
и когда выводила графики получается прямая линия((
но если просто выводит график синуса, то все нормально

@vital792 · 15.01.2014, 13:53

Сообщение от crupps

y1=sin(f*t);

sin(2*pi*f*t)

Новые блоги и статьи Все статьи Все блоги /
Когда логика программы не спасает от человеческих ошибок Maks 18.08.2026 В последнее время всё чаще и чаще сталкиваюсь с таким явлением, как абсолютная невнимательность (или глупость) пользователей. Проявляется это чаще всего на работе в коллективе. Допустим, человек с. . .	Лето уходит kumehtar 17.08.2026	Мысли в слух kumehtar 17.08.2026 Забавно, насколько сейчас стала доступна информация. Например о магии, духовном развитии, медитациях, и других подобных направлениях, ранее зачастую тайных, передаваемых от учителя к ученику. Хотя. . .	Перемещение строк из ТЧ в другой документ с учетом текущего пробега Maks 17.08.2026 Реализация из решения ниже выполнена на примере нетипового документа "Автозапчасти", с ТЧ "Шины". За основу взят алгоритм отсюда: https:/ / www. cyberforum. ru/ blogs/ 359708/ 10838. html Задача: . . .
Саморегулирующийся социальный контракт для сервера cross-section. Hrethgir 14.08.2026 С кодом конечно таких глубоких размышлений пока не было, впрочем я уже привык к алгоритмизации. Суть предмета записи: снова в диалоге с нейросетью (я взял пока себе ник для учётки админа - Rector). . . .	Часы электронные Uhbif79 12.08.2026 Выкладываю программу часов. Программа позволяет: 1. Использовать системное время и дату, 2. Есть возможность вводить время и дату вручную. 3. Реализованы 2 будильника: начало и конец рабочего дня. . . .	Часы с будильником на основе класса QLCDNumber Uhbif79 12.08.2026 Всем добрый день, выкладываю программу часов с будильником на основе класса QLCDNumber. Здесь я пробовал самостоятельно создавал классы, впервые столкнулся с видимостью переменной одного класса из. . .	Установка MinGW GCC 16.2 и CMake 8Observer8 10.08.2026 VK Видео: https:/ / vkvideo. ru/ video-240781534_456239017 YouTube: eY5-5PyI9NM Текстовая версия

@crupps 0 / 0 / 0 Регистрация: 12.01.2014 Сообщений: 14
	15.01.2014, 09:10 [ТС]
	такое и надо было, просто к изначальному сигналу нужно добавить 0.01 герцовую синусойду( что бы было видно, что мы ее вырезаем потом, а сейчас этого просто не заметно, так как ее нет) у меня теперь встает вопрос: как к сигналу y добавить синусоиду с частотой 0.01 Гц?(в итоге должно получится график ЭКГ наложенный на синусоиду(волнистый)) я писала: f=0.01; y1=sin(f*t); а потом пыталась сложить сигналы: Y=y+y1; и когда выводила графики получается прямая линия(( но если просто выводит график синуса, то все нормально 0

Опции темы