1、(0 X。m)使用迭代公式: X(n+1) = (aX(n) +c)modm得到随机数序列X(n)其中周期为50002、 用反函数法,将均匀分布的随机变量变换为具有单边指数分布的随机变量。编写Matlab程序,产生指数分布的随机数。计算并比较理论pdf和从直方图得到的pdf。指数分布随机变量pdf定义为:,为单位阶跃函数。先自行设置取样点数,取a=5;产生均匀分布随机变量,转化为单边指数分布,理论与仿真符合设计题:3、 用Matlab编程分别产生标准正态分布、指定均值方差正态分布、瑞利分布、赖斯分布、中心与非中心2分布的随机数,并画出相应的pdf。y1=normpdf(x,0,1); y2=n
2、ormpdf(x,4,2);瑞丽p1= ncfpdf(x,5,20,10);非中心p= fpdf(x,5,20);中心4、 设输入的随机变量序列X(n)为N=1000独立同分布高斯分布的离散时间序列,均值为0,方差为1,采样间隔0.01s。通过某线性时不变滤波器,输出随机变量序列Y(n)的功率谱密度为:(1) 设计该滤波器(2) 产生随机变量序列Y(n)。X0=0; %设置产生序列的递推公式的初始值:X(0)N=1000; %设置序列的长度rh=0.9; %设置产生序列的递推公式的系数X=zeros(1,N); %定义序列Xw=rand(1,N)-1/2; %产生序列w:在(-1/2,1/2)
3、内均匀分布 %计算序列X的N个样本:X(1),X(2),X(N) X(1)=rh*X0+w(1); for i=2:N X(i)=rh*X(i-1)+w(i);EndX(n)的功率谱密度滤波器的幅度响应附件:实验二 数字基带调制数字通信系统中,基带传输的仿真。用MATLAB编程仿真实现二进制脉冲幅度调制(PAM)数字通信系统的调制过程。要求画出12bit随机输入与对应的已调波形输出。1. 绘出40bit随机输入条件下调制波形形成的眼图。2. 用蒙特卡罗仿真方法计算在信道为加性高斯白噪声时,该系统在不同信噪比下的差错概率。3. 画出该系统的理论误码率(报告中还要写出理论公式),与蒙特卡罗仿真结果
4、比较,是否一致,分析结果。设计题4. 设计FIR根升余弦滤波器,具体指标如下:(1)码片速率为1.28MHz,采样率为4倍码片速率(2)滚降系数0.22,冲激响应序列长度65N_T=8; %冲激响应序列长度为2*N_T*Fs/Fc+1R=0.22 %滚降系数Fc=1.28e+6;Fs=4*Fc; %抽样率为4倍码片速率Tc=1.0e-6/1.28; %码片周期%Num,Den = rcosine(Fc,Fs,sqrt,R);Num=rcosfir(R,N_T,4,Tc,);H,w=freqz(Num,1,1000,wholeH=(H(1:1:501);w=(w(1:Mag=abs(H);db=
5、20*log10(Mag)/max(Mag);pha=angle(H);plot(w/pi,db);grid;axis(0 1 -60 1);xlabel(归一化角频率ylabel(RRC滤波器幅度响应(dB)(1)H,w=freqz(B,A,N)(2)H,w=freqz(B,A,N,whole)(1)中B和A分别为离散系统的系统函数分子、分母多项式的系数向量,返回量H则包含了离散系统频响在 0pi范围内N个频率等分点的值(其中N为正整数),w则包含了范围内N个频率等分点。调用默认的N时,其值是512。(2)中调用格式将计算离散系统在0pi范内的N个频率等分店的频率响应的值。因此,可以先调用f
6、reqz()函数计算系统的频率响应,然后利用abs()和angle()函数及plot()函数,即可绘制出系统在 或 范围内的频响曲线(3)产生一串(-1.1)等概率分布的随机序列,并对该序列进行脉冲成形滤波。实验三 数字频带调制对数字信息的频带传输进行仿真。用MATLAB编程仿真实现二进制相位调制(BPSK)数字通信系统的调制过程。1. 并用蒙特卡罗仿真方法计算在信道为加性高斯白噪声时,该系统在不同信噪比下的差错概率画出该系统的理论误码率,与蒙特卡罗仿真结果比较,是否一致,分析结果。2. QPSK调制,解调与检测的MATLAB仿真,并用蒙特卡罗方法估计该系统在加性高斯白噪声情况下的差错概率。(
7、1) 使用范围在(0,1)内的均匀分布随机数发生器,来产生等概率出现的四 进制符号序列,再将序列映射到对应的信号向量。s11=-j;s10=-1;s00=j;s01=1; %定义QPSK信号:4种可能的取值N=10000; %设置发送数据符号的个数%产生待发送的二进制比特数据流:长度为2Nsignal=rand(1,2*N);qpsk=zeros(1,N); %定义经过调制后的信号序列%产生调制后的信号序列qpskfor i=1: if signal(2*i-1)0.5 if signal(2*i) qpsk(i)=s00; else qpsk(i)=s01; end; else qpsk(i
8、)=s10; else qpsk(i)=s11;end;(2) 利用高斯随机数发生器产生均值为0,方差为N0/2的高斯噪声。NO=(10(SNR_in_DB/10)sgma=sqrt(N0/2); n(1)=gngauss(sgma)(3) 设计检测器,用蒙特卡罗方法估计检测器产生的符号误差。实验四 通信信道建模仿真无线通信信道的仿真实现确定信号的DTFT 谱分析窗对频率分辨率的影响1-11-21-31-42-12-23-1 % Zero padding DFT v=2;dft_vn = fftshift(fft(vn,v*N);figure(3);stem(-v*N/2:v*N/2-1/(v
9、*N/2),abs(dft_vn),.axis(-1 1 0 35);title(DFT spectrum with 64 zeros paddedNormalized digital frequencyv=4;figure(4);DFT spectrum with 3*64 zeros padded% v = 8; dft_vn = fftshift(fft(vn,v*N);figure(5);DFT spectrum with 7*64 zeros padded4-1:产生并绘制10 个高斯-马尔科夫序列样本功率谱.4-2 R=0.54-2 R=0.5功率谱.5实验五 信道衰落的影响与分集接收仿真单径A=0单路径移动台包络幅度-移动距离 单路径移动台包络相位单路径移动台归一化频谱2两径幅度 两径相位两径频谱2 两径R=0.5幅度 两径R=0.5相位两径R=0.5频谱3:3-1 30幅度 3-1 30相位频谱3-1 45幅度 3-1 453-1 90幅度 3-1 90频率3-1 180幅度4-1N=124-1N=2565-1幅度分布N=12 5-1幅度分布N=645-1幅度分布N=2566-1相位分布N=12 6-1相位分布N=646-1相位N=2567-17-1功率分布N=12 7-1功率N=647-1功率分布N=256如有侵权请联系告知删除,感谢你们的配合!
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1