hht参数（采样频率、采样点数）设置问题
请问对这个信号怎么设置采样频率、采样点数等参数 ？如果是个简单信号就没问题，这个复杂的就不会了
fs=50000;
N=10000;
t=0:fs:(N-1)*fs;
% x2=cos(0.5*t+2*cos(0.05*t));
%
% z=x2;%%z包括一个调频信号和另外两个正弦信号
% %%%%%%%%%%%%%%%%%%%%%%%%%% %%瞬时参数
init=205515866;randn('seed',init);
Sig4=rand(1,1000);
% SampFreq=500000;
% t=0:1/SampFreq:0.014;
Sig1=(t>=0&t<=0.014).*(1+sin(2*pi*15000*t)).*cos(2*pi*60000*t+sin(2*pi*15000*t));
Sig2=(t>=0&t<=0.028).*(1+sin(2*pi*20000*t)).*cos(2*pi*150000*t+sin(2*pi*20000*t));
Sig3=(t>=0.1128&t<=0.0084).*cos(2*pi*150000*t.*(1+sin(2*pi*20000*t))).*cos(2*pi*150000*t+sin(2*pi*20000*t));
Sig4=2*randn(size(Sig1+Sig2+Sig3));
Sig=Sig1+Sig2+Sig3+Sig4;
x=hilbert(Sig');
xi=imag(x);
z=Sig';%转置成和xi一致的行列
a=sqrt(z.^2+xi.^2);%瞬时幅度求法。
p=angle(x);%求瞬时相位角
=instfreq(x);
subplot(4,1,1);plot(z);title('原数据');
subplot(4,1,2);plot(a);title('瞬时幅度');
subplot(4,1,3);plot(p);title('瞬时相位');
subplot(4,1,4);plot(fnor);title('瞬时频率');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%边际谱
imf=emd(z);
emd_visu(z,1:length(z),imf,2);
=hhspectrum(imf);
=toimage(A,fa,tt,length(tt));figure;
%=toimage(A,fa);
%E=flipud(E);
for k=1:size(E,1)
bjp(k)=sum(E(k,:))*1/fs;
end
ft=(0:N-3)/N*(fs/2);
plot(ft,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
%%%%%%%%%%%%%%%%%%%%%HHT谱
%=hhspectrum(imf);
%=toimage(A,f);
disp_hhs(E);
页:
[1]