lmfhilbiir.m

来自「这是有关锁相环的一些必用知识.希望能对大家有所帮助哈」· M 代码 · 共 35 行

%希尔伯特变换
%for nn=1:2
clear

fs=40E+6;%采样频率
ts=1/fs;%采样时间间隔
k=10E+10;  %调频斜率
fo=8E+6;%载频
T=40E-6;%脉冲宽度
B=k*T;
n=4000;%抽样点数
t=ts:ts:n*ts;

p=zeros(1,n);
p(300+1:300+T*fs)=ones(1,T*fs);
signal1=1/sqrt(T)*exp(j*2*pi*(fo*t)).*p;

m=300;
mix_cum1=cum4est(signal1,300,1600,60,'unbiased',0,0);
mix_cum1(m+1)=mix_cum1(m+1)/10;
mix_cum1_spect=abs(fft(mix_cum1));
figure
subplot(2,1,1)
plot(abs(mix_cum1))
subplot(2,1,2)
plot(mix_cum1_spect)


[a,b]=hilbiir(ts,5*ts,50);
signal2=filter(a,b,signal1);

figure
plot(signal1)
figure
plot(signal2)