区分总和与千里马

时间:2012-01-13 16:12:08

标签: symbolic-math maxima

我有以下总和:

sum((R[i]-(a*X[i]+b)*t + 1/2*(c*X[i]+d)^2*t)^2/((c*X[i]+d)^2*t), i, 1, N);

我希望区分wrt。一个:

diff(%, a);

但是Maxima(确切地说是wxMaxima)只是打印d / da。我可以吗 使它实际上区分总和(因为N是有限的是 应该分别区分总和中的每个元素吗?

如果我将N设置为某个常数,例如:

sum((R[i]-(a*X[i]+b)*t + 1/2*(c*X[i]+d)^2*t)^2/((c*X[i]+d)^2*t), i, 1, 100);

然后我得到100个元素的明确总和(大约需要2页),和 然后分化工作(但我又得到2页而不是小页 和)。我可以将此结果显示为总和吗?

1 个答案:

答案 0 :(得分:2)

您使用的是哪个版本的Maxima

这是我的千里马会话,你用方程式区分wrt.a而不是替换为N=100

~$ maxima 
Maxima 5.24.0 http://maxima.sourceforge.net
using Lisp SBCL 1.0.51
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) sum((R[i]-(a*X[i]+b)*t + 1/2*(c*X[i]+d)^2*t)^2/((c*X[i]+d)^2*t), i, 1, N);
                                   2
                         (c X  + d)  t
                   N         i                             2
                  ====  (------------- - (a X  + b) t + R )
                  \            2             i           i
                   >    ------------------------------------
                  /                           2                                                                                                
                  ====              (c X  + d)                                                                                                 
                  i = 1                 i                                                                                                      
(%o1)             ------------------------------------------                                                                                   
                                      t                                                                                                        
(%i2) diff(%, a);                                                                                                                              
                                       2                                                                                                       
                             (c X  + d)  t                                                                                                     
                    N            i
                   ====  X  (------------- - (a X  + b) t + R )
                   \      i        2             i           i
(%o2)          - 2  >    --------------------------------------
                   /                            2
                   ====               (c X  + d)
                   i = 1                  i
(%i3) %, N=100;
                                       2
                             (c X  + d)  t
                   100           i
                   ====  X  (------------- - (a X  + b) t + R )
                   \      i        2             i           i
(%o3)          - 2  >    --------------------------------------
                   /                            2
                   ====               (c X  + d)
                   i = 1                  i