




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、 HUNAN UNIVERSITY 通信原理 实验学生姓名 阳铭宇 学生学号 201308030119 专业班级 通信1301 _ _指导老师 余小游 2016 年5 月 5 日实验一 模拟调制Ø 实验预习l 实验目的1. 理解模拟调制的5个基本性质:已调信号的时域表示、已调信号的频域表示、已调信号的带宽、已调信号的功率含量、解调后的信噪比。2. 掌握模拟调制(幅度调制、角调制)的基本原理并且通过仿真加深理解。l 实验原理AM信号时域表达式 式中 m(t) 调制信号,均值为0; A0 常数,表示叠加的直流分量。频域表达式DSB信号时域表达式:无直流分量频域表达式:无载频分量SSB信号
2、双边带信号两个边带中的任意一个都包含了调制信号频谱M(w)的所有频谱成分,因此仅传输其中一个边带即可。这样既节省发送功率,还可节省一半传输频带,这种方式称为单边带调制时域表达式u 保留上边带u 保留下边带频谱表达式FM信号时域表达式频域表达式² 宽带调频:² 窄带调频(NBFM)Ø 实验报告1. DSB-AM调制代码:clear all;close all;clc;t0=0.15;%信号持续时间ts=0.001;%采样周期fc=250;%载波频率snr=20;%信噪比fs=1/ts;%采样频率df=0.3;%频率分辨率t=0:ts:t0;snr_lin=10(sn
3、r/10);%线性信噪比%信息信号m=ones(1,t0/(3*ts),-2*ones(1,t0/(3*ts),zeros(1,t0/(3*ts)+1);c=cos(2*pi*fc.*t);%载波信号u=m.*c;%调制后信号M,m,df1=fftseq(m,ts,df);%傅里叶变换M=M/fs;U,u,df1=fftseq(u,ts,df);%调制后信号的傅里叶变换U=U/fs;C,c,df1=fftseq(c,ts,df); % 载波的傅里叶变换f=0:df1:df1*(length(m)-1)-fs/2; % 频率向量signal_power=spower(u(1:length(t);
4、 % 已调信号功率noise_power=signal_power/snr_lin; % 计算噪声功率noise_std=sqrt(noise_power); % 计算噪声功率的标准差noise=noise_std*randn(1,length(u); % 产生高斯噪声r=u+noise; % 将噪声加到已调信号中R,r,df1=fftseq(r,ts,df); % 混合信号频谱R=R/fs; % 缩放比例signal_power;clf;figure(1);subplot(2,2,1);plot(t,m(1:length(t);xlabel('Time');title(
5、9;The message signal');subplot(2,2,2);plot(t,c(1:length(t);xlabel('Time');title('The carrier');subplot(2,2,3);plot(t,u(1:length(t);xlabel('Time');title('The modulated signal');figure(2);subplot(2,1,1);plot(f,abs(fftshift(M);xlabel('Frequency');title('S
6、pectrum of the message signal');subplot(2,1,2);plot(f,abs(fftshift(U);title('Spectrum of the modulated signal');xlabel('Frequency');figure(3);subplot(2,1,1);plot(t,noise(1:length(t);title('Noise sample'); xlabel('Time');subplot(2,1,2);plot(t,r(1:length(t);title(
7、39;Signal and noise');xlabel('Time');figure(4);subplot(2,1,1);plot(f,abs(fftshift(U);title('Signal spectrum');xlabel('Frequency');subplot(2,1,2);plot(f,abs(fftshift(R); title('Signal and noise spectrum');xlabel('Frequency');实验结果:2. DSB-AM 解调代码:t0=0.15; % 信
8、号持续时间ts=1/1500; % 采样间隔fc=250; % 载波频率fs=1/ts; % 采样频率df=0.3; % 所需的频率分辨率t=0:ts:t0; % 时间向量m=ones(1,t0/(3*ts),-2*ones(1,t0/(3*ts),zeros(1,t0/(3*ts)+1); % 消息信号c=cos(2*pi*fc.*t); % 载波u=m.*c; % 已调信号y=u.*c; % 已调信号与载波相乘M,m,df1=fftseq(m,ts,df); % 傅立叶变换M=M/fs; % 缩放比例 U,u,df1=fftseq(u,ts,df); % 傅立叶变换U=U/fs; % 缩放
9、比例 Y,y,df1=fftseq(y,ts,df); % 傅立叶变换Y=Y/fs; % 缩放比例 f_cutoff=150; % 滤波器的截止频率n_cutoff=floor(150/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);DEM=H.*Y; % 滤波器输出信号的频谱dem=real(ifft(DEM)*fs; % 滤波器输出subplot(
10、3,1,1)plot(f,fftshift(abs(M)title('消息信号的幅度谱')xlabel('频率')subplot(3,1,2)plot(f,fftshift(abs(U)title('已调信号的幅度谱')xlabel('频率')subplot(3,1,3)plot(f,fftshift(abs(Y)title('混频器输出信号的幅度谱')xlabel('频率')figuresubplot(2,1,1)plot(f,fftshift(abs(H)title('低通滤波器频率特性
11、')xlabel('频率')subplot(2,1,2)plot(f,fftshift(abs(DEM)title('解调输出信号的幅度谱')xlabel('频率')figuresubplot(2,1,1)plot(f,fftshift(abs(M)title('消息信号的幅度谱')xlabel('频率')subplot(2,1,2)plot(f,fftshift(abs(DEM)title('解调输出信号的幅度谱')xlabel('频率')figuresubplot(2,1
12、,1)plot(t,m(1:length(t)title('时域消息信号')xlabel('时间')subplot(2,1,2)plot(t,dem(1:length(t)title('解调输出的时域信号')xlabel('时间')实验结果:2. SSB-AM 调制代码:clear all;close all;clc;t0=0.15;%信号持续时间ts=0.001;%采样周期fc=250;%载波频率snr=10;%信噪比fs=1/ts;%采样频率df=0.25;%频率分辨率t=0:ts:t0;snr_lin=10(snr/10);
13、%线性信噪比% the message vectorm=ones(1,t0/(3*ts),-2*ones(1,t0/(3*ts),zeros(1,t0/(3*ts)+1);c=cos(2*pi*fc.*t); % 载波数组udsb=m.*c; % DSB调制信号UDSB,udssb,df1=fftseq(udsb,ts,df); % 傅里叶变换UDSB=UDSB/fs; % 缩放比例f=0:df1:df1*(length(udssb)-1)-fs/2; % 频率数组n2=ceil(fc/df1); % 本地载波% 除去上边带UDSB(n2:length(UDSB)-n2)=zeros(size
14、(UDSB(n2:length(UDSB)-n2);ULSSB=UDSB; % 产生LSSB-AM的频谱M,m,df1=fftseq(m,ts,df); % 傅里叶变换M=M/fs; % 缩放比例u=real(ifft(ULSSB)*fs; % 产生LSSB信号signal_power=spower(udsb(1:length(t)/2; % % 计算信号功率noise_power=signal_power/snr_lin; % 计算噪声功率noise_std=sqrt(noise_power); % 计算噪声功率标准差noise=noise_std*randn(1,length(u); %
15、产生噪声r=u+noise; % 在信号中加入噪声R,r,df1=fftseq(r,ts,df); % 傅里叶变换R=R/fs; % 缩放比例signal_powerclf;figure(1);subplot(2,1,1);plot(t,m(1:length(t);axis(0,0.15,-2.1,2.1);xlabel('Time');title('The message signal');subplot(2,1,2);plot(t,c(1:length(t); xlabel('Time');title('The carrier'
16、;); figure(2);subplot(2,1,1);plot(0:ts:ts*(length(u)-1)/8,u(1:length(u)/8);xlabel('Time');title('The LSSB-AM modulated signal') ;subplot(2,1,2);plot(f,abs(fftshift(ULSSB);xlabel('Frequency');title('Spectrum of the LSSB-AM modulated signal');figure(3);subplot(2,1,1);pl
17、ot(f,abs(fftshift(M);xlabel('Frequency');title('Spectrum of the message signal');subplot(2,1,2);plot(f,abs(fftshift(ULSSB);xlabel('Frequency');title('Spectrum of the LSSB-AM modulated signal');figure(4);subplot(2,1,1);plot(t,noise(1:length(t);title('Noise sample
18、39;) ;xlabel('Time');subplot(2,1,2);plot(t,r(1:length(t);title('Modulated signal and noise');xlabel('Time');figure(5);subplot(2,1,1);plot(f,abs(fftshift(ULSSB);title('Modulated signal spectrum');xlabel('Frequency');subplot(2,1,2);plot(f,abs(fftshift(R);title(&
19、#39;Modulated signal noise spectrum');xlabel('Frequency');4. SSB-AM 解调代码:t0=0.15; % 信号持续时间ts=1/1500; % 采样间隔fc=250; % 载波频率fs=1/ts; % 采样频率df=0.25; % 所需的频率分辨率t=0:ts:t0; % 时间向量m=ones(1,t0/(3*ts),-2*ones(1,t0/(3*ts),zeros(1,t0/(3*ts)+1); % 消息信号c=cos(2*pi*fc.*t); % 载波udsb=m.*c; % 已调的双边带信号UDSB,
20、udsb,df1=fftseq(udsb,ts,df); % 傅立叶变换UDSB=UDSB/fs; % 缩放比例n2=ceil(fc/df1); % 载波在频率向量中的位置UDSB(n2:length(UDSB)-n2)=zeros(size(UDSB(n2:length(UDSB)-n2); % 去掉双边带的上边带ULSSB=UDSB; % 产生下边带频谱M,m,df1=fftseq(m,ts,df); % 消息信号的频谱M=M/fs; f=0:df1:df1*(length(M)-1)-fs/2; % 频率向量u=real(ifft(ULSSB)*fs; % 从频谱中得到下边带信号y=u.
21、*cos(2*pi*fc*0:ts:ts*(length(u)-1);%混频Y,y,df1=fftseq(y,ts,df); % 混频器输出的频谱Y=Y/fs; % 缩放比例f_cutoff=150; % 选择滤波器的截止频率n_cutoff=floor(150/df); % 设计滤波器H=zeros(size(f); H(1:n_cutoff)=4*ones(1,n_cutoff); H(length(f)-n_cutoff+1:length(f)=4*ones(1,n_cutoff); DEM=H.*Y; % 滤波器输出信号的频谱dem=real(ifft(DEM)*fs; % 滤波器输出
22、subplot(3,1,1)plot(f,fftshift(abs(M)title('消息信号的频谱')xlabel('频率')subplot(3,1,2)plot(f,fftshift(abs(ULSSB)title('已调信号的频谱')xlabel('频率')subplot(3,1,3)plot(f,fftshift(abs(Y)title('混频器输出的幅度谱')xlabel('频率')figure实验结果:5. 频率调制(FM)代码:t0=0.15; % 信号持续时间ts=0.0005; %
23、 采样间隔fc=200; % 载波频率fs=1/ts; % 采样频率df=0.25; % 所需的频率分辨率kf=50; % 偏差常数t=0:ts:t0; % 时间向量m=ones(1,t0/(3*ts),-2*ones(1,t0/(3*ts),zeros(1,t0/(3*ts)+1); % 消息信号int_m(1)=0;for i=1:length(t)-1 % m的积分 int_m(i+1)=int_m(i)+m(i)*ts;endM,m,df1=fftseq(m,ts,df); % 傅立叶变换M=M/fs; % 缩放比例f=0:df1:df1*(length(m)-1)-fs/2; % 频
24、率向量u=cos(2*pi*fc*t+2*pi*kf*int_m); % 调频信号U,u,df1=fftseq(u,ts,df); U=U/fs; subplot(2,1,1)plot(t,m(1:length(t)axis(0 0.15 -2.1 2.1)xlabel('时间')title('消息信号波形')subplot(2,1,2)plot(t,u(1:length(t)axis(0 0.15 -2.1 2.1)xlabel('时间')title('已调信号波形')figuresubplot(2,1,1)plot(f,abs
25、(fftshift(M) xlabel('频率')title('消息信号的幅度谱')subplot(2,1,2)plot(f,abs(fftshift(U) title('已调信号的幅度谱')xlabel('频率')subplot(3,1,1)plot(f,fftshift(abs(Y)title('混频器输出的幅度谱')xlabel('频率')subplot(3,1,2)plot(f,fftshift(abs(H)title('低通滤波器')xlabel('频率')
26、subplot(3,1,3)plot(f,fftshift(abs(DEM)title('解调输出的幅度谱')xlabel('频率')figuresubplot(2,1,1)plot(t,m(1:length(t)title('时域的消息信号')xlabel('时间')subplot(2,1,2)plot(t,dem(1:length(t)title('解调输出的时域波形')xlabel('时间')6. FM 解调代码:t0=0.2; % 信号持续时间ts=0.001; % 采样间隔fc=250;
27、% 载波频率snr=20; % 信噪比(dB)fs=1/ts; % 采样频率df=0.25; % 所需的频率分辨率kf=100; % 偏差常数t=-t0/2:ts:t0/2; % 时间向量m=sinc(100*t); % 消息信号int_m(1)=0;for i=1:length(t)-1 int_m(i+1)=int_m(i)+m(i)*ts;endM,m,df1=fftseq(m,ts,df); % 傅立叶变换M=M/fs; % 缩放比例f=0:df1:df1*(length(m)-1)-fs/2; % 频率向量Int_m=int_m(1:length(t);u=cos(2*pi*fc*t
28、 + 2*pi*kf*Int_m); % 调频信号U,u,df1=fftseq(u,ts,df); U=U/fs; v,phase=env_phas(u,ts,250); % 解调,找出u的相位phi=unwrap(phase); % 保存初始相位dem=(1/(2*pi*kf)*(diff(phi)/ts); % 解调器输出、 区分和扩展相位subplot(2,1,1)plot(t,m(1:length(t)xlabel('时间')title('消息信号波形')subplot(2,1,2)plot(t,u(1:length(t)xlabel('时间
29、39;)title('已调信号波形')figuresubplot(2,1,1)plot(f,abs(fftshift(M) xlabel('频率')title('消息信号的幅度谱')subplot(2,1,2)plot(f,abs(fftshift(U) title('已调信号的幅度谱')xlabel('频率')实验结果: 思考题:一 1.调幅和调频在无线电中是最常见的,调幅是使高频载波的频率随信号改变的调制(AM)。其中,载波信号的振幅随着调制信号的某种特征的变换而变化。调频是使载波频率按照调制信号改变的调制(F
30、M)。已调波频率变化的大小由调制信号的大小决定,变化的周期由调制信号的频率决定调频的抗干扰能力强,失真小,但服务半径小、 调频比调幅抗干扰能力强2、调频波比调幅波频带宽 频带宽度与调制系数有关,即:调制系数大,频带宽。调频中常取调频系数大于1,而调幅系数是小于1的(大于1就过调了,会严重失真),所以,调频波的频带宽度比调幅波的频带宽度大得多。因此音质要好,音频的上限可达10KHZ左右3、调频制功率利用率大于调幅制 发射总功率中,边频功率为传送调制信号的有效功率,而边频功率与调制系数有关,调制系数大,边频功率大。由于调频系数mf大于调幅系数ma,所以,调频制的功率利用率比调幅制高。4、由于调频波
31、比调幅波频带宽,所以也带来选择性的问题,所以不适合在射频的低频段使用,往往使用100M左右的波段。在这个波段的特点直线传播,因此传播的半径相对小些二代码:t=0:0.001:2;m=ones(size(t);m=(t>0&t<=0.1).*0.1+(t>=0.1&t<1).*t+(t>=1&t<1.9).*(-t+2)+(t>=1.9&t<2).*0.1;subplot(2,2,1)plot(t,m)xlabel('原信号')subplot(2,2,2)f1,fftm=T2F(t,m);plot(f
32、1,abs(fftm)xlabel('原信号频谱')kf=25fc=20*1000y=cos(2*pi*fc*t+kf*trapz(t,m)subplot(2,2,3)plot(t,y)xlabel('FM信号')f,ffty=T2F(t,y);subplot(2,2,4)plot(f,ffty)xlabel('FM信号频谱')a=min(m);b=max(m);fa=fc+kf*a/(2*pi);fb=fc+kf*b/(2*pi);sf_max=max(abs(ffty);B=sum(abs(ffty).2)*0.1*fc*fc/sf_max.
33、2;B_done=2*B+2*kf;实验结果:三代码:t0=2; %信号持续时间ts=0.001; %采样时间fs=1/ts;fc=50; %载波频率df=0.3; %分辨率t=0:ts:t0; %时间向量m=(t>=0.1&t<1).*t+(t>=1&t<1.9).*(-t+2); %消息信号c=cos(2*pi*fc.*t); %调制载波y=m.*c; %已调信号w=y.*c; %已调信号乘以相干载波%滤波器Rp=0.1; %通带最大衰减Rs=40; %阻带最小衰减Wp=25/100; %通带截止频率Ws=75/100; %阻带起始频率n,Wn=el
34、lipord(Wp,Ws,Rp,Rs);b,a=ellip(n,Rp,Rs,Wn);Sh=filter(b,a,w); %信号经过椭圆低通滤波器%求消息信号的频谱Mk=fft(m,512); %傅里叶变换Mw=2*pi*abs(fftshift(Mk);fw1=(-256:255)/512*1000;%求已调信号的频谱Yk=fft(y,512); %傅里叶变换Yw=2*pi*abs(fftshift(Yk);fw2=(-256:255)/512*1000;%计算信噪比snrs1=1/0.001;snrs2=1/0.01;snrs3=1/0.05;snrs4=1/0.1;snrs5=1/0.3;
35、snr1=10*log10(snrs1);snr2=10*log10(snrs2);snr3=10*log10(snrs3);snr4=10*log10(snrs4);snr5=10*log10(snrs5);%*power_y=sum(y.2)/length(y); %信号功率%*%噪声功率power_noise1=power_y/(10(snr1/10);power_noise2=power_y/(10(snr2/10);power_noise3=power_y/(10(snr3/10);power_noise4=power_y/(10(snr4/10);power_noise5=powe
36、r_y/(10(snr5/10);%*%加噪信号yn1=y+sqrt(power_noise1)*randn(size(y); yn2=y+sqrt(power_noise2)*randn(size(y); yn3=y+sqrt(power_noise3)*randn(size(y); yn4=y+sqrt(power_noise4)*randn(size(y); yn5=y+sqrt(power_noise5)*randn(size(y); %*%解调w1=yn1.*c; %已调信号乘以相干载波Sh1=filter(b,a,w1); %信号经过椭圆低通滤波器w2=yn2.*c;Sh2=filter(b,a,w2);w3=yn3.*c;Sh3=filter(b,a,w3);w4=yn4.*c;Sh4=filter(b,a,w4);w5=yn5.*c;Sh5=filter(b,a,w5);subplot(2,1,1);plot(t,m);xlabel('Time');title('消息信号');grid onsubplot(2,1,2);plot(t,y);xlabel('Time'
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 二零二五年度分手协议书:精神损害赔偿范本
- 二零二五年度男方道歉夫妻关系修复合同
- 二零二五年度地下车库停车位租赁及夜间停车服务合同
- 二零二五年度贷款装修合同附件说明
- 2025年度电影收益权与游戏改编授权合作协议
- 游泳馆2025年度老年人健康免责协议模板
- 2025年度车辆过户责任及风险分担协议书
- 2025年度特色民宿预订定金合同规范文本
- 2025年度跨境电商加盟合作协议
- 2025-2030年发动机起动机项目商业计划书001
- 专利管理制度管理办法
- 拖拉机和联合收割机驾驶人身体条件证明
- 机电控制与可编程序控制器课程设计
- 广联达钢筋输入规则
- 基于ADAMS的悬置刚度仿真指南
- 弯矩二次分配法EXCEL计算
- 分数小数混合运算练习80题
- 童话故事《老鼠搬鸡蛋》.ppt
- 偏差管理和CAPA_王丽丽
- 北京大学数学物理方法经典课件第五章——傅里叶变换
- 消防安全知识壁报-08灭火器作用及使用方法识别1
评论
0/150
提交评论