How can i Calculate this with matlab?

2 次查看（过去 30 天）

显示更早的评论

maziar 2021-12-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1600770-how-can-i-calculate-this-with-matlab

评论： maziar 2021-12-16

采纳的回答： Star Strider

Alpha=0.1

L=1

t=0.05

x=0.2

采纳的回答

Star Strider 2021-12-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1600770-how-can-i-calculate-this-with-matlab#answer_845355

在 MATLAB Online 中打开

Using the Symbolic Math Toolbox —

syms m

alpha = 1;

L = 1;

t = 0.05;

x = 42; % Not Supplied

arg = exp(-(m*pi/L)^2*alpha*t) * (1 - (-1)^m)/(m*pi) * sin(m*pi*x/L)

arg =

S = symsum(arg, m, 1, Inf)

S =

Svpa = vpa(S)

Svpa =

0.0

I am not certain where the

term is supposed to go, so I included it as part ot the exponent argument because it appears to mne that is where it belongs. Change that if I guessed in error.

Experiment to get different results.

11 个评论
显示 9更早的评论隐藏 9更早的评论

Walter Roberson 2021-12-16

在 MATLAB Online 中打开

The m of symsum() is what is known as a "bound" variable. It should not need to be input to the calculation because it is effectively local to the symsum.

Unfortunately, in some releases of the Symbolic Toolbox, matlabFunction() generates bad code for handling bound variables. Sometimes it fails to sym() it into existence to do the calculations with.

Sometimes it can end up using the bound variable outside of the bound context, generating incorrect results, especially if 'optimize' is turned on. (optimize is off when you are generating to an anonymous function but it is on by default if you are generating to a file)

In the simple case where it just forgets to sym() m into existence, the work-around is

Sfcn = matlabFunction(S, 'vars', {x, m} )

and later

Sx = Sfcn(x, m);

In the case where it generates bad code because of optimize to a file, turn optimize off in the matlabFunction options

In the case where it generates bad references to the bound variable even with optimization off, you have to find a different way to calculate the sum, such as what Yongjian Feng suggests.

maziar 2021-12-16

thanks everybody for suggestions and i must say that i am using Matlab R2019b.

请先登录，再进行评论。

更多回答（1 个）

Yongjian Feng 2021-12-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1600770-how-can-i-calculate-this-with-matlab#answer_845255

编辑：Yongjian Feng 2021-12-1

在 MATLAB Online 中打开

If the series is convergent, one way is to set a cutoff value (for example epsilon = 0.00000001). Then when you do the sum from m going up, compute each term. If a term is smaller than epsilon (when m is big enough) then stop. Something like this

epsilon = 0.000000001; % this determines how accurate your result is
result = 0;
done = false;
m = 0;
while ~done
    next_term = 0; % Do you calculation here for this m value!!!!
    result += next_term; % accumulate
    if next_term < epsilon
        % the increment is small enough, stop the while loop
        done = True;
    end
    m += 1; % go to next term
end