|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clear
rand('seed', 0);
T = 0.05; % 仿真时间
R = 500; % 码速率
fd = 1000; % 载波频差
fc = 2000; % 载波频率
fs = 20000; % 采样率
samp = fs/R; % 每个码元上的采样点数
N = T*fs;
n = 1:N;
x = randint(1, R*T, 2);
y = fskmod(x, 2, fd, samp, fs);
y = y .* exp(i*2*pi*fc/fs*n);
y = real(y);
% z = awgn(y, 20, 'measured');
z = y;
imf=emd(z);
emd_visu(z,1:length(z),imf)
[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))
代码如上;
图如下,请大家对输出的图分析下啊 |
|