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

@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 27.03.2026 Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. При выборе "Спецтехники" (Тип Справочник. Спецтехника), заполняется. . .	Сумматор с применением элементов трёх состояний. Hrethgir 26.03.2026 Тут. https:/ / fips. ru/ EGD/ ab3c85c8-836d-4866-871b-c2f0c5d77fbc Первый документ красиво выглядит, но без схемы. Это конечно не даёт никаких плюсов автору, но тем не менее. . . всё может быть. . .	Автозаполнение реквизитов при создании документа Maks 26.03.2026 Программный код из решения ниже размещается в модуле объекта документа, в процедуре "ПриСозданииНаСервере". Алгоритм проверки заполнения реализован для исключения перезаписи значения реквизита,. . .	Команды формы и диалоговое окно Maks 26.03.2026 1. Команда формы "ЗаполнитьЗапчасти". Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. В качестве источника данных. . .
Кому нужен AOT? DevAlt 26.03.2026 Решил сделать простой ланчер Написал заготовку: dotnet new console --aot -o UrlHandler var items = args. Split(":"); var tag = items; var id = items; var executable = args;. . .	Отправка уведомления на почту при создании или изменении элементов справочника Maks 24.03.2026 Программная отправка письма электронной почты на примере типового справочника "Склады" в конфигурации БП3. Перед реализацией необходимо выполнить настройку системной учетной записи электронной. . .	модель ЗдравоСохранения 5. Меньше увольнений- больше дохода! anaschu 24.03.2026 Теперь система здравосохранения уменьшает количество увольнений. 9TO2GP2bpX4 a42b81fb172ffc12ca589c7898261ccb/ https:/ / rutube. ru/ video/ a42b81fb172ffc12ca589c7898261ccb/ Слева синяя линия -. . .	Midnight Chicago Blues kumehtar 24.03.2026 Такой Midnight Chicago Blues, знаешь?. . Когда вечерние улицы становятся ночными, а ты не можешь уснуть. Ты идёшь в любимый старый бар, и бармен наливает тебе виски. Ты смотришь на пролетающие. . .

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

	Квадратурные формулы 05.11.2022, 16:10. Показов 1922. Ответов 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

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

Решение