Y_FFT3D_spectra_2D_peak_line_damped-zhonglu°Øs MacBook Pro.m

%% FFT 3D spectra
%% 2D peak line fft
%%
setlatex

%cd G:\orange_backup % any dir
%cd /Volumes/GROUP_BLACK/orange_backup/
clear all
fclose all;
%Gi=[0.05 0.1 0.2 0.3 0.4 0.5 1 1.25 3]; %All the G

%Gi=[0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0];
Gi=[1.0];
mi=[1.5 1.7 2.0 2.2 2.5];
ri='*';

Gi=0.2
mi=2.0
ri=1177

for i3=Gi
    for i4=mi

        close all
        clearvars -except i3 i4 Gi mi ri
        %sys=computer;
        sysdir
        gstr=num2str(i3,['%10.1','f']);
        mstr=num2str(i4,['%10.1','f']);
        rstr=num2str(ri);
        list0=dir(['g',gstr,'*m',mstr,'*r',rstr]);
        if isempty(list0); continue; end
        f1pmax=2.5;
        %A1din = input('input: digits after decimal point for A_1 legend? (default is 3 digits)','s');
        %if size(A1din)==0;
            A1din= '3';
        %end

        maximumA1=inf;

        A2switch=false; % A2=true or A2dA1=false

        mrswitch=0; %mass ratio is a variable or not

        n1=1;%starting number of files
        disp(['starting number= ',num2str(n1)])

        n0=length(list0);
        s=0;sp=0;
        countl=0;
        clear gamall aryaxismax aryaxismin arxaxismax arxaxismin
        
        for ai=n1:n0

            cd (list0(ai).name)
            list=dir('0*');
            if isempty(list); continue; end
            n=length(list);
            

            str=list0(ai).name;
            str(str=='r')=[' '];
            str(str=='u')=[];
            str(str=='n')=[];
            str(str=='_')=[];
            str(str=='z')=[];
            str(str=='l')=[' '];
            str1=str;
            str(str=='g')=' ';
            str(str=='a')=' ';
            str(str=='m')=' ';
            gam=str2num(str);
            gamall(1:length(gam),ai)=gam;
            A1=gam(2);
%             if A1>maximumA1
%                 continue
%             end
            countl=countl+1;
            A2m=[];f2m=[];f1temp=[]; %initialise A2m and f2m
            clear A2dA1 A2p A2peak f1 f2p f2peak h
            for f1i=1:n
                dirnamef=[list0(ai).name,'/',list(f1i).name,'/'];
                disp(dirnamef);
                f1(f1i)=str2double(list(f1i).name);
                T1=1/f1(f1i);

    %             if A1==0.1 && f(i)==0.4
    %                 stop
    %             else
    %                 continue
    %             end

                cd(list(f1i).name)
                [t,y]=textscanty;
                cd ../
                [t,y]=autocut(f1(f1i),t,y,t,y);
                [f2, A2]=fftzp2(t,y);
                %[~,~,A2, f2]=Contfft(t,y,);
                
                f2(f2>f1pmax)=[];
                A2=A2(1:numel(f2));
                A2m(f1i,1:numel(A2))=A2; f2m(f1i,1:numel(f2))=f2;
                
                %if A2switch==false;A2dA1(f1i)=A2p(f1i)/A1;end
                clear ymax ymin yamp yampds t v02 y;
            end
            %% 3D spectra
            
          
            Lf1=length(f1);
            figure(1)
            %waterfall(f1,[],A2m)
            hold on
            for itemp=1:Lf1
                f1temp(1:length(f2m(itemp,:)))=f1(itemp);
                h(itemp)=plot3(f2m(itemp,:),f1temp,A2m(itemp,:),'color','k');
                clear f1temp
            end
            grid on; %grid minor
            
            xlim([0 2.5]);ylim([0 2.5])
            xlabel('$f_2$','Interpreter','LaTex','FontSize',20);
            ylabel('$f_1$','Interpreter','LaTex','FontSize',20);
            zlabel('$A_2$','Interpreter','LaTex','FontSize',20);
            
            campos([12.5839  -11.2161    0.5939])
            [A2peak, IA2peak]=max(A2m,[],2);
            for itemp=1:Lf1
                f2peak(itemp)=f2m(itemp,IA2peak(itemp));
            end

            hp3=plot3(f2peak,f1,zeros(1,Lf1),'color','r','LineStyle','--');
            %zlim([0 2.5])
            cd ../
            %
            
            temp=['Y_FFT_spectra_fig/',list0(ai).name];
            mkdir(temp)
            cd (temp)
            temp=[list0(ai).name,'_Y_3DFFT'];
            %printpdf(gcf,['displacement_verification_','f1_',num2str(f1),'.eps']);
            printpdf(gcf,[temp,'.eps']);
            saveas(gcf,[temp,'.fig'])
            %% 2D f1-f2
            hold off
            close all
            figure(1)
            hp3=plot(f1,f2peak,'color','r');
            grid on; grid minor;
            %set(gca,'XTick',0:0.2:2.5);
            xlabel('$f_1$','Interpreter','LaTex','FontSize',20);
            ylabel('$f_2$','Interpreter','LaTex','FontSize',20);
            setfigure;

            printpdf(gcf,[list0(ai).name,'_2Df1_f2peak.eps']);
            close all
            cd ../;cd ../;
            
        end


    end
end