作一个矩形脉冲的频谱分析 总出错 :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