Квадратурные формулы

@VictorVAlduin · Регистрация: 15.12.2020

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

Получить таблицу значений интеграла
k = 6-10

Подскажите в пайтоне есть какая то библиотека для реализации этого метода?
Например в матлабе я реализовывал вручную через циклы все это

@VictorVAlduin · 05.11.2022, 18:17 **[ТС]**

Квадратурная формула с весом!

Добавлено через 9 минут
ПОдскажите вот я все реализовал.
А как указать вес?

Python
1
2
3
4
5
6
7
8
import numpy as np
from scipy import integrate
for k in range(6, 11):
    a = 1.0+0.2*k
    print(a)
 
f = lambda x: (np.cos(a*x) / a+x**2)*np.exp(-x**2)
print(integrate.quad(f, -np.inf, np.inf))

Добавлено через 17 минут
Например используя квадратурные формулы с весом я бы воспользовался формулой Чебышева-Лагерра

@Red white socks · 05.11.2022, 21:04

Документацию читаем по функциям, которые используем
https://docs.scipy.org/doc/sci... grate.quad

@VictorVAlduin · 05.11.2022, 21:12 **[ТС]**

Может не тот вес указываю?

Python
1
2
3
4
5
6
for k in range(6, 11):
    a = 1.0 + 0.2 * k
    print(a)
 
f = lambda x1: (np.cos(a * x1) / (a + x1 ** 2)) * np.exp(-x1 ** 2)
print(integrate.quad(f, -np.inf, np.inf, weight='cos', wvar=a))

Python
1
2
    raise ValueError("Cannot integrate with this weight from -Inf to +Inf.")
ValueError: Cannot integrate with this weight from -Inf to +Inf.

@Red white socks · 05.11.2022, 23:05

Интегрирование с весовой функцией поддерживает только один бесконечный предел.
Весовую функцию вы должны убирать из f.
Теперь плохие новости для вас. Та в.функция, что вам дана не поддерживается модулем. Боюсь, как бы вам не пришлось писать квадратурную формулу самостоятельно.
Ну и старые грабли: опять бессмысленный цикл:

Python
1
2
3
for k in range(6, 11):
    a = 1.0+0.2*k
    print(a)

@VictorVAlduin · 05.11.2022, 23:11 **[ТС]**

Ну я попробовал сделал вручную на матлабе
Может попробую сделаю также вручную через циклы на пайтоне

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
xk = [0.263560 1.413403 3.596425 7.085810 12.640801];
Ak = [0.521756 0.398667 0.075942 0.003612 0.000023];
 
k = 6:1:10;
a = 1.0 + 0.2*k
s = []
for i = 1:length(a)
    s1 = 0;
    for j = 1:length(xk)
        s1 = s1 + Ak(j)*cos(a(i)*xk(j))/(a(i)+xk(j))
    end
    s = [s, s1]
end
s

@u235 · 06.11.2022, 06:45

Сообщение от Red white socks

Интегрирование с весовой функцией поддерживает только один бесконечный предел.

Если не ошибаюсь, то подинтегральная функция - четная, а значит можно заменить один бесконечный предел на ноль и удвоить результат...

@Red white socks · 06.11.2022, 07:56

u235, да, в том числе

Новые блоги и статьи Все статьи Все блоги /
Отображение реквизитов в документе по условию и контроль их заполнения Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеСпецтехники", разработанного в конфигурации КА2. Данный документ берёт данные из другого нетипового документа. . .	Фото всей Земли с борта корабля Orion миссии Artemis II kumehtar 04.04.2026 Это первое подобное фото сделанное человеком за 50 лет. Снимок называют новым вариантом легендарной фотографии «The Blue Marble» 1972 года, сделанной с борта корабля «Аполлон-17». Новое фото. . .	Вывод диалогового окна перед закрытием, если документ не проведён Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать программный контроль на предмет проведения документа. . .	Программный контроль заполнения реквизита табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать контроль заполнения реквизита "ПричинаСписания". . .
wmic не является внутренней или внешней командой Maks 02.04.2026 Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/	Программная установка даты и запрет ее изменения Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: при создании документов установить период списания автоматически. . .	Вывод данных в справочнике через динамический список Maks 01.04.2026 Реализация из решения ниже выполнена на примере нетипового справочника "Спецтехника" разработанного в конфигурации КА2. Задача: вывести данные из ТЧ нетипового документа. . .	Программное заполнения текстового поля в реквизите формы документа Maks 01.04.2026 Алгоритм из решения ниже реализован на нетиповом документе "ВыдачаОборудованияНаСпецтехнику" разработанного в конфигурации КА2, в дополнении к предыдущему решению. На форме документа создается. . .

@VictorVAlduin 8 / 5 / 1 Регистрация: 15.12.2020 Сообщений: 349

	Квадратурные формулы 05.11.2022, 16:10. Показов 1926. Ответов 7 Метки нет (Все метки) Получить таблицу значений интеграла k = 6-10 Подскажите в пайтоне есть какая то библиотека для реализации этого метода? Например в матлабе я реализовывал вручную через циклы все это Миниатюры 0

@Red white socks 4523 / 1899 / 336 Регистрация: 18.01.2021 Сообщений: 3,489
	05.11.2022, 21:04
	Документацию читаем по функциям, которые используем https://docs.scipy.org/doc/sci... grate.quad 0

@Red white socks 4523 / 1899 / 336 Регистрация: 18.01.2021 Сообщений: 3,489
	06.11.2022, 07:56
	u235, да, в том числе 0

Опции темы

Квадратурные формулы

Решение