MATLAB:不同长度的向量

MATLAB:不同长度的向量,matlab,octave,Matlab,Octave,我想创建一个MATLAB函数,从另一个目录中的文件导入数据,并将其适配到给定的模型中,但因为数据需要过滤(文件中的不同位置有“抖动”数据,例如,在分析运动开始之前没有任何测量) 因此,包含用于拟合的数据的向量最终具有不同的长度,因此我无法在矩阵中返回它们(例如,下面的函数中的x)。我怎样才能解决这个问题 我有很多数据文件,所以我不想使用“手动”方法。我的功能如下。欢迎大家提出意见和建议 datafit.m function [p, x, y_c, y_func] = datafit(patter

我想创建一个MATLAB函数,从另一个目录中的文件导入数据,并将其适配到给定的模型中,但因为数据需要过滤(文件中的不同位置有“抖动”数据,例如,在分析运动开始之前没有任何测量)

因此,包含用于拟合的数据的向量最终具有不同的长度,因此我无法在矩阵中返回它们(例如,下面的函数中的x)。我怎样才能解决这个问题

我有很多数据文件,所以我不想使用“手动”方法。我的功能如下。欢迎大家提出意见和建议

datafit.m

function [p, x, y_c, y_func] = datafit(pattern, xcol, ycol, xfilter, calib, p_calib,    func, p_0, nhl)

    datafiles = dir(pattern);
    path = fileparts(pattern);
    p = NaN(length(datafiles));
    y_func = [];
    for i = 1:length(datafiles)
        exist(strcat(path, '/', datafiles(i).name));
        filename = datafiles(i).name;
        data = importdata(strcat(path, '/', datafiles(i).name), '\t', nhl);
        filedata = data.data/1e3;
        xdata = filedata(:,xcol);
        ydata = filedata(:,ycol);
        filter = filedata(:,xcol) > xfilter(i);
        x(i,:) = xdata(filter);
        y(i,:) = ydata(filter);
        y_c(i,:) = calib(y(i,:), p_calib);
        error = @(par) sum(power(y_c(i,:) - func(x(i,:), par),2));
        p(i,:) = fminsearch(error, p_0);
        y_func = [y_func; func(x(i,:), p(i,:))];
    end
end

示例数据:

我可以想到两种策略:

  • 我将返回一个数组中的数据,每个单元格存储不同长度的向量。您可以使用与数组相同的方式访问数据,但使用大括号:例如c{1}=[1 2 3],c{2}=[1 2 10 8 5]c{3}=[]
  • 如果使向量具有相同的长度,还可以在读取行时过滤垃圾数据

如果内存不是主要问题,请尝试使用不同的值填充向量,例如NaN或Inf-anything,这些值在基于物理上下文的测量中找不到。在为矩阵(*)分配内存之前,可能需要确定最长的数据集。通过这种方式,您可以使用大小相同的矩阵,并且在以后很容易忽略“空数据”

(*)想法。。。首先根据最大文件的大小分配内存。把它装满,比如楠的

matrix = zeros(length(datafiles), longest_file_line_number) .* NaN;
然后运行你的函数。确定最长连续数据集的长度

new_max = length(xdata(filter));
if new_max > old_max
    old_max = new_max;
end
matrix(i, length(xdata(filter))) = xdata(filter);
在函数返回矩阵之前,相应地裁剪矩阵

matrix = matrix(:, 1:old_max);

也许也会发布一些示例数据。