信抗大作业.

ata=askzerodb; data=(askdata);% 第j行所有信号的希尔伯特变化%求信号序列的希尔伯特变换 hask=hilbert(data);% 求信号序列的希尔伯特变换%求瞬时幅度%a_amplitude=sqrtm(y*y.+hask*hask); a_amplitude=abs(hask);% 瞬时幅度 figure(2) plot(b); title(' 瞬时幅度'); phi_phase=angle(hask); %求相位修正值 c_k(1)=0;% 修正相位初始值 for i=2:n if i+1>n break;end if phi_phase(i)-phi_phase(i-1)>pi c_k(i)=c_k(i-1)-2*pi; elseif phi_phase(i)-phi_phase(i+1)>pi c_k(i)=c_k(i-1)+2*pi; else c_k(i)=c_k(i-1); end end c_k(n)=phi_phase(n); c_k; %去卷叠相位 phi_uw_phase=phi_phase+c_k;% 去卷叠相位%去相位线性分量后真正的瞬时相位 phi_NL_phase=phi_uw_phase-2*pi*(1:n)*fc/fs;% 去相位线性分量后真正的瞬时相位 figure(3) plot(phi_phase); title(' 瞬时相位'); %求瞬时频率 f_frequency=1/(2*pi)*diff(phi_uw_phase);% 瞬时频率 figure(4) plot(f_frequency); title(' 瞬时频率'); 信号截图与各瞬时参量截图如下: 0 50 100 150 -1 -0. 5 0 0. 5 1 0 50 100 150 0 1000 2000 3000图1 :时域波形和频域波形