Question

matlab，随机信号问题

问题一：tau =(-length(x)+1:length(x)-1)/16000 什么意思，怎么理解。
问题二：plot(f(1:length(f)/2),N(1:length(f)/2)) 画出来的图为什么横坐标是到8000，f(1:length(f)/2)什么意思

%输入信号x的产生
t=0:1/16000:0.01;
x1=sin(1000*2*pi*t)+sin(2000*2*pi*t)+sin(3000*2*pi*t);
x=awgn(x1,5,'measured'); %加入高斯白噪声
n=x-x1; %高斯白噪声
(b)输入信号及其噪声的分析
%输入信号x自相关系数
x_arr=xcorr(x);
tau =(-length(x)+1:length(x)-1)/16000;
%输入信号x的频谱和功率谱
x_mag=abs(fft(x,2048));
f=(0:2047)*16000/2048;
x_cm=abs(fft(x_arr,2048));
%画出高斯白噪声n的时域图和频域图
figure(1)
subplot(1,2,1)
plot(t,n)
title('高斯白噪声n')
xlabel('t/s')
ylabel('n(t)')
grid on
subplot(1,2,2)
N=fft(n,2048);
plot(f(1:length(f)/2),N(1:length(f)/2))
title('高斯白噪声n的频谱图')
xlabel('f/Hz')
ylabel('幅值')
grid on

Answer 1

问题一：tau =(-length(x)+1:length(x)-1)/16000 什么意思，怎么理解。
length(x)是序列采样长度， tau的结果是两倍长度的采样时钟序列 －－这个tau在程序中没有用

问题二：plot(f(1:length(f)/2),N(1:length(f)/2)) 画出来的图为什么横坐标是到8000，f(1:length(f)/2)什么意思

因为傅里叶变换的结果对称的，故显示一半即可。因为f=(0:2047)*16000/2048，f( length(f)/2＋1) 正好等于采样频率的一半即 (1023+1)*16000/2048=8000