2012-01-13 20 views
1

区分总和我有以下的总和:与千里马

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。 a:

diff(%, a); 

但是Maxima(wxMaxima是精确的)只是打印d/da。我是否可以将它作为 来区分其总和(因为N有限是 应该分别区分总和中的每个元素)?

如果我N设置一些常量,比如:100元的

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); 

然后我得到明确的总和(约需2页),并 然后分化作品(但我爬不起来,而不是2页一个小的 总和)。我可以将这个结果显示为一个总和吗?

回答

2

你使用哪个版本的Maxima

这里是我的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 
+1

我使用Debian的wxMaxima 11.08.0(Maxima版本:5.24.0)。 – Grzenio 2012-01-15 22:13:49

+0

它看起来像相同的Maxima版本。你有和我一样的输出吗?我希望我的输出帮助。 – 2012-01-15 22:29:23

+1

我实际上设法通过设置simplify_sum:true来设法实现差异化。我想知道其他默认选项对于新手来说是违反直觉的。 – Grzenio 2012-01-16 13:38:41