Ряд Фурье

@heplme · Регистрация: 11.11.2021

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

Добрый день, плохо разбираюсь в языке Python и хочу спросить как исправить ошибку "can't multiply sequence by non-int of type 'Pi'" Вот мой код:

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
import sympy 
import numpy as np
import matplotlib.pyplot as plt 
 
 
t = sympy.Symbol('t', real = 'True', nonzero = 'True') 
T = 1
A = 1
 
a0 = A/2
k = [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20] 
ak = 0
bk= (-A/(k*sympy.pi))*((-1)**k - 1)
 
for i in range(2,20):
    
    res += ak*sympy.cos(k*t) + bk*sympy.sin(k*t)
 
Fx=a0/2 + res
 
currFun = Fx
 
NumberOfSamples = 200
time = np.linspace(0, T, NumberOfSamples)
samples = sympy.lambdify(t, currFun, 'numpy')(time)
 
plt.plot(time, samples)

Программа должна по указанным коэффициентам Фурье (а0, аk, bk) нарисовать график. К это числа от 1 до бесконечности. Я указала от 1 до 20. Вот общая формула на всякий случай:

@Catstail · 20.12.2021, 11:12

У тебя k - это не целое число, а список. Список нельзя умножить на пи.

@Gdez · 20.12.2021, 12:06

heplme, если не поздно, без построения графика похоже так:

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
import numpy as np
 
...
...
T = 1
A = 1
 
a0 = A/2
k = np.arange(0,21)
 
ak = np.zeros(21)
bk= (-A/(k*np.pi))*((-1)**k - 1)
 
NumberOfSamples = 200
t = np.linspace(0, T, NumberOfSamples).reshape(-1,1)
fun = lambda x: np.nansum(ak*np.cos(k*x) + bk*np.sin(k*x), 1)
 
Fx=a0/2 + fun(t)
currFun = Fx
...
...

@FFPowerMan · 20.12.2021, 12:19

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

Python
1
2
for i in range(2,20):
    res += ak*sympy.cos(k*t) + bk*sympy.sin(k*t)

- и где i в этой формуле?

Новые блоги и статьи Все статьи Все блоги /
Инструменты 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 и т. д. Сборка примера Скачайте. . .	Использование SDL3-callbacks вместо функции main() на Android, Desktop и WebAssembly 8Observer8 24.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .

@Catstail Супер-модератор 38169 / 21104 / 4307 Регистрация: 12.02.2012 Сообщений: 34,693 Записей в блоге: 14
	20.12.2021, 11:12
	У тебя k - это не целое число, а список. Список нельзя умножить на пи. 0

Ряд Фурье

Решение