python 如何画双纵坐标柱状图

u = [4026        3967        4005        4073        4139        3961        3947        4118        4169        3991        4007        4000        3984        4089        4112        3980        3982        4019        4124        4021        4024        4062        3973        3967        4053        4096        4054        4059        3901        4014        4166        4108        4082        4048        4025        4086        4065        4102        4012        4123        4008        4098        4140        4017        4061        4101        4063        4048        4127        4019        3926        4068        4043        4033        4017        4033        4067        4008        4010        4090        4087        4012        3958        4127        4148        4178        4102        3920        4113        4048        4042        4096        4106        3954        4063        4138        4112        3968        3988        3957        4007        3920        4132        4050        3999        3992        4136        3878        3977        4158        4112        4041        3995        3961        4034        4106        4037        4115        4068        3996        4006        4057        4049        4098        4106        4044        3982        3998        4102        3971        4068        4007        4027        4019        4088        4086        3968        4106        4054        4077        3962        4034        4023        3968        4093        4122        4097        3959        3983        4027        4071        3995        4036        4108        4119        3966        4041        4016        4119        4103        3953        4027        3996        4053        4124        4119        4069        4035        4102        4082        3998        3997        4101        3958        4203        4057        4179        3974        4103        4085        3999        4070        4080        4072        3995        4111        4020        4016        4068        4019        4216        4094        4070        4010        4055        4061        4051        4090        3903        3989];

[a,b] = hist(u,13);

bar(b,a/sum(a)),xlabel('组距'),ylabel('频率'),title('t=1s时频率直方图与正态分布')

ax1=gca ;

ax2=axes('position',get(ax1,'position'),'yaxislocation','right','color','none');

axis([3891 4203 0 0.007]);

hold on

fplot('normpdf(u,4047.227778,62.52050663)',[3891,4203])

axis off

这是我的解决办法

由于那个坐标混乱，不知道怎么弄，所以用了axis off把坐标去掉了。

不知道其他人有什么解决办法！希望可以一起分享