Advertisement
Guest User

Untitled

a guest
Nov 24th, 2017
71
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.91 KB | None | 0 0
  1. clear all;
  2. close all;
  3. FileName = 'p_app1.mat';
  4. FolderName = 'C:\FuzzyTP';
  5. File = fullfile(FolderName, FileName);
  6. load(File); % not: load('File')
  7.  
  8. FileName = 'p_app2.mat';
  9. File = fullfile(FolderName, FileName);
  10. load(File); % not: load('File')
  11.  
  12. FileName = 'p_app3.mat';
  13. File = fullfile(FolderName, FileName);
  14. load(File); % not: load('File')
  15.  
  16.  
  17. FileName = 'p_test1.mat';
  18. File = fullfile(FolderName, FileName);
  19. load(File); % not: load('File')
  20.  
  21.  
  22. FileName = 'p_test2.mat';
  23. File = fullfile(FolderName, FileName);
  24. load(File); % not: load('File')
  25.  
  26.  
  27. FileName = 'p_test3.mat';
  28. File = fullfile(FolderName, FileName);
  29. load(File); % not: load('File')
  30.  
  31. App1 = [p_app1(1:4,1:10),p_app1(1:4,11:20),p_app2(1:4,1:10),p_app3(1:4,1:10)];
  32. App2 = [p_app2(1:4,1:10),p_app2(1:4,11:20),p_app1(1:4,1:10),p_app3(1:4,1:10)];
  33. App3 = [p_app3(1:4,1:10),p_app3(1:4,11:20),p_app1(1:4,1:10),p_app2(1:4,1:10)];
  34. Test1 = [p_test1(1:4,1:10),p_test2(1:4,1:10),p_test3(1:4,1:10)];
  35. Test2 = [p_test2(1:4,1:10),p_test1(1:4,1:10),p_test3(1:4,1:10)];
  36. Test3 = [p_test3(1:4,1:10),p_test1(1:4,1:10),p_test2(1:4,1:10)];
  37.  
  38.  
  39. for i=1:3
  40. stri=int2str(i);
  41. p_App = eval(['App',stri]);
  42. [nbre_caract,nbre_exemp] = size(p_App);
  43. nero_cache = 10; %nbre de neurones de couche cachée
  44. t=[ones(1,20),(-1*ones(1,20))];
  45. % creation d'un réseau
  46. net_rn = newff(minmax(p_App),[nbre_caract,nero_cache,1]);
  47. net_rn=init(net_rn); %initialisation de RN
  48. net_rn.trainparam.epochs=10;%nbre de epoachs
  49. net_rn.trainparam.goal=0.01;%error
  50. net_rn.trainParam.lr=0.9;%taux d'app
  51. net_rn= train(net_rn,p_App,t); %lancer train
  52. p_Test = eval(['Test',stri]);
  53. seuil = 0.5;%
  54. vect=sim(net_rn,p_Test);%simulation enable R
  55. vect(find(vect>=seuil))=1;
  56. vect(find(vect<seuil))=0;
  57. tr(i)=100*(sum(vect(1:10))+(20-sum(vect(11:30))))/30;
  58. end
  59. TR = mean(tr);
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement