clear;
clc
L=100;
LL=1/L;
P=zeros(1,L);
i=360;
RAve=20; %选择功率判定范围
for l=1:L
RR(l)=LL*l*RAve;
RRl=LL*RAve; %选择判定步长，分为100步
end
for k=1:100
A2=k/100; %A2为接收到的信号幅度值，下A1同义
F=rand(1,i)*2*pi;
C=cos(F);
S=sin(F); %得到360个在（0，360度）上均匀分布的角度，并计算正弦
R=((1+A2*C).*(1+A2*C)+(A2*S).*(A2*S)); %计算接收功率，取A1为1(R=(A1)^2.*((1+A2.*C).*(1+A2.*C)+(A2.*S).*(A2.*S)))
for j=1:i
for l=1:L
if ((R(j)>(RR(l)-RRl)))&(R(j)<RR(l))
P(l)=P(l)+1; %判定并统计落在各个统计区间上的次数
end
end
end
end
M=0.3; %通过控制M的取值来控制仿真图形跨度
for m=1:M*L
PP(m)=P(m);
end
l=1:M*L;
bar(l,PP)

[此贴子已经被作者于2007-5-31 9:03:08编辑过]