随即过程作业文档格式.docx
《随即过程作业文档格式.docx》由会员分享，可在线阅读，更多相关《随即过程作业文档格式.docx（18页珍藏版）》请在冰豆网上搜索。
图1随机过程第1题图
程序源代码：
%本段程序是用来产生在[0,1]上均匀分布的白色序列，并对其进行分布、...
...均值、方差和相关检验。
clearall%清空内存空间
X=rand（1,2000）;
%产生在[0,1]均匀分布的白色序列
X50=X（1:
50）%提取序列的前50个数
EX=mean（X）%求取均值
DX=var（X）%求取方差
subplot（2,1,1）;
hist（X,10）%画出分布检验的直方图
title（'
均匀分布的白色序列的分布检验图'
）;
y=zeros（1,21）;
%y=linspace（-10,10,21）
%以下为求取白色序列的相关函数过程
form=-10:
10
absm=abs（m）;
%取数值m的绝对值
sum=0;
forn=1:
2000-absm
sum=sum+（X（n+absm）-EX）*（X（n）-EX）;
end
y（m+11）=sum/2000;
end
m=（-10:
10）;
subplot（2,1,2）;
plot（m,y）%画出白色序列的相关检验图
均匀分布的白色序列的相关检验图'
）
第二题
用PC机产生N（0,1）分布的正态序列{
（1）打印出前50个数；
（2）分布检验；
（3）均值检验；
（4）方差检验；
（5）计算出相关函数
，i=0,±
（1）正态序列的前50个数：
-0.4326-1.66560.12530.2877-1.1465
1.19091.1892-0.03760.32730.1746
-0.18670.7258-0.58832.1832-0.1364
0.11391.06680.0593-0.0956-0.8323
0.2944-1.33620.71431.6236-0.6918
0.85801.2540-1.5937-1.44100.5711
-0.39990.69000.81560.71191.2902
0.66861.1908-1.2025-0.0198-0.1567
-1.60410.2573-1.05651.4151-0.8051
0.52870.2193-0.9219-2.1707-0.0592
（2）在[0,1]间取五个小段检验分布的实际值与理论值如下表所示，分布检验如图2所示
抽样区间
0~0.2
0.2~0.4
0.4~0.6
0.6~0.8
0.8~1.0
实际个数
156
179
131
134
113
理论个数
158.5194
152.3241
140.6503
124.7954
106.4003
EX=0.0012
DX=0.9783
（5）相关检验如图2所示
图2随机过程第2题图
%本段程序是用来产生在N（0,1）分布的正态序列，并对其进行分布、...
X=randn（1,2000）;
%产生在N（0,1）正态分布序列
hist（X,20）%画出分布检验的直方图
N（0,1）分布正态序列的分布检验图'
%以下为求取N（0,1）正态分布序列的相关函数过程
N（0,1）分布正态序列的相关检验图'
%下面这段程序为计算在区间[0,0.2]、[0.2,0.4]、[0.4,0.6]、[0.6,0.8]、...
...[0.8,1.0]五个区间上的实际分布个数
sum=ones（1,5）;
fori=1:
2000
forj=1:
5
if（X（i）>
=（j-1）/5&
X（i）<
j/5）
sum（j）=sum（j）+1;
shiji=sum
%以下为统计上述五个区间的理论分布个数
fork=1:
sum（k）=quad（'
exp（（-x.^2）/2）/sqrt（2*pi）'
（k-1）/5,k/5）;
lilun=2000*sum
第三题
设
为正态白噪声序列N（0,1），
其中
（1）
（2）
（3）
（4）
（1）EX=0.0061
=16.8989
（3）DX=16.8989
（4）中心自相关图如图3所示
图3随机过程第3题图
clearall%清空内存空间
x=randn（1,2001）;
%N（0,1）正态白噪声序列
X=ones（1,2000）;
fori=1:
X（i）=x（i+1）+4*x（i）;
%
EX=mean（X）
DX=var（X）
E2=DX+EX^2%X^2的均值
form=-10:
plot（m,y）
第四题
（1）EX=3.5126e-004
=1.9082
（3）DX=1.9082
（4）中心自相关图如图4所示
图4随机过程第4题图
x=randn（1,2000）;
X
（1）=x
（1）;
1999
X（i+1）=x（i+1）-0.707*X（i）;
%关系式
Y=X（101:
1000）;
EX=mean（Y）
DX=var（Y）
forn=101:
y（m+11）=sum/1900;
第五题
已知
，
，采样周期
若令
求
其中N取5、10、20
画出
和
并比较。
图5随机过程第5题图
clearall
Y=zeros（3,200）;
%定义一个3行200列的矩阵
200
t=i*pi/100;
fork=1:
3
forn=-5*2^（k-1）:
5*2^（k-1）
if（t-n*pi/2）==0%分式分母为零的情况
Y（k,i）=Y（k,i）+sin（n*pi/2）;
else
Y（k,i）=Y（k,i）+sin（n*pi/2）*sin（t-n*pi/2）/（t-n*pi/2）;
end
t=（1:
200）*pi/100;
f=sin（t）;
%X（t）
x=plot（t,f,'
r-'
t,Y（1,:
）,'
k*'
t,Y（2,:
g-.'
t,Y（3,:
b--'
legend（x,'
原函数'
'
N=5采样曲线'
N=10采样曲线'
N=20采样曲线'
第六题
系统方框图：
（1）列出奥斯特姆表；
（2）判断系统稳定性；
（3）如果稳定，求
？
，其中
。
（1）随机过程第6题图1是Astrom表左侧A（s）部分，图2是Astrom表右侧B（s）部分。
图6随机过程第6题图1
图7随机过程第6题图2
（2）由Astrom表可知，系统不稳定。
clc
formatshort
A=zeros（40,21）;
B=zeros（40,20）;
Ak=zeros（1,20）;
Bk=zeros（1,20）;
21
A（1,i）=2*i-1;
B（1,20）=1;
forj=1:
20%求矩阵A的偶数行
ifmod（j,2）==1;
%奇数位（非矩阵的奇数位）
A（2,j）=A（1,j+1）;
else%偶数位
A（2,j）=0;
end;
end;
B（2,:
）=A（2,1:
20）;
Ak
（1）=A（1,1）/A（2,1）;
Bk
（1）=B（1,1）/B（2,1）;
fori=3:
40
ifmod（i,2）==1;
%矩阵奇数行，下面这小段程序是用来求矩阵奇数行的。
forj=（i+1）/2:
2:
20;
A（i,j）=A（i-2,j）;
B（i,j）=B（i-2,j）;
ifj~=20;
ifj==19
A（i,j+1）=A（i-2,j+1）-A（i-1,j+1）*Ak（（i-1）/2）;
B（i,j+1）=B（i-2,j+1）-Bk（（i-1）/2）*B（i-1,j+1）;
A（i,j+2）=A（1,j+2）;
A（i,j+1）=A（i-2,j+1）-A（i-1,j+1）*Ak（（i-1）/2）;
else
A（i,j+1）=A（1,j+1）;
end
else%下面这小段程序是用来求矩阵偶数行的。
forj=i/2:
20
A（i,j）=A（i-1,j+1）;
B（i,:
）=A（i,1:
Ak（i/2）=A（i-1,i/2）/A（i,i/2）;
Bk（i/2）=B（i-1,i/2）/B（i,i/2）;
sum=0;
sum=sum+Bk（j）^2/Ak（j）;
disp（A）
disp（Ak）
disp（B）
disp（Bk）
I=0.5*sum%系统输出过程方差
第七题
（1）Astrom表如下图所示，上部分是A（s）下部分是B（k）
（2）由上图可知，系统稳定。
=2.06560223580002
2.067
formatlong
A=zeros（14,7）;
B=zeros（14,7）;
Ak=zeros（1,7）;
Bk=zeros（1,7）;
A（1,:
）=[10.50.620.010.01210.000050.00006];
B（1,:
）=[000001-0.55];
7
A（2,i）=A（1,（7-i+1））;
）=A（2,:
Ak
（1）=A（1,7）/A（2,7）;
Bk
（2）=B（1,7）/B（2,7）;
14
（7-（i-1）/2）
A（i,j）=A（i-2,j）-Ak（（i-1）/2）*A（i-2,（（7-（i-3）/2）-j+1））;
B（i,j）=B（i-2,j）-Bk（（i-1）/2）*A（i-2,（（7-（i-3）/2）-j+1））;
else%矩阵偶数行，下面这小段程序是用来求矩阵偶数行的。
（7-（i-2）/2）
A（i,j）=A（i-1,7-（i-2）/2-j+1）;
B（i,j）=A（i,j）;
Ak（i/2）=A（i-1,（7-i/2+1））/A（i,（7-i/2+1））;
Bk（i/2）=B（i-1,（7-i/2+1））/B（i,（7-i/2+1））;
7;
sum=sum+B（（2*j-1）,（7-j+1））^2/B（2*j,（7-j+1））;
I=A（1,1）*sum%系统输出过程方差
升级会员 