- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
《MATLAB编程与系统仿真》课程考核说明讲解.doc
重 庆 交 通 大 学
学 生 综 合 设 计 报 告
课程名称 《MATLAB编程与系统仿真》
学 院 信息学院 年级 2012 专业班 通信工程3班
学生姓名 李春燕 学号 631206040323
成 绩
题目B
生成随机01序列n位(注意:n是一个确定的值,为你学号的后位所组成的整数),对生成的序列用(7,4)线性分组码进行编码,生成矩阵G如下:
假设信息速率为1Kbps,载波频率为10KHz,试对生成的随机01序列进行2ASK调制,并在同一图形界面上画出未调信号,已调波形(前10个比特)及调制之后功率谱图。对上述生成的2ASK信号滤波,请设计滤波器,并画出滤波后的波形请设计滤波器,并画出要求:
(1)将生成随机01序列的过程设计成函数,该函数的输入为01随机序列的位数,输出为01随机序列。
(2)将加线性分组码的过程设计成函数,输入为生成的01序列,输出为编码后的序列;
(3)将2ASK调制过程设计成函数,该函数的输入为01序列,载波频率及信息速率,输出为调制之后的信号数据输为调制之后的信号数据滤波器设计利用FDATool工具进行将滤波器参数用save命令保存到数据文件中。
(6)编写主程序调用前面的函数,完成题目中所要求的功能
一、【程序仿真流程图】
流程图1 生成随机01序列n位
流程图3 对信号进行2ASK调制
【主要步骤】
为方便理解和查看,在伪码后加入了少量主要步骤。
(1)
%用rand函数生成1行n(n=323,主函数中输入)列0~1之间的随机序列,然后用round函数对其四舍五 入得到0,1随 机序列A_order%
A_order=round(rand(1,n))
(2)
%用rem函数对0,1随机序列的长度除以4取余数n1,再用zeros函数对序列进行补零生成A_order1向量%
n2=zeros(1,4-n1);
A_order1=[A_order,n2];
%使之长度能够被4整除,然后再用reshape函数对其行和列进行重组转置为m向量%
Signal=reshape(A_order1,4,[]);
m=Signal;
%最后进行用rem对m*G求余,即得到线性分组码的编码后的信号%
(3)
%2ASK调制即是s(t)*cos(),载波信号为carrier=cos()%
%完整的载波周期个数=由每一位bit占用时间/载波的一个周期所占用时间,即 (1/1000bps)/(1/10000Hz)=10,即是一个bit有10个完整的波形周期%
M=carrierFrequency/bitRate;
%一个bit取N1=16个抽样点,取t从2*pi/N1到2*pi*10,步长为2*pi/N1,则t持续160个抽样点,即 是取前10个bit。则有100个完整的波形,用for循环实现。=(2*pi/N1)t%
t=2*pi/16:2*pi/16:2*pi*10;
carrier=[];
for j=1:M
carrier=[carrier,cos(t)];
end
%将输入信号carrier进行2ASK调制,为s(t)*cos(),s(t)是0,1比特流,当s(t)=1,即为原信 号;当s(t)=,0,则为0,用zeros函数将其变为0,最后一起放入transmittedSignal空矢量中,同 时在同一图形界面上画出输入信号、已调信号、功率谱%
transmittedSignal=[];
for i=1:length(Signal)
if (Signal(i)==1)
transmittedSignal=[transmittedSignal,carrier];
else
transmittedSignal=[transmittedSignal,zeros(1,length(carrier))];
end
end
(4)
%在调制信号上叠加平均功率为sig=3*0.3的噪声,在调制信号每一个点上加上0到1之间的随机数产生噪 声影响%
N2=length(transmittedSignal);
noi
文档评论(0)