Question

我有50个观察值作为具有变量值的行和列。对于每个观察，我使用fsolve来求解未知变量q_iz。我正在尝试使用LOOP。我已经为1 ^st行解决了这个问题，但是当MATLAB进入2 ^nd行时会抛出错误。另外，我想将k的求解值存储在(50,1)维度中。

我提出了一个更简单的代码来表示问题，如下所示：

 global price gnw dataq_izc(c) c

price = ones(50,2);
gnw = ones(50,2,2);
dataq_izc = 0.20*ones(50,1);

k0 = 0.1+ zeros(50,1);
for c = 1:50
options=optimset('Display','notify','MaxFunEvals',10^7,'TolX',1e-8,'TolFun',1e-6,'MaxIter',10^5);
% Solve 
k = fzero(@qiz,k0(c,1),options);
fprintf('Quantity\n');
disp(k);
end;

脚本'qiz'是：

function F = qiz(k)

global price gnw dataq_izc(c) c

q_iz = zeros(50,2,2);
i = 1.;
while i < 3
    z = 1.;
    while z < 3
    q_iz(c,i,z) = k(c,1)*price(c,z)/gnw(c,i,z);
    z = z + 1.;
    end
i = i + 1.;
end;

q_izc = zeros(50);
c = 1.;
while c < 51;
i = 1.;
while i < 3
    z = 1.;
    while z < 3
    q_izc(c) = q_izc(c)  + q_iz(c,i,z);
    z = z + 1.;
    end
i = i + 1.;
end;
c = c + 1.;
end

F = q_izc(c) - dataq_izc(c,1);

当我运行这段代码时，它会在第一行解决，但是当它转到第二行时会抛出错误：

FZERO无法继续，因为用户提供的function_handle ==＆gt; QIZ 因以下错误而失败。试图访问k（2,1）;索引界限因为numel（k）= 1。

Answer 1

您已将标量k0(c,1)传递给您的函数qiz。在该功能中，您尝试访问k(c,1)并获得错误说明。

FZERO无法继续，因为用户提供的function_handle ==＆gt; qiz因以下错误而失败。试图访问k（2,1）;因为数字（k）= 1，所以索引越界。

错误告诉您k是numel(k)=1的标量。这是第一个错误。第二个是在qiz函数中你重新分配c的值。您将c定义为global，然后将其用作第二个循环中的索引。当qiz结束时，c的值现在为51.这超出了k0向量的范围。尽量使用能够避免此错误的变量名称。

这是一个工作（无错误）版本。

function k = Test(  )

global c price gnw dataq_izc

price = ones(50,2);
gnw = ones(50,2,2);
dataq_izc = 0.20*ones(50,1);

k0 = 0.1+ zeros(50,1);
k = size(k0);
for c = 1:50
options=optimset('Display','notify','MaxFunEvals',10^7,'TolX',1e-8,'TolFun',1e-6,'MaxIter',10^5);
% Solve 
k(c) = fzero(@qiz,k0(c),options);
fprintf('Quantity: %f\n', k(c));
end

end


function F = qiz(k)

global c price gnw dataq_izc

q_iz = zeros(50,2,2);
i = 1.;
while i < 3
    z = 1.;
    while z < 3
        q_iz(c,i,z) = k*price(c,z)/gnw(c,i,z);
        z = z + 1.;
    end
    i = i + 1.;
end;

q_izc = zeros(50);
c2 = 1.;
while c2 < 51;
    i = 1.;
    while i < 3
        z = 1.;
        while z < 3
            q_izc(c2) = q_izc(c2) + q_iz(c2,i,z);
            z = z + 1.;
        end
        i = i + 1.;
    end;
    c2 = c2 + 1.;
end

F = q_izc(c) - dataq_izc(c,1);

end

多次求解一个方程

1 个答案: