参数估计：λ=1，计算10000点随机数。 程序如下： n=10000; m=1;
x=rand(n,1); for i=1;n if (x(i)==0) x(i)=0.0001; end end u=log(x); y=-(1/m)*u; e=mean(y); d=var(y); hist(y,50);
title('指数分布直方图');
text(0,-170,strcat('均值是',num2str(e))); text(3,-170,strcat('方差是',num2str(d))); 运行结果如下：
均值和方差为：0.99154和0.97058，而指数分布的数学期望计算式为：1/m＝1，与上面所求出的随机数的平均值0. 99154相当接近，指数分布方差的计算公式为：(1/m)的平方，当m=1时代入计算得1与0. 97058相当符合。
作业三：雷达系统仿真
雷达系统仿真流程为：
目标回波 天 线 高放 混频 中放 相干 地物杂波 系统噪声 A/D A/D 脉 冲 压 缩 缓存 MT I 取模 积累 CFAR
雷达目标回波：
这里的雷达回波由三部分组成即：目标回波、地物杂波、系统噪声。其中目标回波这里采用线性调频信号模拟，地物杂波采用相关对数-正态分布的随机数，系统噪声采用高斯白噪声。回波信号由目标回波（动目标），地物杂波（静目标），及系统高斯白噪声组成。
线性调频信号：x=rect(t/mk)exp(jπkt2) (k=B/mk)
目标回波：y＝rect(t/mk)*exp(j*2*pi*((f1+k*t/2).*t+fd*i*T)) 地物杂波（静目标）：y＝rect(t/mk)*exp(j*2*pi*((f1+k*t/2).*t)) 系统噪声（高斯白噪声）： z＝0.2*randn(1,N)。
1．雷达的回波模型：
回波
目标回波 杂波回波 十 噪声 加噪声的目标回波的图形：
2．放大器：
其中高频放大器的增益为k，带通滤波器的设计采用的是FIR滤波器，其中通频带的带宽取为信号的带宽B的5/3倍,过渡带的宽度为B，采用汉明窗进行加权，数字滤波器的通频带和阻带的值利用w??T（w为数字角频率，?为模拟角频率，T为采样周期）计算。
高放图：
3．混频器：
X（t）
参考源 X 带通滤波器 其中x(t)?cos(2?f1t)，参考源为：cos?2?f2t?，两路信号相乘以后的输出为：11cos[2?(f1?f2)t]?cos[2?(f1?f2)t]22，这里取带通滤波器的通带的中心频率为1cos[2?(f1?f2)t]f1?f2于是通过带通滤波器以后的输出为：2，于是目标回波
经高频放大器以后通过混频器后的时域图形频域以及滤波器的频率响应分别如下图所示：
4．中频放大器的设计：
中频放大器的设计与高频放大器的设计除滤波器的通带范围不同以外其余的均相同，模型框图与高频放大器完全相同，回波经混频以后通过中频放大器后的时域图形频域以及滤波器的频率响应分别如下图所示：