当前位置:首页 >> 信息与通信 >>

矩阵奇异值分解算法


function [ Psvd ] = PsvdSpectrum( M_zhenyuan,SNAP_kuaipai,d_X,thea_0,th ea_1,SNR ) N=2;%???????? C=3e8; f0=10e6; lamda=C/f0; d=d_X*lamda; fs=1000; ts=1/fs; t=(0:SNAP_kuaipai-1)*ts; a=[0:M_zhenyuan-1]';%?ó?????? j=sqrt(-1); s0=exp(j*2*pi*f0*t); s1=exp(j*2*pi*f0*t); %?ó???÷?????? a_theta0=exp(j*2*pi*d/lamda*a*sin(thea_0/180*pi)); a_theta1=exp(j*2*pi*d/lamda*a*sin(thea_1/180*pi)); A=[a_theta0 a_theta1]; S=[s0;s1]; X0=A*S; X=awgn(X0,SNR); Rxx=X*X'/length(t); %%???÷??·??? [V,D]=eig(Rxx); v1=V(:,M_zhenyuan);%×??ó???÷?ò?? q=4; p=M_zhenyuan-q+1; %????qxp???óY1 for k=1:q Y1(k,:)=v1(k:k+p-1).'; end %------------------------------------------------------%?-?????ì??·?????·¨ %------------------------------------------------------[U1,D1,V1]=svd(Y1); i=[1:q]';

Un=U1(:,(N+1):q); theta=-90:0.1:90; for ii=1:length(theta) a_theta=exp(j*2*pi*d/lamda*i*sin(theta(ii)*pi/180)) ; WW=a_theta'*Un*Un'*a_theta; Psvd(ii)=abs(1/WW); end Psvd=10*log10(Psvd/max(Psvd)+eps); end



相关文章:
矩阵的奇异值分解及其应用
特征值和奇 异值在大部分人的印象中,往往是停留在纯粹的数学计算中。而且线性代数或者矩阵论里面,也 很少讲任何跟特征值与奇异值有关的应用背景。 奇异值分解是...
奇异值分解
奇异值分解 - 2矩阵奇异值分解算法原理与分析 本章将会讨论奇异值分解(SVD)的基本概念以及矩阵奇异值分解的 Jacobi 算法原理, 包括 Jacobi 旋转变换的基础知识, ...
并行计算奇异值分解
并行计算奇异值分解--Jacobi 旋转 鉴于矩阵奇异值分解 SVD 在工程领域的广泛应用(如数据压缩、噪声去除、数值分析等等,包括在 NLP 领域的潜在 语义索引 LSI 核心...
基于奇异值分解的图像压缩处理
? 三.基于奇异值分解的图像压缩编码 1.图像的压缩 奇异值分解算法的一个重要的特征是可以降维。设A是 m ? n 型矩阵,表明 A是 n 个 m 维向量, 若A的秩...
奇异值分解
奇异值分解 - 地球物理系反演报告 实验一 奇异值分解计算广义逆 G+ 专姓学 业:地球物理学 名: 号: 指导教师:邵广周 实验一 奇异值分解计算广义逆 G+ ...
矩阵的奇异值分解
§2 矩阵奇异值分解 定义 设 A 是秩为 r 的 m × n 复矩阵, AT A ...?. ?0 2 ? 计算 U = AV1Σ ?1 = 1 ?1 2 ? ? ? ,则 A 的奇异...
对于矩阵奇异值分解的探究
?iui , i ? 1,2...n é 2 1 -2 3 ùêú 据此,可计算计算矩阵 A = ê 0 4 1 1 ú 的奇异值分解 ê 2 5 -2 4 ú ? ? ?2 0 2? ?1 ...
矩阵奇异值分解在图像隐藏中的应用
S ? ?W 3) 将新矩阵 S ' 进行奇异值分解 S ' ? U1S1V1T 4)得到含水印的图像 ? A ? US1V T 3、数据隐藏算法?1? 数字水印技术实际上就是一种...
我所理解的奇异值分解
矩阵奇异值分解计算技术中的应用 [J]. 计算机与现代 化.2006 年.第 6 期:67—68. [2] 基于奇异值分解的彩色图像水印算法* 曾文权,熊祥光,余爱民. 基于...
奇异值分解和图像主分量复原_20180129
奇异值分解和图像主分量复原_20180129 - 奇异值分解(SVD)和图像矩阵的分解测试 ? ·SVD 简单介绍 在很多情况下,数据的绝大部分信息往往集中在很小一部分数据上,...
更多相关标签: