为什么显示错误使用 plot 数据必须为可转换为双精度值的数值、日期时间、持续时间或数组。应该怎么

fs=2000;%采样频率
N=fs*2; %整个图由N个样点构成
tscale=0.01;% X轴显示的时间长度，单位为秒，5个周期
dt=tscale/N; %每个样点间的时间间隔
t=0:dt:(N-1)*dt;
f=500;%生成信号频率500Hz 
p0=1;
ps=p0.*sin(2*pi*f.*t);%%入射信号声压
midu=1000;%%密度
c=1500;%%声速单位m/s
syms sita arfa;%%信号入射水平方位角
%arfa=pi/6;%%信号入射方向俯仰角
vx=1/(midu*c).*sin(2*pi*f.*t)*cos(sita)*sin(arfa);
vy=1/(midu*c).*sin(2*pi*f.*t)*sin(sita)*sin(arfa);
vz=1/(midu*c).*sin(2*pi*f.*t)*cos(arfa);
%【质点振速法】
sita=atan(mean(vy)./mean(vx));
arfa=atan(sqrt(mean(vx).^2+mean(vy).^2)./mean(vz));
figure(1);
plot(t,sita);
title('500Hz正弦信号');
xlabel('t/s');
ylabel('y');

看看你的 plot 的数据是什么类型的，类型不符