exam4_wn.m.txt

ami 信号仿真实验源代码

echo off 
clear all 
close all 
N=1000; 
M=50; 
Rx_av=zeros(1,M+1); 
Sx_av=zeros(1,M+1); 
x=linspace(-25,25,M+1); 
for j=1:10, > 取10次实现到整体平均 
> 产生N 维均匀分布随机变量，在-1/2和1/2之间 
X=rand(1,N)-1/2; 
Rx=Rx_est(X,M); > 本次实现的自相关系数 
Sx=fftshift(abs(fft(Rx))); > 本次实现到功率谱 
Rx_av=Rx_av+Rx; > 子相关系数之和 
Sx_av=Sx_av+Sx; > 功率谱之和 
end; 
Rx_av=Rx_av/10; > 自相关系数的整体平均 
Sx_av=Sx_av/10; > 功率谱到整体平均 
> Plotting comments follow 
subplot(2,1,1) 
plot(x,fftshift(Rx_av)) 
grid on 
title('autocorrelation')>画自相关函数图 
subplot(2,1,2) 
plot(x,fftshift(Sx_av)) >画功率谱 
grid on 
title('spectrum')