Построение фильтра ФНЧ

@lode00 · Регистрация: 23.12.2017

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

Добрый день!
Помогите кто-нибудь реализовать цифровой ФНЧ!
Код вроде бы написал, но это не то, что я хотел получить...
А пытался я реализовать код, который был написан на Питоне

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
clear;clc
Fs = 8000 % Частота дискретизации (Hz)
Fc = 2000 % Частота среза (Hz)
 
N = Fs %Кол-во образцов за 1 с
 
x = range(0, N)
y = range(0, N)
 
f1 = 100 %Hz
f2 = 3500 % Hz
 
N_TO_DISPLAY = 2*Fs/f1
 
% Генерация тестового сигнала
t = 1:N;
x = (sin(2 * pi * f1 * t/ Fs)+sin(2 * pi * f2 * t / Fs)) / 2;
 
% Расчет необходимых значений для фильтра
dt = 1.0/Fs;
tau = 1.0/(2.0*pi*Fc);
alpha = dt / (dt + 2*tau);
 
% Применение фильтра
i=1:N;
y = alpha.*(x(2:end)+x(1:end - 1));
y = y(2:end) + (1-2*alpha)*y(1:end - 1);
% Графики оригинальных и отфильтрованных спектров
plot(x(1:N_TO_DISPLAY),y(1:N_TO_DISPLAY))
% xlabel("Time")
% ylabel("Sample value")

Добавлено через 1 минуту
Вот код на Python:

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
import numpy as np
import matplotlib.pyplot as plot
 
Fs = 8000 # Sampling frequency (Hz)
Fc = 2000 # Cutoff frequency (Hz)
 
N = Fs # Number of samples per 1 s
 
x = range(0, N)
y = range(0, N)
 
f1 = 100 # Hz
f2 = 3500 # Hz
 
N_TO_DISPLAY = 2*Fs/f1
 
# Generating test signal
for t in range(0, N):
    x[t] = (np.sin(2.0*np.pi*f1*t/Fs) + np.sin(2.0*np.pi*f2*t/Fs))/2.0
 
# Calculating necessary values for filter
dt = 1.0/Fs
tau = 1.0/(2.0*np.pi*Fc)
alpha = dt / (dt + 2*tau)
 
# Applying filter
for i in range(1, N):
    y[i] = alpha*(x[i]+x[i-1]) + (1-2*alpha)*y[i-1]
 
# Drawing original and filtered spectrums
plot.clf()
plot.xlabel("Time")
plot.ylabel("Sample value")
plot.plot(x[0:N_TO_DISPLAY], 'r')
plot.plot(y[0:N_TO_DISPLAY], 'g')
plot.draw()

@Alexis · 15.04.2019, 11:19

У Вас ошибка в выводе графиков.
вместо строки 29 напишите

Matlab M
1
plot(x(1:N_TO_DISPLAY)); hold on; plot(y(1:N_TO_DISPLAY)); hold off;

Сначала будет выводится график "x" (обычно синим цветом), а поверх него - график "y" (красным цветом).

В Вашем же варианте выводятся не графики сигналов "x" и "y", а график зависимости "x" от "y".

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .
И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»	SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .