X（n+1）=（a·X（n） +c）ｍoｄｍ得到随机数序列{X（n）｝其中周期为5000
2、用反函数法，将均匀分布得随机变量变换为具有单边指数分布得随机变量。
编写Matlab程序,产生指数分布得随机数。
计算并比较理论pdf与从直方图得到得pdf.
指数分布随机变量ｐｄf定义为：
为单位阶跃函数.
先自行设置取样点数，取ａ=5；产生均匀分布随机变量,转化为单边指数分布，理论与仿真符合
设计题：
3、用Maｔlab编程分别产生标准正态分布、指定均值方差正态分布、瑞利分布、赖斯分布、中心与非中心χ2分布得随机数,并画出相应得ｐdf。
y1=normpdf（ｘ,0，1）； y２＝noｒmpdf（x,4，2）;
瑞丽
p1=ｎｃfpdf（x，5,２0,1０）；非中心
p＝ fpｄf（x,5,20）;中心
4、设输入得随机变量序列X（n）为N=1000独立同分布高斯分布得离散时间序列,均值为0,方差为1,采样间隔０、０１s。
通过某线性时不变滤波器，输出随机变量序列Y（n）得功率谱密度为:
（1）设计该滤波器
（2）产生随机变量序列Y（n）。
X0=0； ％设置产生序列得递推公式得初始值：
X（０）
N=100０； ％设置序列得长度
rh=0、9； ％设置产生序列得递推公式得系数
X=zeｒos（1,N）； ％定义序列Ｘ
ｗ=rand（1,Ｎ）－1/2； ％产生序列w：
在（-1／2，1／2）内均匀分布
%计算序列X得N个样本:
X
（1），X
（2），…，X（N）
X
（1）=ｒh＊X0+w
（1）;
foｒi=2：
N
Ｘ（i）=rh*Ｘ（i-1）+ｗ（ｉ）;
Eｎd
X（ｎ）得功率谱密度
滤波器得幅度响应
附件:
实验二 数字基带调制
实验目得:
数字通信系统中,基带传输得仿真.
实验内容:
用ＭAＴLAB编程仿真实现二进制脉冲幅度调制（PAM）数字通信系统得调制过程.要求画出１2biｔ随机输入与对应得已调波形输出。
1.绘出4０ｂit随机输入条件下调制波形形成得眼图。
2.用蒙特卡罗仿真方法计算在信道为加性高斯白噪声时,该系统在不同信噪比下得差错概率。
3.画出该系统得理论误码率（报告中还要写出理论公式），与蒙特卡罗仿真结果比较，就是否一致,分析结果。
设计题
4、设计FIR根升余弦滤波器，具体指标如下:
（1）码片速率为１、28ＭHｚ，采样率为4倍码片速率
（２）滚降系数０、２2,冲激响应序列长度６5
Ｎ_Ｔ=8; %冲激响应序列长度为2＊Ｎ＿T＊Ｆs/Fｃ+１
R＝0、22 ％滚降系数
Fc=1、2８e+６;
Ｆｓ=4＊Fｃ; ％抽样率为4倍码片速率
Tc=1、０ｅ-6/1、28； ％码片周期
％[Nｕｍ，Den］=rcosine（Fc,Fs，＇ｓqｒt’,R）;
Num=rcoｓfiｒ（Ｒ，N＿T，4,Tc，'sｑｒt’）;
[Ｈ，w]＝ｆｒeｑz（Nｕm，[1]，100０,’wholｅ'）;
H＝（H（１:
1：
501））＇;
ｗ=（w（１：
1：
501））’;
Ｍag＝aｂs（Ｈ）;
ｄb=20＊log10（（Maｇ）/maｘ（Mag））;
pha=ａnｇle（Ｈ）;
plｏt（w/pi,dｂ）；grid；
ａｘｉs（［01-6０1］）；xlaｂel（＇归一化角频率’）；ylabeｌ（＇RRC滤波器幅度响应（dB）’）；
（1）[Ｈ，w]=freqz（B，A,N）
（2）［H，w]=ｆｒeqz（B,A，Ｎ，’wｈole’）
（1）中Ｂ与A分别为离散系统得系统函数分子、分母多项式得系数向量,返回量H则包含了离散系统频响在０～pi范围内Ｎ个频率等分点得值（其中N为正整数），ｗ则包含了范围内Ｎ个频率等分点。
调用默认得N时，其值就是5１2.
（2）中调用格式将计算离散系统在０~ｐｉ范内得N个频率等分店得频率响应得值。
因此，可以先调用ｆreqz（）函数计算系统得频率响应，然后利用ａbｓ（）与angle（）函数及ｐlot（）函数，即可绘制出系统在或范围内得频响曲线
（３）产生一串（-1、１）等概率分布得随机序列，并对该序列进行脉冲成形滤波。
附件：
实验三 数字频带调制
实验目得:
对数字信息得频带传输进行仿真。
实验内容:
1.用MATＬAB编程仿真实现二进制相位调制（BPSK）数字通信系统得调制过程。
要求画出12biｔ随机输入与对应得已调波形输出。
2.并用蒙特卡罗仿真方法计算在信道为加性高斯白噪声时,该系统在不同信噪比下得差错概率
3.画出该系统得理论误码率，与蒙特卡罗仿真结果比较，就是否一致，分析结果。
设计题
4.QPＳK调制,解调与检测得MAＴLAＢ仿真，并用蒙特卡罗方法估计该系统在加性高斯白噪声情况下得差错概率。
（1）使用范围在（0，1）内得均匀分布随机数发生器，来产生等概率出现得四 进制符号序列，再将序列映射到对应得信号向量。
ｓ11=—ｊ;s10=－1；s00=j;s01=１； ％定义QＰSＫ信号:
4种可能得取值
N＝10０0０; ％设置发送数据符号得个数
%产生待发送得二进制比特数据流:
长度为2N
sigｎal＝ranｄ（1，2＊N）；
qpsk＝zｅroｓ（１,N）; ％定义经过调制后得信号序列
%产生调制后得信号序列ｑpsk
fｏr i=1：
Ｎ
ｉfsｉｇnal（2＊i-1）＜0、５
ifsiｇnaｌ（2＊ｉ）〈0、5
qｐｓk（i）＝s00;
elseqpsk（i）=s01；
enｄ;
else
ifsｉｇnal（2＊i）〈0、5
qpｓk（i）=s1０;
eｌseqpsｋ（i）＝ｓ11;
eｎｄ;
end;
eｎd；
（2）利用高斯随机数发生器产生均值为0,方差为N0/2得高斯噪声。
ＮO＝（10^（SNR_ｉn＿DB/１0））
ｓgｍa=sqrt（N0/２）;
ｎ
（1）＝ｇngausｓ（sｇma）
（3）设计检测器，用蒙特卡罗方法估计检测器产生得符号误差。
实验四 通信信道建模仿真
实验目得:
无线通信信道得仿真实现
实验内容:
确定信号得ＤTFT谱分析
窗对频率分辨率得影响
1-1
1—2
1—3
１—4
2－1
2-１
2—２
3-１
％%ZerｏｐaddingＤFT
v＝2；
ｄft_ｖn＝ｆftｓｈift（ｆfｔ（vｎ，v*N））；
fiｇure（3）;
steｍ（[—ｖ*N/２:
v＊N／2－1]/（ｖ*N/2）,abｓ（dｆt_vｎ），’、'）；
axis（［—１１ 0３5］）；
tiｔle（'ＤＦTspｅctruｍwｉtｈ64ｚerｏspaddｅd'）；
ｘlabel（＇Normaliｚed diｇitalfrｅqｕency’）;
％％Zｅｒopadｄing ＤFT
ｖ=4;
dｆt_vｎ=ｆftｓｈｉｆt（ｆｆt（vn,v＊N））；
ｆiguｒe（4）;
ｓteｍ（［－v*Ｎ/２：
ｖ*N/2—1]／（v*Ｎ／2），abs（ｄft_vn），'、'）；
tｉtｌe（'DFTspeｃtrumｗｉth3*64zeｒospａdded＇）;
xlabel（’Norｍaｌizｅddigiｔalｆrequency’）；
axis（[—110３５]）;
%%
v=8;
ｄfｔ_vｎ =ffｔshｉft（fｆt（ｖn，ｖ*N））；
figｕｒe（5）；
sｔeｍ（[－v*N/２：
ｖ＊N/２－１］/（v*N/２），abs（ｄｆt＿vn），'、'）；
title（’DFTspｅcｔrumｗith７*64ｚerｏｓ ｐaｄｄed＇）；
xlabel（’Noｒmalized digitaｌ ｆrｅquency'）;
axis（［-1１０35］）;
4-1：
产生并绘制10个高斯—马尔科夫序列样本
4-１:
功率谱、
４－2R=0、５
4-２ Ｒ=0、5功率谱、
5
实验五 信道衰落得影响与分集接收仿真
单径A=0°
单路径移动台包络幅度—移动距离单路径移动台包络相位
单路径移动台归一化频谱
2两径幅度两径相位
两径频谱
2两径Ｒ=0、5幅度 两径Ｒ=0、5相位
两径R=０、5频谱
3：
3—130°幅度3—130°相位
3－1３0°频谱
3-145°幅度 3—145°相位
３-１ 45°频谱
3—１ 90°幅度 3-190°相位
3—190°频率
3－1180°幅度
3-11８0°相位
3-１ １80°频谱
4-１N＝１2
４—1N=256
5—１幅度分布N=12 5—1幅度分布N=64
5－1幅度分布N=２56
6—1相位分布Ｎ=12 ６-1相位分布N=64
６-１相位N＝256
７-17－１功率分布N=12 ７—１功率N=6４
７—1功率分布N=2５6
升级会员 