matlab周期矩形脉冲幅度谱及频谱,Matlab矩形脉冲频谱分析程序出错

作一个矩形脉冲的频谱分析 总出错 :L 。哪位大侠帮看看，谢谢了。

??? Error using ==> sym.double at 29

DOUBLE cannot convert the input expression into a double array.

If the input expression contains a symbolic variable, use the VPA function instead.

Error in ==> fbsbpp at 16

An(1)=double(vpa(A0,Nn));

function [A_sym,B_sym]= CTFStpshsym

% 采用符号计算求[0,T]内时间函数的三角级数展开系数，并绘制其双边频谱。

% 函数的输出为数值量

% Nn 输出数据的准确位数

% A_sym 第1 元素是直流项，其后元素依次是1,2,3...次谐波cos 项展开系数

% B_sym 第2,3,4,...元素依次是1,2,3...次谐波sin 项展开系数

% T T=m*tao， 信号周期

% Nf 谐波的阶数

% Nn 输出数据的准确位数

% m (m=T/tao)周期与脉冲宽度之比,如m=4,8,16,100 等

% tao 脉宽:tao=T/m

syms t n y

if nargin<3;Nf=input('所需展开的最高谐波次数:Nf=');end

T=input('pleas Input 信号的周期T=');

if nargin<5;Nn=32;e