你的问题出在MATLAB的变量作用域规则上：脚本中定义的变量属于基础工作区，而和脚本同文件的局部函数默认无法直接访问这些工作区变量。下面给你三种可行的解决方案，按推荐程度排序：

1. 将常量作为参数传递（最推荐） #

这是最符合软件工程规范的做法，函数的输入输出清晰，不会依赖模糊的全局状态，也避免了后续维护时的变量冲突问题。你只需要把函数需要用到的所有外部变量（g、K、h、L、D）都作为参数传入函数：

修改后的脚本代码：

% Define constants
K=0.5; % minor loss coefficient (square edged)
g=9.8; % gravity in m/s2
h = 10; % 这里替换成你实际的h值
L = 50; % 替换成你的管道长度值
D = 0.2; % 替换成你的管道直径值
f1= 2; % value obtained from external function
v1=velocity1(f1, g, K, h, L, D);

对应的函数代码：

function v1=velocity1(f, g, K, h, L, D)
    v1= sqrt((2*g*h)/(1+(f*(L./D))+K));
end

2. 使用嵌套函数 #

如果不想传递太多参数，可以把整个脚本改成一个主函数，然后将velocity1作为嵌套函数放在主函数内部。嵌套函数可以直接访问父函数（主函数）中定义的所有变量：

function LvsQ_plot()
    % 所有常量和变量定义在主函数内部
    K=0.5; % minor loss coefficient (square edged)
    g=9.8; % gravity in m/s2
    h = 10; % 替换为实际值
    L = 50; % 替换为实际值
    D = 0.2; % 替换为实际值
    f1= 2; % value obtained from external function
    v1=velocity1(f1);

    % 嵌套函数，可直接访问主函数的变量
    function v1=velocity1(f)
        v1= sqrt((2*g*h)/(1+(f*(L./D))+K));
    end
end

3. 使用全局变量（不推荐） #

虽然这种方法能快速解决问题，但全局变量会让代码的依赖关系变得不清晰，容易出现意外的变量覆盖或修改问题，后续维护难度大。如果一定要用，需要在脚本和函数中都声明全局变量：

脚本代码：

% 声明全局变量
global K g h L D
K=0.5; % minor loss coefficient (square edged)
g=9.8; % gravity in m/s2
h = 10; % 替换为实际值
L = 50; % 替换为实际值
D = 0.2; % 替换为实际值
f1= 2; % value obtained from external function
v1=velocity1(f1);

函数代码：

function v1=velocity1(f)
    % 同样声明全局变量
    global K g h L D
    v1= sqrt((2*g*h)/(1+(f*(L./D))+K));
end

注意：你的代码里h、L、D这些变量也没有在函数中定义，不管用哪种方法，都要确保这些变量能被函数访问到，否则还会出现类似的未定义变量错误。

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