MATLAB绘制矩阵行列式异常：仅得单个y值，如何显示全x对应值？

嘿，我明白你遇到的问题了！你想要对每个x值对应的3×3矩阵计算行列式，但现在只得到一个y值，核心原因是你的代码把所有x对应的矩阵元素拼成了一个大的非方阵，而det()函数只能对单个方阵计算行列式，所以最后只返回了一个结果，完全不是你想要的逐个计算效果。

问题根源拆解 #

当x是一个1×11的向量时，你构造的矩阵[(cos(x)/x)^2,2*(cos(x)/x),2*(cos(x)/x);2*(cos(x)/x),8,8;2*(cos(x)/x),8,12]其实是一个3×33的矩阵（每行是3个1×11的向量直接拼接而成），这显然不是符合det()要求的方阵。MATLAB虽然没报错，但只返回了一个非预期的行列式值，自然无法匹配x的长度。

你真正需要做的是：对每个x中的元素，单独构造一个3×3的矩阵，再计算它的行列式。

两种实用解决方法 #

方法1：使用arrayfun（向量化写法，简洁高效）

arrayfun可以帮你对向量中的每个元素批量应用自定义计算，非常适合这种场景：

x = 0.1:0.1:1; % 注意：去掉x=0，因为cos(0)/0会得到NaN，导致行列式无意义
% 定义单个x值对应的行列式计算逻辑
calc_single_det = @(xi) det([(cos(xi)/xi)^2, 2*(cos(xi)/xi), 2*(cos(xi)/xi);
                            2*(cos(xi)/xi), 8, 8;
                            2*(cos(xi)/xi), 8, 12]);
y = arrayfun(calc_single_det, x); % 对每个x元素计算行列式
figure
plot(x,y)

方法2：使用for循环（直观易懂，适合新手理解）

如果你更喜欢清晰的流程，可以手动遍历每个x元素：

x = 0.1:0.1:1; % 同样去掉x=0避免NaN
y = zeros(size(x)); % 预先分配存储空间，提升代码运行效率
for i = 1:length(x)
    xi = x(i);
    % 构造当前xi对应的3×3矩阵
    current_mat = [(cos(xi)/xi)^2, 2*(cos(xi)/xi), 2*(cos(xi)/xi);
                   2*(cos(xi)/xi), 8, 8;
                   2*(cos(xi)/xi), 8, 12];
    y(i) = det(current_mat); % 计算单个矩阵的行列式并存储
end
figure
plot(x,y)

额外小提示 #

• 如果必须保留x=0的情况，可以单独处理这个点：比如用极限公式计算cos(x)/x在x→0时的值（极限为1），手动设置y(1)对应的行列式值。
• 两种方法最终都会得到与x长度一致的y向量，这样绘图时就能完整显示所有x对应的行列式结果了！

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