Matlab for 循环语句

发布网友 发布时间：2022-04-20 09:49

我来回答

共3个回答

热心网友 时间：2023-07-17 19:23

没仔细研究你那个拟合函数的用法。

不过像这种需要分段处理数据的情况很多，有种在 matlab 里很常用的技巧感觉你可以学学：

假设原始数据（xdata）是一列 100 个数，你需要一次处理 13 个，那么下面这段代码先将这 100 数“折叠”成一个 13 行的矩阵（x）。 我这个例子里是折叠成 13 * 8 的矩阵了。

因为 100 不能整除 13，所以会在最后补零。如果你不需要补零，而是把最后的零头扔掉，那么就把 

ncol = ceil(nxdata/n);

改成 

ncol = fix(nxdata/n);

最后就可以在循环里一列一列的处理数据了。我这个例子只是 y(:,i) = x(:,i) 而已，你的情况应该就是调用你自已写的拟合函数。

最后再用

y = y(:);

把数据还原成一维的。

如果你需要把之前补零的值删掉也可以在还原成一维的时候去掉。

以下是一个小例子，仔细研究一下就能学会这个小技巧了。

xdata = (1:100).';
nxdata = length(xdata);

nrow = 13;  % 假设你需要一次处理 13 个数据
ncol = ceil(nxdata/nrow);
x = zeros(nrow, ncol, 'like', xdata);
size(x)
x(1:nxdata) = xdata(:);

y = zeros(nrow, ncol, 'like', xdata);

for i = 1:ncol
    y(:,i) = x(:,i);
end

y = y(:);

如果你在处理数据之前，能保证折叠之后是整除的，即无需补零。还用我这个小例子，比如需要折叠成 10 行，也可以用 reshape(xdata, 10, [])。[] 是让 matlab 自己算整除后是多少列。如：

x2 = reshape(xdata, 10, []);
size(x2)

热心网友 时间：2023-07-17 19:24

x(:,i)，y(:,i)改成x{1，i}，y{1，i},下面的也一样，你试试

热心网友 时间：2023-07-17 19:24

x=x(:,i);y=y(:,i);n=5;
ovm(:,i)=BcreateFit(x,y,n);
改成：
u=x(:,i);v=y(:,i);n=5;
ovm(:,i)=BcreateFit(u,v,n);
试试