Add files via upload

foryichuanqi · Feb 21, 2024 · 602805c · 602805c
1 parent c380308
commit 602805c
Show file tree

Hide file tree

Showing 44 changed files with 11,358 additions and 0 deletions.
diff --git a/code/comparison experiments/IGBT/DPA/DPA.py b/code/comparison experiments/IGBT/DPA/DPA.py
diff --git a/code/comparison experiments/IGBT/prediction_value_and_error_IGBT_3_4_error_band.py b/code/comparison experiments/IGBT/prediction_value_and_error_IGBT_3_4_error_band.py
diff --git a/code/comparison experiments/IGBT/tafcn2022/LOGS-LSTM-Keras-CMAPSS.txt b/code/comparison experiments/IGBT/tafcn2022/LOGS-LSTM-Keras-CMAPSS.txt
diff --git a/code/comparison experiments/IGBT/tafcn2022/TaFCN.py b/code/comparison experiments/IGBT/tafcn2022/TaFCN.py
diff --git a/code/comparison experiments/cacle/DPA_our_method/DPA.py b/code/comparison experiments/cacle/DPA_our_method/DPA.py
diff --git a/.../comparison experiments/cacle/prediction_value_and_error_cacle_3536_38_error_band_bian.py b/.../comparison experiments/cacle/prediction_value_and_error_cacle_3536_38_error_band_bian.py
diff --git a/code/comparison experiments/cacle/tafcn2022/LOGS-LSTM-Keras-CMAPSS.txt b/code/comparison experiments/cacle/tafcn2022/LOGS-LSTM-Keras-CMAPSS.txt
diff --git a/code/comparison experiments/cacle/tafcn2022/TaFCN.py b/code/comparison experiments/cacle/tafcn2022/TaFCN.py
diff --git a/code/comparison experiments/nasa/DPA_our_method/DPA.py b/code/comparison experiments/nasa/DPA_our_method/DPA.py
diff --git a/code/comparison experiments/nasa/prediction_value_and_error_nasa_6_7 .py b/code/comparison experiments/nasa/prediction_value_and_error_nasa_6_7 .py
diff --git a/code/comparison experiments/nasa/tafcn2022/LOGS-LSTM-Keras-CMAPSS.txt b/code/comparison experiments/nasa/tafcn2022/LOGS-LSTM-Keras-CMAPSS.txt
diff --git a/code/comparison experiments/nasa/tafcn2022/TaFCN.py b/code/comparison experiments/nasa/tafcn2022/TaFCN.py
diff --git a/code/comparison experiments/xiandao/DPA_our_method/DPA.py b/code/comparison experiments/xiandao/DPA_our_method/DPA.py
diff --git a/code/comparison experiments/xiandao/prediction_value_and_error_xiandao_c1_c4_to_c6.py b/code/comparison experiments/xiandao/prediction_value_and_error_xiandao_c1_c4_to_c6.py
diff --git a/code/comparison experiments/xiandao/tafcn2022/LOGS-LSTM-Keras-CMAPSS.txt b/code/comparison experiments/xiandao/tafcn2022/LOGS-LSTM-Keras-CMAPSS.txt
diff --git a/code/comparison experiments/xiandao/tafcn2022/TaFCN.py b/code/comparison experiments/xiandao/tafcn2022/TaFCN.py
diff --git a/code/data_read/read_mat_for_IGBT.py b/code/data_read/read_mat_for_IGBT.py
@@ -0,0 +1,161 @@
+# -*- coding: utf-8 -*-
+"""
+Created on Tue Mar 14 15:25:26 2023
+
+@author: Administrator
+"""
+
+import numpy as np
+
+import scipy as sp
+
+import math
+
+from numpy import matmul as mm
+from math import sqrt,pi,log, exp
+
+import xlrd
+
+import matplotlib.pyplot as plt
+
+import numpy as np
+
+from scipy.optimize import minimize, rosen, rosen_der
+
+from scipy.stats import linregress
+
+from scipy.stats import norm
+
+
+import scipy.io as scio
+
+# path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 2\Device2  1.mat'
+# matdata = scio.loadmat(path)
+# print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+# x=matdata["measurement"][0][0][1][0][0][2][0][0][4][0]
+
+# plt.plot(range(0,len(x),10000),x[range(0,len(x),10000)])
+# plt.show()
+
+
+
+invaribale_index=3    #########dtype=[('dt', 'O'), ('gateSignalVoltage', 'O'), ('gateEmitterVoltage', 'O'), ('collectorEmitterVoltage', 'O'), ('collectorEmitterCurrentSignal', 'O')]),
+
+
+
+
+
+
+
+
+path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 2\Device2  1.mat'
+matdata = scio.loadmat(path)
+# print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+Vces_list=[]
+
+for i in range(len(matdata["measurement"][0][0][1][0])):
+
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][110000])
+    Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][-1])
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][50000])
+
+
+print(Vces_list)
+
+plt.plot(range(len(Vces_list)),Vces_list)
+plt.show()
+
+
+
+
+
+path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 3\Device3  1.mat'
+matdata = scio.loadmat(path)
+# print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+Vces_list=[]
+
+for i in range(len(matdata["measurement"][0][0][1][0])):
+
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][110000])    
+    Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][-1])
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][50000])
+
+
+print(Vces_list)
+
+plt.plot(range(len(Vces_list)),Vces_list)
+plt.show()
+
+
+
+
+path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 4\Device4  1.mat'
+matdata = scio.loadmat(path)
+# print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+Vces_list=[]
+
+for i in range(len(matdata["measurement"][0][0][1][0])):
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][110000])    
+    Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][-1])
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][50000])
+
+
+print(Vces_list)
+
+plt.plot(range(len(Vces_list)),Vces_list)
+plt.show()
+
+
+
+
+
+path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 5\Device5  1.mat'
+matdata = scio.loadmat(path)
+# print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+Vces_list=[]
+
+for i in range(len(matdata["measurement"][0][0][1][0])):
+
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][110000])    
+    Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][-1])
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][50000])
+
+print(Vces_list)
+
+plt.plot(range(len(Vces_list)),Vces_list)
+plt.show()
+
+
+
+
+
+# epoch_index=6
+
+
+# path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 2\Device2  1.mat'
+# matdata = scio.loadmat(path)
+# # print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+# # x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+# x=matdata["measurement"][0][0][1][0][epoch_index][2][0][0][invaribale_index][0]
+
+# print(x)
+# plt.plot(range(0,len(x),10000),x[range(0,len(x),10000)])
+# plt.show()
+
diff --git a/code/data_read/test.py b/code/data_read/test.py
@@ -0,0 +1,170 @@
+# -*- coding: utf-8 -*-
+"""
+Created on Tue Mar 14 17:20:51 2023
+
+@author: Administrator
+"""
+
+# -*- coding: utf-8 -*-
+"""
+Created on Tue Mar 14 15:25:26 2023
+
+@author: Administrator
+"""
+
+import numpy as np
+
+import scipy as sp
+
+import math
+
+from numpy import matmul as mm
+from math import sqrt,pi,log, exp
+
+import xlrd
+
+import matplotlib.pyplot as plt
+
+import numpy as np
+
+from scipy.optimize import minimize, rosen, rosen_der
+
+from scipy.stats import linregress
+
+from scipy.stats import norm
+
+
+import scipy.io as scio
+
+# path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 2\Device2  1.mat'
+# matdata = scio.loadmat(path)
+# print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+# x=matdata["measurement"][0][0][1][0][0][2][0][0][4][0]
+
+# plt.plot(range(0,len(x),10000),x[range(0,len(x),10000)])
+# plt.show()
+
+
+
+invaribale_index=3    #########dtype=[('dt', 'O'), ('gateSignalVoltage', 'O'), ('gateEmitterVoltage', 'O'), ('collectorEmitterVoltage', 'O'), ('collectorEmitterCurrentSignal', 'O')]),
+
+
+
+
+
+
+
+
+path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 2\Device2  1.mat'
+matdata = scio.loadmat(path)
+# print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+Vces_list=[]
+
+for i in range(len(matdata["measurement"][0][0][1][0])):
+
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][110000])
+    Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][-1])
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][50000])
+
+
+print(Vces_list)
+
+plt.plot(range(len(Vces_list)),Vces_list)
+plt.show()
+
+
+
+
+
+path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 3\Device3  1.mat'
+matdata = scio.loadmat(path)
+# print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+Vces_list=[]
+
+for i in range(len(matdata["measurement"][0][0][1][0])):
+
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][110000])    
+    Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][-1])
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][50000])
+
+
+print(Vces_list)
+# print()
+
+plt.plot(range(len(Vces_list)),Vces_list)
+plt.show()
+
+
+
+
+path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 4\Device4  1.mat'
+matdata = scio.loadmat(path)
+# print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+Vces_list=[]
+
+for i in range(len(matdata["measurement"][0][0][1][0])):
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][110000])    
+    Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][-1])
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][50000])
+
+
+print(Vces_list)
+
+plt.plot(range(len(Vces_list)),Vces_list)
+plt.show()
+
+
+
+
+
+path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 5\Device5  1.mat'
+matdata = scio.loadmat(path)
+# print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+Vces_list=[]
+
+for i in range(len(matdata["measurement"][0][0][1][0])):
+
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][110000])    
+    Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][-1])
+    # Vces_list.append(matdata["measurement"][0][0][1][0][i][2][0][0][invaribale_index][0][50000])
+
+
+# print(x)
+
+plt.plot(range(len(Vces_list)),Vces_list)
+plt.show()
+
+
+
+
+
+# epoch_index=9
+
+
+# path = r'F:\桌面11.17\project\RUL_with_highly_small_sample\dataset\IGBTAgingData_04022009\Data\Thermal Overstress Aging with Square Signal at gate and SMU data\Aging Data\Device 2\Device2  1.mat'
+# matdata = scio.loadmat(path)
+# # print(matdata["measurement"][0][0][1][0][0][2][0][0][4])
+# # x=matdata["measurement"][0][0][1][0][0][2][0][0][invaribale_index][0]
+
+# x=matdata["measurement"][0][0][1][0][epoch_index][2][0][0][invaribale_index][0]
+
+# print(x)
+# plt.plot(range(0,len(x),10000),x[range(0,len(x),10000)])
+# plt.show()
+