How to plot a function dependent on no. of iterations, provided it plots till nth iteration where n is the number where the function is converged?

3 次查看(过去 30 天)
Hello All,
I am trying to plot a function vs no. of iterations. In the for loop, i have given the convergence condition, but I am unable to plot for f(x) vs no. iterations (till convergence). Please help.
Its of the form:
for m=1:n
pu(m)=(-1)^m)/((2*m)+1))*(exp(-((2*m)+1)))
e = e+pu;
if (abs(pu)/e)<(10^(-6)) %convergence check
break
end
end
hold all
plot(1:n,pu,'linestyle','none','marker','*')

采纳的回答

VBBV
VBBV 2022-1-28
编辑:VBBV 2022-1-29
n= 10;
e = 0;
for m=1:n
pu(m)=(-1)^m/((2*m)+1)*(exp(-((2*m)+1)))
e = e+pu(m);
if (abs(pu(m))/e)>(10^(-6)) %convergence check
break
end
end
pu = -0.0166
pu = 1×2
-0.0166 0.0013
pu = 1×3
-0.0166 0.0013 -0.0001
pu = 1×4
-0.0166 0.0013 -0.0001 0.0000
pu = 1×5
-0.0166 0.0013 -0.0001 0.0000 -0.0000
pu = 1×6
-0.0166 0.0013 -0.0001 0.0000 -0.0000 0.0000
pu = 1×7
-0.0166 0.0013 -0.0001 0.0000 -0.0000 0.0000 -0.0000
pu = 1×8
-0.0166 0.0013 -0.0001 0.0000 -0.0000 0.0000 -0.0000 0.0000
pu = 1×9
-0.0166 0.0013 -0.0001 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000
pu = 1×10
-0.0166 0.0013 -0.0001 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000
hold all
plot(1:n,pu,'linestyle','-','marker','.')
Check the convergence condition
  7 个评论
VBBV
VBBV 2022-1-29
legend(mrkr{cidx},legendInfo)
Put this line inside the inner for loop. After legendInfo. Everytime the inner for loop exits, it's over written with initial values for mrkr and clr index.
VBBV
VBBV 2022-1-29
编辑:VBBV 2022-1-29
u=zeros(n,npoints); % change also This line in beginning
u(m,cidx)=pu(m);
hold all
plot(1:m,u(1:m,cidx),'linestyle','none','marker',mrkr{cidx},'Color',clr{cidx})
You also need to change these two lines inside the for loop. As you can see the u value is dependent on pu which is again varying with each outer iteration of for loop.

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Graphics Performance 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by