MATLAB 中的分段函数

Question

我今年才开始使用 MATLAB for Uni，作为一个家庭作业测验问题，我被要求创建一个分段函数（我称之为“piecewise_method”），它能够在“x”上建立一个不同的方程" 取决于 "x" 是低于 0、介于 0 和 8 之间还是高于 8。这是我目前编写的代码。

function solution = piecewise_method(x)

% given the value of the input x, the function
% piecewise_method will choose from one of the three specified
% equations to enact upon the given value of x
% and give the user a solution

solution = zeros(size(x));
e = exp(1);

for j = 1:size(x)
    a = x(j);

    if a < 0
        solution(j) = -a.^3 - 2*a.^2 + 3*a;
    elseif (a >= 0) && (a <= 8)
        solution(j) = (12/pi)*sin(pi*a./4);
    else
        solution(j) = ((600.*e.^(a-8))./(7*(14+6.*e.^(a-8))) - 30/7);
    end
end

使用输入运行时...

x = -3:12

它为变量解产生这个结果...

解决方案 =

 0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0

现在这向我表明该数组正在正确创建，但由于某种原因，for 循环没有正常运行，或者没有按预期运行。我尝试从基本级别开始多次重建 for 循环，但是当我开始输入方程式时，它又开始分崩离析，所以我相信我的方程式可能有问题（这就是为什么我到处都放括号的原因，以防万一）。

问题还要求我使用if语句，所以我无法尝试使用其他方法来产生分段方法函数，并且从我的搜索来看，似乎没有很多分段函数中if语句的例子.

如果您能提供任何可以帮助我使用此功能的建议，将不胜感激，谢谢！

附：如果您有任何建议可以在将来改进我的问题，那也很好！

Answer 1

您应该在for 循环中使用length 而不是size。

size 函数的输出是 x 的维度，对于您的 x=-3:12 示例，它返回 size(x)=[1 16]。然后您的for 循环将为j=1:size(x) 运行，即j=1:1，即j=1。

length 的输出是x 的最大数组维度的长度，如here 所列。在您的示例中：length(x) = 16，然后是 j=1:length(x)=1:16。

或者，您可以使用size(x,2)，它将返回x 的第二维的大小，在这种情况下与length 相同。

Answer 2

如 cmets 和 RadioJava's answer 中所述，您需要检查循环。我会使用size(x,2)（x 中的列数）或numel(x)（x 中的元素数）。

for j = 1:numel(x)
    % ...

其他人推荐length(x)，即max(size(x))。我通常会尽量避免这种情况，因为它不明确您要查找矩阵的哪个维度。

除了只是重复信息之外，我想展示您可以通过逻辑索引更有效地做到这一点，并完全消除有问题的循环......

function solution = piecewise_method(x)

% given the value of the input x, the function
% piecewise_method will choose from one of the three specified
% equations to enact upon the given value of x
% and give the user a solution

solution = zeros(size(x));
e = exp(1);

idx = (x<0);
solution(idx) = -x(idx).^3 - 2*x(idx).^2 + 3*x(idx);
idx = (x>=0 & x<=8);
solution(idx) = (12/pi)*sin(pi*x(idx)/4);
idx = (x>8);
solution(idx) = (600*e.^(x(idx)-8))./(7*(14+6*e.^(x(idx)-8))) - 30/7;

请注意，您可以轻松地在数组中包含条件和（匿名）函数并循环它们以使其更加灵活。