采样频率为1000，采样数据点为128，输入信号为S=sin(2*pi*50t)+0.4sin(3*2*pi*50*t)+0.3sin(5*2*pi*50*t)+0.05sin(7*2*pi*50*t)的加窗插值fftmatlab算法

这里提供一种可能的算法：

1. 生成时间序列t，从0到（128-1）/1000秒，步长为1/1000秒。

2. 计算输入信号S在时间序列t上的取值，得到长度为128的信号向量s。

3. 对s进行汉宁窗处理，得到窗口向量w。

4. 将s和w相乘，得到加窗后的信号向量s_w。

5. 对s_w进行插值，插值倍数为M，即插值后的信号点数为M*128。插值方法可以选择线性插值或样条插值。

6. 对插值后的信号向量进行FFT计算，得到频域向量f。

7. 计算f的幅度谱abs_f和相位谱angle_f。

8. 绘制频谱图，横坐标为频率，纵坐标为幅度谱，可以使用matlab中的plot函数或stem函数。

9. 输出幅度谱和相位谱，可以使用matlab中的disp函数或fprintf函数。

注意事项：

1. 在插值之前，可以对信号进行零填充，以增加FFT分辨率。

2. 在绘制频谱图时，需要对横坐标进行换算，从频率序列转换为对应的频率值。