SSB.rar_SSB 希尔伯特_matlab希尔伯特_ssb_ssb调制_单边带

global dt df N T f t close all dt=0.002; %取样时间间隔 N=2^15; %取样点数 T=N*dt; %截短时间 df=1/T; %频域取样间隔 Bs=N*df/2; %系统带宽 f=[-Bs+df/2:df:Bs]; %频域横坐标 t=[-T/2+dt/2:dt:T/2]; figure(1) set(1,'Position',[200,200,600,600]) fs=10; %载波频率 fm=1; %调制信号频率 sm=sin(2*pi*fm*t); %调制信号 Sm=t2f(sm); %调制信号频谱 s=sin(2*pi*fm*t).*cos(2*pi*fs*t)-(-cos(2*pi*fm*t)).*sin(2*pi*fs*t); %单边带上边频调制信号.sin(2*pi*fm*t)的希尔伯特变换为-cos(2*pi*fm*t); S=t2f(s);%已调信号频谱 s0=s.*cos(2*pi*fs*t); %解调 S0=t2f(s0);%解调信号频谱 %lpf=2*sin(2*pi*2*t)./(2*pi*2*t) %LPF=t2f(lpf) Y=LPF(2).*S0 %滤波输出频谱 y=f2t(Y);%滤波后输出的波形 %以下为作图： subplot(3,2,1) plot(t,sm) xlabel('t(ms)'); ylabel('s(t)'); axis([-0.5,0.5,min(sm),max(sm)]); title('输入调制信号波形') subplot(3,2,2) ASm=abs(Sm); plot(f,ASm) xlabel('f(Mhz)') axis([-3,3,min(ASm),max(ASm)]) title('调制输入信号频谱') subplot(3,2,3) plot(t,s) xlabel('t(ms)') ylabel('s') title('SSB上边频已调信号') axis([-0.5,0.5,min(s),max(s)]) subplot(3,2,4) AS=abs(S); xlabel('t(ms)') ylabel('S') title('已调信号频谱') hold on plot(f,AS) axis([-12,12,min(AS),max(AS)]) subplot(3,2,5); y=real(y); plot(t,y) xlabel('t(ms)') ylabel('y') axis([-0.5,0.5,min(y),max(y)]) title('解调输出信号波形') AY=abs(Y); subplot(3,2,6) plot(f,AY) title('解调输出信号频谱') axis([-3,3,min(AY),max(AY)])