0

我想创建一个 MATLAB 函数来从另一个目录中的文件导入数据并将它们拟合到给定模型，但是因为需要过滤数据（文件中的不同位置有“thrash”数据，例如之前没有测量分析的运动开始）。

因此，包含用于拟合的数据的向量最终具有不同的长度，因此我无法在矩阵中返回它们（例如，下面的函数中的 x）。我该如何解决这个问题？

我有很多数据文件，所以我不想使用“手动”方法。我的功能如下。欢迎所有和建议。

数据拟合.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

样本数据： http ://hastebin.com/mokocixeda.md

4

2 回答 2

1

There are two strategies I can think of:

  • I would return the data in a vector of cells instead, where the individual cells store vectors of different lengths. You can access data the same way as arrays, but use curly braces: Say c{1}=[1 2 3], c{2}=[1 2 10 8 5] c{3} = [ ].
  • You can also filter the trash data upon reading a line, if that makes your vectors have the same length.
于 2012-12-01T15:47:13.597 回答
0

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

(*) 想法...首先根据最大文件的大小分配内存。用例如 NaN 填充它

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);
于 2012-12-04T14:47:46.050 回答