象征性的分段表达式转换成一个函数如何处理

16个视图(30天)

显示旧的评论

Vishal Tripathi 2016年3月6日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/271748-how-to-convert-symbolic-piecewise-expression-into-a-function-handle

评论道: 沃尔特·罗伯森 2016年4月6日

我需要一个象征性的分段表达式转换为matlab函数处理。

我试着使用matlabFunction()但它抛出一个symengine错误的不平衡或意想不到的括号或支架。

任何帮助将不胜感激。

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

Vishal Tripathi 2016年3月21日

编辑:Vishal Tripathi 2016年3月26日

                               clc;信谊T;信谊E;信谊v;信谊z;
                              
                               假设(v > 0 & v < 6000);
                              
                               有趣= int (exp (e / T) (8.314 *), T, 300年,T);
                              
                               s = 18 * 10 ^ 3;%的意思
                              
                               fun1 = (1 / (v *√(2 *π)))* exp(-((学位)^ 2)/ (2 * v ^ 2));
                              
                               一个= 1.95 e + 13;%频率因子
                              
                               b = 10/60;%升温速率
                              
                               fun2 = int (exp(有趣(a / b) *) * fun1, E, s 3 * v, s + 3 * v);
                              
                               α1 = (z - fun2);%计算α
                              
                               T = (350 - 400);
                              
                               alphaexp = (0.01134 - 0.04317);
                              
                               alpha2 =(潜艇(α))
                              
                               minfunc = (alpha2-alphaexp) ^ 2
                              
                               error1 = (sum (minfunc))
                              
                               error2 = matlabFunction (error1)
                              
                               错误= @ (x) error2 (x (1), (2))
                              
                               [x, fval] = fminsearch(错误,2000 [1])

当调用matlabFunction,错误显示

错误使用symengine > makeFhandle(第109行)错误:不平衡或意想不到的括号或支架。

错误symengine(60)行

信谊误差/ matlabFunction(第125行)g = symengine (makeFhandle, varnames、身体);

错误i2 =(15)行错误matlabFunction (error1)

抱歉迟到的回复,任何帮助是急需的。

谢谢。

注:-编辑-编辑代码按照下面的讨论。

Vishal Tripathi 2016年3月26日

编辑:Vishal Tripathi 2016年3月26日

我认为集成产生了一个分段函数,因为E或v,而不是因为t .但是我想我可能已经找到了解决我的问题。(通过指定一个v的支架,所以没有形成分段函数)

是的,我想在这两个变量同时最小化。

可悲的是另一个问题,出现在使用fminsearch函数认为,尽管E作为一个未定义的函数/变量,而我已经注明在fun2集成变量E,限制包含v (E甚至不应该存在在最后的误差函数)

fminsearch有问题,还是我做一些错误的地方?

P。S -我已经编辑上面的代码,供您参考。

登录置评。

在回答这个问题。

答案(1)

沃尔特·罗伯森 2016年3月27日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/271748-how-to-convert-symbolic-piecewise-expression-into-a-function-handle answer_215175

如果你写出来与最终集成步骤没有完成,你扩大()条款,结合指数和简化,可以表明,最低必须发生在z是alphaexp的均值。

v = 0和E = 18000附近,积分非常陡峭。v = 0.002486920643673647643722704附近的关键是让你头发少积分是如此巨大,以至于无法计算有效,但头发更让你不可或缺的非常小——例如值的1/10 ^ (10 ^ 13)。我认为这意味着至少是你的上界,v < 6000但之间的区别,以及任何v > 0.002486920643673647643722704将需要计算数以十亿美元计的数字来确定不同。固有的不准确常数如8.314中没有证明这些计算。

4评论
显示3年长的评论隐藏3年长的评论

Vishal Tripathi 2016年3月27日

编辑:Vishal Tripathi 2016年4月6日

道歉的值。这些只是诱饵的主要代码中的错误(有点长)

我会然而感激如果你能解释扩大,合并和简化的步骤。我附加的主要代码,与实际值。

                                  clc;clearvars;清晰的所有;
                                 
                                  信谊T;信谊E;信谊v1;信谊v2;信谊v3;
                                 
                                  假设(v1 > 0 & < 50000);
                                 
                                  假设(v2 > 0 & v2 < 50000);
                                 
                                  假设(v3 > 0 & v3 < 60000);
                                 
                                  b = 10/60;
                                 
                                  有趣= int (exp (e / T) (8.314 *), T, 300年,T);
                                 
                                  s1 = 175.6 * 10 ^ 3;
                                 
                                  fun11 = (1 / (v1 *√(2 *π)))* exp (- ((E-s1) ^ 2) / (2 * v1 ^ 2));
                                 
                                  a1 = 10 ^ 14.52;
                                 
                                  fun12 = int (exp(有趣(a1 / b) *) * fun11, E, s1-3 * v1, s1 + 3 * v1);
                                 
                                  α1 = 1 - fun12;
                                 
                                  s2 = 185.4 * 10 ^ 3;
                                 
                                  fun21 = (1 / (v2 *√(2 *π)))* exp (- ((E-s2) ^ 2) / (2 * v2 ^ 2));
                                 
                                  a2 = 10 ^ 13.64;
                                 
                                  fun22 = int (exp(有趣(a2 / b) *) * fun21, E, s2-3 * v2, s2 + 3 * v2);
                                 
                                  alpha2 = 1 - fun22;
                                 
                                  s3 = 195.4 * 10 ^ 3;
                                 
                                  fun31 = (1 / (v3 *√(2 *π)))* exp (- ((E-s3) ^ 2) / (2 * v3 ^ 2));
                                 
                                  a3 = 10 ^ 13.98;
                                 
                                  fun32 = int (exp(有趣(a3 / b) *) * fun31, E, s3-3 * v3, s3 + 3 * v3);
                                 
                                  alpha3 = 1 - fun32;
                                 
                                  α=(α+ alpha2 + alpha3) / 3
                                 
                                  alphaexp = (0.01134 - 0.04317);% 0.06494 0.08783 0.17053 0.32533 0.49142 0.55575 0.59242 0.6367 0.678 0.71621 0.75124 0.78442 0.81727);
                                 
                                  T = (350 - 400);% T = [350:50:1050];
                                 
                                  minfunc =(潜艇(α)-alphaexp) ^ 2
                                 
                                  error1 =总和(minfunc)
                                 
                                  错误= matlabFunction (error1)
                                 
                                  [xfinal, fval] = fminsearch (@ (x)错误(x (1), (2) x (3)), (4300 3500 32000))
                                 
                                  方差= xfinal
                                 
                                  错误=√fval

沃尔特·罗伯森 2016年4月6日

请不要一个变量名称“错误”干扰使用关键的MATLAB函数误差()和为别人很难读的代码。

请出示的输出error1 matlabFunction调用的输出。

我有一个模糊的记忆,可能有一个错误在int()函数的转换处理。

登录置评。

在回答这个问题。

类别

数学和优化符号数学工具箱数学假设

找到更多的在假设在帮助中心和文件交换

社区寻宝

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

开始狩猎!

象征性的分段表达式转换成一个函数如何处理

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

答案(1)

4评论
显示3年长的评论隐藏3年长的评论

另请参阅

类别

标签

社区寻宝

象征性的分段表达式转换成一个函数如何处理

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

答案(1)

4评论 显示3年长的评论隐藏3年长的评论

另请参阅

类别

标签

社区寻宝

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

4评论
显示3年长的评论隐藏3年长的评论