博客

%内孤立波过浅体，流场时间序列
clear;clc;
    load(mat);
    y=flipud(y); v=-v;%quiver(x,y,u,v);
    C0=sqrt(U.^2+V.^2);
    scale=scalek(k);
    [x,y]=meshgrid(0:interval:xmax,linspace(ystart,yend,ynumber));
    contourf(x,y,C0,100,'linestyle','none');%绘制云图
    colormap gray;%设置颜色格式，还有 jet summer white pink spring等
    c=colorbar('eastoutside');%设置colorbar位置
    caxis([0 0.02]);
    c.Label.String = 'Velocity (m/s)';
    set(gcf,'position',[100,100,1500,900]);
    name=['figurename'];
    title(name);
    set(gca, 'fontsize',fontsize,'fontname','Times New Roman ') ;%设置图片字体字号大小
    xlabel('Time (s)','FontName','Times New Roman','FontSize',fontsize);
    ylabel('Z (m)','FontName','Times New Roman','FontSize',fontsize);
    img=gcf;
    box on
    axis equal;
    print(img, '-dpng', '-r300', ['E:\18.12实验 mat 1229\图片\timeseries_19',num2str(yita(k)),'_300.png']);%储存图像
    close;
end