131. 132. 133. 134. 135. 136. 137. 138. 139. 140. 141. 142. 143. 144. 145. 146. 147. 148. 149. 150. 151. 152. 153. 154. 155. 156. 157. 158. 159. 160. 161. 162. 163. 164. 165. 166. 167. 168. 169. 170. 171. 172. 173. 174. 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))
175. 176. 177. 178. 179. 180. 181. 182. 183. 184. 185. 186. 187. 188. 189. 190. 191. 192. 193. 194. 195. 196. 197. 198. 199. 200. 201. 202. 203. 204. 205. 206. 207. 208. 209. 210. 211. 212. 213. 214. 215. 216. 217. 218. end end
function disp_hhs(im,t,inf)
% DISP_HHS(im,t,inf)
% displays in a new figure the spectrum contained in matrix \% (amplitudes in log). %
% inputs : - im : image matrix (e.g., output of \% - t (optional) : time instants (e.g., output of \% - inf (optional) : -dynamic range in dB (wrt max) % default : inf = -20 %
% utilisation : disp_hhs(im) ; disp_hhs(im,t) ; disp_hhs(im,inf) % 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));
219. im = log10(im/M+1e-300); 220. 221. inf=inf/10; 222. 223. 224. imagesc(t,fliplr((1:size(im,1))/(2*size(im,1))),im,[inf,0]); 225. set(gca,'YDir','normal') 226. xlabel(['time']) 227. ylabel(['normalized frequency']) 228. title('Hilbert-Huang spectrum') 229. function [f,z]=fftfenxi(t,y) 230. L=length(t);N=2^nextpow2(L); 231. ?t默认计算的信号是从0开始的 232. t=linspace(t(1),t(L),N);deta=t(2)-t(1); 233. m=0:N-1; 234. f=1./(N*deta)*m; 235. %下面计算的Y就是x(t)的傅里叶变换数值 236. %Y=exp(i*4*pi*f).*fft(y)%将计算出来的频谱乘以exp(i*4*pi*f)得到
频移后[-2,2]之间的频谱值 237. Y=fft(y); 238. z=sqrt(Y.*conj(Y)); 239. 240. 241. 242.
复制代码
4.总结。
（1）边际谱与傅里叶谱的比较：
意义不同：边际谱从统计意义上表征了整组数据每个频率点的累积幅值
分布，而傅里叶频谱的某一点频率上的幅值表示在整个信号里有一个含有此频率的三角函数组分。
作用不同：边际谱可以处理非平稳信号，如果信号中存在某一频率的能量出现，就表示一定有该频率的振动波出现，也就是说，边际谱能比较准确地反映信号的实际频率成分。而傅里叶变换只能处理平稳信号。
（2） HHT与Hilbert变换的比较：
Hilbert变换只是单纯地求信号的瞬时振幅，频率和相位，有可能出现
没有意义的负频率；HHT变换先将信号进行EMD分解，得到的是各个不同尺度的分量，对每一个分量进行Hilbert变换后得到的是有实际意义的瞬时频率。
PS:运行上面的程序需要装时频工具箱，我仅将用到的emd分解的程序贴到下面。