在符号表达式上使用 fplot 的问题

Question

这是一个 fplot 不绘制任何内容的示例：

a=0.336;
Ta=9.476;
Te=1.208;
Tw=1.498;
eqh=[0.661;0.619;0.568];
ex=[-1.24;-1.346;-1.441];
en=-ex;
ey=[0.376;0.705;0.968];
eqx=[-0.309;-0.357;-0.392];
eh=[1.594;1.583;1.545];
eyqh=[0.642;0.78;0.897];
a0=a*Ta*Te^2;
syms bt Td Ki Kp;
a1=sym([]);a2=sym([]);a3=sym([]);a4=sym([]);exqh=sym([]);

for i=1:3
        Kp=1/bt;
        exqh(i)=en(i)*eqh(i)+eqx(i)*eh(i);
        Ki=1/(bt*Td);
        a1(i)=Ta*Tw*eqh(i)+a*Te^2*en(i)+a*Te^2*ey(i)*Kp;
        a2(i)=a*Te^2*ey(i)*Ki+Ta+Tw*exqh(i)-eyqh(i)*Kp*Tw;
        a3(i)=en(i)+ey(i)*Kp-eyqh(i)*Ki*Tw;
        a4(i)=ey(i)*Ki;
        assume(bt~=0)
        f=@(bt) a1(i)
        fplot(f,[0.01 1],'b')
        hold on
end

这是另一个例子，fplot 工作：

syms y x;
y=@(x) 2/x+6;
z=y;
assume(x~=0)
fplot(z,[-1 1],'b')

我无法理解这两种情况之间的区别。在我看来，它们是相同的。

谁能解释一下为什么上面的例子不起作用，但下面的例子起作用？

Answer 1

fplot 需要一个返回数值的表达式。 f 返回一个 sym（符号表达式），所以它不起作用，而 z 返回一个数字 - 所以 它 可以。这就是所有的区别。

请注意，在工作示例中，您覆盖了第 2^nd 行的 y sym，这意味着您向 fplot 提供了“正确的”function handle，“by错误”。事实上，您需要更少的代码来让该示例工作：

fplot( @(x)2./x+6, [-1 1], 'b');

让顶级代码工作的最简单方法是使用matlabFunction。换句话说，您可以通过更改此行来修复您的代码：

f=@(bt) a1(i)

到这里：

f = matlabFunction(a1(i));