文档介绍：通信原理实验一——实验报告
主要补充的代码如下：
B0=15000;%低通滤波器的带宽
B1=100;%尖峰滤波采用的带宽
hlpf=fir1(100,[B0]/(
通信原理实验一——实验报告
主要补充的代码如下：
B0=15000;%低通滤波器的带宽
B1=100;%尖峰滤波采用的带宽
hlpf=fir1(100,[B0]/(fs/2));%低通滤波器f1=38000;
hbpf=fir1(100,[f1-B0/2,f1+B0/2]/(fs/2));%带通滤波器f2=19000;
hcpf=fir1(100,[f2-B1/2,f2+B1/2]/(fs/2));%带通滤波器-尖峰滤波
s1_t=filter(hlpf,1,s_t);%低通滤波获得低通信
号
s2_t=filter(hbpf,1,s_t);%带通滤波获得带
通信号
s3_t=filter(hcpf,1,s_t);%带通滤波获得尖峰信
号
MusicLen=2^ceil(log2(MusicLen)+2);
df=fs/MusicLen;
fk=(-MusicLen/2:MusicLen/2-1)*df;
f_range=[-60000,60000];
s1_f=Ts*fftshift(fft(s1_t,MusicLen));%低通信号频谱
s2_f=Ts*fftshift(fft(s2_t,MusicLen));%带通信号频谱
s3_f=Ts*fftshift(fft(s3_t,MusicLen));%尖峰信号频谱
s_f=Ts*fftshift(fft(s_t,MusicLen));%原始信号频谱
%f_0=max(s3_f);
figure(1)
subplot(411),plot(fk,abs(s_f));title('原始音乐信号"s_t"的频谱');
subplot(412),plot(fk,abs(s1_f));title('低通滤波后"s1_t"的频谱');
subplot(413),plot(fk,abs(s2_f));title('带通滤波后"s2_t"的频谱');
subplot(414),plot(fk,abs(s3_f));title('尖峰信号的频谱');
s2_t=2.*s2_t.*s3_t.^2-s2_t;
%把
s2_t
进行频谱搬移
s2_t=filter(hlpf,1,s2_t);
%对
s2_t
进行低频滤波，得
到低频成分
s2_f=Ts*fftshift(ff