在matlab中绘制方程

时间:2014-01-30 15:17:47

标签: matlab plot

我正在尝试在matlab中绘制方程式几天,我无法让它看起来正确。这是为了学校,所以我知道最终结果应该是什么样子。我得到类似但不一样的东西。

这是我应该得到的情节:

enter image description here

在图片下面它说我应该使用公式5.8,但我认为用它得到这条曲线是不可能的,因为它只给出了一个结果,并且在图上显示为点,这里是公式5.8:

enter image description here

如果你读到它上面的前几行,我认为很明显有一个打印错误,我应该使用公式5.6(但我不确定)所以这里是公式5.6:

enter image description here

该等式所需的数据也高于5.6。这是我的Matlab代码:

p1=1.581;
p2=-5.534;
p3=0.5523;
om1=1.214;
om2=0.001414;
om3=2.401;
c1=-0.3132;
c2=3.297;
c3=-2.381;
t=0:0.5:5

Ca=0.2132;
Ra=2.275;
V1=(p1*Ca*(sin(om1*t+c1)-om1*Ra*Ca*cos(om1*t+c1)))/(1+(om1^2)*(Ra*Ca)^2)
V2=(p2*Ca*(sin(om2*t+c2)-om2*Ra*Ca*cos(om2*t+c2)))/(1+(om2^2)*(Ra*Ca)^2)
V3=(p3*Ca*(sin(om3*t+c3)-om3*Ra*Ca*cos(om3*t+c3)))/(1+(om3^2)*(Ra*Ca)^2)
V=V1+V2+V3
plot(t,V,'.')

这是我得到的:

enter image description here

因此,在t = 0时代替零,我得到一个负数,同样适用于t = 4:5。我的错误是错误的等式还是我的代码不好?

1 个答案:

答案 0 :(得分:0)

以下是使用您的代码的情节,除了我使用

t=0:0.05:5;

a plot

看起来非常类似于本书的图。我同意指数期限可能是差异接近零的原因。