如何简化函数句柄?

18视图(30天)

显示旧的评论

阿maheshwari 2021年8月19日

0
链接

这个问题直接联系

https://it.mathworks.com/matlabcentral/answers/1437209-how-to-simplify-function-handles

评论道: 沃尔特·罗伯森 2021年8月20日

答:接受沃尔特·罗伯森

Y =

function_handle与价值:

@ (s) (F (s) - q (s)) / P (s)

看,我有三个函数F, Q和P Y是对这三个函数由代数操作。有没有一种方法,得到简化的表达式Y的函数变量(s) ?当我在命令窗口中键入Y,我把结果如上所示。

作为一个简化的例子:

                          F = @ (s);
                         
                          Q = @ (s) s ^ 2;
                         
                          P = @ (s) ^ 3 - 3 *年代^ 2 - 1;
                         
                          Y = @ (s) (F (s) - Q (s)) / P (s);

当我在命令窗口中键入Y,看到的是:

Y =

function_handle与价值:

@ (s) (F (s) - q (s)) / P (s)

我希望看到的是:

Y = (s ^ 2) /(^ 3 - 3 *年代^ 2 - 1)

接受的答案

沃尔特·罗伯森 2021年8月19日

0
链接

直接链接到这个答案

https://it.mathworks.com/matlabcentral/answers/1437209-how-to-simplify-function-handles answer_770914

你用一个标签象征性的这意味着你可能会象征性的工具箱。如果是这样那么

                              信谊年代
                             
                              Y = Y (s)

您可能想要简化()结果。

7评论
显示6年长的评论隐藏6年长的评论

沃尔特·罗伯森 2021年8月19日

F = @ (s);

Q = @ (s) s。^ 2;

P = @ ^ 3 - 3 * s (s)。^ 2 - 1;

Y = @ (s) (F (s) - Q (s)) / P (s);

信谊年代真正的

Y = Y (s)

y =

dY = diff (y,年代)

dY =

关键=解决(dY)

关键=

d2Y = diff (dY)

d2Y =

crit2 =潜艇(d2Y,年代,关键)

crit2 =

definitely_max =总(crit2 < 0,“未知”,“假”);

definitely_min =总(crit2 > 0,“未知”,“真正的”);

not_sure = ~ (definitely_max | definitely_min);

max_locations =关键(definitely_max)

max_locations =

max_values = Y (max_locations)

max_values =

vpa (max_values)

ans =

0.42303605210393856972993261326406

min_locations =关键(definitely_min)

min_locations =

min_values = Y (min_locations)

min_values =

vpa (min_values)

ans =

other_locations =关键(not_sure)

other_locations =空信谊:0-by-1

other_values = Y (other_locations)

other_values =空信谊:0-by-1

沃尔特·罗伯森 2021年8月20日

F = @ (s);

Q = @ (s) s。^ 2;

P = @ ^ 3 - 3 * s (s)。^ 2 - 1;

Y = @ (s) (F (s) - Q (s)) / P (s);

信谊年代真正的

Y = Y (s)

y =

iY = ilaplace (y)

iY =

iY matlabFunction ()

错误使用symengine
代码生成失败由于意想不到的“RootOf”类型的对象。

信谊误差/ matlabFunction > mup2mat(第432行)
res1 = mupadmex (symobj:: generateMATLAB, r.s,另,spa);

信谊误差/ matlabFunction > mup2matcell(第401行)
r = mup2mat (c{1},真的,sparseMat);

信谊误差/ matlabFunction(第172行)
身体= mup2matcell(乐趣,opts.Sparse);

即ilaplace()生成一个根()(真的RootOf)和matlabFunction()无法处理这些。

理论上,这是一个三次多项式,ilaplace可以产生准确的公式而不是使用一个根()或RootOf()占位符……但它没有,没有Mathworks所提供的方法可以安排“表达式。

几周前,我发布了一个企图做这类结构的转换。我将试着挖掘的代码,看看我能不能让它为你工作。

沃尔特·罗伯森 2021年8月20日

下面的公式是输出注意:它的符号形式ilaplace依斯攀根从对称到封闭公式。

后策划调查,试图找出为什么fplot()显示不连续。

根据X / YY值,matlabFunction公式()版本是生产价值与虚构的组件大约e-17那里。放大试图寻找可见无穷大,没有帮助。但采取的一些地方产生非零虚部和评估他们象征性的50位显示单纯规则的输出和很小的虚部;虚部的振幅增加数字的数量减少。

这里正在发生的事情是,matlabFunction()版本是遭受舍入误差(甚至是象征性的版本)。根()的解决ilaplace一对共轭,所以理论上的积极的和消极的虚构的组件应该取消,但在实践中由于浮点圆滑,他们不。

在这个特定的情况下,至少在+ t,你应该能够把真正的()的值。

格式长g

数字(50)

F = @ (s);

Q = @ (s) s。^ 2;

P = @ ^ 3 - 3 * s (s)。^ 2 - 1;

Y = @ (s) (F (s) - Q (s)) / P (s);

信谊年代真正的

Y = Y (s)

y =

iY = ilaplace (y)

iY =

N = @ (expr N) expr (N);

resolve_a_root = @ (X6) n(解决((X6, 1),孩子(X6, 2),“maxdegree”4)、儿童(X6, 3));

resolve_a_summand = @ (X5) mapSymType (X5,“根”,resolve_a_root);

resolve_nth_summand = @ (X2, N) resolve_a_summand(潜艇((X2, 1),孩子(X2, 2), N));

map_symsum = @ (X2)和(arrayfun (resolve_nth_summand, repmat (X2, 1双(儿童(X2, 4)孩子们(X2, 3) + 1)),孩子(X2, 3):儿童(X2, 4)));

expand_symsum = @ (X1) mapSymType (X1,“symsum”,map_symsum);

公式=简化(expand_symsum (iY))

公式=

y = matlabFunction(公式);

fplot (y)[10 1的军医,军医]);

fplot (y,[4.257475 4.25749的军医的军医]);-1.001 ylim ([1]);标题(“变焦fplot数字”)

fplot(公式[4.257475 4.25749的军医的军医]);-1.001 ylim ([1]);标题(“变焦fplot符号”)

X = linspace(4.25749 4.257475的军医,军医,1000);

YY = y (X);

次要情节(2,1,1);情节(X,真正的(YY));标题(“固定数值评价真实”)

次要情节(2,1,2);情节(X,图像放大(YY));标题(“固定数值评价图像放大”)

X5idx =找到(图像放大(YY), 5);

X5 = X (X5idx)

X5 = 1×5

0.000425747501501501 0.000425747504504504 0.000425747513513513 0.000425747516516516 0.00042574751951952

日元= YY (X5idx)

日元=

-1.00085203903026 + 1.38777878078145 e-17i -1.00085203903628 - 1.38777878078145 e-17i -1.00085203905432 - 1.38777878078145 e-17i e-17i -1.00085203906033 + 1.38777878078145 e-17i -1.00085203906635 - 1.38777878078145

信谊t

潜艇(公式,t, X5 (:))

ans =

vpa (ans)

ans =

登录置评。

类别

数学和优化符号数学工具箱数学微积分

找到更多的在微积分在帮助中心和文件交换

下载188bet金宝搏

符号数学工具箱

释放

R2020b

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

翻译的