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を区別したい。あ:

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ページを取得します）。この結果を合計として表示することはできますか?

4

1 に答える 1

2

Maxima のどのバージョンを使用していますか?

これは、数式を wrt.a で微分し、 に置き換えた Maxima の私のセッションです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
于 2012-01-15T12:56:25.573 に答える