Обратная интерполяция

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

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

Нужна помощь! Не успеваю ребята пишу дипломную работу.Буду благодарен за помощь!
Посредством обратной интерполяции определить с заданной точностью ε корень
уравнение, размещенное на отрезке [a;b], с шагом h.
Можно с помощью SYMPY
У меня есть пример кода на матлаб

eps = 1e-5
a = 0.3
b = 0.4
h = 0.02

@u235 · 30.10.2022, 14:15

Сообщение от VictorVAlduin

У меня есть пример кода на матлаб

Пример кода на Matlab нужно выкладывать не скриншотом, а текстом в тегах Matlab.

@VictorVAlduin · 30.10.2022, 14:23 **[ТС]**

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
clear
syms y(x)
y = x^2 - pow(cos(pi*x), 2)
X = 0.3:0.02:0.4
Y = double(subs(y, x, X)
XY = [X; Y]
y0 = 0
n = length(X)
syms Dy(y)
L = sym('L', [1, n])
k = 1:n;
for i = 1:n
    Dy = y -Y;
    L(i) = prod(Dy(k ~= i));
    dy = Y(i) - Y;
    L(i) = L(i) ./ prod(dy( k ~= i));
end
 
f = sum(L.*X);
f = simplify(f)
for i = 1:n
    double(subs(f, y, Y(i)) - X(i))
end
x0 = double(subs(f, y, y0))

@VictorVAlduin · 01.11.2022, 19:59 **[ТС]**

Я попробовал переписать код но столкнулся с некоторыми ошибками.
Нужна помощь!
Еще я знаю что нужно написать
for x in range(0.3,0.4, 0.02)
но я не очень понимаю куда этот цикл втыкнуть
Еще не совсем как перевести ~= на язык пайтона

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
y_x = symbols('y(x)')
y = x ^ 2 - pow(cos(pi * x), 2)
# X = 0.3:0.02:0.4
Y = float(y.subs(y, x, X))
XY = [X, Y]
y0 = 0
n = len(X)
D_y = symbols('Dy(y)')
L = symbols('L', [1, n])
for k in range(1,n):
    for i in range(1, n):
        Dy = y - Y
        L[i] = np.prod(Dy(k ~= i))
        dy = Y[i] - Y
        L[i] = L[i] / np.prod(dy( k ~= i))
 
 
f = sum(L*X);
f = simplify(f)
for i in range(1, n):
    print(float(y.subs(f, y, Y[i]) - X[i]))
 
x0 = float(y.subs(f, y, y0))

@u235 · 01.11.2022, 20:05

Сообщение от VictorVAlduin

Еще не совсем как перевести ~= на язык пайтона

в Python это

Python
1
dy[k != i]

При условии, что dy это не список, а np.array

@VictorVAlduin · 01.11.2022, 23:55 **[ТС]**

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
    symbols('y(x)')
    y = x**2 - np.power(cos(np.pi * x), 2)
    X = np.arange(0.3, 0.4, 0.02)
    Y = float(y.subs(y, x, X))
    XY = [X, Y]
    y0 = 0
    n = len(X)
    symbols('Dy(y)')
    L = symbols('L', [1, n])
    for k in range(1, n):
        for i in range(1, n):
            Dy = y - Y
            L[i] = np.prod(Dy[k != i])
            dy = Y[i] - Y
            L[i] = L[i] / np.prod(dy[k != i])
 
 
    f = sum(L*X);
    f = simplify(f)
    for i in range(1, n):
         print(float(y.subs(f, y, Y[i]) - X[i]))
 
    x0 = float(y.subs(f, y, y0))

Почему происходит следущая ошибка? Почитал воспользовался нужным методом np.power
y = x**2 - np.power(cos(np.pi * x), 2)
TypeError: unsupported operand type(s) for ** or pow(): 'list' and 'int'

Новые блоги и статьи Все статьи Все блоги /
Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек SDL3 и Box2D из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия SDL 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual. . .	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 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»

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

	Обратная интерполяция 30.10.2022, 14:08. Показов 1612. Ответов 5 Метки нет (Все метки) Нужна помощь! Не успеваю ребята пишу дипломную работу.Буду благодарен за помощь! Посредством обратной интерполяции определить с заданной точностью ε корень уравнение, размещенное на отрезке [a;b], с шагом h. Можно с помощью SYMPY У меня есть пример кода на матлаб eps = 1e-5 a = 0.3 b = 0.4 h = 0.02 Миниатюры 0

Обратная интерполяция

Решение