figure(2)
plot(h2,kr,'LineWidth',2)
hold on
plot(thickness,kr1,'r*')
line([1.65,1.65],[1,1.45],'Color','g','LineWidth',2) %垂直线
text(1.65,1.1,' \leftarrow skin depth') %文字标注
xlabel('thickness(mm)')
ylabel('ratio AC to DC ')
legend('Analyze','FEM')
% title('Fre=400Hz')
saveas(gcf,'skinthick.png')