线性代数Python计算：正交向量组的扩张

戌崂石

已于 2024-05-21 18:46:08 修改

阅读量687

点赞数 1

分类专栏：线性代数文章标签：线性代数 python 开发语言

于 2022-06-16 15:38:22 首次发布

本文链接：https://blog.csdn.net/u012958850/article/details/125316905

版权

线性代数专栏收录该内容

30 篇文章 11 订阅

订阅专栏

本文通过Python展示如何利用已知正交向量组扩展欧氏空间的正交基，以向量(1,1,1)和(1,-2,1)为例，使用mySolve函数求解线性方程组。关键步骤包括构造矩阵A和零向量o，解方程并获取非零解a3，最终构建新的正交基。

摘要由CSDN通过智能技术生成

在这里插入图片描述
已知 $\boldsymbol{\alpha}_1=\begin{pmatrix}a_{11}\\a_{12}\\\vdots\\a_{1n}\end{pmatrix},\boldsymbol{\alpha}_2=\begin{pmatrix}a_{21}\\a_{22}\\\vdots\\a_{2n}\end{pmatrix},\cdots,\boldsymbol{\alpha}_k=\begin{pmatrix}a_{k1}\\a_{k2}\\\vdots\\a_{kn}\end{pmatrix}\in$ ℝ $^n$ 两两正交（ $k < n$ ），设矩阵
$\boldsymbol{A}=\begin{pmatrix}a_{11}&a_{12}&\cdots&a_{1n}\\a_{21}&a_{22}&\cdots&a_{2n}\\\vdots&\vdots&\ddots&\vdots\\a_{k1}&a_{k2}&\cdots&a_{kn}\end{pmatrix}$
解方程组 $\boldsymbol{Ax}=\boldsymbol{o}$ ，取一非零解即可将两两正交组扩展为 $\boldsymbol{\alpha}_1,\cdots,\boldsymbol{\alpha}_k,\boldsymbol{\alpha}_{k+1}$ 。
例1 用Python将向量组 $\boldsymbol{\alpha}_1=\begin{pmatrix}1\\1\\1\end{pmatrix},\boldsymbol{\alpha}_2=\begin{pmatrix}1\\-2\\1\end{pmatrix}$ 扩展成欧氏空间ℝ $^3$ 的在一组正交基。
解：这只需调用博文《线性方程组的通解》中定义的mySolve函数解方程组 $\boldsymbol{Ax}=\boldsymbol{o}$ ，其中 $\boldsymbol{A}=\begin{pmatrix}1&1\\1&-2\\1&1\end{pmatrix}$ ， $\boldsymbol{o}=\begin{pmatrix}0\\0\end{pmatrix}$ 。

import numpy as np              #导入numpy
a1=np.array([1,1,1])            #向量a1
a2=np.array([1,-2,1])           #向量a2
A=np.vstack((a1,a2))            #构建矩阵A
o=np.zeros((2,1))               #零向量
X=mySolve(A,o)                  #解方程组
a3=X[:,1]                       #非零解
print(a3)

程序的第2、3行设置已知正交向量a1和a2，第4行用a1，a2作为行向量构造矩阵A，第5行调用numpy的zeros函数构造零向量o。第6行调用博文《线性方程组的通解》定义的mySolve函数，计算解集X。第7行将X中的第2列（下标为1），即基础解系中的第1个非零解置为a3。运行程序，输出

[-1.  0.  1.]

即取 $\boldsymbol{\alpha}_3=\begin{pmatrix}-1\\0\\1\end{pmatrix}$ ， $\boldsymbol{\alpha}_1,\boldsymbol{\alpha}_2,\boldsymbol{\alpha}_3$ 构成欧氏空间ℝ $^3$ 的一个正交基。
写博不易，敬请支持：
如果阅读本文于您有所获，敬请点赞、评论、收藏，谢谢大家的支持！
代码诚可贵，原理价更高。若为AI学，读正版书好。