现有电子表格模型的Matlab源代码

debugcn 发表于 Dev

知道

我正在mat实验室做一个项目，我正在处理大量数据。我需要为电子表格中现有的模型创建一个源代码，其中计算已在电子表格中逐行完成。我在编写电子表格的最后5列的代码时遇到问题，其中第一列的第二行的结果将继续进行使用第五列第一行的结果，依此类推。有关信息，

dh / dx是第一列
GWFLOW是第二栏
GWDROP是第三列
NETGWCHANGE是第4列
REVGW是第五列

Excel文件的屏幕截图，可清晰查看

这是我做的代码：

     %% dh/dx
     dh_dx=zeros(numel(rocky.RAIN),1);
     % for the first row
     dh_dx(1)=(ingwl-effecbase)/meandx;
     %for the remaining dh/dx
        for s=2:numel(HER);
            dh_dx(s)=(REVGW(s-1)-effecbase)/meandx;
        end

     %% GW flow for unit width
     GWFLOW=zeros(numel(rocky.RAIN),1);
     for t=1:numel(HER);
         GWFLOW(t)=K*dh_dx(t)*aqb;
     end

   %% Drop of GW due to GW flow
   GWDROP=zeros(numel(rocky.RAIN),1);
   for u=1:numel(HER);
       GWDROP(u)=(GWFLOW(u)/(aqb*por));
   end

   %% Net Groundwater change
   NETGWCHANGE=zeros(numel(rocky.RAIN),1);
   for v=1:numel(HER);
       NETGWCHANGE(v)=INCGW(v)-GWDROP(v);
   end

   % Revised GW level Estimate
   REVGW=zeros(numel(rocky.RAIN),1);
   %for the first row
   REVGW(1)=ingwl+NETGWCHANGE(1);
   for w=2:numel(HER);
       REVGW(w)=REVGW(w-1)+NETGWCHANGE(w)
   end

先感谢您

用户名

问题在于，您需要在代码中单独计算所有数据，而需要在同一循环中一起计算所有数据：

    % Allocate data
    dh_dx       = zeros(numel(rocky.RAIN),1);
    GWFLOW      = zeros(numel(rocky.RAIN),1);
    GWDROP      = zeros(numel(rocky.RAIN),1);
    NETGWCHANGE = zeros(numel(rocky.RAIN),1);
    REVGW       = zeros(numel(rocky.RAIN),1);

    % Seed the calculation on the first row
    dh_dx(1)       = (ingwl-effecbase)/meandx;
    GWFLOW(1)      = K*dh_dx(1)*aqb;
    GWDROP(1)      = (GWFLOW(1)/(aqb*por));
    NETGWCHANGE(1) = INCGW(1)-GWDROP(1);
    REVGW(1)       = ingwl + NETGWCHANGE(1);

    % Loop trough the other rows and calculate fields
    for s = 2:numel(HER)
            dh_dx(s)       = (REVGW(s-1)-effecbase)/meandx;
            GWFLOW(s)      = K*dh_dx(s)*aqb;
            GWDROP(s)      = GWFLOW(s)/(aqb*por);
            NETGWCHANGE(s) = INCGW(s) - GWDROP(s);
            REVGW(s)       = REVGW(s-1) + NETGWCHANGE(s);
    end;

本文收集自互联网，转载请注明来源。

如有侵权，请联系[email protected] 删除。