各位前辈，程序运行到[E,tt1]=toimage(A,fa,tt,length(tt));报错如下：

??? Error: File: toimage.m Line: 59 Column: 1

This statement is not inside any function.

(It follows the END that terminates the definition of the function "toimage".)

Error in ==> Untitled at 30

[E,tt1]=toimage(A,fa,tt,length(tt));

我感觉到时定义toimage函数时出了问题，但是不知道如何改正，请前辈帮我指点一下，谢谢！程序中的错误我用红色表示了。

全部程序如下：

clear;clc;clf;

N=256;

t=linspace(1,2,N);deta=t(2)-t(1);fs=1/deta;

x=15*sin(2*pi*t*5)+51*sin(2*pi*t*15)+25*sin(2*pi*t*30);

z=x;

c=emd(z);

[imf,ort,nbits]=emd(x,t,[0.05,0.5,0.05]);

L=size(imf,1);

X=0;

for i=1:L;

X=imf(1,:)+X;

end

Y=x-X;

[m,n]=size(c);

for i=1:m;

a=corrcoef(c(i,:),z);

xg(i)=a(1,2);

end

xg;

for i=1:m-1;

mse(i)=mean(c(i,:).^2.2)-mean(c(i,:),2).^2;

end

mmse=sum(mse);

for i=1:m-1;

mse(i)=mean(c(i,:).^2.2)-mean(c(i,:),2).^2;

mseb(i)=mse(i)/mmse*100;

end

hold off

[A,fa,tt]=hhspectrum(c);

[E,tt1]=toimage(A,fa,tt,length(tt));disp_hhs(E,tt1)

figure(4)

for i=1:size(c,1)

faa=fa(i,:);

[FA,TT1]=meshgrid(faa,tt1)

surf(FA,TT1,E)

title('HHT时频谱三维显示')

xlabel('t/s');

ylabel('f/Hz');

zlabel('a/V');

hold on

end

E=flipud(E);

for k=1:size(E,1)

bjp(k)=sum(E(k,:))*1/fs;

end

f=(1:N-2)/N*(fs/2);

[A,f,tt] = hhspectrum(x,t,l,aff)

error(nargchk(1,4,nargin));

if nargin < 2

t=1:size(x,2);

end

if nargin < 3

l=1;

end

if nargin < 4

aff = 0;

end

if min(size(x)) == 1

if size(x,2) == 1

x = x';

if nargin < 2

t = 1:size(x,2);

end

end

Nmodes = 1;

else

Nmodes = size(x,1);

end

lt=length(t);

tt=t((l+1):(lt-l));

for i=1:Nmodes

an(i,:)=hilbert(x(i,:)')';

f(i,:)=instfreq(an(i,:)',tt,l)';

A=abs(an(:,l+1:end-l));

if aff

disprog(i,Nmodes,max(Nmodes,100))

end

end

disp_hhs(im,t,inf)

figure

colormap(bone)

colormap(1-colormap);

if nargin==1

inf=-20;

t = 1:size(im,2);

end

if nargin == 2

if length(t) == 1

inf = t;

t = 1:size(im,2);

else

inf = -20;

end

end

if inf >= 0

error('inf doit etre < 0')

end

M=max(max(im));

im = log10(im/M+1e-300);

inf=inf/10;

imagesc(t,fliplr((1:size(im,1))/(2*size(im,1))),im,[inf,0]);

set(gca,'YDir','normal')

xlabel(['time'])

ylabel(['normalized frequency'])

title('Hilbert-Huang spectrum')

toimage.m程序如下：

function [im,tt] = toimage(A,f,t,splx,sply)

% [im,tt] = TOIMAGE(A,f,t,splx,sply) transforms a spectrum made

% of 1D functions (e.g., output of "spectreh") in an 2D image

%

% inputs :   - A    : amplitudes of modes (1 mode per row of A)

%            - f    : instantaneous frequencies

%            - t    : time instants

%            - splx : number of columns of the output im (time resolution).

%                     If different from length(t), works only for uniform

%                     sampling.

%            - sply : number of rows of the output im (frequency resolution).

% outputs :  - im   : 2D image of the spectrum

%            - tt   : time instants in the image

%

% utilisation : [im,tt] = toimage(A,f);[im,tt] = toimage(A,f,t);[im,tt] = toimage(A,f,sply);

%              [im,tt] = toimage(A,f,splx,sply);[im,tt] = toimage(A,f,t,splx,sply);

DEFSPL = 400;

if nargin < 3

t = 1:size(A,2);

sply = DEFSPL;

splx = length(t);

else

if length(t) == 1

tp = t;

t = 1:size(A,2);

if nargin < 4

sply = tp;

splx = length(t);

else

if nargin > 4

error('too many arguments')

end

sply = splx;

splx = tp;

end

else

lt = length(t);

if nargin < 5

sply = splx;

splx = lt;

end

if nargin < 4

sply = DEFSPL;

splx = lt;

end

if nargin > 5

error('too many arguments')

end

end

end

end

lt=length(t);

im=[];

im(splx,sply) = 0;

for i=1:size(f,1)

for j = 1:lt

ff=floor(f(i,j)*2*(sply-1))+1;

if ff <= sply % in case f(i,j) > 0.5

im(floor(j*(splx-1)/lt)+1,ff)=im(floor(j*(splx-1)/lt)+1,ff)+A(i,j);

end

end

end

for i = 1:splx

tt(i) = mean(t(floor((i-1)*lt/(splx))+1:floor(i*lt/(splx))));

end

im=fliplr(im)';

if issparse(A) || ~isreal(A) || length(size(A)) > 2

error('A argument must be a real matrix')

end

if issparse(f) || ~isreal(f) || length(size(f)) > 2

error('f argument must be a real matrix')

end

if any(size(f)~=size(A))

error('A and f matrices must have the same size')

end

if issparse(t) || ~isreal(t) || ~isvector(t) || length(t)~=size(A,2)

error('t argument must be a vector and its length must be the number of columns in A and f inputs')

end

if ~isscalar(splx) || ~isreal(splx) || splx ~= floor(splx) || splx <= 0

error('splx argument must be a positive integer')

end

if ~isscalar(sply) || ~isreal(sply) || sply ~= floor(sply) || sply <= 0

error('splx argument must be a positive integer')

end

if any(diff(diff(t))) && splx ~= length(t)

warning('toimage:nonuniformtimeinsants','When splx differs from length(t), the function only works for equally spaced time instants. You may consider reformating your data (using e.g. interpolation) before using toimage.')

end

f = min(f,0.5);

f = max(f,0);

indf = round(2*f*(sply-1)+1);%

indt = repmat(round(linspace(1,length(t),splx)),size(A,1),1);%创建1*size(A,1)大小的矩阵，矩阵中的数据为round(linspace(1,length(t),splx))

im = accumarray([indf(,indt(],A(,[sply,splx]);

indt = indt(1,;

tt = t(indt);

ff = (0:sply-1)*0.5/sply+1/(4*sply);

end