【矩阵论】7. 矩阵理论——算子范数

这篇具有很好参考价值的文章主要介绍了【矩阵论】7. 矩阵理论——算子范数。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

矩阵论的所有文章,主要内容参考北航赵迪老师的课件

[注]由于矩阵论对计算机比较重要,所以选修了这门课,但不是专业搞数学的,所以存在很多口语化描述,而且对很多东西理解不是很正确与透彻,欢迎大家指正。我可能间歇性忙,但有空一定会回复修改的。

矩阵论
1. 准备知识——复数域上矩阵,Hermite变换
1.准备知识——复数域上的内积域正交阵
1.准备知识——Hermite阵,二次型,矩阵合同,正定阵,幂0阵,幂等阵,矩阵的秩
2. 矩阵分解——SVD准备知识——奇异值
2. 矩阵分解——SVD
2. 矩阵分解——QR分解
2. 矩阵分解——正定阵分解
2. 矩阵分解——单阵谱分解
2. 矩阵分解——正规分解——正规阵
2. 矩阵分解——正规谱分解
2. 矩阵分解——高低分解
3. 矩阵函数——常见解析函数
3. 矩阵函数——谱公式,幂0与泰勒计算矩阵函数
3. 矩阵函数——矩阵函数求导
4. 矩阵运算——观察法求矩阵特征值特征向量
4. 矩阵运算——张量积
4. 矩阵运算——矩阵拉直
4.矩阵运算——广义逆——加号逆定义性质与特殊矩阵的加号逆
4. 矩阵运算——广义逆——加号逆的计算
4. 矩阵运算——广义逆——加号逆应用
4. 矩阵运算——广义逆——减号逆
5. 线性空间与线性变换——线性空间
5. 线性空间与线性变换——生成子空间
5. 线性空间与线性变换——线性映射与自然基分解,线性变换
6. 正规方程与矩阵方程求解
7. 范数理论——基本概念——向量范数与矩阵范数
7.范数理论——基本概念——矩阵范数生成向量范数&谱范不等式
7. 矩阵理论——算子范数
7.范数理论——范数估计——许尔估计&谱估计
7. 范数理论——非负/正矩阵
8. 常用矩阵总结——秩1矩阵,优阵(单位正交阵),Hermite阵
8. 常用矩阵总结——镜面阵,正定阵
8. 常用矩阵总结——单阵,正规阵,幂0阵,幂等阵,循环阵


【矩阵论】7. 矩阵理论——算子范数

7.2 算子范数

7.2.1 定义

C n C^n Cn 上任一向量范数 ∥ X ∥ V \Vert X\Vert_V XV 都产生一个矩阵范数 ∥ A ∥ = max ⁡ x ≠ 0 { ∥ A X ∥ V ∥ X ∥ V } \Vert A\Vert=\max_{x\neq 0}\limits \{\frac{\Vert AX\Vert_V}{\Vert X\Vert_V}\} A=x=0max{XVAXV} , X ∈ C n X\in C^n XCn ,且有相容关系 ∥ A X ∥ V ≤ ∥ A ∥ ⋅ ∥ X ∥ V \Vert AX\Vert_V\le \Vert A\Vert\cdot \Vert X\Vert_V AXVAXV ,这种范数 ∥ A ∥ \Vert A\Vert A为算子范数

向量范数生成的矩阵范数

  • 若存在常数 M,使 ∀ X ∈ C n \forall X\in C^n XCn ,有 ∥ A X ∥ V ≤ M ∥ X ∥ V \Vert AX\Vert_V\le M\Vert X\Vert_V AXVMXV ,则 ∥ A ∥ ≤ M \Vert A\Vert\le M AM ,即 ∥ X ∥ V \Vert X\Vert_V XV 的算子范数是使上述不等式成立的最小范数

【矩阵论】7. 矩阵理论——算子范数

a. 证明算子范数是范数

【矩阵论】7. 矩阵理论——算子范数

矩阵范数与向量范数相容性:

∥ A ∥ = max ⁡ X ≠ 0 { ∥ A X ∥ V ∥ X ∥ V } ⇒ ∥ A X ∥ V ∥ X ∥ V ≤ max ⁡ X ≠ 0 { ∥ A X ∥ V ∥ X ∥ V } = ∥ A ∥ ⇒ ∥ A X ∥ V ≤ ∥ A ∥ ⋅ ∥ X ∥ V \begin{aligned} &\Vert A\Vert=\max_{X\neq 0}\{\frac{\Vert AX\Vert_V}{\Vert X\Vert_V}\}\Rightarrow \frac{\Vert AX\Vert_V}{\Vert X\Vert_V}\le \max_{X\neq 0}\{\frac{\Vert AX\Vert_V}{\Vert X\Vert_V}\}=\Vert A\Vert\\ &\Rightarrow \Vert AX\Vert_V\le \Vert A\Vert\cdot \Vert X\Vert_V \end{aligned} A=X=0max{XVAXV}XVAXVX=0max{XVAXV}=AAXVAXV

  • 由于 ∥ A X ∥ V \Vert AX\Vert_V AXV 是X各分量的连续函数,故在有界闭集 ∥ X ∥ V \Vert X\Vert_V XV 上可取最大值,因此上述定义是有意义的,即存在 X 0 X_0 X0 使 max ⁡ ∥ X ∥ V = 1 ∥ A X ∥ V = ∥ A X 0 ∥ V \max_{\Vert X\Vert_V=1}\limits \Vert AX\Vert_V=\Vert AX_0\Vert_V XV=1maxAXV=AX0V ,在 X 0 X_0 X0 处取最大值

7.2.2 常见算子范数

A = ( a i j ) ∈ C n . n , X ∈ C n A=\left(a_{ij}\right)\in C^{n.n},X\in C^n A=(aij)Cn.nXCn ,则向量范数 ∥ X ∥ 1 , ∥ X ∥ 2 , ∥ X ∥ ∞ \Vert X\Vert_1,\Vert X\Vert_2,\Vert X\Vert_\infty X1,X2,X 产生的算子范数为:

  • ∥ X ∥ 1 \Vert X\Vert_1 X1 产生 ∥ A ∥ 1 = max ⁡ 1 ≤ j ≤ n ∑ i = 1 n ∣ a i j ∣ \Vert A\Vert_1=\max_{1\le j\le n}\limits \sum_{i=1}\limits^n\vert a_{ij}\vert A1=1jnmaxi=1naij (列范数),且 ∥ A X ∥ 1 ≤ ∥ A ∥ 1 ⋅ ∥ X ∥ 1 \Vert AX\Vert_1\le \Vert A\Vert_1\cdot \Vert X\Vert_1 AX1A1X1
  • ∥ X ∥ 2 \Vert X\Vert_2 X2 产生 ∥ A ∥ 2 = λ 1 , λ 1 为 A H A \Vert A\Vert_2=\sqrt{\lambda_1} ,\lambda_1为A^HA A2=λ1 ,λ1AHA 的最大特征值(谱范数),且 ∥ A X ∥ 2 ≤ ∥ A ∥ 2 ⋅ ∥ X ∥ 2 \Vert AX\Vert_2 \le \Vert A\Vert_2\cdot \Vert X\Vert_2 AX2A2X2
  • ∥ X ∥ ∞ \Vert X\Vert_\infty X 产生 ∥ A ∥ ∞ = max ⁡ i = 1 ∑ j = 1 n ∣ a i j ∣ \Vert A\Vert_\infty =\max_{i=1}\limits\sum_{j=1}\limits^n\vert a_{ij}\vert A=i=1maxj=1naij (行范数),且 ∥ A X ∥ ∞ ≤ ∥ A ∥ ∞ ⋅ ∥ X ∥ ∞ \Vert AX\Vert_\infty\le \Vert A\Vert_\infty\cdot \Vert X\Vert_\infty AXAX
a. 必要条件

任一矩阵范数 ∥ I ∥ ≥ 1 \Vert I\Vert \ge 1 I1

  • 任一算子范数必有 ∥ I ∥ = 1 \Vert I\Vert=1 I=1

    由公式 ∥ A ∥ = max ⁡ x ≠ 0 { ∥ A X ∥ V ∥ X ∥ V } \Vert A\Vert=\max_{x\neq 0}\limits \{\frac{\Vert AX\Vert_V}{\Vert X\Vert_V}\} A=x=0max{XVAXV} ,可知 ∥ I ∥ = max ⁡ X ≠ 0 { ∥ I X ∣ V ∥ X ∥ V } = 1 \Vert I\Vert=\max_{X\neq 0}\limits \{\frac{\Vert IX\vert_V}{\Vert X\Vert_V}\}=1 I=X=0max{XVIXV}=1

  • 若某个 ∥ I ∥ ∗ > 1 \Vert I\Vert_*>1 I>1 ,则 ∥ ∙ ∥ ∗ \Vert \bullet\Vert_* 不是算子范数

eg

【矩阵论】7. 矩阵理论——算子范数

可见 ∥ A ∥ F \Vert A\Vert_F AF ∥ x ∥ 2 \Vert x\Vert_2 x2 是相容的,而 ∥ A ∥ 2 \Vert A\Vert_2 A2 作为 ∥ x ∥ 2 \Vert x\Vert_2 x2 的算子范数是相容的,但与 ∥ A ∥ F \Vert A\Vert_F AF 不同

  • ∥ A ∥ 2 ≤ ∥ A ∥ F \Vert A\Vert_2\le \Vert A\Vert_F A2AF ,即 ∥ X ∥ V \Vert X\Vert_V XV 的算子范数是使 ∥ A ∥ ≤ ∥ A X ∥ V ∥ X ∥ V \Vert A\Vert\le \frac{\Vert AX\Vert_V}{\Vert X\Vert_V} AXVAXV 成立的最小范数

7.2.3 算子范数判断矩阵级数收敛

a. 矩阵级数收敛定义

∑ k = 0 ∞ A k = A 0 + A 1 + ⋯ + A k + ⋯ 收敛于 A    ⟺    lim ⁡ k → ∞ ( A 0 + A 1 + ⋯ + A k ) = A 其中 A k ∈ C n , n , 记为 ∑ k = 0 ∞ A k = A 0 + A 1 + ⋯ + A k = A \begin{aligned} &\sum_{k=0}^{\infty}A_k=A_0+A_1+\cdots+A_k+\cdots收敛于A\iff \lim_{k\rightarrow \infty}(A_0+A_1+\cdots+A_k)=A\\ &其中A_k\in C^{n,n},记为 \sum_{k=0}^\infty A_k=A_0+A_1+\cdots+A_k=A \end{aligned} k=0Ak=A0+A1++Ak+收敛于Aklim(A0+A1++Ak)=A其中AkCn,n,记为k=0Ak=A0+A1++Ak=A

绝对收敛本身必收敛 : 若 ∑ k = 0 ∞ ∥ A k ∥ = ∥ A 0 ∥ + ∥ A 1 ∥ + ⋯ + ∥ A k ∥ + ⋯ \sum_{k=0}\limits^\infty \Vert A_k\Vert=\Vert A_0\Vert+\Vert A_1\Vert+\cdots+\Vert A_k\Vert+\cdots k=0Ak=A0+A1++Ak+ 收敛,则 ∑ k = 0 ∞ A k \sum_{k=0}\limits^\infty A_k k=0Ak 收敛

  • 若某个范数 ∥ A ∥ < 1 \Vert A\Vert<1 A<1 ,则 ∑ k = 0 ∞ A k \sum_{k=0}\limits^{\infty}A^k k=0Ak 绝对收敛

    证明:
    ∥ A ∥ < 1 ⇒ ∥ A k ∥ ≤ ∥ A ∥ k 且 ∑ k = 0 ∞ ∥ A ∥ k = ∥ I ∥ I − ∥ A ∥ ( 绝对收敛 ) \begin{aligned} &\Vert A\Vert<1\Rightarrow \Vert A^k\Vert\le \Vert A\Vert^k且 \sum_{k=0}^{\infty}\Vert A\Vert^k=\frac{\Vert I\Vert}{I-\Vert A\Vert}(绝对收敛) \end{aligned} A<1AkAkk=0Ak=IAI(绝对收敛)

  • ρ ( A ) < 1 ⇒ ∑ k = 0 ∞ A k \rho(A)<1\Rightarrow\sum_{k=0}\limits^{\infty}A^k ρ(A)<1k=0Ak 绝对收敛
    ρ ( A ) < 1 ⇒ 某范数 ∥ A ∥ < 1 ⇒ ∑ k = 0 ∞ A k 绝对收敛 \rho(A)<1\Rightarrow 某范数 \Vert A\Vert<1\Rightarrow \sum_{k=0}^{\infty}A^k绝对收敛 ρ(A)<1某范数A<1k=0Ak绝对收敛

ρ ( A ) < 1 \rho(A)<1 ρ(A)<1 ,则 ∑ k = 0 ∞ A k = ( I − A ) − 1 \sum_{k=0}^{\infty}\limits A^k=(I-A)^{-1} k=0Ak=(IA)1

证明:
( I − A ) ( I + A + ⋯ + A k + ⋯   ) = I ( I + A + ⋯ + A k + ⋯   ) − A ( I + A + ⋯ + A k + ⋯   ) = ( I + A + ⋯ + A k + ⋯   ) − ( A + ⋯ + A k + ⋯   ) = I ∴ ( I − A ) − 1 = ∑ k = 0 ∞ A k = I + A + ⋯ + A k + ⋯ \begin{aligned} &(I-A)(I+A+\cdots+A^k+\cdots)=I(I+A+\cdots+A^k+\cdots)-A(I+A+\cdots+A^k+\cdots)\\ &=(I+A+\cdots+A^k+\cdots)-(A+\cdots+A^k+\cdots)=I\\ &\therefore (I-A)^{-1}=\sum_{k=0}^\infty A^k=I+A+\cdots+A^k+\cdots \end{aligned} (IA)(I+A++Ak+)=I(I+A++Ak+)A(I+A++Ak+)=(I+A++Ak+)(A++Ak+)=I(IA)1=k=0Ak=I+A++Ak+

  • ∥ A ∥ < 1 \Vert A\Vert<1 A<1 ,则 ∑ k = 0 ∞ A k = ( I − A ) − 1 \sum_{k=0}^{\infty}\limits A^k=(I-A)^{-1} k=0Ak=(IA)1
b. 幂级数收敛

设幂级数 ∑ k = 0 ∞ c k x k \sum_{k=0}^\infty\limits c_kx^k k=0ckxk 的收敛半径为 R R R ,则

  • ρ ( A ) < R ⇒ ∑ k = 0 ∞ c k A k \rho(A)<R\Rightarrow \sum_{k=0}^\infty\limits c_kA^k ρ(A)<Rk=0ckAk 绝对收敛
  • 某一范数 ∥ A ∥ < R ⇒ ∑ k = 0 ∞ c k A k \Vert A\Vert<R\Rightarrow \sum_{k=0}^\infty\limits c_kA^k A<Rk=0ckAk 绝对收敛
  • ρ ( A ) > R \rho(A)>R ρ(A)>R ,则 ∑ k = 0 ∞ c k A k \sum_{k=0}^\infty\limits c_kA^k k=0ckAk 发散
  • ρ ( A ) = R \rho(A)=R ρ(A)=R ∑ k = 0 ∞ c k A k \sum_{k=0}^\infty\limits c_kA^k k=0ckAk 都有可能

SP:

  • e x = ∑ k = 0 ∞ 1 k ! x k e^x=\sum_{k=0}^\infty\limits \frac{1}{k!}x^k ex=k=0k!1xk

    s i n x = ∑ k = 0 ∞ ( − 1 ) k x 2 k + 1 ( 2 k + 1 ) ! sinx=\sum_{k=0}^\infty\limits (-1)^k\frac{x^{2k+1}}{(2k+1)!} sinx=k=0(1)k(2k+1)!x2k+1

    c o s x = ∑ k = 0 ∞ ( − 1 ) k x 2 k ( 2 k ) ! cosx=\sum_{k=0}^\infty\limits (-1)^k\frac{x^{2k}}{(2k)!} cosx=k=0(1)k(2k)!x2k 的收敛半径都是 R = + ∞ R=+\infty R=+ ,故任意方阵都满足收敛条件 ρ ( A ) < R \rho(A)<R ρ(A)<R

e A = ∑ k = 0 ∞ 1 k ! A k e^A=\sum_{k=0}^\infty\limits \frac{1}{k!}A^k eA=k=0k!1Ak s i n A = ∑ k = 0 ∞ ( − 1 ) k A 2 k + 1 ( 2 k + 1 ) ! sinA=\sum_{k=0}^\infty\limits (-1)^k\frac{A^{2k+1}}{(2k+1)!} sinA=k=0(1)k(2k+1)!A2k+1 c o s A = ∑ k = 0 ∞ ( − 1 ) k A 2 k ( 2 k ) ! cosA=\sum_{k=0}^\infty\limits (-1)^k\frac{A^{2k}}{(2k)!} cosA=k=0(1)k(2k)!A2k 都绝对收敛


  • l n ( 1 + x ) = ∑ k = 1 ∞ ( − 1 ) k − 1 x k k ln(1+x)=\sum_{k=1}^\infty\limits (-1)^{k-1}\frac{x^k}{k} ln(1+x)=k=1(1)k1kxk 的收敛半径为 r = 1 r=1 r=1

    ( 1 − x ) − 1 = 1 + x + ⋯ + x ⋯ k (1-x)^{-1}=1+x+\cdots+x^k_\cdots (1x)1=1+x++xk 的收敛半径 r = 1 r=1 r=1

    ( 1 − x ) − 2 = 1 + 2 x + 3 x 2 + ⋯ + k x k − 1 + ⋯ (1-x)^{-2}=1+2x+3x^2+\cdots+kx^{k-1}+\cdots (1x)2=1+2x+3x2++kxk1+ 的收敛半径 r = 1 r=1 r=1

l n ( I + A ) = ∑ k = 1 ∞ ( − 1 ) k − 1 A k k ln(I+A)=\sum_{k=1}^\infty\limits (-1)^{k-1}\frac{A^k}{k} ln(I+A)=k=1(1)k1kAk 的绝对收敛条件为 ρ ( A ) < 1 \rho(A)<1 ρ(A)<1 ∥ A ∥ < 1 \Vert A\Vert<1 A<1

( I − A ) − 1 = ∑ k = 0 ∞ A k (I-A)^{-1}=\sum_{k=0}^\infty\limits A^k (IA)1=k=0Ak 的绝对收敛条件为 ρ ( A ) < 1 \rho(A)<1 ρ(A)<1 ∥ A ∥ < 1 \Vert A\Vert<1 A<1

( I − A ) − 2 = ∑ k = 1 ∞ k A k − 1 (I-A)^{-2}=\sum_{k=1}^\infty\limits kA^{k-1} (IA)2=k=1kAk1 的绝对收敛条件是 ρ ( A ) < 1 \rho(A)<1 ρ(A)<1 ∥ A ∥ < 1 \Vert A\Vert<1 A<1文章来源地址https://www.toymoban.com/news/detail-424608.html

到了这里,关于【矩阵论】7. 矩阵理论——算子范数的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 【数值计算方法(黄明游)】解线性代数方程组的迭代法(一):向量、矩阵范数与谱半径【理论到程序】

       注意:速读可直接跳转至“4、知识点总结”及“5、计算例题”部分   当涉及到线性代数和矩阵理论时, 向量、矩阵范数以及谱半径 是非常重要的概念,下面将详细介绍这些内容: a. 定义及性质   考虑一个 n n n 维向量 x x x ,定义一个实值函数 N ( x ) N(x) N ( x ) ,

    2024年01月25日
    浏览(37)
  • 计算方法 | 范数(向量:1范数、2范数、无穷范数;矩阵:行范数、列范数)

    0范数: 向量中非零元素的个数。 1范数: 为绝对值之和。 2范数: 通常意义上的模。 无穷范数:取向量的最大值。 行范数:矩阵中每行绝对值之和的最大值 列范数:矩阵中每列绝对值之和的最大值 范数对于数学的意义?1范数、2范数、无穷范数_yangpan011的博客-CSDN博客_无穷

    2024年02月15日
    浏览(28)
  • 矩阵分析与计算学习记录-向量范数与矩阵范数

    本章知识重点: 向量范数:定义、性质、等价性、分析性质 矩阵范数:定义、算子范数 矩阵范数与向量范数的相容性 矩阵的普半径及应用:普半径、矩阵序列及级数中的应用 矩阵的条件数及应用:矩阵的条件数、误差估计中的应用 范数,是具有“长度”概念的函数。 在线

    2024年02月05日
    浏览(37)
  • L1范数,L2范数,L2,1范数(向量范数、矩阵范数、正则化)

    参考文章如下:https://blog.csdn.net/lqzdreamer/article/details/79676305                          https://blog.csdn.net/lqzdreamer/article/details/79676305         一般常用范数来衡量向量,向量的Lp范数定义为:          Lp范数示意图:         从图中可以看出,p的取值在 [0,1) 之间,范数

    2023年04月09日
    浏览(27)
  • 用python实现矩阵和向量的范数(包括一范数,二范数,无穷范数)

    首先,导入需要用到的库: 创建一个矩阵和一个向量并输出: 计算矩阵的第一范数: 计算矩阵的第二范数: 计算矩阵的无穷范数: 计算向量的第一范数: 计算向量的第二范数: 计算向量的无穷范数: 完整代码如下:

    2024年04月22日
    浏览(24)
  • 证明矩阵二范数的平方等于转置矩阵与原矩阵相乘后的二范数

    定理: 对于任意的矩阵 A ∈ R n × m A in R^{ntimes m} A ∈ R n × m ,有 ∥ A ∥ 2 2 = ∥ A T A ∥ 2 left|Aright|_2^2=left|A^TAright|_2 ∥ A ∥ 2 2 ​ = ​ A T A ​ 2 ​ 证明: 假设矩阵 A T A A^TA A T A 最大特征值为 λ lambda λ ,即 ∥ A ∥ 2 2 = λ left|Aright|_2^2=lambda ∥ A ∥ 2 2 ​ = λ ,设

    2024年02月05日
    浏览(29)
  • 证明转置矩阵的二范数与原矩阵的二范数相等

    定理: 对于任意的矩阵 A ∈ R n × m A in R^{ntimes m} A ∈ R n × m ,有 ∥ A ⊤ ∥ 2 = ∥ A ∥ 2 left|A^topright|_2=left|Aright|_2 ​ A ⊤ ​ 2 ​ = ∥ A ∥ 2 ​ 证明: 根据矩阵二范数的定义, ∥ A ∥ 2 = λ m a x ( A ⊤ A ) left|Aright|_2=sqrt{lambda_{max}(A^top A)} ∥ A ∥ 2 ​ = λ ma x ​ (

    2024年02月06日
    浏览(25)
  • <3>【深度学习 × PyTorch】必会 线性代数 (含详细分析):点积 | 矩阵-向量积 | Hadamard积 | 矩阵乘法 | 范数/矩阵范数

      拍照的意义在于你按下快门的那一刻,万里山河的一瞬间变成了永恒。   🎯作者主页: 追光者♂🔥          🌸个人简介:   💖[1] 计算机专业硕士研究生💖   🌟[2] 2022年度博客之星人工智能领域TOP4🌟   🏅[3] 阿里云社区特邀专家博主🏅   🏆[4] CSDN-人工智能领域

    2024年02月05日
    浏览(41)
  • MATLAB矩阵范数最小值

    % 现在有一个矩阵B是2行×100列, % 想求该矩阵每一列×该列转置后形成的矩阵与一矩阵A(2行×2列)之差的范数, % 100列乘完后,定位出上述100个结果中范数最小值,定位出最小值所在B中哪列, % 及对应B中对应哪个元素 代码如下:

    2023年04月17日
    浏览(22)
  • 3.3 向量与矩阵的范数

      要学习向量与矩阵的范数,我会采取以下几个步骤: 了解基本概念:首先,我会了解向量和矩阵的范数的基本概念和定义,以及它们的性质和特点,这是理解和掌握范数的基础。 学习具体算法:其次,我会学习具体的算法和计算方法,如计算向量的L1、L2、无穷范数,计算

    2023年04月22日
    浏览(23)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包