编写的plot语句报错了。

sharkN · 发表于 2020-11-25 18:25

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

我用一个调用函数来计算一个区间上的函数值并画图。此调用函数的输出是x = [x1;  x2;  x3];我想要把三个值的图像都画出来，编写的plot语句却都错了。错误信息：
错误使用 plot。
数据必须为可转换为双精度值的数值、日期时间、持续时间或数组。

出错 Untitled3 (line 3)
plot(y,c(:,1));

以下是程序，使用的64位R2018。
y=9.9:0.01:10.1;
c=Solve3Polynomial(1,2.*k-20,power(k-10,2)+0.00005,-0.000001);
plot(y,c(:,1));
plot(y,c(:,2));
plot(y,c(:,3));
function x = Solve3Polynomial(a,b,c,d)
% 范盛金. 一元三次方程的新求根公式与新判别法[J]. 海南师范学院学报,1989,2(2):91-98.
syms x1 x2 x3;
A = b.*b - 3.*a.*c; if abs(A) < 1e-14; A = 0;  end
B = b.*c - 9.*a.*d; if abs(B) < 1e-14; B = 0;  end
C = c.*c - 3.*b.*d; if abs(C) < 1e-14; C = 0;  end
DET = B.*B - 4.*A.*C; if abs(DET) < 1e-14;  DET = 0;  end
if (A==0)&(B==0)
x1 = -c./b;    x2 = x1 ; x3 = x1;
end
if DET > 0
Y1 = A.*b + 1.5.*a.*(-B + sqrt(DET));
Y2 = A.*b + 1.5.*a.*(-B - sqrt(DET));
y1 = nthroot(Y1,3);  y2 = nthroot(Y2,3);
x1 = (-b-y1-y2)./(3.*a);
vec1 = (-b + 0.5.*(y1 + y2))./(3.*a);
vec2 = 0.5.*sqrt(3).*(y1 - y2)./(3.*a);
x2 = 0;
x3 = 0;
clear Y1 Y2 y1 y2 vec1 vec2;
end
if DET == 0 & (A ~= 0) & (B ~= 0)
K = (b.*c-9.*a.*d)./(b.*b - 3.*a.*c); K = round(K,14);
x1 = -b./a + K; x2 = -0.5.*K; x3 = x2;
end
if DET < 0
sqA = sqrt(A);
T = (A.*b - 1.5.*a.*B)./(A.*sqA);
theta = acos(T);
csth  = cos(theta./3);
sn3th = sqrt(3).*sin(theta./3);
x1 = (-b - 2.*sqA.*csth)./(3.*a);
x2 = (-b + sqA.*(csth + sn3th))./(3.*a);
x3 = (-b + sqA.*(csth - sn3th))./(3.*a);
clear sqA T theta csth sn3th;
end
x = [x1;  x2;  x3];
end

麻烦大神帮忙指点一下，非常感谢！

yxlk · 发表于 2020-11-25 18:40

提示说，数据必须为可转换为双精度值的数值、日期时间、持续时间或数组。所以你的数据是什么类型？

帐号		自动登录	找回密码
密码			注册

编写的plot语句报错了。

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1