k1.zip_kernel matlab

clc clear [Nx]=size(ax,1); %scaling mean_X =mean(ax); std_X=std(ax); ax=ax-mean_X(ones(Nx,1),:); ax=ax./std_X(ones(Nx,1),:); c=kk; %kernel width % build kernel matrix [N]=size(ax,1); for i=1:N for j=1:N K(i,j)=exp(-norm(ax(i,:)-ax(j,:))^2/c); end end n1=ones(N,N); N1=1/N*n1; Kp=K-N1*K-K*N1+N1*K*N1; [Nt]=size(axt,1); for i=1:Nt for j=1:N Kt(i,j)=exp(-norm(axt(i,:)-ax(j,:))^2/c); end end n2=ones(1,N); N2=1/N*n2; for i=1:Nt Kpt(i,:)=Kt(i,:)-N2*K-Kt(i,:)*N1+N2*K*N1;%中心化 end %calculating the number of PCs [~,s,v]=svd(Kp); lamda=s*(1/(N-1)); P=v; lamda=diag(lamda); % B=length(find(lamda>1e-5)); npc=1; while sum(lamda(1:npc))/sum(lamda(1:B))<0.85 npc=npc+1; end lamda=diag(lamda); for i=1:B P(:,i)=P(:,i)/(sqrt((N-1)*lamda(i,i))); % load matrix P end for i=1:Nx for k=1:B tx(i,k)=P(:,k)'*Kp(i,:)'; txt(i,k)=P(:,k)'*Kpt(i,:)'; end end [Nxt,~]=size(axt); for i=1:Nxt for k=1:B txt(i,k)=P(:,k)'*Kpt(i,:)'; end end Zx=tx'; % whenting signals Zxt=txt';