???使用 ==> 从 sym 转换为双精度时出错是不可能的

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)

当我尝试对p绘制lambda时，我得到"从sym转换到double是不可能的"错误消息。我不是特别擅长Matlab，所以任何帮助都会非常感激!

Answer 1

plot函数只适用于数字输入。lambda和P都是符号表达式(至少在将P改写为积分后的向量之前)，不能直接转换为浮点数。如果尝试double(sym('exp(x)'))，你会得到同样的错误。你有两个选择。第一个是符号工具箱中的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为负值的警告。