cesar 3 weeks ago
parent
commit
45a0cdb787
3 changed files with 14 additions and 11 deletions
  1. BIN
      paper/fig4.png
  2. BIN
      paper/fig5.png
  3. 14
    11
      v3.py

BIN
paper/fig4.png View File


BIN
paper/fig5.png View File


+ 14
- 11
v3.py View File

460
         nrows=NumFeaturesToPlot, ncols=1, figsize=(15, 10), dpi=80, facecolor="w", edgecolor="k",sharex=True
460
         nrows=NumFeaturesToPlot, ncols=1, figsize=(15, 10), dpi=80, facecolor="w", edgecolor="k",sharex=True
461
     )
461
     )
462
     for i in range(NumFeaturesToPlot):
462
     for i in range(NumFeaturesToPlot):
463
+        x=[]
464
+        y=[]
465
+        for k in anomalous_data_indices:
466
+            if (k+TIME_STEPS)<x_test.shape[0]:
467
+                x.append(k+TIME_STEPS)
468
+                y.append(x_test[k+TIME_STEPS,0,indexesToPlot[i]]*stdevs[i]+means[i])
469
+        axes[i].plot(x,y ,color='black',marker='.',linewidth=0,label="Fail detection" )
470
+
471
+
463
         init=0
472
         init=0
464
         end=testRanges[0][1]
473
         end=testRanges[0][1]
465
         axes[i].plot(range(init,end),x_test[testRanges[0][0]:testRanges[0][1],0,indexesToPlot[i]]*stdevs[i]+means[i],label="No fail")
474
         axes[i].plot(range(init,end),x_test[testRanges[0][0]:testRanges[0][1],0,indexesToPlot[i]]*stdevs[i]+means[i],label="No fail")
466
         init=end
475
         init=end
467
         end+=(testRanges[1][1]-testRanges[1][0])
476
         end+=(testRanges[1][1]-testRanges[1][0])
468
         for j in range(1,NumberOfFailures+1):
477
         for j in range(1,NumberOfFailures+1):
469
-            axes[i].plot(range(init,end),x_test[testRanges[j][0]:testRanges[j][1],0,indexesToPlot[i]]*stdevs[i]+means[i],label="Fail type "+str(j), color=colorline[j-1])
478
+            axes[i].plot(range(init,end),x_test[testRanges[j][0]:testRanges[j][1],0,indexesToPlot[i]]*stdevs[i]+means[i],label="Fail type "+str(j), color=colorline[j-1],linewidth=1)
470
             if j<NumberOfFailures:
479
             if j<NumberOfFailures:
471
                 init=end
480
                 init=end
472
                 end+=(testRanges[j+1][1]-testRanges[j+1][0])
481
                 end+=(testRanges[j+1][1]-testRanges[j+1][0])
473
-        x=[]
474
-        y=[]
475
-        for k in anomalous_data_indices:
476
-            if (k+TIME_STEPS)<x_test.shape[0]:
477
-                x.append(k+TIME_STEPS)
478
-                y.append(x_test[k+TIME_STEPS,0,indexesToPlot[i]]*stdevs[i]+means[i])
479
-        axes[i].plot(x,y ,color='grey',marker='.',linewidth=0,label="Fail detection" )
482
+
480
 
483
 
481
         if i==(NumFeatures-1):
484
         if i==(NumFeatures-1):
482
             axes[i].legend(loc='right')
485
             axes[i].legend(loc='right')
621
             if (k+TIME_STEPS)<x_test.shape[0]:
624
             if (k+TIME_STEPS)<x_test.shape[0]:
622
                 x.append(k+TIME_STEPS)
625
                 x.append(k+TIME_STEPS)
623
                 y.append(x_test[k+TIME_STEPS,0,indexesToPlot[i]]*stdevs[i]+means[i])
626
                 y.append(x_test[k+TIME_STEPS,0,indexesToPlot[i]]*stdevs[i]+means[i])
624
-        axes[i].plot(x,y ,color='grey',marker='.',linewidth=0,label="Fail detection" )
627
+        axes[i].plot(x,y ,color='black',marker='.',linewidth=0,label="Fail detection" )
625
 
628
 
626
         if i==(NumFeatures-1):
629
         if i==(NumFeatures-1):
627
             axes[i].legend(loc='right')
630
             axes[i].legend(loc='right')
668
             if (k+TIME_STEPS)<x_test.shape[0]:
671
             if (k+TIME_STEPS)<x_test.shape[0]:
669
                 x.append(k+TIME_STEPS)
672
                 x.append(k+TIME_STEPS)
670
                 y.append(x_test[k+TIME_STEPS,0,indexesToPlot[i]]*stdevs[i]+means[i])
673
                 y.append(x_test[k+TIME_STEPS,0,indexesToPlot[i]]*stdevs[i]+means[i])
671
-        axes[i+1].plot(x,y ,color='grey',marker='.',linewidth=0,label="Fail detection" )
674
+        axes[i+1].plot(x,y ,color='black',marker='.',linewidth=0,label="Fail detection" )
672
         if i==0:
675
         if i==0:
673
             axes[i+1].legend(loc='right')
676
             axes[i+1].legend(loc='right')
674
         s=''
677
         s=''
729
             if (k+TIME_STEPS)<x_test.shape[0]:
732
             if (k+TIME_STEPS)<x_test.shape[0]:
730
                 x.append(k+TIME_STEPS)
733
                 x.append(k+TIME_STEPS)
731
                 y.append(x_test[k+TIME_STEPS,0,indexesToPlot[i]])
734
                 y.append(x_test[k+TIME_STEPS,0,indexesToPlot[i]])
732
-        axes[i].plot(x,y ,color='grey',marker='.',linewidth=0,label="Fail detection" )
735
+        axes[i].plot(x,y ,color='black',marker='.',linewidth=0,label="Fail detection" )
733
 
736
 
734
         if i==0:
737
         if i==0:
735
             axes[i].legend(bbox_to_anchor=(0.9, 0.4))
738
             axes[i].legend(bbox_to_anchor=(0.9, 0.4))

Powered by TurnKey Linux.