How to plot a convergence graph for a double sum?

2 次查看（过去 30 天）
Elizabeth McLean
Elizabeth McLean 2021-11-10
编辑： VBBV 2021-11-11
I'm trying to show converengence for a double sum but I can't figure out how to plot the convergence. I think I need to create an array to store the previous values of the series but I'm not sure how to do that.
This is the series I'm trying to graph:
This is my current code:
T=0; %terms used in series
q_m=0;
M=15000;
N=15000;
x=2.5;
y=3.5;
for n=1:N
for m=1:M
T=T+1;
q_m1 = (400/(m*n*pi^2))*(cos(n*pi*4/10)-cos(n*pi*3/10))*(cos(m*pi*3/10)-cos(m*pi*2/10))*sin(m*pi*x/10)*sin(n*pi*y/10);
q_m = q_m + q_m1;
end
end
q_m
T

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

回答（1 个）

VBBV
VBBV 2021-11-10
编辑：VBBV 2021-11-11
Ran in:
T=0; %terms used in series
M=15;
N=15;
q_m=zeros(M,1);
qq_m=zeros(M,1);
x=2.5;
y=3.5;
for n=1:N
for m=1:M
T=T+1;
q_m1 = (400/(m*n*pi^2))*(cos(n*pi*4/10)-cos(n*pi*3/10))*(cos(m*pi*3/10)-cos(m*pi*2/10))*sin(m*pi*x/10)*sin(n*pi*y/10);
q_m(m) = q_m1;
end;
qq_m(n) = sum(q_m);
end
plot(qq_m)
T
T = 225
You can modify the code as above to plot convergence
  4 个评论
显示 2更早的评论
Star Strider
Star Strider 2021-11-11
@Elizabeth McLean — It is appropriate and appreciated to Accept answers that solve the problem.
Elizabeth McLean
Elizabeth McLean 2021-11-11
Sorry about that, this is my first time asking a question. I've accepted the answer now.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by