0

我写了以下matlab函数

function [t x] = MSD(xo, z, fo, c)
T = 1/fo; t = 0:T/10:10*T; fd = fo * sqrt(1-z*z);
wo = 2*pi*fo; w1 = wo * z; wd = 2*pi * fd;
x = xo * exp(-w1*t) .* cos(wd*t);
grid on;
plot(t,x,c); 

xlabel('time [s]'); 
ylabel('displacement');
s = sprintf('unforced Mass-Spring Damper [damped freq: %.3f Hz]', fd); title(s);
end

当我运行它时，命令窗口中出现以下错误：

MSD(.1, .7,.4, .2)
??? Error using ==> plot
Data must be a single matrix Y or a list of pairs X,Y

Error in ==> MSD at 13
plot(t,x,c);
4

2 回答 2

5

Acorbe 的回答是正确的，但它没有解释为什么。

简而言之，您可以使用该plot命令绘制一个向量x与一个向量y，如下所示：

plot(x, y)

或在相同的轴上绘制多个图形（向量x1对y1、向量x2对y2等...），例如：

plot(x1, y1, x2, y2, x3, y3, ...)

后一种语法要求输入向量的数量是偶数，因为它成对地绘制它们。

我相信您希望显示为与图表c相交的线。x(t)为此，您需要创建一个新函数y(t) = c以便绘制它：

y = c * ones(size(t));  % # Or c * ones(1, length(t)) like Acorbe has shown

这只是创建了一个与 t 大小相同的向量，其中每个元素等于c。只有这样，您才能将它绘制在与 相同的轴上x(t)，如下所示：

plot(t, x, t, y)

或更短的形式：

plot(t, x, t, c * ones(size(t)))


希望这能说明问题！

于 2012-11-08T18:44:20.740 回答
2

因此你的电话应该是

 plot(t,x,t,c*ones(1,length(t)));

简化后，等效于：

 plot(t,x);
 hold on

 c_vect = c * ones(1, length(t));  //you need a vector (constant in this case) 
                                   //to be plotted against t!!
                                   //ones(1, length(t)) will give you
                                   //[1111 ... 1] (as many ones as 
                                   //the entries of t)
 plot(t,c_vect);
于 2012-11-08T18:33:35.370 回答