用值替换变量,但不进行评估

时间:2017-04-05 03:09:39

标签: maxima

假设我有以下表达式:

(%i1) (8*x)*(log(x) / log(10));

(%i2) X^2;

现在,因为我想知道我可以选择什么常量值来使语句%i1 is O(%i2)为真,所以我在循环中评估它们是这样的:

for a:1 thru 10 do print(%i1, "=", ev(%i1, x=a), %i2, "=", ev(%i2, numer, x=a));

输出结果为:

8 x log(x)          2
---------- = 0.0 , x  = 1 
 log(10)

8 x log(x)                        2
---------- = 4.816479930623698 , x  = 4 
 log(10)

8 x log(x)                        2
---------- = 11.45091011327189 , x  = 9 
 log(10)

8 x log(x)                        2
---------- = 19.26591972249479 , x  = 16 
 log(10)

8 x log(x)                        2
---------- = 27.95880017344075 , x  = 25 
 log(10)

8 x log(x)                        2
---------- = 37.35126001841489 , x  = 36 
 log(10)

8 x log(x)                        2
---------- = 47.32549024079837 , x  = 49 
 log(10)

8 x log(x)                        2
---------- = 57.79775916748438 , x  = 64 
 log(10)

8 x log(x)                        2
---------- = 68.70546067963139 , x  = 81 
 log(10)

8 x log(x)           2
---------- = 80.0 , x  = 100 
 log(10)

我想让输出更容易引人注目,例如:

8 1 log(1)          2
---------- = 0.0 , 1  = 1 
 log(10)

8 2 log(2)                        2
---------- = 4.816479930623698 , 2  = 4 
 log(10)

8 3 log(3)                        2
---------- = 11.45091011327189 , 3  = 9 
 log(10)

[snip]

8 10 log(10)           2
---------- = 80.0 , 10  = 100 
 log(10)

如何在不评估表达式的情况下告诉Maxima在循环的每次迭代中用a替换x的值?

我搜索了手册,但我没有找到任何看似相关的内容。

2 个答案:

答案 0 :(得分:1)

使用"空"功能:

(%i1) display2d: false $
(%i2) prefix("") $
(%i3) almost_subst(a, x, e):= subst(""(a), x, e) $
(%i4) almost_subst(10, x, 8*x*log(x)/log(10));
(%o4) (8* 10*log( 10))/log(10)

答案 1 :(得分:1)

Maxima中的许多操作都是通过名为"简化"的过程来执行的,这意味着应用身份来简化"更简单的"表达。例如。 1 + 1简化为2,sin(0)简化为0等

为了获得你想要的效果,我们必须通常禁用简化,以便评估表达式但不简化。但是为了获得数值,我们需要为这些结果实现简化。

这是要做的事情。

(%i16) simp : false $
(%i17) for x in [1,2,3,4,5]
         do print (ev(%i1) = ev(%i1, simp, numer), ev(%i2) = ev(%i2, simp));
       log(1)          2
(8 1) (-------) = 0.0 1  = 1
       log(10)
       log(2)                        2
(8 2) (-------) = 4.816479930623698 2  = 4
       log(10)
       log(3)                       2
(8 3) (-------) = 11.4509101132719 3  = 9
       log(10)
       log(4)                        2
(8 4) (-------) = 19.26591972249479 4  = 16
       log(10)
       log(5)                        2
(8 5) (-------) = 27.95880017344075 5  = 25
       log(10)
(%o17)                               done

请注意,我写了for x in [1, 2, 3, 4, 5] ...而不是for x:1 thru 5 ...。这是因为后者使用算术,这需要简化。尝试两种方式,我认为你会看到差异,我相信这是非常有启发性的。

Nota bene我使用了%i1%i2的相同值。