Моделирование работы цифрового фильтра

@ZamyatinDanya · Регистрация: 17.02.2019

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

Здравствуйте, учусь программировать на Python'e. Ну, как учусь, прохожу, так сказать, экспресс-курс. Есть задание: смоделировать работу цифрового фильтра, то есть получить на выходе сигнал, зная отсчеты цифрового сигнала на входе, а также цифровую импульсную характеристику цепи.
С сигналом на входе, вроде, разобрался. Там ничего сложного нет, ибо надо задать обычный прямоугольный импульс

Кликните здесь для просмотра всего текста

Python
1
2
3
4
5
6
7
8
9
10
U = 100
ar = []
for i in range (U):
    if 5 < i < 20:
        ar.append (1)
    else:
        ar.append (0)
print (ar)
pp.plot(ar)
pp.show()

Дальше требуется создать функцию, производящую расчет выходного сигнала. Параметром функции должен являться список с отсчетами сигнала на входе (x), а возвращаемым значением должен быть список с отсчетами сигнала на выходе (y), равный по длине входному. Данная функция должна использовать в качестве отсчетов импульсной характеристики следующий тьюпл, обеспечивающий задержку сигнала на 10 тактов:
coeff = (0, 0, 0, 0, 0, 0, 0, 0, 0, 1)

Вообще не понимаю, что надо сделать. Нет, понятно, что нужно записать выходную функцию, но у нас пока типа неизвестно, что цепь будет интегрирующей или дифференцирующей. По какому принципу должна формироваться эта функция и какой синтаксис следует использовать?

Добавлено через 1 час 15 минут
Окей, что-то получилось

Кликните здесь для просмотра всего текста

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
37
38
39
40
41
t = int(input("Введите постоянную времени цепи"))
U = t*5
print (U, "Длина выходного сигнала")
a = U*0.2
b = U*0.6
#U = 100 #постоянная времени цепи
ar = [] #Входной сигнал
for i in range (U):
    if a < i < b:
        ar.append (1)
    else:
        ar.append (0)
print (ar)
pp.plot(ar)
pp.show()
 
#Выходной сигнал
coeff = [0, 0, 0, 0, 0, 0, 0, 0, 0, 1]
N = int(input("Введите количество отсчетов импульсной характеристики"))
while N <= t:
    ad = []
    ad = [0]*N
    pp.plot (ad)
    pp.show ()
    break
 
h=[]
for n in range (N):
    h.append (((1/t)*exp(-n/t)))
pp.plot(h)
pp.show()
 
n=0
for n in range (U):
    if n <= a:
        continue
    for i in range (N):
        ad[n] = ad[n]+(ad[n-i]*h[i])
pp.plot(ad)
pp.plot(ar)
pp.show()

Однако выскакивает ошибка IndexError: list index out of range В чем дело?

@Garry Galler · 23.03.2020, 17:14

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

выскакивает ошибка IndexError: list index out of range

Где?

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

В чем дело?

list index out of range Вы перешагнули за границы списка.
Понятия не имею какого и где. Это написано в ошибке, которую вы полностью не привели.

Добавлено через 1 минуту

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

ad = []
ad = [0]*N

Не нужно объявлять пустой список, если вы все равно сразу же его перезаписываете списком с нолями.
Не нужно делать break, если вы не ставите никакого условия для break.

@ZamyatinDanya · 23.03.2020, 17:35 **[ТС]**

Сообщение от Garry Galler

Где?

Traceback (most recent call last):
File "D:/Programms/Qt/Work/zamyatin/sem3/44.py", line 38, in <module>
ad[n] = ad[n]+(ad[n-i]*h[i])
IndexError: list index out of range

@Garry Galler · 23.03.2020, 17:51

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

ad[n] = ad[n]+(ad[n-i]*h[i])

Ну это для вас подсказка. Я все равно ваш код не понимаю.
Вы чего-то вставляете в список ad, который формируется за каким-то ... в цикле while
Причем вставляете всегда в индекс 0, потому что n у вас не инкрементируется.
Причем из n вы еще вычитаете единицу и получается 0 -1 = -1 и затем чой-то еще там умножаете... в общем какой-то полный сумбур с индексами.

Добавлено через 2 минуты
Вы список ad вообще проверяли? Скока там индексов есть? Может он у вас вообще пустой....

@ZamyatinDanya · 23.03.2020, 18:32 **[ТС]**

Сообщение от Garry Galler

Вы список ad вообще проверяли?

Не проверял... Надо будет проверить. Главное, что программа работает

Ну, правда, с некоторыми оговорками

@Garry Galler · 23.03.2020, 18:50

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

Главное, что программа работает

Ну, в программировании не принято называть работающими программы выходящие за границы массивов.
Обычно такой подход заканчивается сгоранием в атмосфере новеньких многомиллионных спутников и падением самолетов в океан.
Если вы не курсе - это называется "ошибка на единицу".

@u235 · 23.03.2020, 20:26

ZamyatinDanya, при линейной свертке двух последовательностей итоговая последовательность будет иметь размер n1+n2+1. Видимо где-то в этом у вас ошибка и из-за этого и происходит IndexError

Добавлено через 3 минуты
29-я строка к чему?

Добавлено через 4 минуты
Где у вас свертка с КИХ фильтром длины 10? coeff?
Такое ощущение, что вы не совсем понимаете как ваш код работает и что должен делать..

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Обработчик клика мыши в браузере ПК и касания экрана в браузере на мобильном устройстве 8Observer8 02.02.2026 Содержание блога Для начала пошагово создадим рабочий пример для подготовки к экспериментам в браузере ПК и в браузере мобильного устройства. Потом напишем обработчик клика мыши и обработчик. . .	Философия технологии iceja 01.02.2026 На мой взгляд у человека в технических проектах остается роль генерального директора. Все остальное нейронки делают уже лучше человека. Они не могут нести предпринимательские риски, не могут. . .	SDL3 для Web (WebAssembly): Вывод текста со шрифтом TTF с помощью SDL3_ttf 8Observer8 01.02.2026 Содержание блога В этой пошаговой инструкции создадим с нуля веб-приложение, которое выводит текст в окне браузера. Запустим на Android на локальном сервере. Загрузим Release на бесплатный. . .	SDL3 для Web (WebAssembly): Сборка C/C++ проекта из консоли 8Observer8 30.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .
SDL3 для Web (WebAssembly): Установка Emscripten SDK (emsdk) и CMake для сборки C и C++ приложений в Wasm 8Observer8 30.01.2026 Содержание блога Для того чтобы скачать Emscripten SDK (emsdk) необходимо сначало скачать и уставить Git: Install for Windows. Следуйте стандартной процедуре установки Git через установщик. . . .	SDL3 для Android: Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 29.01.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами. Версия v3 была полностью переписана на Си, в. . .	Инструменты COM: Сохранение данный из VARIANT в файл и загрузка из файла в VARIANT bedvit 28.01.2026 Сохранение базовых типов COM и массивов (одномерных или двухмерных) любой вложенности (деревья) в файл, с возможностью выбора алгоритмов сжатия и шифрования. Часть библиотеки BedvitCOM Использованы. . .	SDL3 для Android: Загрузка PNG с альфа-каналом с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 28.01.2026 Содержание блога SDL3 имеет собственные средства для загрузки и отображения PNG-файлов с альфа-каналом и базовой работы с ними. В этой инструкции используется функция SDL_LoadPNG(), которая. . .

@u235 5516 / 2869 / 571 Регистрация: 07.11.2019 Сообщений: 4,759
	23.03.2020, 20:26
	ZamyatinDanya, при линейной свертке двух последовательностей итоговая последовательность будет иметь размер n1+n2+1. Видимо где-то в этом у вас ошибка и из-за этого и происходит IndexError Добавлено через 3 минуты 29-я строка к чему? Добавлено через 4 минуты Где у вас свертка с КИХ фильтром длины 10? coeff? Такое ощущение, что вы не совсем понимаете как ваш код работает и что должен делать.. 1

Опции темы