基于Matlab的模拟信号传输系统实验报告

实验七、基于Matlab 的模拟信号传输系统实验报告幅度的调制解调

1，AM 调制原理2，AM 的解调原理

Sam(t

）Sp(t)Sd(t)c(t)=coswc

t

相干解调的一般模型

%041210308AM 调制解调部分代码close all ;

clear;

clc;

……

figure;

yin=awgn(y1,3);

subplot(2,1,1);

plot(t,yin);title('加噪后信号');%%加噪声后的信号

axis([1e-34e-3-55]);

fcutsb=[16000175002250024000];%%接收机前端带通滤波

magsb=[010];

devsb=[0.050.010.05];

[nb,Wnb,betab,ftypeb]=kaiserord (fcutsb,magsb,devsb,fs);

hhb=fir1(nb,Wnb,ftypeb,kaiser(n b+1,betab),'noscale');%带通滤波器；[Hb,f]=freqz(hhb,1,1024,fs);st_pb=fftfilt(hhb,yin);subplot(2,1,2);plot(t,st_pb);title('加噪后频谱');axis([1e-34e-3-55]);

figure;

st_noise=st_pb.*xc;

subplot(2,1,1);

plot(t,st_noise);title('带通滤波后');

axis([1e-34e-3-22]);fcutsl=[250030000];%%混频后低通滤波magsl=[10];devsl=[0.010.05];[nl,Wnl,betal,ftypel]=kaiserord(fcutsl,magsl,devsl,fs);hhl=fir1(nl,Wnl,ftypel,kaiser(n l+1,betal),'noscale');st_pl=fftfilt(hhl,st_noise);subplot(2,1,2);plot(t,st_pl);title('带通滤波后频谱');axis([1e-34e-301]);grid;

LPF

角度的调制解调

FM的解调原理

窄带的相干解调原理框图如下所示：

Ｓfm(t)Si(t)

C（t）

窄带信号的相干解调

FM的非相干解调模型

%041210308FM调制解调部分代码

……w1=0;w2=0;

for m=1:length(t)

w1=mt(m)+w2;

w2=mt(m)+w1;

fi(m)=w1/(2*fs);

end

fi=fi*2*pi/max(abs(fi));

HI=cos(kfm*fi);

HQ=sin(kfm*fi);

yo=A*cos(2*pi*fc*t).*HI-A*sin(2 *pi*fc*t).*HQ;

subplot(2,1,1);

plot(t,yo);title('时域波形'); axis([1e-34e-3-22]);

Yo=fft(yo);

subplot(2,1,2);

plot(f,fftshift(abs(Yo))/1e6); title('频谱图');%%载波频谱

axis([-4e44e401]);

figure;

%%添加噪声

yoo=awgn(yo,30);

subplot(3,1,1);

plot(t,yoo);title('加噪后波形');%%加噪声后的信号

axis([1e-34e-3-22]);

%%带通滤波

KSband=2*(3+1)*f0;

fcutsb=[fc-KSband-2000

fc-KSband fc+KSband fc+KSband+2000];%%接收机前端带通滤波

magsb=[010];

devsb=[0.050.010.05];

[nb,Wnb,betab,ftypeb]=kaiserord (fcutsb,magsb,devsb,fs);

hhb=fir1(nb,Wnb,ftypeb,kaiser(n b+1,betab),'noscale');%带通滤波器；

st_pb=fftfilt(hhb,yoo);

%subplot(3,1,2);

st_pb=st_pb/1e6;

%plot(t,st_pb);

%axis([1e-34e-3-2e-62e-6]);

%figure;

%Yst_pb=fft(st_pb);

%subplot(2,1,1);

%plot(f,fftshift(abs(Yst_pb))); %%载波频谱

BPF LPF微分

%axis([-4e44e401e6]);

%figure('Name','FM解调波形')

%微分器设计

for i=1:length(t)-1%接收信号通过微分器处理

diff_st_pb(i)=(st_pb(i+1)-st_pb (i))/dt;

end

st_noise=abs(hilbert(diff_st_pb) );

subplot(3,1,2);

plot(t,[st_noise*200]);title('

解调后波形');

axis([1e-34e-304]);

%%隔直

%KSbandh=2*(3+1)*f0;

fcutsh=[0.013000];

magsh=[01];

devsh=[0.010.05];

[nh,Wnh,betah,ftypeh]=kaiserord (fcutsh,magsh,devsh,fs);

hhh=fir1(nh,Wnh,ftypeh,kaiser(n h+1,betah),'noscale');

st_out=fftfilt(hhh,st_noise*20) ;

subplot(3,1,3);

plot(t,[st_out0]);title('解调后波形');

axis([1e-34e-3-22]);