AM与FM性能的比较Word格式.docx
《AM与FM性能的比较Word格式.docx》由会员分享,可在线阅读,更多相关《AM与FM性能的比较Word格式.docx(23页珍藏版)》请在冰豆网上搜索。
FM;
调制解调;
MATLAB;
目录
1课程设计目的…………………………………………………………1
2课程设计要求…………………………………………………………1
3相关知识………………………………………………………………1
4课程设计分析…………………………………………………………2
5仿真……………………………………………………………………4
6结果分析………………………………………………………………11
7参考文献………………………………………………………………17
1课程设计目的
(1)分别设计AM,FM,调制与解调通信系统,分别在理信道和非理想信道中运行,并得出仿真结果,分析比较两种种不同的调制方式所存在的不同。
(2)熟悉MATLAB文件中M文件的使用方法,包括函数、原理和方法的应用。
(3)增强在通信原理仿真方面的动手能力与自学能力。
2.课程设计要求
(1)系统建模
(2)确定仿真算法
(3)建立仿真模型
(4)设计仿真程序
(5)运行仿真程序
(6)输出仿真结果并进行分析
3.相关知识
通信系统是为了有效可靠的传输信息,信息由信源发出,以语言、图像、数据为媒体,通过电(光)信号将信息传输,由信宿接收。
通信系统又可分为数字通信与模拟通信。
基于课程设计的要求,下面简要介绍模拟通信系统。
信源是模拟信号,信道中传输的也是模拟信号的系统为模拟通信。
模拟通信系统的模型如图3-1所示。
图3-1模拟通信系统模型
调制器:
使信号与信道相匹配,便于频分复用等。
发滤波器:
滤除调制器输出的无用信号。
收滤波器:
滤除信号频带以外的噪声,一般设N(t)为高斯白噪声,则Ni(t)为窄带白噪声。
4.课程设计分析
4.1AM调制原理
幅度调制是用调制信号去控制高频正弦载波的幅度,使其按调制信号的规律变化的过程。
幅度调制器的一般模型如图4-1所示。
图4-1幅度调制模型
在图2-1中,若假设滤波器为全通网络(
=1),调制信号
叠加直流
后再与载波相乘,则输出的信号就是常规双边带(AM)调幅.AM调制器模型如图4-2所示
图4-2AM调制模型
AM信号波形的包络与输入基带信号
成正比,故用包络检波的方法很容易恢复原始调制信号。
但为了保证包络检波时不发生失真,必须满足
,否则将出现过调幅现象而带来失真。
AM信号的频谱是由载频分量和上、下两个边带组成(通常称频谱中画斜线的部分为上边带,不画斜线的部分为下边带)。
上边带的频谱与原调制信号的频谱结构相同,下边带是上边带的镜像。
显然,无论是上边带还是下边带,都含有原调制信号的完整信息。
故AM信号是带有载波的双边带信号,它的带宽信号带宽的两倍。
4.2FM调制原理
频率调制的一般表达式[4.1]为:
[4.1]
FM和PM非常相似,如果预先不知道调制信号的具体形式,则无法判断已调信号是调频信号还是调相信号。
[4.2]
[4.3]
图(2-3)所示的产生调频信号的方法称为直接调频法,图(4.4)所示的产生调频信号的方法称为间接调频法。
由于实际相位调制器的调节范围不可能超出
,因而间接调频的方法仅适用于相位偏移和频率偏移不大的窄带调制情形,而直接调频则适用于宽带调制情形。
根据调制后载波瞬时相位偏移的大小,可将频率调制分为宽带调频(WBFM)与窄带调频(NBFM)。
宽带与窄带调制的区分并无严格的界限,但通常认为由调频所引起的最大瞬时相位偏移远小于30°
时,
[4.4]
称为窄带调频。
否则,称为宽带调频。
为方便起见,无妨假设正弦载波的振幅A=1,则由式(4.1)调频信号的一般表达式,得
=
[4.5]
通过化解,利用傅立叶变化公式可得NBFM信号的频域表达式:
(4.6)
在NBFM中,由于下边频为负,因而合成矢量不与载波同相,而是存在相位偏移
,当最大相位偏移满足式[4.2]时,合成矢量的幅度基本不变,这样就形成了FM信号。
图4-1NBFM信号频谱
5.仿真
5.1.MATLAB源程序代码
‘am’
fm=100;
fc=500;
fs=5000;
Am=1;
A=2;
N=512;
K=N-1;
n=0:
N-1;
t=(0:
1/fs:
K/fs);
yt=Am*cos(2*pi*fm*t);
figure
(1)
subplot(1,1,1),plot(t,yt),title('
频率为3000的调制信号f1的时时域波'
);
y0=A+yt;
y2=y0.*cos(2*pi*fc*n/fs);
y3=fft(y2,N);
%fft变换
q1=(0:
N/2-1)*fs/N;
mx1=abs(y3(1:
N/2));
figure
(2)
subplot(2,1,1);
plot(t,y2);
title('
已调信号的时时域波'
subplot(2,1,2);
plot(q1,mx1);
f1已调信号的频谱'
%绘图
yc=cos(2*pi*fc*t);
figure(3)
subplot(2,1,1),plot(t,yc),title('
载波fc时域波形'
)
yc1=Am*cos(2*pi*fc*n/fs);
y3=fft(yc1,N);
q=(0:
mx=abs(y3(1:
subplot(2,1,2),plot(q,mx),title('
载波fc频谱'
y4=0.01*randn(1,length(t));
%用RANDN产生高斯分布序列
w=y4.^2;
%噪声功率
figure(4)
plot(t,y4);
高斯白噪声时域波形'
y5=fft(y4,N);
q2=(0:
mx2=abs(y5(1:
subplot(2,1,2),plot(q2,mx2),title('
高斯白噪声频域波形'
y6=y2+y4;
figure(5)
subplot(2,1,1),plot(t,y6),title('
叠加后的调制信号时域波形'
q3=q1;
mx3=mx1+mx2;
subplot(2,1,2),plot(q3,mx3),title('
叠加后的调制信号频谱波形'
)%调制
yv=y6.*yc;
%乘以载波进行解调
Ws=yv.^2;
p1=fc-fm;
[k,Wn,beta,ftype]=kaiserord([p1fc],[10],[0.050.01],fs);
%Fir数字低通滤波
window=kaiser(k+1,beta);
%使用kaiser窗函数
b=fir1(k,Wn,ftype,window,'
noscale'
%使用标准频率响应的加窗设计函数
yt=filter(b,1,yv);
yssdb=yt.*2-2;
figure(6)
subplot(2,1,1),plot(t,yssdb),title('
经过低通已调信号的时域波形采样'
y9=fft(yssdb,N);
mx=abs(y9(1:
经过低通已调信号频域波形'
)%解调
ro=y9-yt;
W=(yt.^2).*(1/2);
R=W/w
r=W/ro
G=r/R
‘fm’
dt=0.001;
%设定时间步长
t=0:
dt:
1.5;
%产生时间向量
am=5;
%设定调制信号幅度
fm=5;
%设定调制信号频率
mt=am*cos(2*pi*fm*t);
%生成调制信号
fc=50;
%设定载波频率
ct=cos(2*pi*fc*t);
%生成载波
kf=10;
%设定调频指数
int_mt
(1)=0;
fori=1:
length(t)-1
int_mt(i+1)=int_mt(i)+mt(i)*dt;
%求信号m(t)的积分
end%调制,产生已调信号
sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt);
%调制信号
添加高斯白噪声
sn1=10;
%设定信躁比(小信噪比)
sn2=30;
%设定信躁比(大信噪比)
sn=0;
%设定信躁比(无信噪比)
db=am^2/(2*(10^(sn/10)));
%计算对应的高斯白躁声的方差
n=sqrt(db)*randn(size(t));
%生成高斯白躁声
nsfm=n+sfm;
%生成含高斯白躁声的已调信号(信号通%过信道传输)
FM解调
length(t)-1%接受信号通过微分器处理
diff_nsfm(i)=(nsfm(i+1)-nsfm(i))./dt;
end
diff_nsfmn=abs(hilbert(diff_nsfm));
%hilbert变换,求绝对值得到瞬时幅度(包络检波)
zero=(max(diff_nsfmn)-min(diff_nsfmn))/2;
diff_nsfmn1=diff_nsfmn-zero;
时域到频域转换
ts=0.001;
%抽样间隔
fs=1/ts;
%抽样频率
df=0.25;
%所需的频率分辨率,用在求傅里叶变换时,它表示FFT的最小频率间隔
对调制信号m(t)求傅里叶变换
m=am*cos(2*pi*fm*t);
%原调信号
ifnargin==2
n1=0;
else
n1=fs/df;
n2=length(m);
n=2^(max(nextpow2(n1),nextpow2(n2)));
M=fft(m,n);
m=[m,zeros(1,n-n2)];
df1=fs/n;
%以上程序是对调制后的信号u求傅里变换
M=M/fs;
%缩放,便于在频铺图上整体观察
f=[0:
df1:
df1*(length(m)-1)]-fs/2;
%时间向量对应的频率向量
对已调信号u求傅里变换
n2=length(sfm);
U=fft(sfm,n);
u=[sfm,zeros(1,n-n2)];
%以上是对已调信号u求傅里变换
U=U/fs;
%缩放
disp('
按任意键可以看到原调制信号、载波信号和已调信号的曲线'
pause
subplot(3,1,1);
plot(t,mt);
%绘制调制信号的时域图
xlabel('
时间t'
调制信号的时域图'
subplot(3,1,2);
plot(t,ct);
%绘制载波的时域图
载波的时域图'
subplot(3,1,3);
plot(t,sfm);
%绘制已调信号的时域图
已调信号的时域图'
%******************************************
按任意键可以看到原调制信号和已调信号在频域内的图形'
%************figure
(2)*********************
subplot(2,1,1)
plot(f,abs(fftshift(M)))%fftshift:
将FFT中的DC分量移到频谱中心
频率f'
原调制信号的频谱图'
subplot(2,1,2)
plot(f,abs(fftshift(U)))
已调信号的频谱图'
按任意键可以看到原调制信号、无噪声条件下已调信号和解调信号的曲线'
%绘制已调信号的时域图
无噪声条件下已调信号的时域图'
nsfm=sfm;
%绘制无噪声条件下解调信号的时域图
plot((1:
length(diff_nsfmn1))./1000,diff_nsfmn1./400,'
r'
无噪声条件下解调信号的时域图'
%*****************************************
按任意键可以看到原调制信号、小信噪比高斯白噪声条件下已调信号和解调信号已调信号的曲线'
db1=am^2/(2*(10^(sn1/10)));
%计算对应的小信噪比高斯白躁声的方差
n1=sqrt(db1)*randn(size(t));
nsfm1=n1+sfm;
%生成含高斯白躁声的已调信号(信号通
%过信道传输)
diff_nsfm1(i)=(nsfm1(i+1)-nsfm1(i))./dt;
diff_nsfmn1=abs(hilbert(diff_nsfm1));
diff_nsfmn1=diff_nsfmn1-zero;
plot(1:
length(diff_nsfm),diff_nsfm);
%绘制含小信噪比高斯白噪声已调信号的时域图
含小信噪比高斯白噪声已调信号的时域图'
%绘制含小信噪比高斯白噪声解调信号的时域图
含小信噪比高斯白噪声解调信号的时域图'
按任意键可以看到原调制信号、大信噪比高斯白噪声条件下已调信号和解调信号已调信号的曲线'
db1=am^2/(2*(10^(sn2/10)));
%计算对应的大信噪比高斯白躁声的方差
%生成含高斯白躁声的已调信号(信号通过信道传输)
%hilbert变换,求绝对值得到瞬时幅度(包%络检波)
length(diff_nsfm1),diff_nsfm1);
%绘制含大信噪比高斯白噪声已调信号
%的时域图
含大信噪比高斯白噪声已调信号的时域图'
%绘制含大信噪比高斯白噪声解调信号%的时域图
含大信噪比高斯白噪声解调信号的时域图'
6结果分析
AM的调制解调信号波形
图6-1频率为3000的调制信号f1的时时域波
图6-2已调信号的时时域波
图6-3载波fc时域波形
图6-4高斯白噪声时域和频域波形
图6-5叠加后的调制信号时域和频谱波形
图6-6经过低通已调信号时域和频域波形
FM的调制解调信号波形
图6-7FM调制信号
图6-8原调制信号与已调信号频谱图
图6-9调制解调信号时域图
图6-10调制信号与较小高斯白噪声加入后时域图
图6-11调制信号与较大高斯白噪声加入后时域图
经过本次设计比对,可以得出以下结论
AM就是指调幅制广播,也就是指中波及短波广播。
这种无线电波的载频幅度,是随着调制的音频信号而变,优点是接收机的构造简单,能接收远距离的短波段广播及国内普遍采用的中波广播,每个电台的带宽是9KHz(双边带,音频带宽是4500Hz)。
但是由于干扰信号主要是对无线电波的幅度起作用,所以抗干扰能力差、由于带宽较窄,音质较差,适合语音广播。
FM指调频制广播。
载频信号的频率,是随调制的音频信号而变,幅度不变。
优点是频带宽,每个电台达上百KHz,音质好;
干扰信号对广播信号幅度的影响,可以用限幅器予以削除,所以抗干扰能力很强,几乎无杂音。
由于频带很宽,整个中波频段也容不下几个电台,所以只能用超短波进行广播。
超短波段信号只能以直线方式传播,只能为本地听众服务,本地能接收的FM电台数量有限。
7.参考文献
[1]樊昌信,曹丽娜.通信原理(第6版)[M].国防工业出版社,2012.3
[2]孙详,徐流美,吴清等.MATLAB7.0基础教程[M].清华大学出版社,2006.5
[3]韩利竹,王华等.MATLAB电子仿真与应用(第2版)[M].国防工业出版,2003.6
[4]郭文彬,桑林.通信原理-基于Matlab的计算机仿真[M].北京邮电大学出版社,2006.5
[5]曹志刚,钱亚生.现代通信原理[M].清华大学出版社,2002.3
[6]郭仕剑等.MATLAB7.x数字信号处理[M].人民邮电出版社,2006.10
[7]张辉,曹丽娜编著.通信原理学习指导[M].西安电子科技大学,2003.6