PIV实验结果流速矢量图、流线图

clear;clc;
for k=1:n
    str1=['file position'];
    str2=['file2 position'];
    str3=['file3 position'];
    ...
    str=eval(['str',num2str(k)]);
    load([str,'filename']);
    [m n]=size(y);
    y=flipud(y);
    v=-v;
    x=x-x(1,1);
    y=y-y(m,n);
    x1=x(:,1);
    y1=y(1,:);
    [xtemp ytemp]=meshgrid([min(x(1,:)):0.001:max(x(1,:))] ,[min(y(:,1)):0.001:max(y(:,1))]);
    V=sqrt(u.^2+v.^2);
    VV=interp2(x,y,V,xtemp,ytemp);
    set(0,'DefaultFigureVisible', 'on');
    figure
    colormap jet
    caxis([0 0.1])
    c.Label.String = 'Velocity (m/s)';%定义标题
    hold on 
    quiver(x(1:4:end,1:4:end),y(1:4:end,1:4:end),u(1:4:end,1:4:end),v(1:4:end,1:4:end),'b','MaxHeadSize',3,'autoScalefactor',1.7,'autoscale','on');%绘制速度矢量图，定义绘制密度，定义箭头大小，定义箭身大小
    set(gcf,'position',[100,100,1500,900]);%设置图像位置和分辨率
    startx=x(1,90);%设置流线起始位置
    starty=y(35,1);
    h=streamline(x(1:end,1:end),y(1:end,1:end),u(1:end,1:end),v(1:end,1:end),startx,starty);%绘制流速流线图
    set(h, 'Color', [1 0 0]); % 设置流线颜色
    box on
    axis([0 0.78 0.03 0.45]);
    xlabel('X（m）');
    ylabel('Z（m）');
end