民航大学通原理课程实验报告.docxVIP

通信原理实验 通信原理 Matlab实验报告-DSB信号的产生与相干解调 指导教师： 学号：姓名： 实验环境：Matlab 2009 Cracked by Wangyiwei 程序设计思路：模块化设计： DSB信号的产生部分 信号的调制与相干解调部分 理想低通滤波器设计部分 子程序：求傅里叶变换 程序部分及程序相关注释说明： 主程序部分： %**************************初始数据设计段*********************091143324 echo off close all clc t0=0.2; %信号的持续时间,用来定义时间向量 ts=0.0008; %抽样间隔 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); %载波信号函数 %************************信号的调制与解调段********************091143324 u=m.*c; %已调信号 k=u.*c; %相干解调时，DSB信号与同频通向的载波信号相乘 [M,m,df1]=fftseq(m,ts,df); %对调制信号m(t)求傅里叶变换 M=M/fs; %缩放 [U,u,df1]=fftseq(u,ts,df); U=U/fs; [K,k,df1]=fftseq(k,ts,df); K=K/fs; %*********************理想低通滤波器设计段*********************091143324 f_cutoff=100; %滤波器的截止频率 n_cutoff=floor(100/df1); % 设计的滤波器的带宽 f=[0:df1:df1*(length(y)-1)]-fs/2; %频率分量 H=zeros(size(f)); H(1:n_cutoff)=2*ones(1,n_cutoff); H(length(f)-n_cutoff+1:length(f))=2*ones(1,n_cutoff);%滤波器的频带函数 AN=H.*K %相干解调后的信号K，与设计的100HZ滤波器相乘，保留包络 an=real(ifft(AN))*fs; %求出an时域信号 %**************************程序结果输出段**********************091143324 figure(1) subplot(3,1,1) plot(f,fftshift(abs(M)))%fftshift：将FFT中的DC分量移到频谱中心 title(原调制信号的频谱图) xlabel(频率f) subplot(3,1,2) plot(f,fftshift(abs(U))) title(已调信号的频谱图) xlabel(频率f) subplot(3,1,3) plot(f,fftshift(abs(H))) title(下图：低通滤波器的 频谱图 (理想f100hz低通滤波器）) 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(下图：DSB信号与同频同相信号相乘后的 时域图) 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) subplot(3,1,2) plot(t,u(1:length(t))) title(下图：调制后信号时域图) xlabel(时间t) subplot(3,1,3) plot(t,an(1:length(t))) title(下图：解调后信号的时域图) xlabel(时间t) %*************************