Obliczenie całki w pythonie

Hej, chciałem policzyć całkę :

jako wynik chciałbym otrzymać wykres. Z tego co się orientuje to wynik całki to pojedyncza liczba, wiec czy żeby mieć wykres to trzeba liczyć np metodą trapezów całkę miedzy dwoma sąsiednimi punktami?
Mam sygnał i odpowiadające mu próbki czasu nagrane z urządzenia.

for n in range(len(czas)):
    napiecie.append(skala * (zmierzone_napiecie[n] * math.cos(2 * math.pi * 50 * czas[n])))

gdzie : skala to 2/N, zmierzone_napiecie[n] to u(t) we wzorze a czas[n] to kolejne próbki czasu , 50 to częstotliwość w Hz

Co jeszcze należy zrobić z tym wzorem w pythonie , by odpowiadał tej całce?

Co tak naprawdę chcesz osiągnąć? Wskazałeś nielogiczność swojego rozumowania - całka daje wynik liczbowy, zaś wykres można zrobić dla funkcji. Więc skup się na tym. No chyba że chcesz zrobić wykres funkcji wartości całki od przyjętych granic funkcji na przykład.

elwis napisał(a):

Co tak naprawdę chcesz osiągnąć? Wskazałeś nielogiczność swojego rozumowania - całka daje wynik liczbowy, zaś wykres można zrobić dla funkcji. Więc skup się na tym. No chyba że chcesz zrobić wykres funkcji wartości całki od przyjętych granic funkcji na przykład.

Tak, mam na myśli wykres funkcji wartości całki dla kolejnych chwil czasowych (oś x).

mam na myśli wykres funkcji wartości całki dla kolejnych chwil czasowych (oś x).

Jak Chcesz zrobić wykres całki od czasu, przecież czas jest w całce zmienną związaną (w całce liczonej analitycznie, to będzie granica przy t -> inf). Możesz zrobić wykres całki w zależności od funkci u, granic całkowania, itd.

lion137 napisał(a):

mam na myśli wykres funkcji wartości całki dla kolejnych chwil czasowych (oś x).

Jak Chcesz zrobić wykres całki od czasu, przecież czas jest w całce zmienną związaną (w całce liczonej analitycznie, to będzie granica przy t -> inf). Możesz zrobić wykres całki w zależności od funkci u, granic całkowania, itd.

W zasadzie to przedział chciałem określić jako np liczba n i to jest np 200, więc od 0 do 200 i z krokiem takim jaki jest na osi x

Użyj SymPy: https://docs.sympy.org/latest/tutorial/calculus.html#integrals

Jesli masz probki funkcj to pewnie jakas interpolacja wchodzi gre (bo probki => masz zbior dyskretny) i potem dopiero calka

Ta całka, którą zilustrowałeś na załączonym na obrazku definiuje funkcję U w zależności od t_calc. Po całkowaniu względem czasu uzyskasz funkcję, której argumentem jest t_calc.
Chcąc zrobić wykres tej funkcji potrzebujesz wyliczyć wartości dla różnych t_calc. Może łatwiej to zobaczyć jak podstawisz sobie x=t_calc .

A nie chcesz przypadkiem policzyć po prostu numpy.cumsum ?