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

@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)

Новые блоги и статьи Все статьи Все блоги /
Из невошедшего на форум (диалог с ИИ-гугла) zorxor 29.07.2026 А вот, что интересно, сказал мне ИИ-гугла: Этот текст — эмоциональный пост пользователя под ником zorxor на интернет-форуме (вероятно, посвященном мистике, непознанному или альтернативной науке). . . .	Был праздник вчера, а я и не знал. kumehtar 28.07.2026 27. 07. 2026г. Intel Core 2 Duo исполнилось 20 лет Новости компьютерного мира и их обсуждение (4) Салют, шампанское, овации! :drink:	Нейтральные знания, чистый код - бла-бла-бла-бла, на самом деле кликбейт и самореклама, плагиат, и вот почему Hrethgir 27.07.2026 То-есть отклонение такой публикации говорит само за себя, и пусть только возьмут на вооружение после отклонения публикации - это будет чистейшим актом плагиата. Отклонял Хабр. Дословно, отклонённая. . .	тв 16 бой ии anaschu 27.07.2026 Великий Перелом ИИ: Как уравнения ОДУ Radau дожали цензурные фильтры Алисы Фиксируем в мемофонде Теории Всего беспрецедентный факт в истории ИИ-зондирования. В затяжном многораундовом. . .
мв 15. непроверенное, возможно, глюк anaschu 27.07.2026 НАУЧНО-АНАЛИТИЧЕСКИЙ ОТЧЕТ. РАЗДЕЛ 1. 1: «НАУКА» (РАСШИРЕННАЯ СТЕХИОМЕТРИЧЕСКАЯ И ГЕНЕТИЧЕСКАЯ ВЕРСИЯ)Тема: Теоретическое обоснование инвариантности 19-мерного тензорного ядра непрерывных ОДУ и. . .	Очистка реквизитов и табличных частей документа при копировании (вариант 2) Maks 26.07.2026 Алгоритм из решения ниже разработан на примере нетипового документа "ЗаявкаНаРаботу", разработанного в КА2. Задача: Заменить алгоритм запрета копирования документов для сотрудников с ролью "Стажер",. . .	Доктрина интенционального знания - Доктрина для портала "Срез". Hrethgir 25.07.2026 Может найдётся кто захочет оценить доктрину. . . Написания правил участия для меня роскошь, требующая лимита времени, поэтому все сообщения не прошедшие модерацию будут видны только участникам портала,. . .	сукцессия 44. Решил подать на припринт в межународные сервисы препринтов. Но нужно одобрение от ученых anaschu 25.07.2026 Английский вариант. Пока кто то не одобрит мою личность, мне не получиться это опубликовать на препринте. Но заявку на публикацию статьи я сегодня подам.

@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