Matlab布莱克曼窗相频和幅频特性

我想实现图二那样的结果，该怎么修改图一的代码呢？还有plot里的参数改怎么填写呢？