int ir5_ia(){ scltracefa(§, ƒ, ∅); scltracefc("scl-%s (%s) : %s\n", sclver(), sclsec(), temps_char(0)); // http://recherche.telecom-bretagne.eu/iono/ir5_publication.php scl::string_c datname="ir5.dat", datpath="./don/sul/"+datname; int ir5an, ir5mois; double cir5, cjj; char ir5op; // statut observation/prevision std::vector<double> ir5a, ir5o, ir5p, jja, jjo, jjp; // Jour Julien 1er jan 1950 0[h] double jj_dis = date_jj(1950, 1, 1, 0, 0, 0, 1); FILE *fid = fopen(datpath.c_str(), "rt"); while(!feof(fid)){ fscanf(fid, "%d %d %lf %c", &ir5an, &ir5mois, &cir5, &ir5op); fgetc(fid); // on passe la fin de ligne //scltracefc("%d;%d;%lf;'%c'\n", ir5an, ir5mois, ir5val, ir5op); cjj = date_jj(ir5an, ir5mois, 1, 0, 0, 0, 1) - jj_dis; if(ir5op=='O'){ ir5o.push_back(cir5); jjo.push_back(cjj); } else if(ir5op=='P'){ ir5p.push_back(cir5); jjp.push_back(cjj); } } fclose(fid); fid=∅; //--------------------------------------------------------------------- datname.clear(); datname="ir5_archive.dat"; datpath.clear(); datpath="./don/sul/"+datname; fid = fopen(datpath.c_str(), "rt"); while(!feof(fid)){ fscanf(fid, "%d %d %lf %c", &ir5an, &ir5mois, &cir5, &ir5op); fgetc(fid); // on passe la fin de ligne //scltracefc("%d;%d;%lf;'%c'\n", ir5an, ir5mois, ir5val, ir5op); cjj = date_jj(ir5an, ir5mois, 1, 0, 0, 0, 1) - jj_dis; ir5a.push_back(cir5); jja.push_back(cjj); } fclose(fid); fid=∅; int o5n=(int)ir5o.size(), o5i, p5n=(int)ir5p.size(), a5n=(int)ir5a.size(), a5i; double ir5e=-∞_rd; cir5=max<double>(o5n, &ir5o[0]); if(ir5e<cir5){ ir5e = cir5; } cir5=max<double>(p5n, &ir5p[0]); if(ir5e<cir5){ ir5e = cir5; } cir5=max<double>(a5n, &ir5a[0]); if(ir5e<cir5){ ir5e = cir5; } double ir5r = (ir5e-0)*10/100/2; int xn=o5n+a5n, xi; double *x=new double[xn], *t=new double[xn]; for(a5i=0,xi=a5n-1;a5i<a5n;a5i++,xi--){ x[xi]=ir5a[a5i]/ir5e; t[xi]=jja[a5i]; } for(o5i=0,xi=a5n ;o5i<o5n;o5i++,xi++){ x[xi]=ir5o[o5i]/ir5e; t[xi]=jjo[o5i]; } int an=40, cn=23, en=6; int chn=3; int nrn[3], nlt[3]; nrn[0]=an; nrn[1]=cn; nrn[2]=en; nlt[0]=SCL_NLTYPE_LINR; nlt[1]=SCL_NLTYPE_SIGM; nlt[2]=SCL_NLTYPE_SIGM; double η = 0.015; int sbln=(xn-an)/en, sbli; double **xr=new double*[sbln], *cx; double **zr=new double*[sbln], *cz; double **yr_0=new double*[sbln], **yr_1=new double*[sbln], **yr_4=new double*[sbln], *cy; double *y_0=new double[xn], *y_1=new double[xn], *y_4=new double[xn]; for(sbli=0;sbli<sbln;sbli++){ cx=&x[sbli*en]; xr[sbli] =cx; cz=&x[an+sbli*en]; zr[sbli] =cz; cy=&y_0[an+sbli*en]; yr_0[sbli]=cy; cy=&y_1[an+sbli*en]; yr_1[sbli]=cy; cy=&y_4[an+sbli*en]; yr_4[sbli]=cy; } double **xj=newmat<double>(sbln, an, 0); double **zj=newmat<double>(sbln, en, 0); double **yj=newmat<double>(sbln, en, 0); // Poids synaptiques double ***wr = new double**[chn]; wr[0] = newmat<double>(an, 1); wr[1] = newmat<double>(cn, an); wr[2] = newmat<double>(en, cn); double ***wj = new double**[chn]; wj[0] = newmat<double>(an, 1); wj[1] = newmat<double>(cn, an); wj[2] = newmat<double>(en, cn); int itrn=1; tic(); nevronx(chn, &nrn[0], &nlt[0], sbln, xr, xj, zr, zj, itrn, ε_rd, η, 0, yr_0, yj, wr, wj, ∅); scltracefc("x%d : %lf [ms]\n", itrn, tac()); scl::string_c legend = "Itération x"+num_str<int>(itrn); itrn=10; tic(); nevronx(chn, &nrn[0], &nlt[0], sbln, xr, xj, zr, zj, itrn, ε_rd, η, 0, yr_1, yj, wr, wj, ∅); scltracefc("x%d : %lf [ms]\n", itrn, tac()); legend += ";Itération x"+num_str<int>(itrn); itrn=10000; double *eq = new double[itrn]; tic(); nevronx(chn, &nrn[0], &nlt[0], sbln, xr, xj, zr, zj, itrn, ε_rd, η, 0, yr_4, yj, wr, wj, eq); scltracefc("x%d : %lf [ms]\n", itrn, tac()); legend += ";Itération x"+num_str<int>(itrn); for(xi=0;xi<xn;xi++){ y_0[xi]*=ir5e; y_1[xi]*=ir5e; y_4[xi]*=ir5e; } //--------------------------------------------------------------------- scl::string_c imgname=scl::string_c(ƒ)+".svg", imgpath="./srt/"+imgname; graf_dis(imgpath.c_str(), o5n, &jjo[0], &ir5o[0], "curve", p5n, &jjp[0], &ir5p[0], "curve", a5n, &jja[0], &ir5a[0], "curve", 1, &jjo[o5n-1], &ir5o[o5n-1], "curve", 1, &jjp[p5n-1], &ir5p[p5n-1], "grid", "on", "subtitlex", "https://recherche.telecom-bretagne.eu/iono/ir5_publication.php", "linespec", "-l;-r;-l;h+;h+", "dataxis", "x", "legend", "Observation;Prédiction", "legcorn", "ula", "title", "Évolution de l'indice IR5", "axmax", (jjp[p5n-1]+11*12*30.5), "aymin", 0-ir5r, "aymax", ir5e+ir5r, "markersize", 20, "linespeca", "-a", ∅); //--------------------------------------------------------------------- imgname.clear(); imgname=scl::string_c(ƒ)+"-drnr-y.svg"; imgpath.clear(); imgpath="./srt/"+imgname; graf_dis(imgpath.c_str(), xn, t, y_0, "curve", xn, t, y_1, "curve", xn, t, y_4, "grid", "on", "linespec", "-c;-m;-u", "dataxis", "x", "legend", legend.c_str(), "legcorn", "ula", "title", "Dernières sorties d'apprentissage Sigmoïde", "axmax", (jjp[p5n-1]+11*12*30.5), "aymin", 0-ir5r, "aymax", ir5e+ir5r, "linespeca", "-a", ∅); //--------------------------------------------------------------------- imgname.clear(); imgname=scl::string_c(ƒ)+"-eq.svg"; imgpath.clear(); imgpath="./srt/"+imgname; graf_dis(imgpath.c_str(), itrn, ∅, eq, "grid", "on", "linespec", "-u", "title", "Erreur quadratique d'apprentissage Sigmoïde", "xlabel", "Numéro d'itération", "logaxis", "y", "linespeca", "-a", ∅); //--------------------------------------------------------------------- //[>LIBR. MEMOIRE] //--------------------------------------------------------------------- imgpath.clear(); imgname.clear(); legend.clear(); deletemat<double>(an, &wr[0]); deletemat<double>(cn, &wr[1]); deletemat<double>(en, &wr[2]); delete[] wr; wr=∅; deletemat<double>(an, &wj[0]); deletemat<double>(cn, &wj[1]); deletemat<double>(en, &wj[2]); delete[] wj; wj=∅; deletemat<double>(sbln, &yj); deletemat<double>(sbln, &zj); deletemat<double>(sbln, &xj); delete[] y_4; y_4=∅; delete[] y_1; y_1=∅; delete[] y_0; y_0=∅; delete[] yr_4; yr_4=∅; delete[] yr_1; yr_1=∅; delete[] yr_0; yr_0=∅; delete[] zr; zr=∅; delete[] xr; xr=∅; delete[] t; t=∅; delete[] x; x=∅; jjp.clear(); jjo.clear(); jja.clear(); ir5p.clear(); ir5o.clear(); ir5a.clear(); datpath.clear(); datname.clear(); delete[] eq; eq=∅; //--------------------------------------------------------------------- //[<LIBR. MEMOIRE] //--------------------------------------------------------------------- scltracefe(§, ƒ, ∅); return 0; } /* ƒ décorée par le 04-03-2025 23:32:14 */

[>..\xpl\src\scl.ia.cpp.ir5_ia] scl-25.03 (gwin64) : 04-03-2025 23:32:14 [<..\src\sclronx.c.nevronx>] Iteration 0/1 Error 0.346211 x1 : 4.967600 [ms] [<..\src\sclronx.c.nevronx>] Iteration 0/10 Error 0.346211 [<..\src\sclronx.c.nevronx>] Iteration 1/10 Error 0.279188 [<..\src\sclronx.c.nevronx>] Iteration 2/10 Error 0.275200 [<..\src\sclronx.c.nevronx>] Iteration 3/10 Error 0.268745 [<..\src\sclronx.c.nevronx>] Iteration 4/10 Error 0.257574 [<..\src\sclronx.c.nevronx>] Iteration 5/10 Error 0.237437 [<..\src\sclronx.c.nevronx>] Iteration 6/10 Error 0.204412 [<..\src\sclronx.c.nevronx>] Iteration 7/10 Error 0.173927 [<..\src\sclronx.c.nevronx>] Iteration 8/10 Error 0.168145 [<..\src\sclronx.c.nevronx>] Iteration 9/10 Error 0.166023 x10 : 49.496400 [ms] x10000 : 53044.374100 [ms] <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< << END OF DISLIN / VERSION 11.3.2 << << Date : 04.03.2025 Time : 23:33:07 Pageformat: DA4L << << Vectors : 4174 Warnings: 0 Fileformat: SVG << << Metafile: ./srt/ir5_ia.svg << <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< << END OF DISLIN / VERSION 11.3.2 << << Date : 04.03.2025 Time : 23:33:07 Pageformat: DA4L << << Vectors : 10647 Warnings: 0 Fileformat: SVG << << Metafile: ./srt/ir5_ia-drnr-y.svg << <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< << END OF DISLIN / VERSION 11.3.2 << << Date : 04.03.2025 Time : 23:33:07 Pageformat: DA4L << << Vectors : 11939 Warnings: 0 Fileformat: SVG << << Metafile: ./srt/ir5_ia-eq.svg << <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< [<..\xpl\src\scl.ia.cpp.ir5_ia]