testfunction.m

粒子群算法

X1=[-4:0.1:0]';
X2=[0:0.1:4]';
X3=[-5:0.1:-1]';
X4=[-4:0.1:0]';
X5=[0:0.1:4]';
X6=[1:0.1:5]';
X7=[-4:0.1:0]';
X8=[0:0.1:4]';
X9=[0.5:0.1:4.5]';
X10=[0:0.1:4]';

% X1=rand(40,1)*5;
% X=[X1,X1,X1,X1,X1,X1,X1,X1,X1,X1]';
X=[X1,X2,X3,X4,X5,X6,X7,X8,X9,X10]';
Y=sum(X.^2+10*cos(2*pi*X)+10);
% Y=-20*exp(-0.2*(sum(X.^2)/10).^0.5)-exp(sum(cos(2*pi*X))/10)+22.71828182;
% Y=sum(X.^2/4000)-prod(cos(X./sqrt(10)+1));

Y=Y';
% N=Y+rand(size(Y,1),1)/4.*Y;
N=Y+randn(size(Y,1),1)*10;
% y1=testPNN(N',X',X');  
[y2,t2]=testPSVM(N,X',X');
[y3,t3]=testKNN(N',X',X');  
[y4,t4]=testKSVM(N,X',X');

%  E1=abs(Y-y1)./Y;
 E2=abs(Y-y2)./Y;
 E3=abs(Y-y3)./Y;
 E4=abs(Y-y4)./Y;
% meanerror1=mean(E1)
% maxerror1=max(E1)
% stderror1=std(E1)       

meanerror2=mean(E2)
maxerror2=max(E2)
stderror2=std(E2)  
t2
meanerror3=mean(E3)
maxerror3=max(E3)
stderror3=std(E3)  
t3
meanerror4=mean(E4)
maxerror4=max(E4)
stderror4=std(E4)  
t4



% subplot(411);
% plot(X1,y1,'b','linewidth',2); 
%  hold on  
%  plot(X1,Y,'r--','linewidth',2);
%  hold on   
 
 subplot(311);
 ylabel('PCA-SVM\bf%');
 plot(X1,y2,'b','linewidth',2); 
 hold on  
 plot(X1,Y,'r--','linewidth',2);
 hold on   
  
 subplot(312);
 ylabel('KPCA-NN\bf%');
plot(X1,y3,'b','linewidth',2); 
 hold on  
 plot(X1,Y,'r--','linewidth',2);
 hold on 
 
  subplot(313);
  ylabel('KPCA-SVM\bf%');
plot(X1,y4,'b','linewidth',2); 
 hold on  
 plot(X1,Y,'r--','linewidth',2);
 hold on