【无标题】

最新推荐文章于 2024-07-08 08:27:00 发布

weixin_38899721

最新推荐文章于 2024-07-08 08:27:00 发布

阅读量48

点赞数 1

文章标签：信号处理

本文链接：https://blog.csdn.net/weixin_38899721/article/details/130297249

版权

KK接收机
用示波器采到的信号包含SSBI
利用kramers-kronig关系去掉SSBI

KK接收机算法流程


idx = 0;
x_sync = sync_two_signals(wave_resample,waveform_fc1,idx); % % % sync_two_signals.m

[fscale_sync,fre_sync] = pfft(fs1,wave_resample);

[fscale_sync,fre_sync] = pfft(fs1,x_sync);
fig_sync = figure()
plot(fscale_sync,fre_sync)
axis([-1e9 5e9 -90 -30])
title('接收同步截取信号图')

%% KK receiver
x_filter = fft(x_sync);
figure();
plot(abs(x_filter))
x_filter(end/10+1:end*9/10) = x_filter(end/10+1:end*9/10)/1;
figure();
plot(10*log10(abs(x_filter)))
x_filter = ifft(x_filter);

x_filter = x_filter-0*mean(x_filter)+0.2*1;


rec_up = resample(double(x_filter),10e9,10e9);

signal_sqrt = sqrt(rec_up);
signal_log = log(signal_sqrt);
signal_log(isinf(signal_log))=0;
signal_temp = exp(imag(hilbert(signal_log))*1j).*signal_sqrt;
signal_temp_temp = signal_temp+conj(signal_temp);
[fscale_kk,fre_kk] = pfft(fs1*1,signal_temp_temp);
fig_spectrum_kk = figure();
plot(fscale_kk,fre_kk)
axis([-1e9 5e9 -90 -30])



rec_down = resample(signal_temp,10e9,10e9);
signal_recover = (rec_down-0*mean(rec_down));

signal_kk = signal_recover+conj(signal_recover);



[fscale_kk,fre_kk] = pfft(fs1*1,signal_kk);


fig_spectrum_kk = figure();
plot(fscale_kk,fre_kk)
axis([-1e9 5e9 -90 -30])
title('kk-SSBI消除后信号图')