matlab 模拟调相,基于matlab的模拟通信系统调相实现.doc

基于matlab的模拟通信系统调相实现.doc

青 岛 农 业 大 学 理学与信息科学学院 通 信 原 理 课 程 报 告 设 计 题 目 调相实现 学生专业班级 电子信息工程 学生姓名(学号) 指 导 教 师 完 成 时 间 2009年6月12日 实 习(设计)地点 信息楼112 2011年6月12日基于 Matlab 的模拟通信系统——调相实现 一、课程设计目的 本次课程设计是对通信原理课程理论教学和实验教学的综合和总结。通过 这次课程设计，使同学认识和理解通信系统，掌握信号是怎样经过发端处理、 被送入信道、然后在接收端还原。 要求学生掌握通信原理的基本知识，运用所学的通信仿真的方法实现某种 传输系统。能够根据设计任务的具体要求，掌握软件设计、调试的具体方法、 步骤和技巧。对一个实际课题的软件设计有基本了解，拓展知识面，激发在此 领域中继续学习和研究的兴趣，为学习后续课程做准备。 二、课程设计内容 这个课程设计，我主要完成了以下工作： 1、在信号源产生了一个初始的阶跃信号。 2、将阶跃信号加载到余弦载波上 3、进行 FFT(快速傅立叶变换) ，从频域的角度来研究调制波。 4、模拟自然界的噪声，给调制信号加上白噪声。 5、由于白噪声属于高频信号，所以在接收端通过低通滤波器滤除噪声，并还原 出初始信号。 三、设计原理 1、初始信号和载波信号通过调幅形成已调波在空间传输。 2、通过 Matlab 中的白噪声函数在信号中引入白噪声。 3、信号从时域转换到频域的公式(傅里叶变换和逆变换)4、解调：正弦波幅度解调 标准调幅信号的解调可以不用拍频振荡器。调幅信号中的载波实际上起了 拍频振荡波的作用，利用非线性元件实现频率变换，经低通滤波即得到与调幅 信号包络成对应关系的输出。这种方法属于非相干解调。单边带信号的解调需要一个频率和相位与被抑制载波完全一致的正弦振荡 波。使这个由接收机复原的载波和单边带信号相乘，即可实现解调。这种方式 称为同步检波，也称为相干解调。 四、实现方法 1、发生模块 (1)产生阶跃信号 代码： clear all; t0=0.15; ts=0.0001; fc=250; a=0.85;%调制指数 snr=20;%SNR in db %m(t) t=0:ts:2*t0/3; t1=t0/3;%在 t1 处跳变 m1=stepfun(t,t1); m=-3*m1+1; %m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)];\ figure(1); subplot(311); plot(t,m(1:length(t))) title( m(t)信号波形 ); axis([0 2/3*t0 -2.5 1.5])(2)产生载波信号 代码： %c(t) c=cos(2*pi*fc.*t); figure(1); subplot(312); plot(t,c(1:length(t))) title( c(t)载波 ) axis([0 2/3*t0 -1.5 1.5]) (3)将阶跃信号与载波信号相乘，得到调制信号。 代码： %u(t) u=(2+a*m).*c; figure(1); subplot(313); plot(t,u(1:length(t))) title( u(t)调制波 ) axis([0 2/3*t0 -3.5 3.5])(4)将初始信号从时域转换为频域，从频域来对初始信号进行分析 代码： %m(f) mf=fft(m); figure(2); subplot(211); plot(t,mf); title( m(f)信号频谱 ); axis([0 0.01 -2.5 5]) (5)将载波从时域转换为频域，从频域来对载波信号进行分析 代码： %u(f) uf=fft(u); figure(2); subplot(212) plot(t,uf) title( u(f)调制频谱 ); %signal_power dt=0.01; signal_power=sum(u.*u)*dt/(length(t)*dt)2、传输模块 模拟在传输过程中的噪声干扰，对已调信号加上白噪声。 代码： %noise_power snr_lin=10^(snr/10); noise_power=signal_power/snr_lin %u_n(t)=u+n u_n = awgn(u,snr); figure(3); subplot(211); plot(t,u_n); title( u+n 波形 ); %u_n(f) u_nf=fft(u_n); figure(3); subplot(212); plot(t,u_nf) title( u+n 频谱 );3、接收模块 在接受模块中，我们要完成两部分操作 (1)去除噪声 由于白噪声属于高频噪声，所以在这里，我们选择低通滤波器来去除噪声 代码： %滤波前数据 t = 0:0.01:10; f=5*sin(2*pi*20*t)+5*sin(2*pi*60*t); %滤波器设计 Fpass = 1; % Passband Frequency Fstop = 50; % Stopband Frequency Apass = 1; % Passband Ripple (dB) Astop = 60; % Stopband Attenuation (dB) Fs = 2000; % Sampling Frequency h = fdesign.lowpass( fp,fst,ap,ast , Fpass, Fstop, Apass, Astop, Fs); Hd = design(h, equiripple , . MinOrder , any , . StopbandShape , flat ); %f2滤波后的数据 f2 = filter(Hd,f); %对比图 plot(t,f, b ,t,f2, r ); (2)将原信号与载波分离，还原出原始信号 代码： %示范课题 clear all; t0=0.15; ts=0.0001; fc=250; a=1;%调制指数 snr=10;%SNR in db %m(t) t=0:ts:2*t0/3; t1=t0/3;%在 t1 处跳变 m1=stepfun(t,t1); m=-2*m1+1; figure(4);subplot(311); plot(t,m(1:length(t))) title( m(t)信号波形 ); axis([0 2/3*t0 -1.5 1.5])