0

我有以下脚本:

font_size       = 15
font_size_ticks = 14

### create figure ###
fig, (ax1, ax2) = plt.subplots(nrows =1, ncols =2, figsize=(20, 10),sharex=False,sharey=False)
fig.subplots_adjust(hspace = .2, wspace=.15)

### For min, max x axis limits ###
### create cmap ###
n = 5
cmap = pl.cm.rainbow
col = cmap(np.linspace(0,1,n))
maxim1 = []
maxim2 = []
for  i in range(5):
    x =10**10**np.random.rand(100);    y =10**10**np.random.rand(100);
    x1 =10**10**np.random.rand(100);    y1 =10**10**np.random.rand(100);
    maxim1.append(max(x))                                            # subtruct (-small, small=resolution)
    maxim2.append(max(x1))

    ax1.plot(np.sort(y),1/np.sort(x),'.',ls=':',ms=4,linewidth=1.4,color = col[i])
    ax2.plot(np.sort(y1),1/np.sort(x1),'.',ls=':',ms=4,linewidth=1.4,color = col[i])

  
    ax1.set_yscale('log');                        ax2.set_yscale('log')
    ax1.set_xscale('log');                        ax2.set_xscale('log')


    ax1.xaxis.grid(True, "minor", linewidth=.25); ax2.xaxis.grid(True, "minor", linewidth=.25);
    ax1.yaxis.grid(True, "minor", linewidth=.25); ax2.yaxis.grid(True, "minor", linewidth=.25);
    ax1.grid();                                   ax2.grid();

    ax1.set_xlabel('$Waiting \ time \ [s]$', fontsize = font_size);     ax2.set_xlabel('$Duration \ [s]$', fontsize = font_size)
    ax1.set_ylabel('$PDF$', fontsize = font_size);                      ax2.set_ylabel('$PDF$', fontsize = font_size)

    #### Set tick size ####
    ax1.tick_params(axis='x',labelsize=font_size_ticks);               ax2.tick_params(axis='x',labelsize=font_size_ticks);
    ax1.tick_params(axis='y',labelsize=font_size_ticks);               ax2.tick_params(axis='y',labelsize=font_size_ticks);   

    #### Set fig dimensions ####

    ax1.legend( frameon=False,loc=3, fontsize = font_size); ax2.legend( frameon=False,loc=3, fontsize = font_size);

    max1 = max(np.array(maxim1));                 max2 = max(np.array(maxim2))
    ax1.set_xlim([None,max1+1e-3*max1]);          ax2.set_xlim([None,max2+ 2*1e-1*max2])

我试图让网格线和小刻度出现,但这就是我得到的:

在此处输入图像描述

我已经尝试了很多东西,但无法让它发挥作用!关于如何完成这项工作的任何想法?

4

0 回答 0