通信原理课程设计报告
学 院 专 业
班 级 学 号
学生姓名 指导教师
课程成绩 完成日期 2016年X月X日
课程设计成绩评定
学 院 专 业
班 级 学 号
学生姓名 指导教师
课程成绩 完成日期 2016年X月X日
指导教师学生课程设计中评价
评分项目
优
良
中
格
格
课程设计中创造性成果
学生掌握课程容程度
课程设计完成情况
课程设计动手力
文字表达
学态度
规范求
课程设计文质量
指导教师课程设计评定意见
综合成绩 指导教师签字 2016年1月8日
课程设计务书
城南学院 通信工程专业
课程名称
通信原理课程设计
时间
20152016学年第学期17~19周
学生姓名
指导老师
题 目
基MATLABASK调制解调实现
容:
利MATLAB集成环境M文件编写程序实现ASK调制解调求采样频率360HZ绘制出解调前时域频域波形叠加噪声时解调前时频波形根运行结果波形分析该解调程正确性信道信号传输影响
求:
1)熟悉MATLAB中M文件方法掌握ASK调制解调原理基础编写出ASK调制解调程序
2)绘制出ASK信号解调前时域频域中波形观察解调前频谱变化加深ASK信号解调原理理解
3)信号叠加噪声进行解调绘制出解调前信号时频波形改变噪声功率进行解调分析噪声信号传输造成影响
4)老师指导求独立完成课程设计全部容求编写课程设计学年文正确阐述分析设计实验结果
应提交文件:
(1)课程设计学年文
(2)课程设计附件
基MATLABASK调制解调实现
学生姓名: 指导老师:
摘 MATLAB美国MathWorks公司生产科学工程计算专门设计交互式型软件课程设计容利MATLAB集成环境M文件编写程序实现ASK调制解调求采样频率360HZ绘制出解调前时域频域波形叠加噪声时解调前时频波形根运行结果波形分析该解调程正确性信道信号传输影响目熟悉MATLAB中M文件方法掌握ASK调制解调原理基础编写出2ASK调制解调程序绘制出ASK信号解调前时域频域中波形观察解调前频谱变化信号叠加噪声变化终着输入信号噪声增加增误码越严重结加深ASK信号解调原理理解
关键词 ASK调制解调时域谱频域谱高斯白噪声信噪
1 引 言
通信原理通信工程专业门重专业课通信工程专业续专业课基础掌握通信原理课程知识仅坚实专业基础提高处理通信系统问题力素质通课程设计ASK振幅键控调制解调进步理解数字通信基础理助加深通信原理理解
11课程设计目
通设计基MATLABASK调制解调实现深入理解掌握二进制ASK调制解调噪声信号传输影响[1]
通信原理理知识基础加深ASK调制解调设计原理实现方法理解通信信号波形频谱深刻认识仅加强课知识解涉MATLAB编程语言软件基操作常识[2]
掌握调制解调函数应增强动手实践力
12课程设计求
根题目查阅关资料掌握数字带通调制技术通信原理
学会MATLAB中M文件编写方法应MATLAB软件中调制解调函数掌握ASK调制解调原理根原理编写出ASK调制解调程序
绘制出ASK信号解调前时域频域中波形观察解调前频谱变化理解ASK信号解调原理
二进制基带信号叠加噪声解调绘制出解调前信号时频波形然改变噪声功率进行解调记录分析分析噪声信号传输造成影响
根求独立完成课程设计学年文正确阐述分析设计结果出结
13课程设计步骤
产生数字基带信号绘制时域谱频域谱
设置载波频率绘制时域谱频域谱
信号进行数字调制绘制时域谱频域谱
已调信号进行解调绘制时域谱频域谱
已调信号加入高斯噪声绘制时域谱频域谱
加噪声信号进行解调绘制时域谱频域谱
已调信号加入高斯噪声绘制时域谱频域谱
加噪声信号进行解调绘制时域谱频域谱
较信噪时误码率
2 ASK调制解调原理
21 ASK调制原理
ASK指振幅键控方式种调制方式根信号调节正弦波幅度幅度键控通法器开关电路实现载波数字信号10控制通断信号1状态载波接通时传输信道载波出现信号0状态载波关断时传输信道载波传送接收端根载波原出数字信号10二进制幅度键控信号频带宽度二进制基带信号宽度两倍幅移键控法(ASK)载波幅度着调制信号变化简单形式载波二进制调制信号控制通断 时称作开关键控法(OOK)二进制幅度键控记作2ASK2ASK利代表数字信息01基带矩形脉键控连续载波载波时断时续输出载波输出时表示发送1载波输出时表示发送02ASK信号表示
式中载波角频率s(t)单极性NRZ矩形脉序列
中g(t)持续时间高度1矩形脉常称门函数二进制数字
2ASKOOK信号产生方法通常两种:模拟调制(相器法)键控法模拟幅度调制方法法器实现相应调制图21图22:
法器
图21模拟相法
图22键控开关法
22 ASK解调原理
2ASKOOK信号两种基解调方法:非相干解调(包络检波法)相干解调(步检测法)课程设计求相干解调图23:
带通滤波器
相器
低通滤波器
抽样判决器
定时脉
输出
图23相干解调
3 仿真实现程
31 ASK信号产生
图31二进制基带信号时域谱频域谱
先源程序创建M文件定义路径编写语句xceil(rand(1100000)05) 生成段机二进制基带信号中rand产生机矩阵ceil取整函数编写语句FFT1fft(x128) FFT1abs(FFT1)机序列进行傅里叶变换取绝值时域谱频域谱图31
32载波信号波形
图32载波信号时域谱频域谱
调制解调系统中载波信号频率般信号源频率信号源频率12 Hz载波频率设置36 Hz编写正弦函数carrycos(2*pi*Fc*t)进行频域转换FFT2fft(carry256) FFT2abs(FFT2)时域谱频域谱图32
33 ASK调制解调实现
编写数字带通调制函数ydmod(xFcFdFs'ask'2)进行2ASK调制中x输入信号Fc载波频率Fd码速率Fs采样频率FsFd必须正整数进行频域转换FFT3fft(y256) FFT3abs(FFT3)时域谱频域谱图33:
图33已调信号时域谱频域谱
图34解调信号解调时域谱频域谱
编写数字带通调制函数zddemod(yFcFdFs'ask'2)进行2ASK解调中x输入信号Fc载波频率Fd码速率Fs采样频率进行频域转换FFT4fft(z64) FFT4abs(FFT4)调制解调信道中没加入噪声调制前信号(原始信号)调制信号相时域谱频域谱图34
34叠加噪声ASK调制解调
图35叠加噪声调制信号时域谱频域谱
图36叠加噪声调制信号时域谱频域谱
信道中噪声叠加信号噪声始终存通常称加性噪声加性干扰已调信号y分叠加高斯噪声(信噪6)噪声(信噪2)编写Ynt1awgn(y6)Ynt2awgn(y2)进行频域转换时域谱频域谱分图35图36
图37叠加噪声解调信号时域谱频域谱
图38叠加噪声解调信号时域谱频域谱
噪声Ynt1噪声信号Ynt2分解调编写数字带通调制函数z1ddemod(Ynt1FcFdFs'ask'2)z2ddemod(Ynt2FcFdFs'ask'2)进行频域转换时域谱频域谱分图37图38
35误差分析
图39叠加叠加噪声噪声误差数误差率
图310误差数误差率
误码率衡量数字通信系统性重指标信道高斯白噪声干扰二进制2ASK数字调制误码率取决解调器输入信噪编写误码率函数[brPe1]symerr(xz1)[brPe2]symerr(xz2)中br符号误差数Pe1符号误差率绘制图39
调函数semilogy(SNRPe)绘制信噪误码率关系曲线图310出:噪声时(误码率零)相较信噪较时噪声误码率低反信噪较时噪声误码率高
4 遇问题解决办法
次课程设计运MATLAB集成环境M文件实现ASK调制解调编写程中遇种问题通探索老师学帮助解决总结:
产生二进制机序列时没01情况显示三角波解答办法:采ceil取整函数01序列stairs方波画图函数方波显示出
序列0时调制波形反振幅零解决办法:调制函数基础编写循环判断语句序列周期作次循环序列0时调制函数赋值零反波形消失
绘制信噪误码率关系曲线时错误解决办法:调数字调制函数dmod加入噪声Ynt3awgn(ySNR(i))调制函数基础编写循环判断语句显示已调波形
5 结束语
课程设计培养学生综合运学知识发现提出分析解决实际问题锻炼实践力重环节学生实际工作力具体训练考察程次课程设计中信号ASK调制解调程深刻理解掌握ASK设计设计思路设计流程更容易理解通信原理提高通信系统分析力解决力机会通信原理进行模拟仿真专业说尤重三周时间基MatlabASK调制解调课程设计收获基掌握Matlab中函数期间利通信原理MATLAB相关书籍网络资源老师学帮助完成务意识数字信号处理重性时次课程设计中遇许困难时发现解决问题断改进方法实际操作基础理准断完善断创新提高效率事通信类职业提供宝贵验设计课题达举反三步增强实践动手力受益匪浅
该设计然期完成功基满足处理需力限果意方需做进步改善
参考文献
[1] 樊昌信曹丽娜 通信原理 北京:国防工业出版社2012
[2] 黄文梅熊桂林杨勇 信号分析处理—MATLAB语言应 长沙:国防科技学出版社2000
[3] 唐宏岳恒立郑雪峰 MATLAB电子信息类课程中应 北京:电子工业出版社20068
[4] 邓华 MATLAB通信仿真应实例详解 民邮电出版社2003
附录:ASK调制解调程序清单
程序名称:ASKm
程序功:ASK解调解调
程序作者:庄维
修改时间:201618
程序代码
clc 清命令工作窗里容
clear 清存空间变量
产生二进制机序列
xceil(rand(1100000)05) 产生二进制机序列取x 整数
figure(1) 窗口1包含时域谱频域谱
subplot(211) 分块图函数subplot图形窗口分成2块子窗口第1图
stairs(x) 第2图
xlabel('时间 t') x轴标注
ylabel('序列值') y轴标注
title('二进制机序列') 添加图标题
axis([1 21 1 2]) 控制坐标轴范围
grid on 图中添加栅格
机序列进行频谱分析
FFT1fft(x128) 机序列进行傅里叶变换
FFT1abs(FFT1) 傅里叶变换取绝值
figure(1)
subplot(212) 第2图
plot(FFT1)
xlabel('频率 f')
ylabel('幅度 FFT1')
title('机序列频谱')
axis([0 128 0 50])
grid on
载波信号
t1360136020 载波时间范围
Fc36 载波频率
carrycos(2*pi*Fc*t) 正弦载波信号
figure(2) 窗口2包含时域谱频域谱
subplot(211)
plot(carry)
xlabel('时间 t')
ylabel('幅度 carry')
title('载波信号')
axis([1 600 2 2])
grid on
载波信号进行频谱分析
FFT2fft(carry256) 载波信号进行傅里叶变换
FFT2abs(FFT2) 傅里叶变换取绝值
figure(2)
subplot(212)
plot(FFT2)
xlabel('频率 f')
ylabel('幅度 FFT2')
title('载波信号频谱')
axis([0 256 0 100])
grid on
ASK调制
Fd12 Fd码速率Fs采样频率
Fs360
ydmod(xFcFdFs'ask'2)调数字带通调制函数dmod进行2ASK调制
for i120
if x(i)0
yy(30*(i1)+130*i)0
else
yy(30*(i1)+130*i)y(30*(i1)+130*i)
end
end
20机码元进行判码元0该码元周期调制信号零
figure(3)
subplot(211)
plot(yy)
xlabel('时间 t')
ylabel('幅度 y')
title('已调信号')
axis([1 600 2 2])
grid on
已调信号进行频谱分析
FFT3fft(y256) 已调信号进行傅里叶变换
FFT3abs(FFT3) 傅里叶变换取绝值
figure(3)
subplot(212)
plot(FFT3)
xlabel('频率 f')
ylabel('幅度 FFT3')
title('已调信号频谱')
axis([0 256 0 50])
grid on
ASK解调
zddemod(yFcFdFs'ask'2) 调数字带通调制函数dmod进行2ASK解调
figure(4) 傅里叶变换取绝值
subplot(211)
stairs(z)
xlabel('时间 t')
ylabel('幅度 z')
title('解调信号')
axis([1 21 1 2])
grid on
解调信号进行频谱分析
FFT4fft(z64) 解调信号进行傅里叶变换
FFT4abs(FFT4) 傅里叶变换取绝值
figure(4)
subplot(212)
plot(FFT4)
xlabel('频率 f')
ylabel('幅度 FFT4')
title('解调信号频谱')
axis([0 64 0 50])
grid on
加入高斯噪声SNR6
Ynt1awgn(y6) 加入高斯噪声信噪6
figure(5)
subplot(211)
plot(Ynt1)
xlabel('时间 t')
ylabel('幅度 Ynt1')
title('加噪声信号')
axis([1 600 2 2])
grid on
加噪声信号进行频谱分析
FFT5fft(Ynt1256) 加入噪声调制信号进行傅里叶变换
FFT5abs(FFT5) 傅里叶变换取绝值
figure(5)
subplot(212)
plot(FFT5)
xlabel('频率 f')
ylabel('幅度 FFT5')
title('加噪声信号频谱')
axis([0 256 0 50])
grid on
ASK加噪声信号解调误码率
z1ddemod(Ynt1FcFdFs'ask'2)
调数字带通调制函数dmod加噪声信号进行解调
[brPe1]symerr(xz1)
解调加噪声信号误码分析br符号误差数Pe1符号误差率
figure(6)
subplot(211)
stairs(z1)
xlabel('时间 t')
ylabel('幅度 z1')
title('加噪声解调信号')
axis([1 21 1 2])
grid on
加噪声解调信号进行频谱分析
FFT6fft(z164) 加入噪声解调信号进行傅里叶变换
FFT6abs(FFT6) 傅里叶变换取绝值
figure(6)
subplot(212)
plot(FFT6)
xlabel('频率 f')
ylabel('幅度 FFT6')
title('加噪声解调信号频谱')
axis([0 64 0 50])
grid on
加入高斯噪声SNR2
Ynt2awgn(y3) 加入高斯噪声信噪2
figure(7)
subplot(211)
plot(Ynt2)
xlabel('时间 t')
ylabel('幅度 Ynt2')
title('加噪声信号')
axis([1 600 2 2])
grid on
加噪声信号进行频谱分析
FFT7fft(Ynt2256) 加入噪声调制信号进行傅里叶变换
FFT7abs(FFT7) 傅里叶变换取绝值
figure(7)
subplot(212)
plot(FFT7)
xlabel('频率 f')
ylabel('幅度 FFT5')
title('加噪声信号频谱')
axis([0 256 0 50])
grid on
ASK加噪声信号解调误码率
z2ddemod(Ynt2FcFdFs'ask'2)
调数字带通调制函数dmod加噪声信号进行解调
[brPe2]symerr(xz2)
解调加噪声信号误码分析br符号误差数Pe1符号误差率
figure(8)
subplot(211)
stairs(z2)
xlabel('时间 t')
ylabel('幅度 z2')
title('加噪声解调信号')
axis([1 21 1 2])
grid on
加噪声解调信号进行频谱分析
FFT8fft(z264) 加入噪声解调信号进行傅里叶变换
FFT8abs(FFT8) 傅里叶变换取绝值
figure(8)
subplot(212)
plot(FFT8)
xlabel('频率 f')
ylabel('幅度 FFT6')
title('加噪声解调信号频谱')
axis([0 64 0 50])
grid on
误码分析
SNR1010
for i1length(SNR)
Ynt3awgn(ySNR(i)) 加入高斯噪声信噪10dB10dB
Zddemod(Ynt3FcFdFs'ask'2)
调数字带通解调函数ddemod加噪声信号进行解调
[br Pe(i)]symerr(xZ)
解调加噪声信号误码分析br符号误差数Pe(i)符号误差率
end
figure(9)
semilogy(SNRPe) 调semilogy函数绘制信噪误码率关系曲线
xlabel('信噪 SNR(rdB)')
ylabel('误码率 Pe')
title('信噪误码率关系')
axis([10 10 0 1])
grid on
文档香网(httpswwwxiangdangnet)户传
《香当网》用户分享的内容,不代表《香当网》观点或立场,请自行判断内容的真实性和可靠性!
该内容是文档的文本内容,更好的格式请下载文档