フィルターのクリア

Analytical solution of ilaplace

2 ビュー (過去 30 日間)
john
john 2014 年 4 月 19 日
コメント済み: john 2014 年 4 月 19 日
Hello,
numerical solution of ilaplace is done by:
s=ilaplace(result(i,1));
s=feval(symengine, 'float', s);
syms t 'real'
s=vpa(sqrt(2)*real(s),4);
but how can I solve ilaplace analytical? for example:
-((2^(1/2)*C*L2*U*sin((pi*phiU)/180)*s^3)/2 + ((2^(1/2)*C*R3*U*sin((pi*phiU)/180))/2 + (2^(1/2)*C*L2*U*omega*cos((pi*phiU)/180))/2)*s^2 + ((2^(1/2)*U*sin((pi*phiU)/180))/2 + (2^(1/2)*C*R3*U*omega*cos((pi*phiU)/180))/2)*s + (2^(1/2)*U*omega*cos((pi*phiU)/180))/2)/((omega^2 + s^2)*(R1 + R3 + L1*s - C*M12^2*s^3 + C*L1*L2*s^3 + C*L2*R1*s^2 + C*L1*R3*s^2 + C*L2*R3*s^2 + 2*C*M12*R3*s^2 + C*R1*R3*s))

回答 (1 件)

Walter Roberson
Walter Roberson 2014 年 4 月 19 日
ilapace(expression)
is the analytic form. It might involve the symbolic sum over the roots of an expression, perhaps using a RootOf() object. If the polynomial in the RootOf() is degree higher than 4 then you are unlikely to be able to get an analytical solution. If it is degree 2, 3, or 4 then you can in theory extract the analytical solutions, but for anything higher than degree 2 they might be terribly long.
Unfortunately at the moment I do not know how to use MuPAD to request all of the roots of a RootOf() object which is a polynomial of degree 2, 3, o4 4. Possibly indexing the RootOf() might do it. In Maple I would use the allvalues() operation, which is not defined for MuPAD. I cannot test this as I do not have MuPAD.
The analytic transform for the expression you give, with the roots all expanded, runs to pages and pages. I have my system trying to simplify() it, but my expectation is that I will run out of memory.
  1 件のコメント
john
john 2014 年 4 月 19 日
sorry....what does it mean: runs to pages and pages ? I will run out of memory?
thank you

サインインしてコメントする。

カテゴリ

Help Center および File ExchangeCalculus についてさらに検索

タグ

Community Treasure Hunt

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

Start Hunting!

Translated by