fminsearch 在存储数据时循环覆盖

Question

我在 matlab 中有下面提到的代码。我想将所有计算出的 162 行和 4 列写入一个 excel 文件。

当我在代码中使用 xlswrite 时，我只得到一行和 4 列，因为 P 的值在每个迭代步骤中都被覆盖。

如果我在 for 循环中使用另一个循环，执行时间会急剧增加。请帮助将 P 的值至少写入一个数组，我稍后可以将其写入 excel 文件（当我尝试"在赋值 A（I） = B 中，B 和 I 中的元素数量必须相同"错误出现。

请帮忙

function FitSMC_BC
clc
% Parameters: P(1)=theta_S; P(2)=theta_r; P(3)=psib; P(4)=lamda; 
smcdata=xlsread('asimdata');
nn=length(smcdata)-1;
for i=1:nn
    psi=smcdata(:,1);
    thetaObs=smcdata(:,i+1);
    %Make an initial guess:
    Pini=[0.5 0.1 -1 1.5];
    P=fminsearch(@ObFun,Pini,[],psi,thetaObs);
    disp(['result',num2str(i),': P=',num2str(P)]);
    theta=Gettheta(P,psi);
end
function OF=ObFun(P,psi,thetaObs)
theta=Gettheta(P,psi);
OF=sqrt(mean((theta - thetaObs).^2));
function theta=Gettheta(P,psi)
SoilPars.theta_S=P(1);
SoilPars.theta_r=P(2);
SoilPars.psib=P(3);
SoilPars.lamda=P(4);
[theta]=thetaFun(psi,SoilPars);
function [theta]=thetaFun(psi,SoilPars)
theta_S=SoilPars.theta_S;
theta_r=SoilPars.theta_r;
psib=SoilPars.psib;
lamda=SoilPars.lamda;
theta=theta_r+((theta_S-theta_r)*((psib./psi).^lamda));
theta(psi>psib)=theta_S;

Answer 1

你可以用

P(i,:) = fminsearch(@ObFun,Pini,[],psi,thetaObs);

P将每个计算（4 个元素向量）存储在新行中。

您也可以在for循环之前初始化P P = nan(nn, 4);

然后使用 xlswrite 在 Excel 文件中写入P。

Answer 2

我没有深入研究过你的代码，但据我所知，你有两个选择：

1. 创建一个矩阵P并在整个矩阵上使用xlswrite。在我看来，这是最合理的方法。
2. 在循环中使用文件交换中的xlswrite1。这将稍微增加执行时间，但不会像使用常规xlswrite那样多，因为它专门设计用于循环内部。它之所以如此快，是因为它只打开和关闭一次 Excel 文件，而常规xlswrite每次调用函数时都会打开和关闭它。

Answer 3

你似乎知道如何使用索引，所以我不确定你为什么只是做这样的事情：

P = zeros(size(smcdata,1),nn)
for i=1:nn
    ...
    P(:,i) = fminsearch(@ObFun,Pini,[],psi,thetaObs);
    disp(['result',num2str(i),': P=',num2str(P(:,i))]);
    theta = Gettheta(P(:,i),psi); % Why is this here? Are you writing it to file too?
end
xlswrite('My_FileName.xls',P);

或者，您可以在循环的每次迭代中调用xlswrite（可能更慢），并使用如下所示的内容附加新数据：

for i=1:nn
    ...
    P = fminsearch(@ObFun,Pini,[],psi,thetaObs);
    disp(['result',num2str(i),': P=',num2str(P)]);
    theta = Gettheta(P,psi); % Why is this here? Are you writing it to file too?
    xlswrite('My_FileName.xls',P,1,['A' int2str((i-1)*size(P,2)+1)]);
end

当然，您的代码不可运行，因此您必须调试任何其他小错误。此外，由于smcdata似乎是矩阵而不是向量，因此您应该小心使用它length。您可能应该使用 size .