2012-01-13 42 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