Вычисление интегралов в Python

@Dionisius · Регистрация: 19.10.2013

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

Начал писать функцию для алгоритмов автокорреляционного и спектрального анализа, требуется использовать разложение в ряд Фурье по следующей формуле:
u(t) = a0 + E(m=1..inf)(a[m]*cos(Wm*t)+b[m]*sin(Wm*t)

При вычислении интеграла выводится следующая ошибка:

File "C:/Users/M.video/Desktop/autocor.py", line 40, in <module>
a = float(input((2/D)*integrate.quad(f1,Tr, Tl)))

TypeError: can't multiply sequence by non-int of type 'float'

Как с этим бороться?

Фрагмент кода дан ниже:

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
import numpy as np
import math
import matplotlib.pyplot as plt
import mpl_toolkits.mplot3d.axes3d as p3
import matplotlib.animation as animation
from scipy import integrate
import scipy.special as special
 
from numpy.fft import rfft, rfftfreq
from numpy.random import uniform
from math import sin, cos, pi, sqrt
 
def show(fun, arg):
 
    print(type(fun), ':', fun)
    print('arg =',arg,'=> fun(arg) =', fun(arg))
        
a=0.01
b=1
 
Tr=2
Tl=1
D=Tr-Tl
h=0.01
n=D/h
m=100
Vm = m/D
Tm = D/m
sum =0
a0=(1/D)*(Tl-Tr)
Am=0
 
for j in range (100):
    f1 = lambda j: u(t)*cos((2*pi*j)/D)
    f2 = lambda j: u(t)* sin((2*pi*j)/D)
    a = float(input((2/D)*integrate.quad(f1,Tr, Tl)))
    b = float(input((2/D)*integrate.quad(f2,Tr, Tl)))
    sum+=(a*cos(2*pi*j/D)+b*sin(2*pi*j/D))
    Am+=sqrt(a*a+b*b)
U=a0+sum

@grgdvo · 27.01.2019, 13:59

См. мануал. Функция quad возвращает несколько результатов (последовательность). Вам надо как-то так написать

Python
1
2
a, aerr = integrate.quad(...)
b, berr = integrate.quad(...)

@Dionisius · 27.01.2019, 15:18 **[ТС]**

С интегралами разобрался, теперь появились сомнения о вычислении суммы квадратов двух чисел (сумма под квадратным корнем):

Python
1
sqsum = float(sqrt(pow(a,2)+pow(b,2)))

Компилятор выдает следующую ошибку:
File "C:/Users/M.video/Desktop/autocor.py", line 53, in <module>
sqsum = float(sqrt(pow(a,2)+pow(b,2)))

TypeError: unsupported operand type(s) for ** or pow(): 'tuple' and 'int'

@grgdvo · 27.01.2019, 17:23

У вас осталась такая же ошибка.
"a" содержит набор значений вместо одного значения

Добавлено через 7 минут
И да, увидел не сразу... А зачем вам input там внутри?? Что вы пытаетесь ввести с командной строки?? Да еще неправильно!

Новые блоги и статьи Все статьи Все блоги /
Инструменты COM: Сохранение данный из VARIANT в файл и загрузка из файла в VARIANT bedvit 28.01.2026 Сохранение базовых типов COM и массивов (одномерных или двухмерных) любой вложенности (деревья) в файл, с возможностью выбора алгоритмов сжатия и шифрования. Часть библиотеки BedvitCOM Использованы. . .	Киев стоит - украинская песня zorxor 28.01.2026 wfWdiRqdTxc О Господи, Вечный, Ты . . . Я помоги, Бесконечный. . . Я прошу Ты. . . Я погибаю, спаси. . . Я прошу Тебя Вечный. . .	Загрузка 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 и т. д. Сборка примера Скачайте. . .

@grgdvo 1263 / 977 / 384 Регистрация: 02.09.2012 Сообщений: 3,021
	27.01.2019, 17:23
	У вас осталась такая же ошибка. "a" содержит набор значений вместо одного значения Добавлено через 7 минут И да, увидел не сразу... А зачем вам input там внутри?? Что вы пытаетесь ввести с командной строки?? Да еще неправильно! 0

Вычисление интегралов в Python

Решение