私はsymsum関数の助けを調べました、そしてあなたは本当に良い例を持っています、これを試してください:
syms x;
syms k real;
symsum(x^k/sym('k!'), k, 0, inf)
このコマンドは系列を評価し、実際にはに評価します。ご覧のとおり、「k」に依存するシリーズの用語を指定する必要があります。次に、symsumコマンドで、「k」を0からinfまで合計することを指定する必要があります。
したがって、たとえば、次のようにすることができます。
syms x;
syms k real;
ak = (-1)^k*x^(2*k+1)/sym('(2*k+1)!');
sum_ak = symsum(ak, k, 0, inf); % gives back sin(x)
dak = diff(ak,x);
sum_dak = symsum(dak, k, 0, inf); % should give back cos(x), but does not
A5 = symsum(ak, k, 0, 5); % add only the first values of the series
DA5 = symsum(dak, k, 0, 5); % add the derivated terms of the series
複数のシンボリック変数ukを宣言し、それらを合計することができます。
syms x;
syms k real;
n = 5;
for i = 0:n
eval(['syms u',num2str(i),' real;']);
end
A = cell(1,n);
for i=1:n
A{i} = u0;
for j=1:i
eval(['A{i} = A{i} + u',num2str(j),';']);
end
end
A{3} % check the value of A{i}
お役に立てれば、