mysort.m

工程计算MATLAB code to calculate the reorthogonalized sine tapers input: N = the length of the time se

function [y, in]=mysort(x)

%%MATLAB function 
%%input:  x  = vector 
%%output: y  = sorted in decreasing order x
%%        in = the index of the original x values; x(in)=y
 
[a, i1] = sort(x);
n=length(x);
for i=1:n
        in(i) = i1(n-i+1);
        y(i) = a(n-i+1);
end