通信建模程序.docVIP

实验六 高斯随机变量直方图程序： subplot(2,2,1) %图形标号为2行2列中的第1个 x = randn(1,100); %产生一个长度为100的白噪声序列 hist(x,20) %画出序列x的20个直方图 ylabel(Ni); xlabel((a)) %横纵坐标名称 subplot(2,2,2) %图形标号为2行2列中的第一个 x = randn(1,100); %产生一个长度为100的白噪声序列 hist(x,5) %画出序列x的5个直方图 ylabel(Ni); xlabel((b)) %横纵坐标名称 subplot(2,2,3) %图形标号为2行2列中的第2个 x = randn(1,1000); %产生一个长度为1000的白噪声序列 hist(x,50) %画出序列x的50个直方图 ylabel(Ni); xlabel((c)) %横纵坐标名称 subplot(2,2,4) %图形标号为2行2列中的第3个 x = randn(1,100000); %产生一个长度为100000的白噪声序列 hist(x,50) %画出序列x的50个直方图 ylabel(Ni); xlabel((d)) %横纵坐标名称 实验图形： 白噪声通过切比雪夫滤波器： settle = 100; % 定义变量并初始化 fs = 1000; % 抽样频率 N = 50000; % 定义数据长度 f = (0:(N-1))*fs/N; % 频率范围 [b,a] = cheby1(5,5,0.1); % 通过切比雪夫滤波器 NN = N+settle; % 变化范围 in = randn(1,NN); % 产生白噪声 out = filter(b,a,in); % 滤波器输出 out = out((settle+1):NN); % 截取初始采样值 window = hanning(N); % 加汉宁窗 winout = out.*window; % 窗滤波器输出 fout = abs(fft(winout,N)).^2; % 功率谱 U = sum(window.*window); % 窗能量 f1out = fout/U; % 功率谱密度 psd1 = 10*log10(abs(f1out)); % 功率谱密度的log值 subplot(2,1,1) %图形标号为2行1列中的第1个 plot(f(1:5000),psd1(1:5000)) %画图 grid; axis([0 100 -70 10]); %加网格及图形范围 xlabel(Frequency, Hz) ylabel(PSD) %坐标名称 K = 25; % 分段数 M = N/K; % 分段长度 fK = (0:(M-1))*fs/M; % 频率范围 d = zeros(1,M); % 初始化向量 psdk = zeros(1,M); %初始化向量 window = hanning(M); % 加汉宁窗 U = sum(window.*window); % 窗能量 for k=1:K for j=1:M index = (k-1)*M+j; d(j) = out(index); end