如何修复矩阵维数的一致性

时间:2013-05-16 20:01:42

标签: matlab matrix matlab-guide

我是MATLAB的初学者。我创建了一个代码,用于读取类型为JPG(相同文件类型)但每个文件大小不同的数据集文件。我做了什么我随机读取了10个文件并将这些文件的内容放入一个名为'fr'的矩阵(带有单元格数组的类型)中(如下所示)。这个'fr'包含10列,这意味着10个文件被读取。我提取(读取)的文件内容是二进制文件字节(0-255)但是,我提取了每个文件的内容并将它们放入一个名为'out1'的变量中,但问题是如何使用不同的文件访问这些内容每个文件的大小,我已经通过使用函数(cellfun)解决了它,我通过添加零填充空白单元格。 最终,矩阵工作,但我现在面临的问题是如何从读取的整个文件中减去这个矩阵'out1'?我坚持这一部分。这个错误表明矩阵维度必须达成一致:如最后所述。

这是输出:     ???使用==>时出错减去     矩阵维度必须一致。

Error in ==> PCATEST2 at 50
B = (out1 - repmat(AMean,[n 1])) / repmat(AStd,[n 1]);

EDU>> whos
  Name                Size              Bytes  Class     Attributes

  AMean               1x10                 80  double              
  AStd                1x10                 80  double              
  ans                 1x1                   8  double              
  f                   1x57                114  char                
  fid                 1x1                   8  double              
  files              50x1               32870  struct              
  fr                  1x10            6715472  cell                
  i                   1x1                   8  double              
  j                   1x1                   8  double              
  m                   1x1                   8  double              
  maxLength           1x1                   8  double              
  n                   1x1                   8  double              
  out1           101077x10            8086160  double              
  ridx               50x1                 400  double              
  st1                 1x10               1460  cell

这段代码:

f ='/Users/nsa/Documents/MATLAB/jpg-data1/';
files = dir(fullfile(f,'*.jpg'));
ridx = randi(numel(files),size(files)); %choose files from f randomly

 for i = 1:10  %randomly pick 10 files
 st1(i) = {files(ridx(i)).name}; 

  for j= i:length(st1)

    fid = fopen(fullfile(f, st1{i}),'rt');

    fr{j} = fread(fullfile(fid));

    fclose(fid); 

 end

end
maxLength=max(cellfun(@(fr)numel(fr),fr));
out1 = cell2mat(cellfun(@(fr)cat(1,fr,zeros(maxLength-length(fr),1)),fr,'UniformOutput',false));

[n m] = size(fr);

AMean = cellfun(@mean,fr);

AStd = cellfun(@std,fr);

B = (out1 - repmat(AMean,[n 1])) / repmat(AStd,[n 1]);


??? Error using ==> minus
Matrix dimensions must agree.

Error in ==> TEST2 at 50
B = (out1 - repmat(AMean,[n 1])) / repmat(AStd,[n 1]);

任何人都可以帮助我并指导我修复并克服此错误吗?

1 个答案:

答案 0 :(得分:1)

你应该这样做:

[n m] = size(out1); 

因为fr只有1行而out1只有一行。

但是,您甚至可以使用bsxfun()

跳过重新匹配
B = bsxfun(@rdivide, bsxfun(@minus,out1,AMean), AStd);