画出来的图不对，不知怎么改了

这是文献中的图

我写的代码上这样的，就是画不对，不知是逻辑出错了，还是单位没有统一
clc;clear;
f=20000; % 20kHz
A=12; % 12um
n=20002pi/60; % r/min
fia=90; % 90度
R=40; % 4mm
v=50; % 5mm/s
a=1; % 1mm
w0=npi/30;
w=2pif;
t=linspace(0,2pi,200);
t0=30a/(npiR);
x=Rsin(w0*(t-t0))+Asin(wt)+v*(t-t0);
y=Asin(w.t+fia);
z=R-Rcos(w0.(t-t0));
plot3(x,y,z)