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

10 次查看（过去 30 天）

显示更早的评论

Cheung Ka Ho 2017-7-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/347116-solving-a-sum-of-series-of-exponential-function-with-a-sum-of-series-of-cosine-function-inside

回答： mohammed alzubaidy 2021-5-16

采纳的回答： Torsten

在 MATLAB Online 中打开

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.

Could anyone please help me on this error?

3 个评论
显示 1更早的评论隐藏 1更早的评论

Torsten 2017-7-5

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(...) ?

Please clarify.

Best wishes

Torsten.

Cheung Ka Ho 2017-7-5

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.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Torsten 2017-7-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/347116-solving-a-sum-of-series-of-exponential-function-with-a-sum-of-series-of-cosine-function-inside#answer_272956

编辑：Torsten 2017-7-5

在 MATLAB Online 中打开

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 个评论
显示 -1更早的评论隐藏 -1更早的评论

Cheung Ka Ho 2017-7-10

Thank you!

请先登录，再进行评论。

更多回答（2 个）

Matthew Taliaferro 2017-7-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/347116-solving-a-sum-of-series-of-exponential-function-with-a-sum-of-series-of-cosine-function-inside#answer_272695

编辑：Matthew Taliaferro 2017-7-3

在 MATLAB Online 中打开

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 个评论
显示 -1更早的评论隐藏 -1更早的评论

Cheung Ka Ho 2017-7-4

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?

请先登录，再进行评论。