making a function work
11 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Richard
el 16 de Feb. de 2012
I've written the following function which imports .txt files into matlab. The .txt files can either be recorded at hourly intervals or 4minute intervals where depending on the initial resolution the script will calculate the hourly or daily averages.
function [Daily, Hourly] = calc_avg(pathName)
TopFolder = pathName;
dirListing = dir(fullfile(TopFolder,'*.txt'));%Lists the folders in the directory specified
%by pathName.
for i = 1:length(dirListing);
fileToRead1{i} = (dirListing(i).name);
%lists all of the .txt files in the TopFolder
end
cell_all = arrayfun(@(i1)importdata(fullfile(pathName,dirListing(i1).name)),...
(1:length(dirListing))','un',0);
%apply function to each element of the array.
d = cat(2,cell_all{:});
%concatenate arrays along each column (i.e. 2)
%find the length of the dataset, which will provide information on the
%amount of averaging required.
if length(d) == 365,...
error('error: daily averages already calculated');
elseif length(d) == 8760;
daily = squeeze(mean(reshape(d,24,size(d,1)/24,[])));
else length(d) == 131400;
hourly = squeeze(mean(reshape(d,15,size(d,1)/15,[])));
daily = squeeze(mean(reshape(d,360,size(d,1)/360,[])));
end
%find which averages have been calculated:
A = exist('hourly','var');
%if A == 1 it means that hourly values had to be calculated therefore
%the data if of high resolution (minutes).
if A == 1;
hourly = mat2cell(hourly,size(hourly,1),cellfun('size',cell_all,2)).';
daily = mat2cell(daily,size(daily,1),cellfun('size',cell_all,2)).';
elseif A == 0;
daily = mat2cell(daily,size(daily,1),cellfun('size',cell_all,2)).';
end
%create cell in the same format as 'cell_all' where cellfun applies the
%same function to each cell in a cell array. 'size' is used to create
%the same format.
for i=1:length(dirListing);
[~,name{i}] = fileparts(fileToRead1{i});
%obtain the name of each of the .txt files (dirListing)
end
%Generate a structure for the averages calculated.
if A == 1;
for i=1:length(dirListing);
Daily.(genvarname(name{i})) = daily{i};
Hourly.(genvarname(name{i})) = hourly{i};
end
elseif A == 0;
for i=1:length(dirListing);
Daily.(genvarname(name{i})) = daily{i};
end
end
The script works fine if I simply run it as a script i.e. avoid using the function and just type the path Name into the second line.
What am I missing which would make this work as a function?
0 comentarios
Respuesta aceptada
Más respuestas (1)
Honglei Chen
el 16 de Feb. de 2012
Maybe the path is messed up, what is 'ThePathName'? I'd suggest to put calc_avg on the path, then cd to the directory where the files are located, and call the following
calc_avg(pwd)
to see if it fixes the issue.
3 comentarios
Honglei Chen
el 16 de Feb. de 2012
As Jan mentioned in his comment, in this case Hourly is not defined. Therefore you need to put something default for Hourly if A is 0. Maybe add something like Hourly = [] or Hourly = nan;
Ver también
Categorías
Más información sobre Logical en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!