#3:
#6: (the "show that if" condition follows from the property, the question is likely misstated)
Thanks, I actually wanted to get rid of the earlier condition that f(x)≥x for all x, and I did that.
Answer to question 1.
Answer to question 2.
Answer to question 3.
For question 2
you haven't proven f is continuous
For question 3 you say
f is a contraction map because f is differentiable and ... ∀x:|f′(x)|<1
I would think proving this is part of what is asked for.
#6:
#7:
#8:
#9:
#10:
Q1
Q2
Q3
Ex 6:
Ex 7:
Ex 8:
Ex1
Let x0:=x, let D:=d(x0,x1), and let k∈N. Then d(xk,xk+1)=d(f(xk−1,f(xk))≤qd(xk−1,xk)⋯≤qkD. For each ϵ∈R+, we find an n∈N such that qn<ϵD, then d(xn,xn+1)≤qnD<ϵ. This proves that (xk)k∈N is a Cauchy-sequence, which (because (X,d) is complete) means it converges to some point x∗∈X.
Furthermore, given a position n∈N, we have
d(xn,x∗)≤∞∑k=nd(xk,xk+1)≤∞∑k=nqkD=qnD∞∑k=0qk=qnD11−q=D⋅qn1−q<D⋅(q1−q)n.
Ex2
Given a any sequence (xk)k∈N in X, it converges to some point x∗, and it's easy to see that x∗ is a fixed point of f. Let y∗ be a fixed point of X. Then, d(f(x∗),f(y∗))=d(x∗,y∗), hence x=y. (Otherwise, this contradicts the fact that f is a contraction.)
Ex3
Choose X⊊R2 as X:=R×N+. Then X is complete because, given any Cauchy sequence (xk)k∈N, it's easy to prove that there is an n∈N such that all but finitely many xk are in the subspace R×{n}. However, the map f:X→X given by f(x,n):=f(x+1n,n+1) has no fixed points since it moves each point by at least 1. (And it's straight-forward to verify that f is a weak contraction.)
#1
d(fn(x),fn+1(x))≤qnd(x,f(x)) - can show by induction.
∀m>n,d(fn(x),fm(x))≤qn+qn+1+…+qm−1≤qn1−q→n→∞0
Therefore, fn(x) is a Cauchy sequence, and since (X, d) is complete, it must have a limit in X. Suppose y=limn→∞fn(x) . Then d(y,fn+1(x))≤qd(y,fn(x)) , therefore d(y,fn(x))≤qnd(x,y)
#2
Suppose y=limn→∞fn(x) . Let's show that y is a fixed point. Indeed, for any n, d(fn(x),f(y))≤qd(fn+1(x),y) , and if we take the limit in both sides, we get d(y,f(y))≤qd(y,y)=0 .
Let's show uniqueness: suppose x and y are fixed points, then d(x,y)=d(f(x),f(y))≤qd(x,y) , therefore d(x,y) = 0.
#3
X=[1,+∞) , f(x) = x + 1/x.
#4
Suppose f(x)=ϵ||x||2+h(x) , where h is some convex function and ϵ<1/2. Take x,y∈Rn. Since h is convex on segment [x,y], its directional derivative is nondecreasing. Its directional derivative is a projection of gradient of g on the [x,y] line. Therefore, we have ⟨∇h(y),y−x⟩≥⟨∇h(x),y−x⟩ , or ⟨∇h(y)−∇h(x),y−x⟩≥0 . Hence,
||x−∇f(x)−y+∇f(y)||=||x−2ϵx−∇h(x)−y−2ϵy+∇h(y)||≤(1−2ϵ)||x−y||
Therefore, g is a contraction mapping, and from problem 1 it follows that the gradient descent converges exponentially quickly.
#5
Suppose A is an NxN positive matrix, and e is its minimal entry. (Then e < 1/N). Then we can write A = eJ + (1 - Ne)Q, where J is a matrix whose entries are all 1, and Q is a matrix whose entries are all nonnegative and the sum of each column is 1 (because the sum of each column is 1 in A and Ne in J). Suppose x and y are probability distributions, i.e. N-dimensional vectors with nonnegative entries whose sum is 1. Then ||Ax−Ay||1=||eJ(x−y)+(1−Ne)Q(x−y)||1=(1−Ne)||Q(x−y)||
Denote x+=max(x−y,0) , x−=min(x−y,0) (pointwise max/min). Then x−y=x+−x− ,||Q(x−y)||1=||Qx+−Qx−||1≤||Qx+||1+||Qx−||1=||x+||1+||x−||1=||x−y||1 ,
so ||Ax−Ay||1≤(1−Ne)||x−y||1 . The space of all probability distributions with metric induced by ||.||1- norm is a compact subset of Rn, so it is a complete metrics space, therefore, the sequence An(x) converges to a unique fixed point.
#6
Let us assume x≤f(x) (the proof for x≥f(x) is the same). Then, from monotonicity of f, x≤f(x)≤f(f(x))≤… is an ascending chain. This sequence cannot have more that |P| distinct elements, so an element of this sequence is going to repeat: fm(x)=fn(x),m<n,m<|P| . Then all the inequalities in fm(x)≤fm+1(x)≤…≤fn(x) must be equalities, so f(fm(x))=fm(x) , fm(x) is a fixed point.
I have fond memories of the contraction mapping theorem, because it was the first fixed point theorem I ever learned.
1:
Let c = d(x, f(x)). Then d(f(x), f^2(x)) <= q d(x, f(x)) = qc, and in general d(f^n(x), f^{n+1}(x)) <= q^n c.
Now, d(f^n(x), f^{m}(x)) <= c * \sum_{n <= i <= m} q^i which is a geometric series. The sum of this is then c q^n (1 - q^{m-n})/(1 - q) <= the full series (which converges since q < 1) = c q^n/(1 - q)
This can be made arbitrarily small, so the sequence is Cauchy, so since we are in a complete metric space it converges.
The distance to the limit point is at most that value c q^n/(1 - q), because all terms after the nth iterate are within that of the nth iterate. This gives us the exponential convergence.
2:
First, there's at most one fixed point for a contraction mapping. For if we had two, then we could take d(x,y) = d(f(x), f(y)) <= q d(x,y) < d(x,y) for the two fixed points x,y. But then, we get a contradiction - unless x = y, in which case q d(x,y) = d(x,y) = 0.
For existence: the limit of the iterates of any x is our fixed point. This is because f(lim f^n(x)) = lim f(f^n(x)) = lim f^{n+1}(x) = lim f^n(x), where the first equality is from continuity of f
3:
Take the space of infinite binary sequences, where the distance we assign is 1/N where N is the earliest index for which two sequences differ (or 0 if they are the same always). This is a complete metric space, actually (as you can easily check). But the right shift operator is continuous (since if you differ no earlier than the Nth spot, then the image of you differs by no earlier than the (N+1)th spot), and has no fixed point, despite the fact that it also is a weak contraction (by virtue of pushing any differences to the right). This is because N/(N+1) isn't bounded by a constant.
4: TODO
Take two points x, y. We'd like to show that there's an epsilon such that ||(x - eps grad(f)(x)) - (y - eps grad(f)(y))|| is bounded by ||x - y|| times some constant less than 1.
We have ||x - y + eps (grad(f)(y) - grad(f)(x))|| =
5: TODO
6:
WLOG suppose x <= f(x). Then f(x) <= f^2(x), and we can repeat upwards. For n > |P|, by the pigeonhole principle the nth iterate is equal to some earlier iterate, so there's n and m such that m < n f^n(x) = f^m(x). Thus all the iterates between these two are also equal. But then f(f^m(x)) = f^m(x) , so f^m(x) is a fixed point, but this is the same as f^n(x).
7:
f^n(x) is the same as max {f(f^m(x)) : m < n}. So let's define f^alpha(x) as sup {f(f^beta(x)) : beta < alpha}
WLOG suppose again that x <= f(x). Now if alpha is a limit ordinal than f(f^beta(x)) for any beta < alpha will just be equal to f^gamma(x) with gamma = beta + 1 < alpha. Thus this is really the sup of all previous iterates, and thus larger than it. For successor ordinals, we need to iterate the previous iterate once more - but we can use monotonicity here to get that we are larger than the previous iterate.
Thus we see that f^beta(x) <= f^alpha(x) for all beta < alpha.
Suppose we have some alpha such that |alpha| > |L|. Then if this sequence is made of all unique elements, we would be able to monotonically inject alpha into L by sending each smaller ordinal to the associated iterate. This totally ordered chain of iterates has an ordinal type, which we just mapped alpha to an initial element of, and so we have |alpha| <= |L| as ordinals, a contradiction (if instead we were to take them as cardinals, we'd still have a contradiction, because we can ordinal inject the least ordinal with the same cardinality with alpha, which then gives us an ordinal injection to L). Thus we have some repeat where f^beta(x) = f^alpha(x) for beta < alpha. But then all the iterates between these two are also equal.
So f(f^(beta)(x)) = f^beta(x), so f^beta(x) is a fixed point of f, f^beta(x) = f^alpha(x) so f^alpha(x) is a fixed point of f.
8:
The question is asking whether the set of fixed points has a sup and an inf.
Suppose we have some subset S of fixed points. I hope that the sup of these fixed points in the original lattice is actually a fixed point.
For any element x, the supremum s is larger than x, and so f(x) = x <= f(s). Thus f(s) is an upper bound for the subset, and thus s <= f(s). But then for some alpha, f^alpha(s) is a fixed point of f that's an upper bound. So now I hope f^alpha(s) is our supremum in the fixed point lattice.
Suppose there was an upper bound that was a fixed point, u. Since s is a supremum in the encompassing lattice L, we have s <= u. Thus f(s) <= f(u) = u. In general, f^n(s) <= f^n(u) = u. We can also take a sup of them to get f^beta(s) <= u for limit ordinals beta. Therefore f^alpha(s) <= u. But then f^alpha(s) is smaller than all upper bounds that are fixed points, and so f^alpha(s) is a supremum in the fix point lattice
Likewise for infimums.
Incidentally: The least element l satisfies l <= f(l), so there's an alpha such that the alpha iterate is a fixed point. If x is a fixed point, then since l <= x for all x, we have f^alpha(l) <= f^alpha(x) = x. But then f^alpha(l) is the least fixed point.
9: TODO
P(A) is clearly a lattice. To take sups and infs, we can just take arbitrary unions and intersections. (I think technically you need the axiom of choice, but whatever). Given a function f: A -> B, we have that S <= S' <= A implies f^img(S) <= f^img(S') (we don't need injectivity for this), so the induced f is monotone.
Since f,g are injective, we know that the induced functions on the powerset lattices are also injective - that is, that f^img(S) != f^img(S') if S != S'. Also, for injective functions, f(X - Y) = f(X) - f(Y)
...
10:
Using problem 9, there is a set A' and B' such that f(A') = B' and A - A' = g(B - B')
Now, we'd like to biject A' with B' and A - A' with B - B'. To do this, we can take the restriction of f on A' (restricting the codomain to B' as well), and likewise take the restriction of g. The restrictions are bijections, because they are surjections and we already knew they were injections.
Then we can take the union of the relation associated to the restriction of f and the inverse of the relation associated with the restriction of g. This is then a function A -> B, and it's a bijection.
(The second half made me realize how much more comfortable I am with abstract exercises than with regular Analysis à la Ex4.)
Ex6
If f(x)≤x, then all fn(x) are comparable to each other: we have
f(x)≤x⟹f(f(x))≤f(x)⟹f3(x)≤f2(x)
and so on. Furthermore, if n∈N is such that fn(x)≠fn−1(x), then fk(x)≠fk−1(x) for all k∈[n] as well (verify by looking at the contrapositive). Consequently (set n:=|P|), if fn(x) were not a fixed point of f, then fn+1(x)<fn(x), and hence {x,f(x),f2(x),...,fn+1(x)}⊆P, which means P would have |P|+2 elements.
If x≤f(x), we get x≤f(x)≤f2(x) and so on, leading to the same argument.
Ex7
Wlog, assume x≤f(x). Set f0(x):=x. Given any non-limit ordinal β, we find a predecessor α and set fβ(x):=f(fα(x)). Given any limit ordinal ω, we set fω(x):=sup{fα(x)|α∈ω}.
Suppose this doesn't define fα(x) for all ordinals α. Then, there is some smallest ordinal α∗ such that fα∗(x) is not defined. This immediately yields a contradiction (regardless of whether α∗ is a limit ordinal or not).
We want this construction to have the properties that f(x)=x⟹fβ(x)=x and that x≤y⟹fβ(x)≤fβ(y). Thus, let x,y∈L and β be an ordinal. If β has a predecessor, the check for both properties are easy. If not, then fβ(x)=sup{fα(x)|α∈β} and fβ(y)=sup{fα(y)|α∈β}. Then, for the first property, note that the upper-bound of a one-element set is just the element itself. For the second, note that each element in the first set is smaller than some element in the second set, so fβ(y) is an upper-bound for the first set, which implies that fβ(x)≤fβ(y) since fβ(x) is the lowest upper-bound.
Now, given an ordinal α, our construction defines a function ϕ:α→L. If f(fα(x))≠fα(x), then the chain doesn't become stationary at any earlier point either (to verify, take a smallest α such that [f(fα(x))≠fα(x) but the chain is stationary for smaller ordinals] and derive a contradiction), and hence ϕ is injective, proving that α≤|L|. (This is the generalized version of the argument from Ex6.)
Ex8
Let f:L→L be monotonic and let L′ be the set of fixed points of f. Then L′ inherits the partial order from L; what needs doing is verify the least upper-bound property. So let X⊆L′. Then, X has a least upper-bound u in L.
Let ω be some ordinal with ω>|L|. From the previous exercise, we know that f(fω(u))=fω(u). Choose the smallest α such that f(fα(u))=fα(u). Then, fα(u)∈L′ and x≤u≤fα(u), hence fα(u) is an upper-bound of X.
It remains to show that it is the least upper-bound. Thus, let u′∈L′ be another upper-bound of x. Then, u≤u′ in L, hence fα(u)≤fα(u′)=u′ (apply Ex7).
Ex9
A least upper-bound is obtained via ⋃ on all sets, and the greatest lower-bound via ⋂. (Easy checks.) Given any function f:A→B, we trivially have X⊆Y⟹f(X)⊆f(Y); injectivity is not needed.
We define
We need to verify that g(B−f(A′))=A−A′, then A′ and f(A′) are the desired sets.
x∈A(k−1)−A(k)=g(B−f(A(j)))⊆g(B−f(A′))
Ex10
Let A′ be the set constructed in Ex9. Then, we can define a bijection ϕ:A→B via
ϕ:x↦{f(x)x∈A′g−1(x)x∉A′
Ex5 (this is super ugly but I don't think it's worth polishing and it does work. All important ideas are in the first third of the proof, the rest just inelegantly resolves the details.)
We define our metric space as (X,d) where X:=[0,1]d is the set of probability distributions, and d(x,y)=∑dj=1|xj−yj|. Let x,y∈X and let Δ:=x−y, then d(Ax,Ay) can be computed as
d∑i=1|(Ax−Ay)i|=d∑i=1|(AΔ)i|=d∑i=1|d∑k=1ak,iδk|≤d∑i=1d∑k=1|ai,kδk|=d∑i=1δi
where the last step holds because multiplying a vector with the state-transition matrix leaves the sum of entries unchanged. (Reasonably easy to verify using that each column of A sums up to 1.)
If x≠y, then Δ has at least one negative entry and the inequality is strict. In that case, let k=argmaxi∈{1,...,d}δi and ℓ=argmini∈{1,...,d}δi. In particular, we have δk>0>δℓ. Note that, when two numbers a,b∈R have different sign, then |a+b|=||a|−|b|| and thus |a|+|b|−|a+b|=min(2|a|,2|b|). Therefore, the amount that gets canceled out is at least
d∑i=1|ai,kδk|+|ai,ℓδℓ|−|ai,kδk+ai,ℓδℓ|=d∑i=12min(|ai,kδk|,|ai,ℓδℓ|)
Let a′ be the smallest entry in A, then we can lower-bound the above as
d∑i=1a′2min(|δk|,|δℓ|)=2a′dmin(|δk|,|δℓ|)
Wlog, let |δk|>|δℓ|. Let K be the sum of all postive entires of Δ, then ∑di=1|δi|=2K, so the term we want to lower-bound is 2a′dδℓ2K=a′dδℓK. The sum of the negative entries is −K, which means that the one with largest norm among them has norm at least 1d−1|K|. Thus, the relative decrease is at least a′dKd−1K=a′dd−1>a′
Then, d(x,y)−d(A(x),A(y))d(x,y)≥a′, hence d(A(x),A(y)d(x,y)≤1−a′. This proves that A is a contraction; apply Banach's theorem.
This is the third of three sets of fixed point exercises. The first post in this sequence is here, giving context.
Note: Questions 1-5 form a coherent sequence and questions 6-10 form a separate coherent sequence. You can jump between the sequences.
Let (X,d) be a complete metric space. A function f:X→X is called a contraction if there exists a q<1 such that for all x,y∈X, d(f(x),f(y))≤q⋅d(x,y). Show that if f is a contraction, then for any x, the sequence {xn=fn(x0)} converges. Show further that it converges exponentially quickly (i.e. the distance between the nth term and the limit point is bounded above by c⋅an for some a<1)
(Banach contraction mapping theorem) Show that if (X,d) is a complete metric space and f is a contraction, then f has a unique fixed point.
If we only require that d(f(x),f(y))<d(x,y) for all x≠y, then we say f is a weak contraction. Find a complete metric space (X,d) and a weak contraction f:X→X with no fixed points.
A function f:Rn→R is convex if f(tx+(1−t)y)≤tf(x)+(1−t)f(y), for all t∈[0,1] and x,y∈Rn. A function f is strongly convex if you can subtract a positive parabaloid from it and it is still convex. (i.e. f is strongly convex if x↦f(x)−ε||x||2 is convex for some ε>0.) Let f be a strongly convex smooth function from Rn to R, and suppose that the magnitude of the second derivative ∥∇2f∥ is bounded. Show that there exists an ε>0 such that the function g:Rn→Rn given by x↦x−ε(∇f)(x) is a contraction. Conclude that gradient descent with a sufficiently small constant step size converges exponentially quickly on a strongly convex smooth function.
A finite stationary Markov chain is a finite set S of states, along with probabilistic rule A:S→ΔS for transitioning between the states, where ΔS represents the space of probability distributions on S. Note that the transition rule has no memory, and depends only on the previous state. If for any pair of states s,t∈ΔS, the probability of passing from s to t in one step is positive, then the Markov chain (S,A) is ergodic. Given an ergodic finite stationary Markov chain, use the Banach contraction mapping theorem to show that there is a unique distribution over states which is fixed under application of transition rule. Show that, starting from any state s, the limit distribution limn→∞An(s) exists and is equal to the stationary distribution.
A function f from a partially ordered set to another partially ordered set is called monotonic if x≤y implies that f(x)≤f(y). Given a partially ordered set (P,≤) with finitely many elements, and a monotonic function from P to itself, show that if f(x)≥x or f(x)≤x, then fn(x) is a fixed point of f for all n>|P|.
A complete lattice (L,≤) is a partially ordered set in which each subset of elements has a least upper bound and greatest lower bound. Under the same hypotheses as the previous exercise, extend the notion of fn(x) for natural numbers n to fα(x) for ordinals α, and show that fα(x) is a fixed point of f for all x∈X with f(x)≤x or f(x)≥x and all |α|>|L| (|A|≤|B| means there is an injection from A to B, and |A|>|B| means there is no such injection).
(Knaster-Tarski fixed point theorem) Show that the set of fixed points of a monotonic function on a complete lattice themselves form a complete lattice. (Note that since the empty set is always a subset, a complete lattice must be nonempty.)
Show that for any set A, (P(A),⊆) forms a complete lattice, and that any injective function from A to B defines a monotonic function from (P(A),⊆) to (P(B),⊆). Given injections f:A→B and g:B→A, construct a subset A′ of A and a subset of B′ of B such that B′=f(A′) and A−A′=g(B−B′).
(Cantor–Schröder–Bernstein theorem) Given sets A and B, show that if |A|≤|B| and |A|≥|B|, then |A|=|B|. (|A|≤|B| means there is an injection from A to B, and |A|=|B| means there is a bijection)
Please use the spoilers feature - the symbol '>' followed by '!' followed by space -in your comments to hide all solutions, partial solutions, and other discussions of the math. The comments will be moderated strictly to hide spoilers!
I recommend putting all the object level points in spoilers and including metadata outside of the spoilers, like so: "I think I've solved problem #5, here's my solution <spoilers>" or "I'd like help with problem #3, here's what I understand <spoilers>" so that people can choose what to read.
Tomorrow's AI Alignment Forum Sequences post will be "Approval-directed agents: overview" by Paul Christiano in the sequence Iterated Amplification.
The next post in this sequence will be released on Saturday 24th November, and will be 'Fixed Point Discussion'.