Я создал функцию Matlab, и я хотел бы преобразовать ее в Python для использования с моим веб-приложением.От Matlab до Python - Решите функцию
Я преобразовал (.m файл в .py-файл) почти все, используя OMPC. Однако я не могу заставить функцию solve()
работать (я использую библиотеку sympy).
Это Matlab линия:
SBC = solve(sqrt((xa-x)^(2)+(ya-y)^(2))-sqrt((xb-x)^(2)+(yb-y)^(2))-D12==0,sqrt((xa-x)^(2)+(ya-y)^(2))-sqrt((xc-x)^(2)+(yc-y)^(2))-D13==0,[x,y]);
И это Python строки, где x
и y
символы (с x = Symbol('x')
и y = Symbol('y')
):
sbc = solve(
sqrt((xa - x) ** (2) + (ya - y) ** (2))
- sqrt((xb - x) ** (2) + (yb - y) ** (2))
- D12 == 0,
sqrt((xa - x) ** (2) + (ya - y) ** (2))
- sqrt((xc - x) ** (2) + (yc - y) ** (2))
- D13 == 0,
[x, y]
)
С помощью этого кода Python, я получив False
вместо результата (отлично работает с кодом Matlab).
Я что-то упустил?
EDIT:
И с этим, я получаю []
:
# -*- coding: utf-8 -*-
from sympy import *
def alg(xa=None, ya=None, za=None, Ta=None, xb=None, yb=None, zb=None, Tb=None, xc=None, yc=None, zc=None, Tc=None, xd=None, yd=None, zd=None, Td=None, RSSIA=None, RSSIB=None, RSSIC=None, RSSID=None, txPower=None, n=None):
n = 2
c = 3 * 10 ** 8
TOA12 = Ta - Tb
TOA13 = Ta - Tc
TOA14 = Ta - Td
D12 = TOA12 * c
D13 = TOA13 * c
D14 = TOA14 * c
x, y = symbols('x y')
eqs = [sqrt((xa - x) ** (2) + (ya - y) ** (2)) - sqrt((xb - x) ** (2) + (yb - y) ** (2)) - D12,
sqrt((xa - x) ** (2) + (ya - y) ** (2)) - sqrt((xc - x) ** (2) + (yc - y) ** (2)) - D13]
print solve(eqs, [x, y])
alg(3,1,0,21.8898790015,4,6,0,21.8898790005,10,4,0,21.88987900009,9,0.5,0,21.889879000105,23.9,23.85,23.9,23.95,24,1)
О, я вижу. Спасибо за это. Я только что увидел, что это было написано в первых строках части Solvers в документе Sympy doc :(Однако у меня все еще есть проблема, так как результатом является '[]', а не '{x: ..., y: ...} '. Любая идея? – Anthony
Не могли бы вы отредактировать вопрос и добавить код, который дает вам' [] '?! Как написано, я не получил результата либо потому, что моя сессия Python убивается через пару минут, вы пытаетесь решать слишком сложные уравнения ...;) – Cleb
Готово! :) Моя сессия не была убита, но для «решения» требуется 5 секунд. – Anthony