泛函分析作业:2025-09-15
泛函分析作业,练习 Newton 迭代、压缩映射、不动点定理及积分方程的可解性。
本页目录5 节
1.1.2 #
(Newton 法) 设 是定义在 上的二次连续可微的实值函数, 使得 . 求证: 存在 的邻域 , 使得 , 迭代序列
是收敛的, 并且
证明
令 。由 和 ,可在 附近选取闭区间 ,使 在 上不为零,且
这样的区间存在,是因为 连续且 。 由中值定理和 ,对 有
因此 把完备度量空间 映入自身,并且是压缩映射。 Banach 不动点定理给出唯一不动点 ,以及任意 的迭代收敛性。 取 即得题目结论。 这里在闭区间 上应用定理,再得到邻域内的结论;开区间本身不必完备。
下面保留 Mworks 程序与对应结果图,程序未在本轮重新运行。
clear()
clc()
f(x)=exp.(x).-2 #求此函数的零点
df(x)=exp.(x) #此函数的导数
#绘制f的函数图像
t=0:0.1:1.1
plot(t,f(t),"k")
hold("on")
#绘制直线y=0
plot(t,zeros(size(t)),"r")
x0=0 #初始猜测
N=3 #总迭代次数
x=zeros(1,N+1)
x[1]=x0;
for ii=1:N
x[ii+1]=x[ii]-f(x[ii])/df(x[ii])
text(x[ii],f(x[ii]),string(ii-1),color="black",fontsize=14)
plot(x[ii:ii+1],f(x[ii:ii+1]),"b")
plot(x[ii+1],f(x[ii+1]),"r*")
end
text(x[N+1],f(x[N+1]),string(N),color="black",fontsize=14)
saveas(gcf(),"1d1d2.eps")
[x;
abs.(x.-log(2))] #零点为log(2)

1.1.4 #
设 是度量空间上的压缩映射, 求证: 是连续的.
证明
记度量空间为 .
设 是该度量空间上的任意收敛列, 且收敛至 . 那么对于 , 由于 是压缩映射, 故存在 满足 , 于是 , 根据夹逼准则可知 , 所以 是连续的.
1.1.5 #
设 是压缩映射, 求证: 也是压缩映射, 并说明逆命题不一定成立.
证明
记度量空间为 ,设 的压缩系数为 。 逐次应用 ,得到
又 仍将 映入自身,且 ,故它是压缩映射。
逆命题不成立。在 上取通常距离,令
因为 ,它是到自身的映射。 在靠近 时大于 ,所以 不是压缩映射。 但 ,从而
由中值定理可知 是压缩映射。
构造的关键是链式法则 。 更一般地,令 ,取整数 与 ,可要求
上例取 , 是 与 的中点。
1.1.6 #
设 是 中的有界闭集, 映射 满足: , 求证: 在 中存在唯一的不动点.
本题按通常欧氏距离理解 ,并假定 非空。若 是任意度量,仅“闭且有界”并不保证下面用到的紧性;相应的一般形式应假设 是非空紧度量空间。
证明
由 和三角不等式可得
从而 是连续的.
由于 是 中的有界闭集, 所以存在最小值点 , 使得 是最小的. 如果 那么 就是不动点.
反之, 如果 , 那么 , 从而 更小, 与 是最小值点矛盾.
而不动点的唯一性和压缩映像原理相同, 假设有两个不动点, 那么 , 矛盾.
综上, 在 中存在唯一不动点.
注: 在度量空间中的连续要转变思维用度量去衡量, 即在取 时, 是对 的限制. 注意区分度量在欧式距离中的连续, 和度量空间中的连续.
1.1.7 #
对于积分方程
其中 为给定函数, 为常数, , 求证: 存在唯一解 .
证明
在完备空间 上取一致范数 。 令 、,则原方程等价于
定义 。 它将 映入自身;对任意 ,差 为常数,故
压缩系数是 ,而不是可能为负的 。 由 Banach 不动点定理, 存在且唯一;乘以 是 上的双射,故 也存在且唯一。
还可以直接核对该解。对变换后的方程积分,得到
题设 保证分母不为零;代回原式即可验证。
下面保留 Mworks 程序与对应结果图,程序未在本轮重新运行。
clear()
clc()
n=20; #离散节点总数
t=1/(2*n):1/n:1-1/(2*n) #居中节点
iter=5 #总迭代次数
x=zeros(n,iter+1)
plot(t,x[:,1],"b*",linewidth=4) #绘制初始猜测(蓝色点)
hold("on")
text(t[end],x[end,1],string(0),color="black",fontsize=14)
lambda=0.5
y=(1-lambda).*exp.(t)
for ii=1:iter
x[:,ii+1]=lambda/n*exp.(t .- transpose(t))*x[:,ii]+y
color_value=2*ii/iter-1 #取值范围调整为-1 到 1, 为使线段颜色可以rgb插值渐变
if color_value>0
plot(t,x[:,ii+1],color=[color_value,1-color_value,0]) #绘制此步结果
else
plot(t,x[:,ii+1],color=[0,1+color_value,-color_value]) #绘制此步结果
end
text(t[end],x[end,ii+1],string(ii),color="black",fontsize=14)
end
x_true=exp.(t)
plot(t,x_true,"r*",linewidth=4) #绘制精确结果(红色点)
saveas(gcf(),"1d1d7.eps")
[x x_true]

讨论
评论
正在加载评论…