-
2021-04-28 03:12:41
标准正交基概述
2011~2012 学年春季学期课程论文
课程名称:线性代数与几何
论文题目: 探索欧氏空间的标准正交基
作者姓名: XXX 学 号:XXX
成 绩:
论文评语:
评阅人:
评阅日期:
注:后附课程论文的正文
探索欧氏空间的标准正交基
XXX XXX XXX
摘要:本文首先给出欧式空间标准正交基的定义以及性质,然后给出将欧氏空间
任意一组线性无关的向量组进行正交化和单位化的过程,最后给出两组标准正交
基之间的关系。
关键词:标准正交基 正交化 正交矩阵
正文:
一.引言:标准正交基对欧氏空间有特别的重要性.在欧氏空间中研究线性变换和二次型时,
都需要通过标准正交基之间的过渡来简化线性变换和二次型矩阵.因此研究欧氏空间的标准
正交基很有意义.
二.定义:
V V
定义1:设 是欧氏空间.如果, V 满足条件() 0 ,就称 正交,记为 .
V
中由两两正交的非零向量组成的向量组 , 称为正交向量组.如果 的基M 是正交向量
1 k
组,就称M 为正交基.由两两正交的单位向量组成的基称为标准正交基.
, V 1,i 1, n
定义2:I= 是欧氏空间 的一个基,且 ,又 ,
1 n , 0, i j i
i j
V
称I 为 的一个标准正交基.
三.性质:
V
性质1. , (1)为 的一个基,则以下条件等价:
1 n
1 (1)是标准正交基
1 i j
2 ,
i j
0 i j
3 V, , , ,
1 1 2 2 n n
n
4 , V, x , x , y , y , 则, x y
1 n 1 n i i
i 1
证明: 当 时, , 1,
1 2 i j i i i
当i j 时由正交向量的定义易知 , 0 ,命题成立.
更多相关内容 -
标准正交基应用于通信原理课程的教学研究
2020-08-08 05:50:04本文在完成了“用信息工具改造线性代数课程”的教学改革基础上,针对通信原理课程教学中存在的问题,以调制内容为例提出了将线性代数中的标准正交基概念引入和将MATIAB用于教学的解决方案。 -
matlab基础学习(五)之映射和正交基
2021-04-20 09:26:04matlab基础学习(5)之映射和正交基一、映射在matlab中,矢量y正交映射到矢量x上的定义如下:Px(y) = /(x^2)*x对于长度为N的列矢量y映射到长度为N的列矢量x上的计算方法如下:yx = (x' * y) * (x' * x) ^ (-1) * x更...matlab基础学习(5)之映射和正交基
一、映射
在matlab中,矢量y正交映射到矢量x上的定义如下:
Px(y) = /(x^2)*x
对于长度为N的列矢量y映射到长度为N的列矢量x上的计算方法如下:
yx = (x' * y) * (x' * x) ^ (-1) * x
更一般的, 长度为N的列矢量y映射到M维的子空间的 N x M的矩阵x,计算方法如下:
yX = X * (X' * X)^(-1) * X' * y
正交映射,像任何的有限维线性操作符,可通过矩阵表示,对于N x M矩阵,计算方法为:
PX = X * (X' * X)^(-1) * X'
PX被称为映射矩阵。
子空间映射是一个例子,其中矩阵的功率线性代数记号非常重要。
代码实现如下:
>> X = [[1;2;3],[1;0;1]]
X =
1 1
2 0
3 1
>> PX = X*(X'*X)^(-1) * X'
PX =
0.66666666666667 -0.33333333333333 0.33333333333333
-0.33333333333333 0.66666666666667 0.33333333333333
0.33333333333333 0.33333333333333 0.66666666666667
>> y = [2; 4;6]
y =
2
4
6
>> yX= PX* y
yX =
2.00000000000000
4.00000000000000
6.00000000000000
>>
二、正交基
在matlab中,函数orth()对于给定的矢量集空间计算正交基。
在matlab中输入help orth查看orth的使用:
>> help orth
ORTH Orthogonalization.
Q = ORTH(A) is an orthonormal basis for the range of A.
That is, Q'*Q = I, the columns of Q span the same space as
the columns of A, and the number of columns of Q is the
rank of A.
Class support for input A:
float: double, single
See also svd, rank, null.
Reference page in Help browser
doc orth
>>
以下通过对N=3的线性独立的基利用orth()进行正交化。
>> v1 = [1;2;3]
v1 =
1
2
3
>> v2 = [1;-2;3]
v2 =
1
-2
3
>> v1' * v2
ans =
6
>> V=[v1,v2]
V =
1 1
2 -2
3 3
>> W= orth(V)
W =
-0.31622776601684 0.00000000000000
0.00000000000000 -1.00000000000000
-0.94868329805051 0.00000000000000
>> W= orth(V)
W =
-0.31622776601684 0.00000000000000
0.00000000000000 -1.00000000000000
-0.94868329805051 0.00000000000000
>> W1=W(:,1)
W1 =
-0.31622776601684
0.00000000000000
-0.94868329805051
>> W2=W(:,2)
W2 =
0.00000000000000
-1.00000000000000
0.00000000000000
>> W1' * W2
ans =
-3.157438761082832e-016
>> W1'*W1
ans =
1.00000000000000
>> W2'*W2
ans =
1
>> W'*W
ans =
1.00000000000000 -0.00000000000000
-0.00000000000000 1.00000000000000
>> x = 2 * v1 - 3 * v2
x =
-1
10
-3
>> c1 = x' * W1
c1 =
3.16227766016838
>> c2 = x' * W2
c2 =
-10
>> xw = c1 * W1 + c2 * W2
xw =
-1.00000000000000
10.00000000000000
-3.00000000000000
>> error = x - xw
error =
1.0e-014 *
0.13322676295502
0
0.08881784197001
>> norm(error)
ans =
1.601186416994689e-015
>> x = [1; 0; 0]
x =
1
0
0
>> ca = x' * W1
ca =
-0.31622776601684
>> c1 = x' * W1
c1 =
-0.31622776601684
>> c2 = x' * W2
c2 =
1.387778780781446e-016
>> xw = c1 * W1 + c2 * W2
xw =
0.10000000000000
-0.00000000000000
0.30000000000000
>> error = x - xw
error =
0.90000000000000
0.00000000000000
-0.30000000000000
>> norm(error)
ans =
0.94868329805051
>> W* error
??? Error using ==> mtimes
Inner matrix dimensions must agree.
>> W'* error
ans =
1.0e-015 *
-0.16653345369377
-0.09984698057522
>>
-
稀疏零空间和正交:稀疏矩阵的零空间和正交基-matlab开发
2021-06-01 14:58:50使用带行置换的 QR 分解计算稀疏矩阵的 NULL 空间和 ORTHOGONAL 基的两个简单函数。 对于 FULL 矩阵,Matlab 库存函数 NULL 和 ORTH ...) 开始,QR 分解可用于稀疏矩阵,可用于估计正交基,而无需在 FULL 中转换矩阵。 -
关于标准正交基
2013-05-06 00:02:36主要介绍标准正交基的定义,性质,定理及定理的证明 -
matlab orth函数 求矩阵的正交基
2017-03-14 10:54:10返回矩阵A正交基。 B列与A列具有相同空间。 B列向量正交向量满足B'*B = eye(rank(A)) B的列数是A的秩。 举例: >> a2 a2 = 1 4 7 2 5 8 3 6 9 >> orth(a2) ans = -0.4797 0.7767 ...
B = orth(A)
返回矩阵A正交基。B列与A列具有相同空间。B列向量正交向量满足B'*B = eye(rank(A))B的列数是A的秩。
举例:>> a2
a2 =
1 4 7 2 5 8 3 6 9
>> orth(a2)
ans =
-0.4797 0.7767 -0.5724 0.0757 -0.6651 -0.6253
>> ans'*ans
ans =
1.0000 -0.0000 -0.0000 1.0000
>> -
Matlab构建标准正交小波基
2021-11-18 16:09:15matlab 构建小波基矩阵做压缩感知构建小波基矩阵时没有像构建cos余弦基那么方便,可以直接采用'dct(N)' 获得,笔者在不断的搜寻下汇总了3个方法(已知长度,展开水平,和小波类型生成小波基矩阵),供大家参考,
1. 压缩感知稀疏基之离散小波变换_彬彬有礼的专栏-CSDN博客
function [ ww ] = dwtmtx( N,wtype,wlev )
%DWTMTX Discrete wavelet transform matrix
% This function generates the transform matrix ww according to input
% parameters N,wtype,wlev .
%Detailed explanation goes here
% N is the dimension of ww
% wtype is the wavelet type
% wlev is the number of decomposition level
%NOTE: The extension mode must be Periodization('per')
[h,g]= wfilters(wtype,'d'); %Decomposition low&high pass filter
L=length(h); %Filter length
h_1 = fliplr(h); %Flip matrix left to right
g_1 = fliplr(g);
loop_max = log2(N);
loop_min = double(int8(log2(L)))+1;
if wlev>loop_max-loop_min+1
fprintf('\nWaring: wlev is too big\n');
fprintf('The biggest wlev is %d\n',loop_max-loop_min+1);
wlev = loop_max-loop_min+1;
end
ww=1;
for loop = loop_max-wlev+1:loop_max
Nii = 2^loop;
p1_0 = [h_1 zeros(1,Nii-L)];
p2_0 = [g_1 zeros(1,Nii-L)];
p1 = zeros(Nii/2,Nii);
p2 = zeros(Nii/2,Nii);
for ii=1:Nii/2
p1(ii,:)=circshift(p1_0',2*(ii-1)+1-(L-1)+L/2-1)';
p2(ii,:)=circshift(p2_0',2*(ii-1)+1-(L-1)+L/2-1)';
end
w1=[p1;p2];
mm=2^loop_max-length(w1);
w=[w1,zeros(length(w1),mm);zeros(mm,length(w1)),eye(mm,mm)];
ww=ww*w;
clear p1;clear p2;
end%The end!!!
end虽然看不太懂,但是使用的时候发现这个与我看的一篇文献不太符合
2 https://sparse-lex.readthedocs.io/en/latest/book/wavelets/wavelet_toolbox/daubechies.html
结合以上网址的介绍写了这个
function PsiT = dwtmtx_R2( N,dbtype,wlev )
[LoD,HiD,LoR,HiR] = wfilters(dbtype);
old_dwt_mode = dwtmode('status','nodisp');
dwtmode('per');
PsiT = zeros(N, N);
for i=1:N
unit_vec = zeros(N, 1);
unit_vec(i) = 1;
[coefficients, levels] = wavedec(unit_vec, wlev, LoD,HiD);
PsiT(:, i) = coefficients;
end
end这个的计算结果使我可以重现文献中的结果
3.wavelab 850 包的使用
How to get started with wavelets and Wavelab - CCM
n=256;
qmf= MakeONFilter('Daubechies',4) ;
XI=eye(n);W=zeros(n);for i=1:n,W(:,i)=FWT_PO(XI(:,i),1,qmf);end
反正这个结果不是我想要的
本人并不懂其中原理,目的只是想要构建压缩感知的基函数矩阵,以上如有疏漏还请见谅
-
标准正交基应用于通信原理课程的教学研究
2021-01-30 21:43:24基于整合教学内容的目的,本文探讨了将众多调制方式(包括模拟调制和数字调制)统一为正交调制模型的教学方案。为了在满足教学要求的前提下降低教学难度,提出了用MATLAB进行仿真演示的教学方法。教学实践证明,该... -
Matlab使用-orth函数(正交基)
2019-06-20 15:31:53Test B = orth(A) ...返回矩阵A正交基。 B列与A列具有相同空间。 B列向量正交向量满足B'*B = eye(rank(A)) B的列数是A的秩。 举例: >> a2 a2 = 1 ... -
(标准)正交基,正交投影,正交分解定理,最佳逼近定理,格拉姆-施密特方法求正交基(手算+MATLAB),...
2021-02-21 11:23:54正交(orthogonality)、正交集(orthogonal set)、单位正交集 定义:如果向量uuu和向量vvv是相互正交的,则有:u⋅v=0u\cdot v=0u⋅v=0。 因为对于零向量,都有0T⋅v=00^T\cdot v=00T⋅v=0,所以零向量与RnR^nRn中... -
matlab 判断信号正交
2020-03-24 21:47:55>> dot(y,yi) ans = 2.4980e-15 >> sum(y.*yi) ans = 2.4980e-15 >> dot([1,0],[0,1]) ans = 0 -
关于PCA的理解
2021-04-21 22:37:04最普通的一组正交基是标准正交基,实验采样的结果通常可以看作是在标准正交基下表示的。因为一般的观测者都是习惯于取摄像机的屏幕坐标,即向上和向右的方向作为观测的基准。也就是说,标准正交基表现了数据观测的... -
2001_10小波变换及其MATLAB工具的应用_10505158.pdf
2021-04-24 10:50:552001_10小波变换及其MATLAB工具的应用_10505158[General Information]书名=2001.10小波变换及其MATLAB工具的应用作者=郑治真等编著页数=156SS号出版日期=2001年10月第1版封面页书名页版权页前言页目录页第一部分 ... -
cgrscho:经典的 Gram-Schmidt 正交化程序。-matlab开发
2021-05-31 03:04:00这个过程从任何一组 N 个线性无关的向量构建一个标准正交基。 显然,通过跳过归一化步骤,我们也可以简单地形成一个正交基。 这个过程的关键要素是每个新的标准正交基向量是通过减去下一个线性无关向量到集合中到... -
ortha.m:相对于矩阵 A 的正交归一化-matlab开发
2021-05-30 06:37:53ORTHA 正交归一化相对于矩阵 A Q=ortha(A,X) Q=ortha('Afunc',... 计算 X 范围的标准正交基 Q,相对于使用正定和自伴随矩阵 A 的标量积。 也就是说,Q'*A*Q = I,Q 的列与 X 的列跨越相同的空间,并且秩(Q)=秩(X)。 -
正交表的构造方法及Matlab实现
2021-04-30 05:18:28* 基金项目: 国家自然科学基金( 11571094) 正交表的构造方法及 Matlab 实现* 河南师范大学数学与信息科学学院( 453007) 庞善起 鹿姗姗 正交试验设计是研究多因素水平的试验设计方法,它是一种高效率、快速、经济的... -
201208-PyTorch求解矩阵正交基
2020-12-09 06:16:20参考 MATLAB Link orth is obtained from U in the singular value decomposition, [U,S] = svd(A,‘econ’). If r = rank(A), the first r columns of U form an orthonormal basis for the range of A. 步骤... -
在线计算专题(12):矩阵的特征值、特征向量、正交变换与二次型与常见矩阵分解
2021-04-24 16:37:305、正交变换与二次型的标准化 例 构建一个正交变换将下列二次型化为标准形 二次型的矩阵为 计算矩阵的特征值与特征向量,参考输入表达式为eigensystem {{1,1,0},{1,0,-1},{0,-1,1}} 计算得到特征值与对应的特征向量... -
Matlab中矩阵的归一化和正交化函数
2016-07-28 17:11:16返回矩阵A的正交基,B的列与A的列具有相同的空间,B的列向量是 正交向量 ,满足B'*B = eye(rank(A)),B的列数是A的秩。 应用举例 A=[4 0 0;0 3 1;0 1 3]; B=orth(A) Q=B'*B B = 0 1.0000 0 -0.7071 0 -0.7071... -
Matlab实现IQ基带生成器
2021-04-19 01:20:03IQ Baseband BuilderThe MATLAB program IQ Baseband Builder (IQBB) allows you to generate ideal as well as non-ideal IQ baseband signals using a GUI. When building the signals the user has control over ... -
matlab如何做此代数题
2021-04-20 02:58:392013-08-11二次型化为标准型的问题二次型f(X1,2014考研数学:暑期把握两个字“啃”和“钻”到目前为止,大家应该已经把教材看完一遍,并做了一些辅助习题,基础相对比较扎实了,接下来进入强化阶段。这个阶段需要... -
正交采样Matlab代码-Parameterized_ROM_Interp:局部降阶基插值应用于具有随机系数场的抛物线扩散方程
2021-06-08 11:11:16局部降阶基插值应用于具有随机系数场的抛物线扩散方程 在这项工作中,在具有对数正态系数场的参数化抛物线扩散方程上描述和数值探索了局部降阶基础 (RB) 插值方法。 RB 及其相应的降阶模型 (ROM) 对于参数 λ ∈ R^p... -
matlab如何将矩阵标准化 命令是什么?
2021-04-21 12:51:45tril 抽取下三角阵 triu 抽取上三角阵 conj 共轭矩阵 companion 伴随矩阵 det 行列式的值 norm 矩阵或向量范数 nnz 矩阵中非零元素的个数 null 清空向量或矩阵中的某个元素 orth 正交基 rank 矩阵秩 trace 矩阵迹 ... -
数值方法解出的近似旋转矩阵转化为正交的标准旋转矩阵 旋转矩阵标准化 正交矩阵
2020-10-06 20:51:34在求解涉及旋转的优化问题过程中,通过解线性方程组或者采用优化方法得到的旋转矩阵的数值解,往往没有考虑到正交约束,会得到一个存在误差的近似旋转矩阵,也就是说这个矩阵并不是一个严格正交的标准正交矩阵,不... -
matlba 正交基
2017-03-01 11:07:02matlba 正交基 -
数字信号处理Day2-小波基与规范正交化
2014-05-08 19:09:04我们可以定义出4096个基,分别是某一位是0其他是1,在这种情况下,如果我们传输图片,那么就相当于传输原始数据 假设传到一半,网络坏了。 于是,我们得到 我们可以计算原图像和这图像的差距 error = I - I_...