function maximize/minimize in Smath Studio

function maximize/minimize in Smath Studio - Help! - Messages

#1 Posted: 4/28/2021 4:47:54 PM
Denis

Denis

0 likes in 14 posts.

Group: User

Hi, guys
I need help for my graduation work, i need to use maximize function in Smath Studio, but i have no idea how to do it .
Is there a plugin for maximize / minimize functions in Smath Studio?
#2 Posted: 4/28/2021 10:58:15 PM
overlord

overlord

552 likes in 1332 posts.

Group: Moderator

Wrote

Hi, guys
I need help for my graduation work, i need to use maximize function in Smath Studio, but i have no idea how to do it .
Is there a plugin for maximize / minimize functions in Smath Studio?


You can write your own functions.
maximize and minimize are peak points of slope.
Take derivative of equation, solve it.
Check the examples below.

Regards

maximize_minimize.sm (16 KiB) downloaded 83 time(s).

maximize_minimize with different solvers.sm (28 KiB) downloaded 83 time(s).

maximize_minimize.png
1 users liked this post
Alvaro Diaz Falconi 4/29/2021 12:29:00 AM
#3 Posted: 4/29/2021 9:00:50 AM
Denis

Denis

0 likes in 14 posts.

Group: User

I have a program in MathCad,the problem is that the marshmallow recipe is being optimized here,

I need a maximization function to find the maximum value of each ingredient.

sbalansir. + organoleptika zefira - f-cii Ketle.xls (21 KiB) downloaded 26 time(s).

the exel file contains the values ​​of each marshmallow component that are used in the task.

sbalansir. + organoleptika zefira - f-cii Ketle.xmcd (250 KiB) downloaded 31 time(s).

The whole problem is that the function K(x) is six-dimensional and there is no idea how to find its derivatives.

I need to do the same in Smath Studio, but there are no Given / Maximize / Minimize functions in Smath,

I don't know how to implement this.

Here is everything that I could implement in Smath from the Mathcad file.

diplomnayaRabota.sm (33 KiB) downloaded 32 time(s).

p.s If something is not clear, ask, I will try to answer.
Sorry for my bad eng
#4 Posted: 4/29/2021 9:22:32 AM
Alvaro Diaz Falconi

Alvaro Diaz Falconi

992 likes in 1674 posts.

Group: User

Wrote



p.s If something is not clear, ask, I will try to answer.



Hi, can you upload a pdf instead the mathcad xmcd file, please?

Best regards.
Alvaro.
#5 Posted: 4/29/2021 9:56:41 AM
Denis

Denis

0 likes in 14 posts.

Group: User

Wrote

Wrote



p.s If something is not clear, ask, I will try to answer.



Hi, can you upload a pdf instead the mathcad xmcd file, please?

Best regards.
Alvaro.




everything is written in Russian, if you have any questions, I will explain.

Mathcad - sbalansir. + organoleptika zefira - f-cii Ketle.pdf (245 KiB) downloaded 46 time(s).
#6 Posted: 4/29/2021 10:18:43 AM
Denis

Denis

0 likes in 14 posts.

Group: User

Wrote

Wrote

Hi, guys
I need help for my graduation work, i need to use maximize function in Smath Studio, but i have no idea how to do it .
Is there a plugin for maximize / minimize functions in Smath Studio?


You can write your own functions.
maximize and minimize are peak points of slope.
Take derivative of equation, solve it.
Check the examples below.

Regards

maximize_minimize.sm (16 KiB) downloaded 83 time(s).

maximize_minimize with different solvers.sm (28 KiB) downloaded 83 time(s).

maximize_minimize.png



Unfortunately, this does not work with my function, because it is more complex than yours.
#7 Posted: 4/29/2021 11:44:59 AM
gallyon1

gallyon1

0 likes in 6 posts.

Group: User

Why are denominators 0? B2j≠B1j. Am i wrong?
#8 Posted: 4/29/2021 12:16:06 PM
overlord

overlord

552 likes in 1332 posts.

Group: Moderator

Wrote

Unfortunately, this does not work with my function, because it is more complex than yours.


Here is the first solution for your graphics. Rest are in the example.
This is an ugly and fast solution compiled in a haste.
I am pretty sure someone in the forum will come with a more classy method.

Regards

diplomnayaRabota_rev.sm (40 KiB) downloaded 36 time(s).

2021-04-29_18-12.png
1 users liked this post
Denis 4/29/2021 1:03:00 PM
#9 Posted: 4/29/2021 12:38:37 PM
Denis

Denis

0 likes in 14 posts.

Group: User

Wrote

Wrote

Wrote

Unfortunately, this does not work with my function, because it is more complex than yours.


Here is the first solution for your graphics. Rest are in the example.
This is an ugly and fast solution compiled in a haste.
I am pretty sure someone in the forum will come with a more classy method.

Regards

diplomnayaRabota_rev.sm (40 KiB) downloaded 23 time(s).

2021-04-29_18-12.png



Thanks for your feedback , but I need to calculate maximize from the function K(x),
not from y(j, x).

Check this pdf Mathcad file:
Mathcad - sbalansir. + organoleptika zefira - f-cii Ketle.pdf (245 KiB) downloaded 37 time(s).

After the function K (x), the vector x is given as the initial approximation
and after that the Given and Miximize functions are used for the K(x) function.
I need to get vector X in Smath.
p.s sorry for my bad eng



#10 Posted: 4/30/2021 6:31:24 AM
Denis

Denis

0 likes in 14 posts.

Group: User

Wrote

Graphs are doctored and more, but can't sanity more
Cheers ...Jean

diplomnayaRabota.sm (45 KiB) downloaded 31 time(s).



Thanks for your corrections , but the main question was how to get this vector X without the Given and Maximize functions,

because it's obtained from these functions: "X: = Maximize (K, x) ?

1.png
#11 Posted: 4/30/2021 7:24:17 AM
overlord

overlord

552 likes in 1332 posts.

Group: Moderator

Wrote

Thanks for your corrections , but the main question was how to get this vector X without the Given and Maximize functions,

because it's obtained from these functions: "X: = Maximize (K, x) ?


That is a total different solver. Smath don't have it integrated.
Maple software has it as Maximize() and Minimize() but plugin doesn't parse it under linux.
At least I couldn't do it with maple plugin.

Maxima can solve that kind of equations with maximize_lp() and minimize_lp().
You have to load(simplex) before using this two functions.
Maybe maxima plugin could work but I can't try it, need a windows user to try it.
I use maxima under terminal, maxima plugin doesn't work with linux.

Regards.
#12 Posted: 4/30/2021 9:03:36 AM
Denis

Denis

0 likes in 14 posts.

Group: User

Wrote

Wrote

but the main question was how to get this vector X without the Given and Maximize functions,
because it's obtained from these functions: "X: = Maximize (K, x) ?


The bloc solver Given ... Find, Minerr, Maximize, Minimize
is Mathsoft proprietary specific to Mathcad. No such equivalent Smath.
Given/Minerr for fitting model function to data set is done otherwise
in Smath via Conjugate Gradient ... same or even better.
If there is another way of doing in Smath, expect Designers to help.
Sorry Collab, can't help directly on that ... Jean.





I tried to do this for a simple function through the "roots" function, and in the end I found the maximum,

finding the maximum for a simple function.sm (12 KiB) downloaded 59 time(s).

but there are no ideas how to do this for K (X), since it depends on 6 arguments of the vector x.

Is it possible to somehow implement the same algorithm in my task?
#13 Posted: 5/1/2021 7:42:10 AM
Denis

Denis

0 likes in 14 posts.

Group: User

Wrote

Wrote

Is it possible to somehow implement the same algorithm in my task ?


The Ancestor of the savant solve is the TKsolve.
The designers of TKsolve are the Founders of Mathsoft/Mathcad.
No Mathcad gurus succeeded to crack the bloc, probably crypted deep.



It's a pity that this method is not solved.
#14 Posted: 5/1/2021 4:57:32 PM
Alvaro Diaz Falconi

Alvaro Diaz Falconi

992 likes in 1674 posts.

Group: User

Hi. There are no minimizers in SMath, with the probable exception of a library in Maxima that can be read and executed, as Overlod point.

However, a kind of minimizer can be implemented as shown in the attached, and incidentally, obtain a solution (just a little bit) better than Mathcad's.

diplomnayaRabota_nleqsol.sm (83 KiB) downloaded 85 time(s).

Best regards.
Alvaro.
1 users liked this post
Denis 5/1/2021 8:09:00 PM
#15 Posted: 5/3/2021 8:00:05 AM
Denis

Denis

0 likes in 14 posts.

Group: User

Wrote

Hi. There are no minimizers in SMath, with the probable exception of a library in Maxima that can be read and executed, as Overlod point.

However, a kind of minimizer can be implemented as shown in the attached, and incidentally, obtain a solution (just a little bit) better than Mathcad's.

diplomnayaRabota_nleqsol.sm (83 KiB) downloaded 85 time(s).

Best regards.
Alvaro.



Thank you so much for your work. I appreciate it.
#16 Posted: 5/4/2021 8:14:35 AM
Denis

Denis

0 likes in 14 posts.

Group: User

Wrote

Hi. There are no minimizers in SMath, with the probable exception of a library in Maxima that can be read and executed, as Overlod point.

However, a kind of minimizer can be implemented as shown in the attached, and incidentally, obtain a solution (just a little bit) better than Mathcad's.

diplomnayaRabota_nleqsol.sm (83 KiB) downloaded 85 time(s).

Best regards.
Alvaro.



Hi, I didn't really understand how you made such a function for the previous maximize problem,
I wanted to ask you to help me solve a simpler problem on this topic.

Milk shake optimization.sm (17 KiB) downloaded 44 time(s).

pdf Mathcad file:

Mathcad - mag.2020 Aleksandr Baranov Optimizacija molochnogo koktejjlja.pdf (99 KiB) downloaded 61 time(s).

In this program, I need to repeat everything in the same way as in the screenshots of MathCad, i attached them to Smath.

I would be very grateful if you could help me.

If you have any questions, ask, I will answer them.
#17 Posted: 5/4/2021 11:41:17 AM
Denis

Denis

0 likes in 14 posts.

Group: User

Wrote

Wrote

In this program, I need to repeat everything in the same way as in the screenshots of MathCad, I attached them to Smath


Horrible screen shot, nothing to see for understanding.
Seemingly, you have a system to solve ... something like Butler.

Solve Given_Find [CO2_6179 Butler].sm (162 KiB) downloaded 37 time(s).



I'll try to make a pdf Mathcad file.
Some formulas simply dont fit into the full sheet and go beyond the boundaries.
#18 Posted: 5/4/2021 1:33:46 PM
Denis

Denis

0 likes in 14 posts.

Group: User

Wrote

Wrote

Some formulas simply don't fit into the full sheet and go beyond the boundaries.


Your Mathcad document is surely badly dressed.



like it turned out to make pdf
Mathcad - mag.2020 Aleksandr Baranov Optimizacija molochnogo koktejjlja.pdf (99 KiB) downloaded 61 time(s).
#19 Posted: 5/4/2021 5:08:24 PM
Alvaro Diaz Falconi

Alvaro Diaz Falconi

992 likes in 1674 posts.

Group: User

Wrote


Hi, I didn't really understand how you made such a function for the previous maximize problem,



Hi. Point is: how to solve that in SMath? Because SMath haven't maximizer, and what you have near to that is al_nleqsolve, you can read that in this way: How to use a traditional solver for seek a maximum with constraints?

So, you have a function K(x) with x element of R^6, a guess xo and the guess value for the maximum K(xo). First, write Check(x), for restrict x for the constraints. Next, because numerical solvers for roots needs the evaluation of the change in the target value, assume that you are looking for the zero of f(x) = K(x) - (1+p)*Ko, where p is some percent increment over Ko. Why? Because you guess that xo it's a good solution, so, Ko = K(xo) "must" to be a good guess for the maximum. The Maximize ad hoc function it's just the implementation for that: here the trick is use K(Check(x)) instead K(x). You can try to change the p-value for get the best answer that you can.

Wrote



In this program, I need to repeat everything in the same way as in the screenshots of MathCad, i attached them to Smath.

I would be very grateful if you could help me.

If you have any questions, ask, I will answer them.



The Mathcad's screenshots are perfect. Here you have a little difference: the restrictions involving another function Kc(x), but it's very similar, and need to write CheckX for x and CheckKc for Kc(x). The trick for CheckKc(x): No one, but three:

First, return some information to the solver about how the things are going. Instead of return 0 or 1, return 0 for Kc is in the good interval, and +/- normi(x) if not.
Second, include CheckKc(CheckX(x)) as second row in the target function for try to get zeros on both components.
Third, include CheckKc in the exit condition for the main loop.

Milk shake optimization nleqsol.sm (1 MiB) downloaded 60 time(s).


Hope that's helps.
Best regards.

Alvaro.
1 users liked this post
Denis 5/4/2021 6:05:00 PM
#20 Posted: 5/4/2021 7:21:08 PM
Davide Carpi

Davide Carpi

1416 likes in 2873 posts.

Group: Moderator

You might try even these

GoldenSectionSearch.sm (21 KiB) downloaded 232 time(s).
GradientDescent.sm (16 KiB) downloaded 191 time(s).
NCGM.sm (19 KiB) downloaded 172 time(s).
LevenbergMarquardt.sm (45 KiB) downloaded 205 time(s).
NewtonMethod.sm (18 KiB) downloaded 236 time(s).
NelderMead.sm (33 KiB) downloaded 187 time(s).
GaussNewton.sm (20 KiB) downloaded 201 time(s).
If you like my plugins please consider to support the program buying a license; for personal contributions to me: paypal.me/dcprojects
  • New Posts New Posts
  • No New Posts No New Posts