在另一个分段函数分段函数作为条件

8视图(30天)

显示旧的评论

克里斯蒂娜其中Gisolo 2020年11月28日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/666998-piecewise-function-as-condition-in-another-piecewise-function

评论道: 克里斯蒂娜其中Gisolo2020年11月29日

答:接受沃尔特·罗伯森

大家好,

我用符号数学工具箱来定义一个matlabFunction用于仿真软件。万博1manbetx

我有问题在评估的代码如下:

                         x =符号(“x”,(N - 1));
                        
                         L = 100;
                        
                         z (1, - 1) = (2) - x (1);%δx (N)
                        
                         z(1, - 1) =分段(z (1) < 0 L + z (1) z (1) > = 0, z (1));
                        
                         V(1,1) =分段(z (1) < = 3 0 3 < z (1) < 5, z (1), 10);
                        
                         matlabFunction (V,“文件”,“函数”,“var”,{z});

原始代码更复杂,但这是一个简化版本。

我得到以下错误消息:

                         错误使用symengine
                        
                         不能评估布尔。
                        
                         错误信谊/ mupadmexnout(第1057行)
                        
                         = mupadmex (fcn, args {:});
                        
                         错误信谊/ matlabFunction >优化(第468行)
                        
                         [tvalues f, tnames] = mupadmexnout (“symobj:: optimizeWithIntermediates”f: {});
                        
                         错误信谊/ matlabFunction > writeMATLAB(第443行)
                        
                         [f tvalues tnames] =优化(f, optim);
                        
                         错误信谊/ matlabFunction(第183行)
                        
                         g = writeMATLAB (varnames乐趣、文件,输出、主体选择。优化选择。稀疏,opts.Comments);
                        
                         错误在主(第94行)
                        
                         matlabFunction (V,“文件”,“函数”,“var”,{z});

(第94行)是指的什么地方我叫matlabFunction。

我不确定,但我想改变些什么,我认为给的问题是,在第二分段函数使用的评估另一个分段函数的输出,因为在第二分段函数使用z(1)从第一个分段。

你认为它可能是问题吗?因为如果我交换订单不能两个分段函数。不幸的是我想评估z(1),然后修改它之前评估V (1)。

z(1)是一个象征性的变量,但也许并不喜欢它,因为从1日分段z(1)是一个分段函数,是不好的用它来评估另一个分段的条件。

我试图使用一个替代第二分段,实现饱和函数,但我不知道如果有另一种方式,也许更聪明。你怎么认为?

提前谢谢你

0评论
显示1年长的评论藏1年长的评论

登录置评。

在回答这个问题。

接受的答案

沃尔特·罗伯森 2020年11月28日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/666998-piecewise-function-as-condition-in-another-piecewise-function answer_558853

matlabFunction (V,“文件”,“函数”,“var”,{x});

为我工作。你不能使用z var matlabFunction因为你需要独立变量的名称列表。

警告:在最近的版本中,优化了matlabFunction写入文件中默认情况下坏了。之前是固定我建议关闭优化:

matlabFunction (V,“文件”,“函数”,“var”{x},“优化”,0);

注意:当你使用一个表达式,包括分段matlabFunction()生成的代码不矢量化的任何变量出现在分段()条件测试。

5个评论
显示4年长的评论隐藏4年长的评论

克里斯蒂娜其中Gisolo 2020年11月28日

非常感谢你,

你是对的,其实我原来的代码中有一个独立的变量输入matlabFunction,就像下面的:

                                  %符号变量宣言
                                 
                                  x =符号(“x”,(N - 1));
                                 
                                  z =符号(“z”,(N - 1));
                                 
                                  Vopt =符号(“v”,(N - 1));
                                 
                                  % %植物的定义
                                 
                                  b = 20;% (1 / s)
                                 
                                  L = 100;% [m]
                                 
                                  z (N - 1) = (1) - x (N);%δx (N)
                                 
                                  Vopt (N - 1) =分段(z (N, 1) < b, 0 b < = z (N, 1) < L 0 z (N, 1) > = L, 0);
                                 
                                  为我= 1:n - 1,
                                 
                                  z (1) = x (i + 1)(我);%δx(我)
                                 
                                  Vopt(我,1)=分段(z(我,1)< b, 0 b < = z(我,1)< L 0 z(我,1)> = L, 0);
                                 
                                  结束
                                 
                                  matlabFunction (Vopt“文件”,“OptimalVelocity”,“var”{x},“优化”,0);

现在,这个简单的代码工作!这是由于优化,我相信,所以非常感谢你的帮助!现在我要做我的更复杂的问题和希望与这个提示。

我不明白你是什么意思:

注意:当你使用一个表达式,包括分段matlabFunction(),那么生成的代码不是矢量化分段()中出现的任何变量条件测试。

沃尔特·罗伯森 2020年11月29日

假设你有类似

                                   信谊x y
                                  
                                   f = matlabFunction (x。^ 2 + y ^ 2,“var”,(x, y))

然后输出

@ (x, y) x ^ 2 + y ^ 2

请注意,您可以通过兼容数组的大小,所以你可以使用

f(兰特(3、5),兰德(3、5))

即使你使用

                                   信谊x y
                                  
                                   f = matlabFunction (x。^ 2 + y ^ 2,“var”{(x, y)})

你出去

@ (in) in (: 1)。^ 2 + in1 (:, 2) ^ 2

注意到你可以通过在多个行(有两列),它会将每一行作为一个x, y对:

2)f(兰德(15日)

所以主要是使用matlabFunction时,生成的函数处理得到矢量化(后模式建立的“var”参数)

然而,如果表达式包含一个int()不能被删除掉,或者如果表达式包含分段(),那么生成的函数处理不矢量化。例如,

                                   信谊x y
                                  
                                   f = matlabFunction (int (√sin (x) - x ^ 2), x,, y),“var”,{y})

回馈

@ (y)积分(@ (x) sqrt (sin (x) - x ^ 2), 0.0, y)

这只会被评估为标量参数y。

如果你使用int(),那么您至少可以使用matlabFunction返回一个匿名函数。当你使用分段(),那么您将得到一个错误,除非你告诉它写一个文件使用“文件”选项。生成的代码可能看起来像

                                   1函数V =函数(in)
                                  
                                   2%的功能
                                  
                                   3% V =函数(in)
                                  
                                   4
                                  
                                   5%这个函数生成的符号数学工具箱8.6版。
                                  
                                   6% 11月28 - - 2020 15:49:10
                                  
                                   7
                                  
                                   8 x1 = in (1:);
                                  
                                   9 x2 = in (2:);
                                  
                                   10 t2 = x1;
                                  
                                   11 t3 x2 = t2 +;
                                  
                                   12如果((t3 < = -9.7 e + 1) | ((0.0 < = t3) & (t3 < = 3.0)))
                                  
                                   13 V = 0.0;
                                  
                                   14elseif((t3 < -9.5 e + 1) & (-9.7 e + 1 < t3))
                                  
                                   15 V = t3 + 1.0 e + 2;
                                  
                                   16elseif((3.0 < t3) & (t3 < 5.0))
                                  
                                   17 V = t3;
                                  
                                   18其他的
                                  
                                   19 V = 1.0 e + 1;
                                  
                                   20.结束

注意在12和14和16行,代码假设t3是一个标量,代码不矢量化的输入。

克里斯蒂娜其中Gisolo 2020年11月29日

谢谢你,你很善良。

登录置评。

类别

数学和优化符号数学工具箱数学数论

找到更多的在数论在帮助中心和文件交换

社区寻宝

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

开始狩猎!

在另一个分段函数分段函数作为条件

0评论
显示1年长的评论藏1年长的评论

接受的答案

5个评论
显示4年长的评论隐藏4年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

在另一个分段函数分段函数作为条件

0评论 显示1年长的评论藏1年长的评论

接受的答案

5个评论 显示4年长的评论隐藏4年长的评论

更多的答案(0)

另请参阅

类别

标签

社区寻宝

0评论
显示1年长的评论藏1年长的评论

5个评论
显示4年长的评论隐藏4年长的评论