sin関数での級数の和について

sin関数での級数の和についてです。

sin関数でk=1の時、k=2の時、... 、k=N-1の時の級数の和をそれぞれ出したいのですが

処理が終わらず出来ません。

この場合どのようにsin関数で級数の和を別々で出すにはどのようなプログラムを

打てばいいでしょうか。

以下がプログラムです。

qpskMod = comm.QPSKModulator;

qpskDemod = comm.QPSKDemodulator;

constDiagram = comm.ConstellationDiagram;

ofdmMod = comm.OFDMModulator;

ofdmDemod = comm.OFDMDemodulator;

modDim = info(ofdmMod);

nFrames = 10000;

e=exp(1);

N=128;

myu=0;

nakasiki_1=(0:N-1)';

nakasiki_2=(0:N-1)';

e_2=(0:N-1);

nSymbolsPerFrame = modDim.DataInputSize(1);

data = randi([0 3],N,1);

X_l = qpskMod(data);

p=pi;

e_1=e^((1i*p*myu*(N-1))/N);

for k=1:N

syms m

nakasiki_1(k)=symsum(sin(p*(m-k+myu)),m,1,N-1);

nakasiki_2(k)=symsum(sin(p*(m-k+myu)/N),m,1,N-1);

e_2(k)=symsum(e.^(1i*p*(m-k)*(N-1)/N),m,1,N-1);

end

[Naoya]

こちら、

nakasiki_1

の計算部分について、もう少し問題を簡略化してみますと以下のように表現できると思います。

syms m
k = 1;
myu = 0;
N = 128;
p = pi;
nakasiki_1=symsum(sin(p*(m-k+myu)),m,1,N-1)

これは、

sin((m-1)*pi) の式を m を 1 ～ 127 まで代入してそれぞれを加算することになります

m = 1 の時、 sin(0*pi) = 0

m = 2 の時、 sin(1*pi) = 0

m = 3 の時、 sin(2*pi) = 0

...

m = 127 の時、 sin(126*pi) = 0

sin関数は、 pi の周期で その解が 0 となりますので、その総和も 0の結果として得られているかもしれません。