MATLAB随机信号分析与处理

2021在学中————
课本：

Q：用频率的方法估计投掷均匀硬币出现正面的概率。

clear all;
clf;
clc;
N=[10,50,100,1000,10000,100000];
for i=1:length(N)
    P(i)=sum(round(rand(N(i),1)))/N(i);
end
plot(P,'-*');

例1.11

正态概率密度函数normpdf()
正态概率分布函数normcdf()

clc;clear all
%画出服从均值为1，方差为4的正态概率密度和概率分布曲线
X=-6:0.01:7;
Y=normpdf(X,1,2);
subplot(1,2,1);
axis on;
plot(X,Y);
axis square;
title('正态概率密度函数');
Y=normcdf(X,1,2)
subplot(1,2,2);
plot(X,Y);
title('正态概率分布函数');
axis square;

误差函数erf(x) 补余误差函数erfc(x)

例1.12

%画参数为1的瑞利分布概率密度图
x=0:0.1:3;
p=raylpdf(x,1);
plot(x,p);

例1.13

syms x y A;                     %？？？
f=A*exp(-x-y);
C=int(int(f,x,0,inf),y,0,inf);  %计算分布函数的最大值，最大值C=1，由此可得C=A=1
P1=int(int(f,y,0,1-x),x,0,1);   %计算P{0<X+Y<1}
fx=int(f,y,0,inf);
fy=int(f,x,0,inf);
syms u v;
f=c*exp(-u-v);
F=int(int(f,u,0,y),v,0,x);      %求分布函数的表达式

lamda=11;
M=32768;
x(1)=19;
for n=1:200
    x(n+1)=(mod(lamda*x(n)+11117,M));
end
plot(x/M);
xlabel('n');
ylabel('x(n)');
axis([0 200 0 1])

X(t)称为半二元传输信号

N=200;
ind=find(rand(N,1)>0.5);
z(1:N)=1;
z(ind)=-1;
stairs(1:25,z(1:25));
axis([0 25 -1.5 1.5]);
xlabel('t/s(假定T=1s)');
ylabel('X(t)','FontSize',[12]);

2.6典型的随机过程：白噪声

x=randn(500,1);%产生平稳白噪声
plot(x);
xlabel('n');
ylabel('x(n)');

平稳白噪声500个样本点的波形

例2.29 产生一组随机序列，并画出他们的直方图

x=-2.9:0.1:2.9;
y=normrnd(0,1,1000,1);
hist(y,x);