TVAR处理信号问题，求高手帮忙看看程序

网上下载的TVAR工具包，基于模型分解理论，经过整理后程序如下
signal=load('10.54.txt');
x=signal-mean(signal);
N=length(x);
freqs=200;
t=(N-1)/freqs;
T=0:1/freqs:t;
subplot(2,2,1)
plot(T,signal)
xlabel('时间（s）')
ylabel('幅值')
pvals=[1,30];
ndel=[2,2];
mindel=[1,1];
[popt,delopt,likp,mt,Ct,st,nt] = tvar_lik(x,pvals,ndel,mindel);
delopt;
popt;
likp;
%choose and fit AR(8)
p=popt;
m0=[1,zeros(p-1,1)']';
c0=eye(p)*10;
s0=st;
n0=nt;
del=delopt;
[m1,c,s1,n,e]=tvar(x,p,del,m0,c0,s0,n0);
m=m1;
s=s1;
%compute decomposition,wavelengths etc,then plot
[waves,mods,decomp,nr,nc]=tvar_decomp(x,m);
subplot(2,2,2)
decomp_plot(x,decomp,p)
%sample posterior N times at selected points and plot trajectories of first
%k(at %most)frequencies
%k=4;
%times=T;
%[freq]=tvar_sim(m,c,n,times,k,N);
%sim_plot(times,freq,waves,4,0);
%plot frequencies
%sim_plot(times,freq,waves,4,1);
%plot wavelengths,plot the tvar spectrum over time and freq ranges specified
[times,omega,spec]=tvar_spec(m,s,T,0:200);
contourf(T,0:freqs,abs(spec)); view(75,-30);
现存在问题是已经能够对信号进行分解，但是不能得到其时频谱图，望高手能够指教一下！