Compute_Slepianstuff.m
6e883b19735amaster

Compute_Slepianstuff.m
View Options

			function [V,Zeta,Mu,Xi] = Compute_Slepianstuff(MOI,S,L)

			[V,Zeta] = eig((MOI+MOI')/2);
			Zeta = diag(Zeta);
			[Zeta,idx] = sort(Zeta,'descend');
			V = V(:,idx);
			Mu = (V'SV)./(V'*V);
			Mu = diag(Mu);
			Xi = (V'L^(0.5)SL^(0.5)V)./(V'*V);
			Xi = diag(Xi);

			for v = 1:size(MOI,1)

			[~,idx_norm]=max(abs(V(:,v)));
			if sign(V(idx_norm(1),v))<0
			V(:,v)=-V(:,v);
			end
			end
			end