# Solving a sum of series of exponential function with a sum of series of cosine function inside

조회 수: 127(최근 30일)
Cheung Ka Ho 2 Jul 2017
댓글: Cheung Ka Ho 10 Jul 2017
Here is the equation I'm going to solve.
where epsilon is a variable.
I wrote the following code.
r = 1:(2^h)-1;
v = 1:h-1
T = exp(x*h/1000)+sum(exp((x/1000)*(cos(2*pi*r/(2^h))+cos(2*pi*r*(2^v)/(2^h)))));
I can run the code when h=2; however, when h becomes greater than or equal to 3, an error that input must be a scalar and a square matrix shows up.

#### 댓글 수: 3

Cheung Ka Ho 4 Jul 2017
After I changed the code into 2.^v, another main problem comes out, which is inner matrix dimensions must agree. I just want to calculate the summation. How can I rewrite the code to ignore the dimensions of r and 2.^v?
Torsten 5 Jul 2017
Is "epsilon" in your equation a given function ?
Is "e" in your equation the usual Euler-number ?
A closing round parenthesis is missing in your picture.
What does e{...} mean ? Does it mean exp(...) ?
Best wishes
Torsten.
Cheung Ka Ho 5 Jul 2017
epsilon is only a variable while e is the usual Euler-number. I mistakenly type it as e{...}. It should be e^{...}. Thank you for your help.

로그인 to comment.

### 채택된 답변

Torsten 5 Jul 2017
편집: Torsten 5 Jul 2017
h=...;
epsilon=...;
r=1:2^h-1;
v=1:h-1;
result=sum(exp(epsilon*sum(cos(pi/2^(h-1)*(2.^v).'*r))));
Best wishes
Torsten.

#### 댓글 수: 1

Cheung Ka Ho 10 Jul 2017
Thank you!

로그인 to comment.

### 추가 답변(1개)

Matthew Taliaferro 2 Jul 2017
편집: Matthew Taliaferro 3 Jul 2017
You cannot raise things to a power unless they are scalar or square (like the warning said). If you want to square each element, the notation is a little different.
h = 1:10
h_square = h.^2 % as opposed to h^2, which won't work
You also cannot divide something element by element unless it is a scalar.
r = 1:10; h = 1:10;
r_over_h = r./(h.^2); % as opposed to r/(h^2), which won't do what you think it does

#### 댓글 수: 1

Cheung Ka Ho 4 Jul 2017
Thank you for your help. It solves the above error. However, another problem that inner matrix dimensions must agree comes out. How can I resolve this error?

로그인 to comment.

이 질문에 답변하려면 로그인을(를) 수행하십시오.

Translated by