testawgn.m

来自「这是一个wcdma下行功率控制仿真的源代码」· M 代码 · 共 35 行

%x = complex_awgn(0,0.36);

%t = 1:30;
%subplot(2,2,1);
%plot(t,10*log10(abs(x)),'--b');
%y1=sqrt(-log(u1)).*cos(2*pi*u2);
%y2=sqrt(-log(u1)).*sin(2*pi*u2);
clear;
clc;
N0 = 0.64;
m = 0;
randn('state',sum(100*clock));
u1 = m + sqrt(N0/2)* randn(1,2560);
randn('state',sum(100*clock));
u2 = m + sqrt(N0/2)* randn(1,2560);
ua = complex(u1,u2);
ub = complex(u1,u2);
uc = complex(u1,u2);
ud = complex(u1,u2);
u = ua+ub+uc+ud;
%u = u./abs(u);
n = mean(abs(ua).^2) 
n = mean(abs(u).^2) %均方根值
%m = mean(u1)
%t = mean(u2)
pdfc = normpdf(u2,m,sqrt(N0/2));
%subplot(1,2,1);
plot(u2,pdfc,'b.');
ylabel('f(xi)');
xlabel('xi');

%pdfc = normpdf(u2,m,sqrt(N0/2));
%subplot(1,2,2);
%plot(u2,pdfc,'b.');