5.随机过程通过线性系统 - 随机信号分析实验报告(2)

5.随机过程通过线性系统 - 随机信号分析实验报告

%输入信号的频谱和功率谱 x_mag=abs(fft(x,2048)); f=(0:2047)*16000/2048; x_cm=abs(fft(x_arr,2048));

%画输入信号的时域,频谱图和频谱图 subplot 221 plot(t,x)

title('输入信号x') xlabel('t/s'); ylabel('x(t)'); grid on; subplot 222 plot(tau,x_arr)

title('输入信号x的自相关系数') xlabel('\tau/s')

ylabel('R_x_1(\tau)') subplot 223

plot(f(1:length(f)/2),x_mag(1:length(f)/2)) title('输入信号x的频谱') xlabel('f/Hz') ylabel('幅值') grid on; subplot 224

plot(f(1:length(f)/2),x_cm(1:length(f)/2)) title('输入信号x的功率谱') xlabel('f/Hz')

ylabel('S_x_1（f)')

%---------------高通滤波器的频谱和相频特性-----------------------------% [B,A]=butter(20,4500/(16000/2),'high'); figure(2)

freqz(B,A,2048);

title('高通滤波器的频率特性曲线') grid on

%---------------输入信号通过高通滤波器后的信号y1-----------------------% y1=filter(B,A,x);

%信号y1的均值，方差，均方值和自相关系数 y1_mean = mean(y1) y1_var=var(y1)

y1_st=y1_var+y1_mean^2 y1_arr=xcorr(y1);

%信号y1的频谱和功率谱 y1_mag=abs(fft(y1,2048)); y1_cm=abs(fft(y1_arr,2048));

%画出信号y1的时域图，频谱图和功率谱图 figure(3) subplot 221