【MATLAB】几种特殊矩阵,Hilbert矩阵,Toeplitz矩阵,Vandermonde矩阵......

news/2024/7/5 3:34:30

1、Hadamard矩阵

Hadamard矩阵是由+1和-1元素构成的且满足Hn*Hn’=nI(这里Hn’为Hn的转置,I为单位方阵)n阶方阵。

>> hadamard(4)ans =1     1     1     11    -1     1    -11     1    -1    -11    -1    -1     1

2、Hankel矩阵:是指每一条副对角线上的元素都相等的方阵。

>> x=[1 2 3 4 5 6 7]x =1     2     3     4     5     6     7>> hankel(x)ans =1     2     3     4     5     6     72     3     4     5     6     7     03     4     5     6     7     0     04     5     6     7     0     0     05     6     7     0     0     0     06     7     0     0     0     0     07     0     0     0     0     0     0

hankel(x,y):

返回一个m x n的Hankel矩阵,它的第一列向量为x,最后一行为向量y

>> x=[1 2 3 4 5 6 7]x =1     2     3     4     5     6     7>> y=[7 6 5 4 3 2 2]y =7     6     5     4     3     2     2>> hankel(x,y)ans =1     2     3     4     5     6     72     3     4     5     6     7     63     4     5     6     7     6     54     5     6     7     6     5     45     6     7     6     5     4     36     7     6     5     4     3     27     6     5     4     3     2     2

3、Pascal矩阵:

Pascal矩阵的第一行元素和第一列元素都为1,其余位置处的元素是该元素的左边元素加起上一行对应位置相加而得,如元素Ai,j=Ai,j-1+Ai-1,j。Ai,j表示第i行,第j列位置上的元素

>> pascal(4)ans =1     1     1     11     2     3     41     3     6    101     4    10    20>> pascal(5)ans =1     1     1     1     11     2     3     4     51     3     6    10    151     4    10    20    351     5    15    35    70

4、Rosser矩阵

>> rosserans =611   196  -192   407    -8   -52   -49    29196   899   113  -192   -71   -43    -8   -44-192   113   899   196    61    49     8    52407  -192   196   611     8    44    59   -23-8   -71    61     8   411  -599   208   208-52   -43    49    44  -599   411   208   208-49    -8     8    59   208   208    99  -91129   -44    52   -23   208   208  -911    99

5、Wilkinson矩阵:Wilkinson特征值测试矩阵

>> wilkinson(3)ans =1     1     01     0     10     1     1>> wilkinson(5)ans =2     1     0     0     01     1     1     0     00     1     0     1     00     0     1     1     10     0     0     1     2

6、Hilbert矩阵

也成为H阵,其元素,由于它是一个条件数差的矩阵,所以将它用来作为试验矩阵。

hilb(n):用于生成n x n的Hilbert矩阵

invhilb(n):用于生成n x n的Hilbert矩阵的逆矩阵整数矩阵。

>> a=hilb(3)a =1.0000    0.5000    0.33330.5000    0.3333    0.25000.3333    0.2500    0.2000>> b=invhilb(3)b =9   -36    30-36   192  -18030  -180   180

可以看出Hilbert矩阵和他的逆矩阵都是对称矩阵。

7、Toeplitz矩阵

他由两个向量定义,一个行向量和一个列向量。对称的Toeplitz矩阵由单一向量来定义。

toeplizt(k,r):用于生成非对称Toeplitz矩阵,第一列为k,第一行为r,其余元素等于其左上角元素

toeplitz(c):向量c生成一个对称的Toeplitz矩阵

>> toeplitz(2:4,2:2:8)ans =2     4     6     83     2     4     64     3     2     4>> toeplitz(2:8)ans =2     3     4     5     6     7     83     2     3     4     5     6     74     3     2     3     4     5     65     4     3     2     3     4     56     5     4     3     2     3     47     6     5     4     3     2     38     7     6     5     4     3     2

8、0~1间均匀分布的随机矩阵

r=rand(n):产生nxn的0~1间均匀分布的随机矩阵

r=rand(m,n):产生mxn的0~1间均匀分布的随机矩阵,r=rand(m,n,p.....):与前者类似、

r=rand(size(A)):产生与A相同的0~1间均匀分布矩阵

>> r=rand(3)r =0.7094    0.6797    0.11900.7547    0.6551    0.49840.2760    0.1626    0.9597>> r=rand(3,4)r =0.3404    0.7513    0.6991    0.54720.5853    0.2551    0.8909    0.13860.2238    0.5060    0.9593    0.1493>> a=[1 2;3 4]a =1     23     4>> r=rand(size(a))r =0.2575    0.25430.8407    0.8143

9、标准正态分布随机矩阵

常用randn()函数产生均值为零、方差为1的随机矩阵

调用方式如下:

r=rand(n)

r=rand(m,n)

r=rand(size(A))

r=randn([m,n,......])

>> b=randn(3)b =-0.1924   -1.4023   -0.17740.8886   -1.4224   -0.1961-0.7648    0.4882    1.4193>> b=randn(3,4)b =0.2916   -0.8045   -0.2437   -1.14800.1978    0.6966    0.2157    0.10491.5877    0.8351   -1.1658    0.7223>> b=randn([3,4])b =2.5855   -0.0825   -1.7947    0.1001-0.6669   -1.9330    0.8404   -0.54450.1873   -0.4390   -0.8880    0.3035

10、Vandermonde矩阵

A=vander(V):生成Vandermonde矩阵,矩阵的列是向量v的幂,

>> vander([1 2 3 4])ans =1     1     1     18     4     2     127     9     3     164    16     4     1>> vander([1 ;2 ;3 ;4])ans =1     1     1     18     4     2     127     9     3     164    16     4     1

11、魔方矩阵:每行,每列,两对角线上的元素和相等

>> a=magic(5)a =17    24     1     8    1523     5     7    14    164     6    13    20    2210    12    19    21     311    18    25     2     9


http://lihuaxi.xjx100.cn/news/252661.html

相关文章

微信自研NLP大规模语言模型WeLM:零/少样本即可完成多种NLP任务

近日,微信AI推出自研NLP大规模语言模型WeLM ,该模型是一个尺寸合理的中文模型,能够在零样本以及少样本的情境下完成包多语言任务在内的多种NLP任务。同时,微信AI团队也提供了WeLM的体验网页和API接口,感兴趣的用户可前…

[转帖]什么是光纤的波长?看看有哪些是你不知道的!

什么是光纤的波长?看看有哪些是你不知道的! FShttps://www.feisu.com/bbs/e-1640.html2017-07-01 00:00:001084我们平时最熟悉的光当然是我们肉眼所能看见的光。我们的眼睛对波长在400nm的紫光到700nm的红光很敏 感。但对于携带玻璃纤维的光纤&#xff0…

【MATLAB】矩阵分析之向量和矩阵的范数运算

本片借鉴于 https://blog.csdn.net/u013534498/article/details/52674008 https://blog.csdn.net/left_la/article/details/9159949 向量范数当p1时,即为各个向量的元素绝对值之和 >> norm(x,1)ans 21>> xx 1 2 3 4 5 6>> no…

SQL Server 2008中的Pivot和UnPivot

SQL Server 2008中SQL应用系列--目录索引 今天给新成员讲解PIVOT 和 UNPIVOT示例,顺便整理了一下其用法。这是自SQL Server 2005起提供的新功能。 官方示例:http://msdn.microsoft.com/zh-cn/library/ms177410%28vsql.105%29.aspx 首先看PIVOT示例&#…

C语言存储类关键字

1、static这个关键字有三种用法:(1)第一种是用来修饰局部变量,使之成为静态局部变量;静态局部变量存储在数据段/bss段中,作用域是代码块作用域,生命周期是程序生命周期,链接属性是无…

BZOJ 1096: [ZJOI2007]仓库建设

传送门 斜率优化DP入门题 显然如果在一个位置 i 建一个仓库,且上一个仓库位置为 j 那么从 j1到 i 的物品显然都要存在 i 仓库是最优的 设 $f [ i ]$ 表示在第 i 个工厂建设仓库时,工厂 1 到 i 的物品都转移好的最小花费 考虑上一个仓库的位置 j 设工厂 i…

成为探路者,成就探路者!亚马逊云科技中国峰会精彩回顾

点击上方入口立即【自由构建 探索无限】一起共赴年度科技盛宴!点击阅读原文进入官方小程序观看主题演讲精彩回放前沿技术,大咖云集更多精彩不容错过

【组队学习】【28期】数据采集从入门到精通

数据采集从入门到精通 论坛版块: http://datawhale.club/c/team-learning/38-category/38 开源内容: https://github.com/datawhalechina/team-learning-program/tree/master/CollectData 学习目标 随着数字化的不断推进,数据采集在数据…