コード例 #1
0
ファイル: operator_data.py プロジェクト: kankan1987/algorithm
def loaddata(file, delimeter):
    #加载文件,以delimeter进行分割成行和列
    data = np.loadtxt(file, delimiter=delimeter)
    #输出处理后数据的行和列
    #print('Dimensions: ',data.shape)
    #输出第2行到第6行
    #print(data[1:6,:])
    return (data)
コード例 #2
0
ファイル: plot_gen.py プロジェクト: sharadmv/handwriting-dtw
def main(argv):
    letters = ["gen/"+x+"-0" for x in process_args(argv)]
    to_plot = []

    plot = plt.subplot(3,1,1)
    for letter in letters:
        data = np.loadtxt(letter,delimiter=" ");
        print data
        plot.plot(data, label=letter[4:-2])

    plot.legend(shadow=True,fancybox=True)
    plt.draw()
    plt.show()
コード例 #3
0
    def evaluate_design(self, x_denorm):

        with open('_eval/ACMConfig.h', 'r') as f:
            new_line = []
            for line in f.readlines():
                if '#define CURRENT_LOOP_BANDWIDTH' in line:
                    new_line.append(
                        f'#define CURRENT_LOOP_BANDWIDTH (2*M_PI*{x_denorm[0]:g})\n'
                    )  # 15
                elif '#define DELTA_THE_DAMPING_FACTOR' in line:
                    new_line.append(
                        f'#define DELTA_THE_DAMPING_FACTOR ({x_denorm[1]:g})\n'
                    )  # 0.08
                else:
                    new_line.append(line)
        with open('_eval/ACMConfig.h', 'w') as f:
            f.writelines(new_line)

        if os.path.exists(self.work_dir + '/_eval/main.exe'):
            os.remove(self.work_dir + '/_eval/main.exe')
        subprocess.run([
            "gcc", "_eval/main.c", "_eval/controller.c", "_eval/observer.c",
            "-o", "_eval/main"
        ])
        while not os.path.exists(self.work_dir + '/_eval/main.exe'):
            time.sleep(0.1)
            print('sleep for main.exe')

        # subprocess.run( [self.work_dir+'/_eval/main.exe'] )
        os.system('cd _eval && main')
        while not os.path.exists(self.work_dir + '/_eval/pi_opti.dat'):
            time.sleep(0.1)
            print('sleep for .dat')

        speed_profile = np.loadtxt('_eval/pi_opti.dat', skiprows=1)
        os.remove(self.work_dir + '/_eval/pi_opti.dat')

        TS = 0.000025  # ACMConfig.h
        print(max(speed_profile), 'rpm')
        print(min(speed_profile, key=lambda x: abs(x - 90)), 'rpm')  # [rpm]
        print(np.abs(speed_profile - 90).argmin() * TS * 1000, 'ms')  # [ms]

        rise_time = np.abs(speed_profile - 90).argmin() * TS * 1000

        over_shoot = max(speed_profile)
        if over_shoot <= 100:
            over_shoot = 100

        # modify this function to call matlab and eMach to evaluate the design with free variables as in x_denorm
        return over_shoot, rise_time
コード例 #4
0
def main(argv):
    args = process_args(argv)
    #letters = ["train/sharad-new/"+x for x in args]
    letters = [""+x for x in args]
    #letters = ["data_train/"+x for x in args]
    #letters = ["train/lei-1354774953143/"+x for x in args]
    to_plot = []

    plot = plt.subplot(1,1,1)
    for i in range(len(letters)):
        data = np.loadtxt(letters[i],delimiter=" ");
        x = [point[0] for point in data]
        y = [point[1] for point in data]
        plot.plot(x,y, label=args[i],linewidth=8)

    #plot.legend(shadow=True,fancybox=True,loc='upper center',bbox_to_anchor=(0.5,-0.1),ncol=4)
    plot.set_aspect('1.0')
    #subplots_adjust(bottom=0.25)   # <--
    plt.draw()
    plt.show()