用matlab求频率响应

对一个产生的离散序列作傅立叶变换，分析其频率响应特性（给出幅频与相频特性曲线）-A discrete sequence generated for the Fourier transform to analyze the frequency response characteristics (frequency and amplitude given phase frequency response curve)

close all
N=input（‘请输入FFT的变换区间长度N:\n‘）;   %  N=16 32 64
n=0:N-1;
x1=2/pi*sin（n*2*pi/N）+2/pi*（sin（3*n*2*pi/N）/3）+2/pi*（sin（5*n*2*pi/N）/5）+2/pi*（sin（7*n*2*pi/N）/7）+2/pi*（sin（9*n*2*pi/N）/9）+2/pi*（sin（11*n*2*pi/N）/11）+2/pi*（sin（13*n*2*pi/N）/13）+2/pi*（sin（15*n*2*pi/N）/15）+2/pi*（sin（17*n*2*pi/N）/17）+2/pi*（sin（19*n*2*pi/N）/19）+2/pi*（sin（21*n*2*pi/N）/21）;
subplot（221）;%输出x1图像
plot（nx1）;