???使用 ==> 时出错绘图 - 无法从 sym 转换为 double

Question

我的代码有一些问题。如下：

lambdaz = 1.2;
n = 24;
mu = 0.00055e9;

lambda = sym('lambda','clear');
W = (((2.*mu)./n.^2)).*((lambda.^n)+(lambdaz.^n)+((lambda.^-n).*(lambdaz.^-n))-3);
dW_dlambda = diff(W, lambda);

W2=(((2.*mu)./n.^2).*(lambda.^n))+(((2.*mu)./n.^2).*(lambdaz.^n))+(((2.*mu)./n.^2).*((lambda.^-n).*(lambdaz.^-n)))-(3.*((2.*mu)./n.^2))
dW2_dlambda=diff(W2,lambda)

x=((((lambda.^2).*(lambdaz))-1).^-1).*(dW_dlambda);

x2=((((lambda.^2).*(lambdaz))-1).^-1).*(dW2_dlambda)

P2 = int(x2,lambda)

P=int(x,lambda);
P=(0:1000:26700)

plot(lambda,P)

当我尝试绘制 lambda 与 PI 的关系图时，出现“不可能从 sym 转换为 double”错误消息。我对 Matlab 不是特别擅长，所以我们将不胜感激地收到任何帮助！

Answer 1

plot 函数仅适用于数字输入。 lambda 和 P 都是符号表达式（至少在积分后通过将其设置为等于向量来覆盖 P 之前），无法直接转换到 到 浮点数。如果尝试类似 double(sym('exp(x)')) 之类的操作，您会收到相同的错误。你有两个选择。第一个是 Symbolic 中的 ezplot 函数工具箱：

...
P = int(x,lambda);
ezplot(P,[-5 5]); % Plot's P from lambda = -5 to lambda = 5

或者您可以使用 subs 函数：

...
P = int(x,lambda);
lambda = -5:0.01:5;
plot(lambda,real(subs(P,'lambda',lambda)))
axis([lambda(1) lambda(end) -1e15 1e15])

我使用 real 来抑制针对 lambda 负值的警告。