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

@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?
Такое ощущение, что вы не совсем понимаете как ваш код работает и что должен делать..

Новые блоги и статьи Все статьи Все блоги /
Подключение Box2D v3 к SDL3 для Android: физика и отрисовка коллайдеров 8Observer8 29.01.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами. Версия v3 была полностью переписана на Си, в. . .	Инструменты COM: Сохранение данный из VARIANT в файл и загрузка из файла в VARIANT bedvit 28.01.2026 Сохранение базовых типов COM и массивов (одномерных или двухмерных) любой вложенности (деревья) в файл, с возможностью выбора алгоритмов сжатия и шифрования. Часть библиотеки BedvitCOM Использованы. . .	Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 28.01.2026 Содержание блога SDL3 имеет собственные средства для загрузки и отображения PNG-файлов с альфа-каналом и базовой работы с ними. В этой инструкции используется функция SDL_LoadPNG(), которая. . .	Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL3_image 8Observer8 27.01.2026 Содержание блога SDL3_image - это библиотека для загрузки и работы с изображениями. Эта пошаговая инструкция покажет, как загрузить и вывести на экран смартфона картинку с альфа-каналом, то есть с. . .
Влияние грибов на сукцессию anaschu 26.01.2026 Бифуркационные изменения массы гриба происходят тогда, когда мы уменьшаем массу компоста в 10 раз, а скорость прироста биомассы уменьшаем в три раза. Скорость прироста биомассы может уменьшаться за. . .	Воспроизведение звукового файла с помощью SDL3_mixer при касании экрана Android 8Observer8 26.01.2026 Содержание блога SDL3_mixer - это библиотека я для воспроизведения аудио. В отличие от инструкции по добавлению текста код по проигрыванию звука уже содержится в шаблоне примера. Нужно только. . .	Установка Android SDK, NDK, JDK, CMake и т.д. 8Observer8 25.01.2026 Содержание блога Перейдите по ссылке: https:/ / developer. android. com/ studio и в самом низу страницы кликните по архиву "commandlinetools-win-xxxxxx_latest. zip" Извлеките архив и вы увидите. . .	Вывод текста со шрифтом TTF на Android с помощью библиотеки SDL3_ttf 8Observer8 25.01.2026 Содержание блога Если у вас не установлены Android SDK, NDK, JDK, и т. д. то сделайте это по следующей инструкции: Установка Android SDK, NDK, JDK, CMake и т. д. Сборка примера Скачайте. . .

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