Вычисление интегралов в 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 там внутри?? Что вы пытаетесь ввести с командной строки?? Да еще неправильно!

Новые блоги и статьи Все статьи Все блоги /
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. Пошагово создадим проект для загрузки изображения. . .

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

Опции темы

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

Решение