下面来看看所谓的“边际谱”。“边际谱”就是在Hilbert谱的时频平面上，各频率点振幅在时间总体上的累积，也就是频率相同、总时长上所有振幅的叠加。

（1）脉搏MB291712的边际谱：

    先看看MB291712由包括残余项在内的Hilbert谱计算得到的边际谱：

[E_cimf_MB291712,t_imf_MB291712,Cenf_imf_MB291712]=toimage(A_imf_MB291712,f_imf_MB291712);

%E_cimf_MB291712表示由各IMF分量及残余项计算所得到的Hilbert谱

bjp_cimf_MB291712=zeros(1,size(E_cimf_MB291712,1));

for k=1:size(E_cimf_MB291712,1);

bjp_cimf_MB291712(k)=sum(E_cimf_MB291712(k,:))/size(E_cimf_MB291712,2);

end

plot(Cenf_imf_MB291712,20*log10(bjp_cimf_MB291712))

    运行，得:

    图16-1 bjp_cimf_MB291712

    再看看MB291712由只有IMF分量但不包括残余项在内的Hilbert谱计算得到的边际谱：

E_imf_MB291712,t_imf_MB291712,Cenf_imf_MB291712]=toimage(A_imf_MB291712(1:13,:),f_imf_MB291712(1:13,:));

%E_imf_MB291712表示由各IMF分量但不包括残余项计算所得到的Hilbert谱

bjp_imf_MB291712=zeros(1,size(E_imf_MB291712,1));

for k=1:size(E_imf_MB291712,1);

bjp_imf_MB291712(k)=sum(E_imf_MB291712(k,:))/size(E_imf_MB291712,2);

end

plot(Cenf_imf_MB291712,20*log10(bjp_imf_MB291712))

    运行，得： 

    图16-2 bjp_imf_MB291712

    比起图16-1来，仅仅在0频率（或包括紧挨0频率的低频处），谱线降低了很多，其余地方完全相同。

（2）体温TW486012的边际谱：

[E_cimf_TW486012,t_imf_TW486012,Cenf_imf_TW486012]=toimage(A_imf_TW486012,f_imf_TW486012);

%E_cimf_TW486012表示由TW486012各IMF分量及残余项计算所得到的Hilbert谱

bjp_cimf_TW486012=zeros(1,size(E_cimf_TW486012,1));

for k=1:size(E_cimf_TW486012,1);

bjp_cimf_TW486012(k)=sum(E_cimf_TW486012(k,:))/size(E_cimf_TW486012,2);

end

plot(Cenf_imf_TW486012,20*log10(bjp_cimf_TW486012))

    运行，得： 

    图16-3 bjp_cimf_TW486012

     （3）收缩压GY291712的边际谱：

[E_cimf_GY291712,t_imf_GY291712,Cenf_imf_GY291712]=toimage(A_imf_GY291712,f_imf_GY291712);

%E_cimf_GY291712表示由GY291712各IMF分量及残余项计算所得到的Hilbert谱

bjp_cimf_GY291712=zeros(1,size(E_cimf_GY291712,1));

for k=1:size(E_cimf_GY291712,1);

bjp_cimf_GY291712(k)=sum(E_cimf_GY291712(k,:))/size(E_cimf_GY291712,2);

end

plot(Cenf_imf_GY291712,20*log10(bjp_cimf_GY291712)) 

    运行，得：

    图16-4 bjp_cimf_GY291712

（4）舒张压DY291712的边际谱： 

[E_cimf_DY291712,t_imf_DY291712,Cenf_imf_DY291712]=toimage(A_imf_DY291712,f_imf_DY291712);

%E_cimf_DY291712表示由DY291712各IMF分量及残余项计算所得到的Hilbert谱

bjp_cimf_DY291712=zeros(1,size(E_cimf_DY291712,1));

for k=1:size(E_cimf_DY291712,1);

bjp_cimf_DY291712(k)=sum(E_cimf_DY291712(k,:))/size(E_cimf_DY291712,2);

end

plot(Cenf_imf_DY291712,20*log10(bjp_cimf_DY291712))

    运行，得： 

    图16-5 bjp_cimf_DY291712

（5）均压JY291712的边际谱：

[E_cimf_JY291712,t_imf_JY291712,Cenf_imf_JY291712]=toimage(A_imf_JY291712,f_imf_JY291712);

%E_cimf_JY291712表示由JY291712各IMF分量及残余项计算所得到的Hilbert谱

bjp_cimf_JY291712=zeros(1,size(E_cimf_JY291712,1));

for k=1:size(E_cimf_JY291712,1);

bjp_cimf_JY291712(k)=sum(E_cimf_JY291712(k,:))/size(E_cimf_JY291712,2);

end

plot(Cenf_imf_JY291712,20*log10(bjp_cimf_JY291712))

    运行，得：

    图16-6 bjp_cimf_JY291712

（6）差压CY291712的边际谱： 

[E_cimf_CY291712,t_imf_CY291712,Cenf_imf_CY291712]=toimage(A_imf_CY291712,f_imf_CY291712);

%E_cimf_CY291712表示由CY291712各IMF分量及残余项计算所得到的Hilbert谱

bjp_cimf_CY291712=zeros(1,size(E_cimf_CY291712,1));

for k=1:size(E_cimf_CY291712,1);

bjp_cimf_CY291712(k)=sum(E_cimf_CY291712(k,:))/size(E_cimf_CY291712,2);

end

plot(Cenf_imf_CY291712,20*log10(bjp_cimf_CY291712))

    运行，得：

    图16-7 bjp_cimf_CY291712

    每条曲线形状都大同小异，我还是耐着性子把它们做完了。真无聊……不过我也真伟大啊……把曲线放大还是各不相同的。它的性质其实就相当于各信号的功率谱，只不过计算方法特别而已。姑存勿论……

（本文首发于：http://blog.sina.com.cn/s/blog_6ad0d3de0100zk0f.html

首发时间：2011-09-17 11:22:04）