实现周期信号 s=sin(2πsin(2πt))的频谱仿真的matlab代码是什么

以下是实现周期信号 s=sin(2πsin(2πt))的频谱仿真的matlab代码：

t=linspace(-0.5,0.5,1000);
s=sin(2*pi*sin(2*pi*t));
subplot(2,1,1);
plot(t,s);
title('时域信号');
xlabel('时间t');
ylabel('幅值');
y=fft(s);
f=linspace(-500,500,1000);
subplot(2,1,2);
plot(f,abs(fftshift(y)));
title('频域信号');
xlabel('频率f');
ylabel('幅值');