c8_52.m

来自「经典书籍《通信系统仿真原理与无线应用》chap8部分课后题仿真源代码」· M 代码 · 共 39 行

clear
a=0.8;
m1=0;
m2=1;
sigma1=1;
sigma2=1;

subplot(2,2,1)
x = randn(1,100); 
n1 = (a/(sqrt(2*pi)*sigma1))*exp(((x-m1).^2)/(2*sigma1.^2));
n2 = ((1-a)/(sqrt(2*pi)*sigma2))*exp(((x-m2).^2)/(2*sigma2.^2));
y = n1 + n2;
hist(y,20)
ylabel('Ni'); xlabel('(a)')

subplot(2,2,2)
x = randn(1,100); 
n1 = (a/(sqrt(2*pi)*sigma1))*exp(((x-m1).^2)/(2*sigma1.^2));
n2 = ((1-a)/(sqrt(2*pi)*sigma2))*exp(((x-m2).^2)/(2*sigma2.^2));
y = n1 + n2;
hist(y,5)
ylabel('Ni'); xlabel('(b)')

subplot(2,2,3)
x = randn(1,1000); 
n1 = (a/(sqrt(2*pi)*sigma1))*exp(((x-m1).^2)/(2*sigma1.^2));
n2 = ((1-a)/(sqrt(2*pi)*sigma2))*exp(((x-m2).^2)/(2*sigma2.^2));
y = n1 + n2;
hist(y,50)
ylabel('Ni'); xlabel('(c)')

subplot(2,2,4)
x = randn(1,100000);
n1 = (a/(sqrt(2*pi)*sigma1))*exp(((x-m1).^2)/(2*sigma1.^2));
n2 = ((1-a)/(sqrt(2*pi)*sigma2))*exp(((x-m2).^2)/(2*sigma2.^2));
y = n1 + n2;
hist(y,50)
ylabel('Ni'); xlabel('(d)')