在 MATLAB 中调用 mldivide 时出现矩阵维数错误

Question

运行代码时出现此错误：

错误使用 ==> mldivide 矩阵维度必须一致。

这是我的代码：

%make the plots of phase and group velocity vs discreteness of the grid
c=1;

a=input('Please enter the ratio cdt/dx : ')

figure(1)
R=2:40;
plot(R,phase_vel(R,a)/c)
xlabel('R=l/dx')
ylabel('u_phase/c')

%figure(2)
%plot(group_vel(R,a),R,0,40)
%xlabel('R=l/dx')
%ylabel('u_group/c')

这是我的功能：

function phase_velocity = phase_vel(R,a)
    %numerical phase velocity of the discrete wave 
    c=1;
    phase_velocity=(2*pi*c)/(R*knum(R,a));
end

function group_velocity =group_vel(R,a )
    %numerical group velocity of the discrete wave
    c=1;
    group_velocity=(a*sin(knum(R,a)))/(sin(2*pi*a/R))
end

function knumber = knum(R,a)
    %This is the k wave number
    knumber=acos((1/a)^2*(cos(2*pi*a/R)-1)+1);
end

我该如何解决这个错误？

编辑：我用过 .每个方程中的运算符，我改变了 R=4:40 的限制

Answer 1

如果您的目标是将公式应用于向量 R 中的每个单独的值，那么您应该使用 element-wise arithmetic operators .*、./ 和 .^ 而不是矩阵运算符*、/ 和^。

您的错误可能发生在第一次调用您的函数knum 时，特别是当您尝试计算2*pi*a/R 时。由于2*pi*a 是单个标量值，因此在尝试使用行向量R 执行matrix right division / 时会出错。 真的奇怪的是错误信息：

??? Error using ==> mldivide
Matrix dimensions must agree.

这意味着您正在使用矩阵 left 除法运算符\，而您显然不是。我在 MATLAB R2010b 中对此进行了测试，我的消息中出现了相同的错误函数名称。我认为这可能只是错误消息中的一个错字，我已经给 MATLAB 人员写了一张便条，让他们查看并清除它。

Answer 2

我没有符号数学工具箱，但您的问题似乎是您正在使用plot，这是一个可以处理数字数组并将符号计算结果提供给它的函数。查看 Matlab 帮助，其中主题 Creating Plots of Symbolic Functions 建议使用 ezplot()。或者，您需要评估某些输入值的符号表达式以创建 plot 可以处理的数字数组 - 但您不能使用 double() ，因为它不知道要插入哪些数字到您的变量中.