Перевести код с матлаб на пайтон (обратная интерполяция)

01.11.2022, 20:03. Показов 417. Ответов 0

Метки нет (Все метки)

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

Нужна помощь в переводе кода с матлаб на пайтон

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 - (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))

Вот скажите как записать

Matlab M
1
XY = [X; Y]

это на пайтоне?
Не понимаю вообще что тут значит это?Вектора?
Также само эта строчка

Matlab M
1
L = sym('L', [1, n])

И что значит ~=
В остальном все ясно

Добавлено через 38 минут
Я попробовал переписать код но столкнулся с некоторыми ошибками.
Нужна помощь!
Еще я знаю что нужно написать
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))

Programming

Эксперт

39485 / 9562 / 3019

Регистрация: 12.04.2006

Сообщений: 41,671

Блог

01.11.2022, 20:03

Ответы с готовыми решениями:

Перевести код из пайтон в с++
Требуется перевести код из Python в С++ N = int(input()) a, b = 0, 0 def isint(s): try: int(s) ...

Нужно перевести код с С++ на Пайтон
//Нужна помощь именно с переводом функции, что стоит перед main E_K(...) #include <stdio.h> #include <iostream> #include...

Интерполяция заданной функции второй формулой Ньютона (обратная интерполяция)
Программа накапливает ошибку, да и первые значения совсем уже "приближённые". Явно что-то упустил. Посмотрите со стороны. #include...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.

inter-admin

Эксперт

29715 / 6470 / 2152

Регистрация: 06.03.2009

Сообщений: 28,500

Блог

01.11.2022, 20:03

Помогаю со студенческими работами здесь

Перевести код C++ на Матлаб
#include <iostream> using namespace std; int main() { int n, m; cin >> n >> m; int x; for (int i = 0; i <=...

Перевести код с языка С на матлаб
#include <stdio.h> #define INT_MIN -2147483647 #define INT_MAX 2147483647 int main() { int a = {{ 1,2,3,4,5 }, { 2, 3, 4,...

Перевести код с питона на матлаб
Кто может помочь переписать код с питона на матлаб , буду очень благодарен n = int(input()) x = list(map(float,input().split())) y =...

Помогите, пожалуйста, перевести код с языка С в матлаб)
#include<math.h> #include<stdio.h> #define N 10 int main(){ float x = { 0.0, 0.2, 0.4, 0.6, 0.8, 1.2, 1.4, 1.6, 1.8, ...

Обратная интерполяция
Нужна помощь! Не успеваю ребята пишу дипломную работу.Буду благодарен за помощь! Посредством обратной интерполяции определить с заданной...

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:

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

Наверх