python应用——用python求范数

应用统计的作业实在是手算不出来，于是顺应一下时代需求用python解题，稍微梳理一下用python求范数的内容。

先看一下诱导范数的定义：

由此可以看到，对于一般的矩阵，列和范数和行和范数的求解比较简单，主要是算子范数/谱范数的求解需要先求特征值较为复杂。下面以一道简单的题由此梳理了一下求第二范数的三种思路。

1.直接用numpy的linalg.norm函数求范数

用np.linalg.norm()函数可以直接求解

a=[[1,2,0],[-1,2,-1],[0,1,1]]#写入矩阵
A=np.array(a)#转换成np.array格式
a1=np.linalg.norm(A,ord=1)
a2=np.linalg.norm(A,ord=2)
a3=np.linalg.norm(A,ord=np.Inf)  #分别计算范数
print("第一范数为",a1)
print("第二范数为",round(a2,2))
print("无穷范数为",a3)

得到结果如下图：

2.算得特征多项式对一元三次方程求解得特征值，对最大特征值求平方根得第二范数

计算出转置矩阵与原矩阵乘积后得特征多项式，行列式计算得到一元三次方程求解特征值

#法二：第二范数：先求最大特征值再求#先得到A的转置矩阵与A相乘的矩阵d
b=np.transpose(A)
B=np.array(b)
d=np.matmul(B,A)#根据特征多项式得一元三次方程求解
import sympy as sp # 导入sympy包
x=sp.Symbol('x')
f=x**3-13*(x**2)+38*x-25   #见前面运算过程
x=sp.solve(f)  #得到的解析解比较复杂，故后续转换成浮点数#取四位小数输出
for i in range(0,3):x[i]=round(x[i].evalf(),4) #求出表达式的浮点数print(x[i])
maxX=x[0]
for i in x:   ## 求最大值if i > maxX:maxX= i
print("第二范数为",round(maxX**0.5,4))

得到结果如下图：

3.直接求解特征值与特征向量，对最大特征值求平方根得第二范数

用np.linalg.eig()函数求解特征值和特征向量之后求解。

#直接求解特征值与特征向量
w,v=np.linalg.eig(d)
print('特征值：{}\n特征向量：{}'.format(w,v))
#求最大特征值的算数平方根
W=w.tolist()
maxW= W[0]
for i in W:   ## 求最大值if i > maxW:maxW= i
print("第二范数为",round(maxW**0.5,4))

得到结果如下图：

python应用——用python求范数相关推荐

【LKH算法体验】Python调用LKH算法求TSP问题
[LKH算法体验]Python调用LKH算法求TSP问题一.LKH算法简介 Keld Helsgaun 是丹麦罗斯基勒大学计算机科学专业的名誉副教授. 他于 1973 年在哥本哈根大学获得DIK ...
python中Scipy模块求取积分
python中Scipy模块求取积分的方法: SciPy下实现求函数的积分的函数的基本使用,积分,高等数学里有大量的讲述,基本意思就是求曲线下面积之和. 其中rn可认为是偏差,一般可以忽略不计,wi可 ...
Algorithm：C++/python语言实现之求旋转数组最小值、求零子数组、求最长公共子序列和最长公共子串、求LCS与字符串编辑距离
Algorithm:C++/python语言实现之求旋转数组最小值.求零子数组.求最长公共子序列和最长公共子串.求LCS与字符串编辑距离目录一.求旋转数组最小值 1.分析问题 2.解决思路二.求 ...
python高阶函数求导_Python---高阶函数
# 高阶函数 # 变量可以指向函数 # 以Python内置的求绝对值的函数abs()为例,调用该函数用以下代码 print(abs(-10)) # 只写abs # 输出为: print(abs) # ...
python用for循环求10的因数_python for循环练习（初级）
for循环练习1 for i in range(4):print(i) D:\尚硅谷Python\venv\Scripts\python.exe D:/尚硅谷Python/test.py 0 1 2 ...
廖雪峰讲python高阶函数求导公式_高阶函数 - 廖雪峰 Python 2.7 中文教程
高阶函数英文叫Higher-order function.什么是高阶函数?我们以实际代码为例子,一步一步深入概念. 变量可以指向函数以Python内置的求绝对值的函数abs()为例,调用该函数用以下 ...
python numpy 多项式函数求导求根
python numpy 多项式函数求导求根 """求出多项式的导函数与根 """import numpy as np import m ...
Python实现相空间重构求关联维数——GP算法、自相关法求时间延迟tau、最近邻算法求嵌入维数m
Python实现相空间重构求关联维数--GP算法.自相关法求时间延迟tau.最近邻算法求嵌入维数m GP算法: 若有一维时间序列为{x1,x2,-,xn},对其进行相空间重构得到高维相空间的一系列向量 ...
python内置库求复数的辐角_皮肤与美容—医学专家如是说_中国大学MOOC(慕课)_章节测验答案...
皮肤与美容-医学专家如是说_中国大学MOOC(慕课)_章节测验答案更多相关问题某国开展对外贸易,一般来说可以使为防止麻醉后引起呕吐或窒息,一般手术前禁食.禁饮的时间是 [单选题]Windows ...

python应用——用python求范数

python应用——用python求范数相关推荐

最新文章

热门文章