在 MATLAB 中进行数值计算时，sqrt() 和 exp() 是两个最基础、最常用的函数。它们分别用于计算平方根和自然指数，几乎出现在所有涉及数学建模、工程仿真或数据分析的脚本中。掌握它们的正确用法，能避免常见错误并提升代码效率。
1. 使用 sqrt() 计算平方根
调用 sqrt(x) 函数即可计算任意非负实数 x 的平方根。若输入为负数，MATLAB 不会报错，而是返回复数结果。
- 打开 MATLAB 命令窗口或新建一个
.m脚本文件。 - 输入
y = sqrt(16)，按回车执行。结果为y = 4。 - 尝试 负数输入：
z = sqrt(-9)，返回z = 0.0000 + 3.0000i，表示虚部为 3 的复数。 - 批量处理 数组：
A = [4, 9, 16]; B = sqrt(A)，得到B = [2, 3, 4]。
注意：
sqrt()对矩阵元素逐个运算（即“逐元素运算”），不是矩阵平方根（后者需用sqrtm()）。
2. 使用 exp() 计算自然指数
调用 exp(x) 函数可计算 $ e^x $，其中 $ e \approx 2.71828 $ 是自然对数的底数。
- 计算 单值：
result = exp(1)，返回2.7183（即 $ e $ 的近似值）。 - 验证 欧拉公式：
exp(1i * pi) + 1，结果接近0（因浮点精度误差，实际为极小复数）。 - 处理向量：
t = 0:0.1:1; y = exp(-t)，生成从 1 衰减到约 0.3679 的指数衰减序列。 - 避免溢出：当
x > 709时，exp(x)返回Inf（无穷大），因为超出双精度浮点数范围。
3. 常见组合用法与性能提示
许多科学计算场景需要同时使用这两个函数，例如高斯函数、概率密度或热传导模型。
示例：计算标准正态分布的概率密度
标准正态分布公式为：
$$
f(x) = \frac{1}{\sqrt{2\pi}} e^{-x^2/2}
$$
实现步骤如下：
- 定义 自变量向量：
x = -3:0.01:3; - 计算 分母部分：
denom = sqrt(2 * pi); - 计算 指数部分：
exponent = exp(-x.^2 / 2); - 合成 密度函数：
pdf = exponent / denom;
注意：x.^2 中的点号 . 表示逐元素平方，这是 MATLAB 向量化运算的关键。
4. 错误排查与最佳实践
| 常见问题 | 原因 | 解决方法 |
|---|---|---|
sqrt(-1) 返回复数而非报错 |
MATLAB 默认支持复数 | 若仅需实数结果，先检查输入是否非负：if x >= 0, y = sqrt(x); end |
exp(1000) 返回 Inf |
超出浮点数上限 | 改用对数域计算，如比较 log(a) > log(b) 而非直接计算 exp(a) > exp(b) |
| 数组运算结果维度不符 | 忘记使用点运算符（.） |
确保对数组使用 .^, .*, ./ 等逐元素操作符 |
5. 性能对比：向量化 vs 循环
对于大规模数据，避免用 for 循环逐个调用 sqrt() 或 exp()。
% 低效方式（不推荐）
n = 1e6;
x = rand(n,1);
y = zeros(n,1);
for i = 1:n
y(i) = sqrt(x(i));
end
% 高效方式（推荐）
y = sqrt(x);MATLAB 内置函数高度优化，向量化代码通常快数十倍，且代码更简洁。
6. 扩展：与其他函数的协同
- 与
log()配合：exp(log(x))在x > 0时恒等于x，可用于数值稳定性处理。 - 与
.^结合：x.^0.5等价于sqrt(x)，但前者支持负数的分数幂（可能产生复数），而sqrt()更语义清晰。 - 在绘图中使用：绘制指数增长曲线：
plot(t, exp(0.5*t))；绘制平方根函数：fplot(@(x) sqrt(x), [0, 10])。
记住：sqrt() 和 exp() 是 MATLAB 数值计算的基石。只要确保输入合理、善用向量化、警惕复数与溢出，就能高效完成绝大多数科学计算任务。
暂无评论，快来抢沙发吧！