Algoritmo das Projeções Sucessivas Para Seleção de ... - PPGQ

Anexos
clc
disp('- Seleção de variáveis para classificação via APS -')
disp(' ')
disp('Escolha uma opção e tecle ')
disp('1 - Gerar cadeias de variáveis via APS')
disp('2 - Carregar cadeias de variaveis previamente geradas')
opcao = input('Opção: ');
if (opcao == 1) % Gerar
disp(' ')
disp('O APS sera utilizado para gerar cadeias de variáveis.')
disp(' ')
disp('Entre com o nome de um arquivo para salvar as cadeias')
filename = input('(entre apóstrofos e sem extensão): ');
% No MATLAB, a inicializacao previa de vetores torna a execução mais rápida
L = zeros(N2,Nlambdas);
Xcaln = zeros(Ntrain,Nlambdas);
X = zeros(Ntrain,Nlambdas);
% Auto-Escalonamento das colunas de Xpooled
opcao = input('Deseja fazer auto-escalonamento para fins de selecao ? (Sim: 1 - Nao: 0)');
if (opcao == 1)
disp('Sera feito auto-escalonamento!')
for i=1:Nlambdas
x = Xpooled(:,i);
Xcaln(:,i)=x/std(x);
end
else
disp('Nao sera feito auto-escalonamento!')
Xcaln = Xpooled;
end
espera = input('Tecle algo para continuar...');
% Aplicacao do APS
clc
normas = sum(Xcaln.^2); % Norma ao quadrado de cada coluna de Xcaln
normamax = max(normas); % Maior norma
t0 = clock; % Comando usado para estimar quanto tempo levará a otimização
for i=1:Nlambdas
X = Xcaln;
X(:,i) = X(:,i)*2*normamax/normas(i); % forca o qr a comecar pela i-esima coluna de Xcaln
[dummy1,dummy2,ordem] = qr(X,0); % ordem contem os comprimentos de onda selecionados a
partir do i-esimo
L(:,i) = ordem(1:N2)'; % Guarda apenas N2 comprimentos de onda
% Apresenta na Tela uma estimativa do tempo restante para o fim da otimização
if (rem(i,10) == 0)
tempo = etime(clock,t0);
home
disp(['Aplicando SPA. Tempo restante = ',num2str(round((Nlambdasi)*tempo/10)),'s',blanks(10)]);
t0 = clock;
end
118