江苏大学通信仿真软件课程设计报告

《江苏大学通信仿真软件课程设计报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、计算机学院通信0902通信仿真软件课程设计报告1、随机生成10个0～10之间的随机数，赋值为变量a，利用MATLAB程序类型中的“函数文件”（function）编写函数，以变量a为输入函数，对这10个数进行排序（升序或者降序），并输出结果到变量b中。操作步骤：1．主程序生成随机数2．子程序调用主程序中的随机数进行升序排列3．结果返回给主程序程序代码：%主程序a=10*rand(1,10);%产生10个0~10之间的随机数b=cmp(a)%调用子程序，将随机数排序后，结果返回给主程序中的变量b%子程序cmp.mfunctiony=cmp(x);form=1:1:9fork=(m+1):1

2、:10ifx(m)>x(k)x1=x(m);x(m)=x(k);x(k)=x1;endendendy=x;运行结果：b=Columns1through60.98411.07571.15822.38402.85542.8629Columns7through105.13836.47178.49569.12391、利用sin和cos函数，生成一个频率为100Hz的余弦波在10个周期内的信号（载波），再生成一个频率为10Hz的正弦波信号（调制信号），将两个信号相乘后，画出波形。（即模拟信号的调制过程）程序代码：t=0:0.001:0.1;y1=sin(200*pi*t);y2=cos(20*p

3、i*t);y3=y1.*y2;subplot(311)plot(t,y1)xlabel('t')ylabel('sin(200*pi*t)')subplot(312)plot(t,y2)xlabel('t')ylabel('cos(20*pi*t)')subplot(313)plot(t,y3);xlabel('t');ylabel('sin(200*pi*t).*cos(20*pi*t)');1、假设某数字通信系统收发信息速率为1kbps，发送端对数字信息进行ASK调制后，使用模拟线路进行传输，其中，载波频率为4kHz，数字“1”对应有载波，数字“0”对应无载波，接收端接收到信号后使

4、用载波信号为模板进行相关解调。考虑信道中存在不同的加性高斯白噪声的情况下，分析此时该系统的误码性能。操作步骤：1、根据假设，系统信息速率为1kbps，考察1s内系统的误码性能，即需要随机生成1000个“0”，“1”数据，每一位信息所占时间片大小为1ms。（round(rand(1,1000))）2、由于系统载波频率为4kHz，生成载波信号Ca=sin(2*pi*4000*t)，进行AM调制，若信息为“1”，载波幅度为1，若信息为“0”，载波幅度为0。例如，若第n位数据为“1”，则在第n位数据对应的时间片（（n-1）ms<=t<=nms）内为正弦载波信号，若每一个时间片内取100个采样点

5、，则对应信号为sin(2*pi*4000*t)，t=n-1+1/100:1/100:n；反之，为0。例如下图为10个离散数据：0，0，0，0，1，0，1，0，0，0，经过调制后的波形。1、信号调制后送入信道，假设该信道存在较大的加性高斯白噪声，使用awgn函数生成在接收端接收到的叠加了噪声的信号。（Rx_m=awgn(Tx_m,-5,0)，其中Tx_m为调制后的信号，Rx_m为在接收端接收的信号，-5表示信噪比为-5dB，0表示信号功率为0dBw）2、在接收端用载波信号为模板进行相关解调。（即将1个时间片（1ms）内的接收信号（100个）与载波信号逐点相乘后累加，与阈值（载波信号逐点相

6、乘后累加值的一半）相比，若大于阈值，则判定接收的信号为“1”，反之为“0”，阈值大小设为载波与载波相关值的一半）3、将接收到的数据与发送数据进行比较，找出，并计算发生错误的位的个数，计算此时的误码率。4、修改awgn函数中的信噪比，重复以上操作，分析误码率变换的规律，并说明原因。程序代码：N=1000;%1000个数据fc=4000;%载波频率fs=100;%每位取100个采样点ts=1/fs;%每个采样点占0.01mst0=1;%整个时间片段是1stc=t0/N;%每个比特是1mst=[0:tc:t0-tc];%整个时间段分为1000份ct=[0:ts/N:tc-ts/N];%Tim

7、eaxisofonebitdurationca=sin(2*pi*fc*ct);%载波信号test=ca.^2;limitation=sum(test);snr=-5;%信噪比bercnt=0;%起始误码个数为0xn=round(rand(1,1000));%1000个随机数据%每个数据进行模拟调制，Tx_m为调制后的信号，Rx_m为加噪后接收端的信号form1=1:Nif(xn(m1)==1)Tx_m((m1-1)*fs+1:(m1*fs))