Question

我在MuPAD中有poly（一个变量，几个参数）。我想以某种方式将它提取到Matlab。

MuPAD中poly的一个例子：

poly((-7/(2*k^2))*X^2 + ((7*(b + k))/(2*k^2) + (7*(b + 2*k))/(2*k^2) + 4/k)*X + (- (b + k)*((7*(b + 2*k))/(2*k^2) + 4/k) + 1), [X])

我想在Matlab中得到它：

x*((7*(b + k))/(2*k^2) + (7*(b + 2*k))/(2*k^2) + 4/k) - (7*x^2)/(2*k^2) - (b + k)*((7*(b + 2*k))/(2*k^2) + 4/k) + 1

有没有方便的方法呢？附： Matlab R2015a

Answer 1

我相信你想将域DOM_POLY的MuPAD对象转换为一般的符号表达式。您可以使用expr功能。我假设你已经在Matlab内部工作了（来自MuPAD，它只是expr(p1)）：

syms b k X
p1 = feval(symengine,'poly',(-7/(2*k^2))*X^2 + ((7*(b + k))/(2*k^2) + (7*(b + 2*k))/(2*k^2) + 4/k)*X + (- (b + k)*((7*(b + 2*k))/(2*k^2) + 4/k) + 1),'[X]')
p2 = feval(symengine,'expr',p1)

返回

p2 =

X*((7*b + 7*k)/(2*k^2) + (7*b + 14*k)/(2*k^2) + 4/k) - (b + k)*((7*b + 14*k)/(2*k^2) + 4/k) - (7*X^2)/(2*k^2) + 1

您还可以使用matlabFunction将表单（上面的p1或p2）转换为矢量化浮点函数句柄，例如：

f = matlabFunction(p1,'Vars',{b k X})

返回

f = 

    @(b,k,X)-(b+k).*(1.0./k.^2.*(b.*7.0+k.*1.4e1).*(1.0./2.0)+4.0./k)+X.*(1.0./k.^2.*(b.*7.0+k.*7.0).*(1.0./2.0)+1.0./k.^2.*(b.*7.0+k.*1.4e1).*(1.0./2.0)+4.0./k)-X.^2.*1.0./k.^2.*(7.0./2.0)+1.0

Answer 2

您可以使用generate::MATLAB函数为任何表达式创建MATLAB代码。要打印函数调用（不带引号），可以使用AVPlayer(URL: URL_WHERE_TO_SAVE)。一个小例子：

print

返回

p := -7*X^2 + 8*X + 1
print(Unquoted, generate::MATLAB(p))

这是该函数的正确MATLAB语法。

从MuPAD中提取多项式到Matlab

2 个答案: