vector_analyse.m

过滤距阵,算数平均值过滤误差值,文件读取.数组操作等.

clear all; close all;
for i=0:1:799
    head = '20090323_wang0';
    speichername='20090323_';
    x_pix=1000;
    y_pix=1016;
    x0_pix=502;
    y0_pix=519;
    r_i=70;
    r_a=120;
    einheit=0.268200989;
    x0=x0_pix*einheit;
    y0=y0_pix*einheit;
    Omeg0=8.057;
    corr=1;
    
    if i<10
        datenname=[head,'000',num2str(i)];
    elseif i<100
        datenname=[head,'00',num2str(i)];
    elseif i<1000
        datenname=[head,'0',num2str(i)];
    else
        datenname=[head,num2str(i)];
    end;
        
    [x,y,u,v] = uvFilter(datenname,r_i,r_a,x0,y0,Omeg0);
    
    uv_betrag=abs(sqrt(u.^2+v.^2));
    a=sum(uv_betrag)/length(uv_betrag);
    b=abs(uv_betrag-a);
    c=0.3*sqrt((sum(b))^2/(length(b)-1));
    d=find(b>c);
    uv=[x,y,u,v];
    P = [];
    Q=[P,d];
    uv(Q,:)=[];    
    x=uv(:,1); y=uv(:,2); u=uv(:,3); v=uv(:,4);
    
    MM=[x,y,u,v];
    mat_format = [speichername,num2str(i),'.mat'];
    save(mat_format,'x','y','u','v')
    asc_format = [speichername,num2str(i),'.asc'];
    dlmwrite(asc_format,MM,'delimiter','\t',...
        'precision',6,'newline', 'pc') % '%.6f'
%     txt_format = [speichername,num2str(i),'.txt'];
%     save(txt_format,'MM','-ascii')

%     figure;
%     quiver(x,y,u,v,1.5),
%     title(['PIV velocity','(Bild Nummer : '...
%         num2str(i),')  Zenter: x0=' num2str(x0),'  y0=' num2str(y0)]),
%     hold on;    
%     %inner cylinder
%     xzyl=[-r_i:0.75:+r_i];
%     yzyl=sqrt(r_i^2-xzyl.^2);
%     area(xzyl,yzyl,0,'facecolor','black');
%     yzyl=-sqrt(r_i^2-xzyl.^2);
%     area(xzyl,yzyl,0,'facecolor','black');
%     %outer cylinder
%     xzyl=[-r_a-10:0.75:r_a+10];
%     yzyl=sqrt(r_a^2-xzyl.^2);
%     area(xzyl,yzyl,r_a+10,'facecolor','black');
%     yzyl=-sqrt(r_a^2-xzyl.^2);
%     area(xzyl,yzyl,-r_a-10,'facecolor','black');    
%     DataAspectRatio=[1 1 1];
%     hold on;
%     plot(0,0,'ro','Markersize',10),
%     axis([-160 160 -r_a-10 r_a+10]);
%     
%     saveas(gcf,[speichername,num2str(i),'.jpg'],'jpg');
%     F(i+1)=getframe(gcf);

    fclose('all');
    clear all;
    close all;
end;
%movie2avi(F,[speichername,'movie.avi'],'fps',4);