Creating a Cell to save n*n matrices

Hi All Basically I have a structure which I am carrying out a sensitivity analysis on. I have 7 accelerogram data sets which I am looping through. I am also changing the stiffness and damping parameters for the dampers used and I am calculating the peak displacements.

So I have saved the peak displacements as an n*n cell for the corresponding stiffness and damping values. Is there a way to save this matrix within another matrix or cell so that I have 7 sets of these matricies? One set for each Record which I am using?

I have tried 3D matrices but it is not quite what I need- unless I have implemented it wrong. I am also struggling with where to place the cell so it saves the data correctly.

Hope you can help and Thanks in advance!

str='E:\01 LOUGHBOROUGH UNI WORK\2017-2018 4th Year\FINAL YEAR\DISSERTATION\SERENA SOLANKI DISS\CODING\california earthquake record';
textfiles=dir([str '*.txt']);

ft=[];
Peak_disp=[];

t_interp1=[];

for num=1:1:7
    textfilename=['Record_', num2str(num), '.txt']; %read Record_ numbers 1-7
    %open textfile 'rt' is for reading a text file only
    FileID=fopen(textfilename,'r');

      acceleration_data=0.01*fscanf(FileID,'%f');
      fclose('all');
      ft{num}=m*acceleration_data.'; %force data is saved in a cell

      kiso_range=0.05:0.05:0.9; %changing the base isolator stiffness of the structure
      iso_range=size(kiso_range,2);
      ziso_range=1:4:73;
      ziso_r=size(ziso_range,2);

      for index=1:iso_range
          disp(index);

          for cc=1:ziso_r

              M=[m 0;0 m];
              K= [k    -2*kiso        
                  -2*kiso  2*kiso+k  ];
              kiso=kiso_range(index)*k/2;
              zeta=0.05;
              zita=[0.05;(zeta*cc);0.05;(zeta*cc);0.05;(zeta*cc)];
              [om,Phi,C]=modan(M,K,zita);
              tau=[1  1  1  1  1  1]';

              b=[0;(-1)]; %force vector
              A=[a  b, c  d];

              dt    = 0.02; %seconds
              NPTS = size(ft{num},2);
              tf    = (NPTS-1)*dt;
              t = 0:dt:tf;

              dt_interp = 0.005;
              t_interp1{num} = dt_interp:dt_interp:tf;
              t_interp=dt_interp:dt_interp:tf;
              ft_interp = F_interp(t, ft{num}, t_interp);

              NPTS_interp = size(t_interp,2);

              yt = zeros(12,NPTS_interp);

              for i=1:(NPTS_interp-1) % Runge-Kutta loop, i = incremental time step      
                  g1 = A*yt(:,i)+b*ft_interp(i);
                  g2 = A*(yt(:,i)+ 0.5*g1*dt_interp)+0.5*b*(ft_interp(i)+ft_interp(i+1));
                  g3 = A*(yt(:,i)+ 0.5*g2*dt_interp)+0.5*b*(ft_interp(i)+ft_interp(i+1));
                  g4 = A*(yt(:,i)+g3*dt_interp)+b*ft_interp(i+1);

                  ynew = yt(:,i)+(g1+2*g2+2*g3+g4)*dt_interp/6;
                  yt(:,i+1)=ynew; %storage
              end
              % yt= displacement, the  indicates the floor number 

              Peak_mat(index,cc)=max(abs(yt(2,:)-yt(1,:))*0.01); 

          end
      end
      Peak_disp{num,index,cc}=max(abs(yt(2,:)-yt(1,:))*0.01); % saving peak displacements for each data set 
  end

[Bob Thompson]

Since I didn't write the code I'm going to need a bit of help with your variable names. If I am interpreting what you're asking correctly then you want a 1x7 (or 7x1) cell array where the cells contain the nxn matrices of values you calculated in the inner loop?

To the best of my knowledge cells can contain any class of items, so you should be able to just define the cell contents as the matrix using indexing.

Peak_disp(num) = {Peak_mat}; % Parentheses around num indicates that we are
                             % defining the entirety of num element. Curly
                             % braces around Peak_mat defines the num element
                             % as a cell containing the Peak_mat matrix.

You shouldn't need to use 'index' or 'cc' in your indexing of Peak_disp, as these values are a fixed value from loops that have already ended.

On a side note, you can loop 'num' using num = 1:7. MATLAB defaults to increasing by 1 when no specific interval is defined.

[serena solanki]

Hi Bob

Thanks for the tip, after a little trial and error I found that

Peak_disp{num}=Peak_mat; 

gave me the answer I needed.

I wanted to create a 1*7 cell which contained the n*n matrices which refer to cc and index- which are for the values calculated in the inner loop.

Again thanks for the help