0

我正在研究最大值的代码,并在其中定义了一个函数。

Ucr(y):=offset+c1[21]+2*sum(realpart(c1[m+20])*cos(2*%pi*(m-1)*y/dp) - imagpart(c1[m+20])*sin(2*%pi*(m-1)*y/dp),m,2,nmax+1);

但是,如果我尝试评估诸如 Ucr(1) 之类的东西,我没有得到答案的数字,我会得到一个很长的方程式,其中包含所有变量和数字。我是 maxima 的新手,不知道我怎么能让它煮成一个漂亮的表情。

我将不胜感激任何帮助。

最好的,本

4

2 回答 2

0

尝试环绕float(trigreduce(rectform( ... )))你的函数体。

于 2015-01-09T18:14:28.737 回答
0

您可以给变量数值

kill(nmax, dp, offset, c1);
numerval(nmax, 1);
numerval(dp, 1);
numerval(offset, 1);
/* fill a list with dummy values */
c1: makelist(1/m+m*%i, m, 1, ev(nmax, numer)+22);
Ucr(y):=offset+c1[21]+2*sum(realpart(c1[m+20])*cos(2*%pi*(m-1)*y/dp) - imagpart(c1[m+20])*sin(2*%pi*(m-1)*y/dp),m,2,nmax+1);

ev(Ucr(1), numer);

或者只是在最后结合数值

kill(nmax, dp, offset, c1);
Ucr(y):=offset+c1[21]+2*sum(realpart(c1[m+20])*cos(2*%pi*(m-1)*y/dp) - imagpart(c1[m+20])*sin(2*%pi*(m-1)*y/dp),m,2,nmax+1);
env: [dp=1, offset=1, nmax=1, c1=makelist(1/m+m*%i, m, 1, 23)];
ev(Ucr(1), env);
于 2013-03-29T09:21:41.940 回答