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

实验七、基于Matlab的模拟信号传输系统实验报告幅度的调制解调1，AM调制原理2，AM的解调原理Sam(t）Sp(t)Sd(t)c(t)=coswct相干解调的一般模型%041210308AM调制解调部分代码closeall;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(nb+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(nl+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;form=1:length(t)w1=mt(m)+w2;w2=mt(m)+w1;fi(m)=w1/(2*fs);endfi=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-2000fc-KSbandfc+KSbandfc+KSband+2000];%%接收机前端带通滤波magsb=[010];devsb=[0.050.010.05];[nb,Wnb,betab,ftypeb]=kaiserord(fcutsb,magsb,devsb,fs);hhb=fir1(nb,Wnb,ftypeb,kaiser(nb+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)));%%载波频谱BPFLPF微分%axis([-4e44e401e6]);%figure('Name','FM解调波形')%微分器设计fori=1:length(t)-1%接收信号通过微分器处理diff_st_pb(i)=(st_pb(i+1)-st_pb(i))/dt;endst_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(nh+1,betah),'noscale');st_out=fftfilt(hhh,st_noise*20);subplot(3,1,3);plot(t,[st_out0]);title('解调后波形');axis([1e-34e-3-22]);