Matlab for 循环语句

Question

就是采用分段拟合拟合的方法，每次都用5阶傅立叶进行函数拟合，一共86400个数据，每次取2800个数据，需要循环30次，即调用30次，这循环语句如何写。附上函数代码：
function ovm=BcreateFit(x,y,n);
xdata=x(:);
ydata=y(:);
ft=fittype(['fourier',num2str(n)']);
options=fitoptions(ft);
opstions.Display='Off'
[fitresult,gof]=fit(xdata,ydata,ft,options);
ovm=fitresult(xdata);
end
自己写的代码是这样的：
for i=1:30
x(:,i)=xdata(2800*(i-1)+1:2800*i);y(:,i)=ydata(2800*(i-1)+1:2800*i);
x=x(:,i);y=y(:,i);n=5;
ovm(:,i)=BcreateFit(x,y,n);
plot(x(:,i),y(:,i),'r');
end
但只调用了一次，即i=1.求大神指导，深表感谢！

Answer 1

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

不过像这种需要分段处理数据的情况很多，有种在 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)

Answer 2

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

Answer 3

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);
试试本回答被网友采纳