implementation of iterative function

2015 3 19
1 个回答
回答已采纳
更新时间：2015 3 19
5 次查看（30 天）

0 个投票

How can we implement the iterative function F(x)= [((1-2*q)*(x-q))/(p-q)]+q we have to iterate it for some values of i and generate the values of xi where xi=F(xi-1). initially let.

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

 采纳的回答

Torsten
Torsten 2015-3-19

0 个投票

n=...;
p=...;
q=...;
x0=...;
f=@(x)(1-2*q)*(x-q)/(p-q)+q;
for ii=1:n
xnew=f(xold);
xold=xnew;
end
Best wishes
Torsten.

8 个评论
显示 6更早的评论 隐藏 6更早的评论

shivani
shivani 2015-3-19
here xold=x0??? and we have to store all the values of xnew into a vector.
shivani
shivani 2015-3-19
with the help of this getting inf value of x in each iteration.
Torsten
Torsten 2015-3-19
ok, replace
x0=...;
by
xold=...;
And don't use p=q -:)
Best wishes
Torsten.
shivani
shivani 2015-3-19
not using p=q still getting inf values. Im=imread(); [n m]=size(A); L=m*n; AVector = A(:); P=A; x0=0.27;p=0.4;q=0.1; xold=x0; f=@(x)(1-2*q)*(x-q)/(p-q)+q; for i=1:n; xnew=f(xold); xold=xnew; end for j = 1:numel(P) i = P(j); output(j) =xold end
Stephen23
Stephen23 2015-3-19
编辑：Stephen23 2015-3-19
@shivani: please learn to format your code correctly to make it readable: use the {} Code button that is above the textbox. The Preview pane is pretty handy too!
shivani
shivani 2015-3-19
not taking p=q still getting inf values.
Im=imread();
[n m]=size(A);
L=m*n;
AVector = A(:);
P=A;
x0=0.27;p=0.4;q=0.1;
xold=x0;
f=@(x)(1-2*q)*(x-q)/(p-q)+q;
for i=1:n;
xnew=f(xold);
xold=xnew;
end
for xi = 1:numel(P);
i = P(xi);
output(xi) =xold;
end
Torsten
Torsten 2015-3-19
Im=imread();
[n m]=size(A);
L=m*n;
x0=0.27;
p=0.4;
q=0.1;
xold=x0;
output=zeros(L);
output(1)=xold;
f=@(x)(1-2*q)*(x-q)/(p-q)+q;
for ii=1:L-1
xnew=f(xold);
xold=xnew;
output(ii+1)=xold;
end
Note that xold becomes bigger and bigger during the Iteration such that you may get xold=Infinity if L is large enough.
Best wishes
Torsten.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Inverse Gaussian Distribution 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by