资源描述：
《民航大学通信原理课程实验报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、通信原理实验通信原理Matlab实验报告-DSB信号的产生与相干解调指导教师：学号：姓名：8/9通信原理实验实验环境：Matlab2009CrackedbyWangyiwei程序设计思路：模块化设计：第一、DSB信号的产生部分第二、信号的调制与相干解调部分第三、理想低通滤波器设计部分第四、子程序：求傅里叶变换程序部分及程序相关注释说明：主程序部分：%**************************初始数据设计段*********************091143324echooffcloseallclct0=0.2;%信号的持续时间,用来定义时间向量ts=0.0
2、008;%抽样间隔fs=1/ts;%抽样频率fc=300;%载波频率t=[-t0/2:ts:t0/2];%时间向量df=0.25;%所需的频率分辨率m=sin(80*t);%调制信号c=3*cos(2*pi*fc.*t);%载波信号函数%************************信号的调制与解调段********************091143324u=m.*c;%已调信号k=u.*c;%相干解调时，DSB信号与同频通向的载波信号相乘[M,m,df1]=fftseq(m,ts,df);%对调制信号m(t)求傅里叶变换M=M/fs;%缩放[U,u,df1]=f
3、ftseq(u,ts,df);U=U/fs;[K,k,df1]=fftseq(k,ts,df);K=K/fs;%*********************理想低通滤波器设计段*********************091143324f_cutoff=100;%滤波器的截止频率n_cutoff=floor(100/df1);%设计的滤波器的带宽f=[0:df1:df1*(length(y)-1)]-fs/2;%频率分量8/9通信原理实验H=zeros(size(f));H(1:n_cutoff)=2*ones(1,n_cutoff);H(length(f)-n_cut
4、off+1:length(f))=2*ones(1,n_cutoff);%滤波器的频带函数AN=H.*K%相干解调后的信号K，与设计的100HZ滤波器相乘，保留包络an=real(ifft(AN))*fs;%求出an时域信号%**************************程序结果输出段**********************091143324figure(1)subplot(3,1,1)plot(f,fftshift(abs(M)))%fftshift：将FFT中的DC分量移到频谱中心title('原调制信号的频谱图')xlabel('频率f')subplo
5、t(3,1,2)plot(f,fftshift(abs(U)))title('已调信号的频谱图')xlabel('频率f')subplot(3,1,3)plot(f,fftshift(abs(H)))title('下图：低通滤波器的频谱图(理想f<100hz低通滤波器）')xlabel('频率f')figure(2)subplot(3,1,1)plot(f,fftshift(abs(K)))title('下图：DSB信号与同频同相信号相乘后的频谱')xlabel('频率f')subplot(3,1,2)plot(t,k(1:length(t)))title('下图：D
6、SB信号与同频同相信号相乘后的时域图')xlabel('时间t')subplot(3,1,3)plot(f,fftshift(abs(AN)))title('下图：DSB信号与同频同相信号相乘后，经过滤波H后的频谱')xlabel('频率f')figure(3)subplot(3,1,1)plot(t,m(1:length(t)))title('下图：原调制信号的时域图')xlabel('时间t')8/9通信原理实验subplot(3,1,2)plot(t,u(1:length(t)))title('下图：调制后信号时域图')xlabel('时间t')subplot(
7、3,1,3)plot(t,an(1:length(t)))title('下图：解调后信号的时域图')xlabel('时间t')%******************************************王乙伟*091143324*2011.12.14*子程序：function[M,m,df]=fftseq(m,ts,df)fs=1/ts;ifnargin==2%nargin为输入参量的个数n1=0;elsen1=fs/df;endn2=length(m);n=2^(max(nextpow2(n1),nextpow2(n2)));%nextp
