Scilab中基于向量变量定义分段函数的问题求助

问题分析 #

• 第一个函数的问题：if语句直接用向量和标量比较时，Scilab仅判断向量的第一个元素，无法对整个向量逐元素判断，导致结果不符合预期。
• 第二个函数的问题：函数定义不能嵌套在for循环内部，且t(i)的索引逻辑错误，混淆了函数参数t和循环变量i。

正确的向量化实现 #

Scilab支持向量化操作，无需循环即可高效处理整个向量。以下是修正后的simula_vg函数：

function vg=simula_vg(t,Tcg,Tfg,Ag)
    // 初始化结果向量，默认值为0
    vg = zeros(size(t));
    // 找到满足区间条件的元素索引
    idx = (t >= Tcg) & (t <= Tfg);
    // 对符合条件的元素计算函数值
    vg(idx) = Ag^2 * (1 - cos(2*%pi*(t(idx)-Tcg)/(Tfg-Tcg)));
endfunction

代码解释 #

1. 先用zeros(size(t))初始化结果向量vg，确保长度和输入t一致，默认值为0，对应t<Tcg或t>Tfg的情况。
2. 用逻辑运算(t >= Tcg) & (t <= Tfg)生成布尔索引向量idx，标记符合区间条件的位置。
3. 对vg中idx为真的位置，代入公式计算对应值，Scilab会自动对向量的对应元素逐运算。

调用示例 #

直接传入向量t和参数即可：

Tc=0;     
Tm=1;   
Tf=1800;  
t=(Tc:1:Tf)';
Ag=2;    
Tcg=200;    
Tfg=400;

vg_result = simula_vg(t, Tcg, Tfg, Ag);

通用写法（适用于其他分段函数） #

对于类似的分段向量函数，都可以遵循以下模板：

function output=func_name(input_vec, param1, param2, ...)
    // 初始化结果向量
    output = zeros(size(input_vec));
    
    // 定义各分段的逻辑索引
    idx1 = 条件1;
    idx2 = 条件2;
    // ... 更多分段
    
    // 为各分段赋值
    output(idx1) = 表达式1;
    output(idx2) = 表达式2;
    // ... 更多分段赋值
endfunction

比如你需要定义的其他涉及Ar、Tcr、Tfr的函数，直接替换模板中的条件和表达式即可。

内容的提问来源于stack exchange，提问作者q19a