matlabFunction on symbolic experssion and the Error using symengine

8 次查看（过去 30 天）
Amin NZH
Amin NZH 2020-3-29
评论： Amin NZH 2020-3-30
Hello
I have a mix double integral and partial differential expression and I want to solve it by symbolic math in matlab.
My code below written to solve the equation. But matlabFunction() does not work on it and gives me an error.
I would be appreciated if anybody could help me. I tried the numerical way to solve this expression but I could not find the right result. The only way that I think work is using symbolic math. Any idea and help is appreciated.
Here is my code:
syms t z
w=1.11;
T=1.2;
zita=0.5;
a=1;
b=1;
d=0;
%%
teta=exp(-((1-t/T)*(t/T))^(-w));
phi=int(teta)/int(teta,0,T);
Y1=[1;1]+[1;2]*phi;
Y2=[2;1]+[1;1]*phi;
eta0=Y1+(z-zita)*Y2;
F=[(1/a)*(diff(eta0(1,1),t)-eta0(2,1)*eta0(1,1)); ...
(1/b)*(diff(eta0(2,1),t)-d*eta0(2,1))];
int1=int(F,z)-int(F,z,0,zita);
etaAux=int(int1,z)-int(int1,z,0,zita);
eta=eta0+etaAux;
etaD=diff(eta,z);
fun=matlabFunction(etaD);
After running the code, I got this error:
Error using symengine
Unable to evaluate to Boolean.
Error in sym/matlabFunction>mup2mat (line 404)
res = mupadmex('symobj::generateMATLAB',r.s,ano,spa,0);
Error in sym/matlabFunction>mup2matcell (line 374)
r = mup2mat(c{1},true,sparseMat);
Error in sym/matlabFunction (line 188)
body = mup2matcell(funs, opts.Sparse);
Could you please help me with this issue?
etaD is 2by1 and matlabFunction() works on etaD(2,1) but does not work on etaD(1,1) and gives me the same error as above.
  1 个评论
Walter Roberson
Walter Roberson 2020-3-30
It has something to do with the square of an integral. It doesn't matter if it is numerator or denominator but it does matter that it is squared

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Walter Roberson
Walter Roberson 2020-3-30
It appears to be a bug. I might look into the internal code when I am more awake, but that code is difficult to read and debug.
In the meantime, there is a workaround. It takes a bit of time to compute it, but it gets there.
The workaround is instead of
matlabFunction(etaD)
use
matlabFunction(rewrite(etaD, 'sincos'))
You will have to wait a bit for it to compute the rewrite()
Do not simplify() the results of the rewrite() or else you get back to code that has the problem
  4 个评论
显示 2更早的评论
Walter Roberson
Walter Roberson 2020-3-30
编辑：Walter Roberson 2020-3-30
Yes, you can matlabFunction(vpa(etaD))
I will file a bug report with Mathworks.
Amin NZH
Amin NZH 2020-3-30
IT WORKS!
Thank you so much for your help and time.
I appreciate your kind guidance.
Best regards

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Code Performance 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by