|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
以下是我的程序:
x1=input('input x1(n):');
x2=input('input x2(n):');
N1=length(x1);n1=0:N1-1;subplot(421);stem(n1,x1);title('x1(n)');grid on;
N2=length(x2);n2=0:N2-1;subplot(422);stem(n2,x2);title('x2(n)');grid on;
% 应用FFT求相关
% 补零使两序列长度为N=N1+N2-1
x3=[x1,zeros(1,N2-1)];
x4=[x2,zeros(1,N1-1)];
N=length(x3);
% 分别求补零后的2N点FFT
X3=fft(x3,2*N);
X4=fft(x4,2*N);
% 应用相关定理
r12=ifftshift(ifft(conj(X3).*X4));
m=-N:N-1;
subplot(412);stem(m,r12/N);title('FFT');grid on;
% 应用xcorr函数进行验证
[r121,lags]=xcorr(x1,x2,'biased');
subplot(413);stem(lags,r121);title('xcorr');grid on;
问题1:应用FFT求信号相关和xcorr函数求的不同,不知道错误出在哪儿了,请大家帮忙!
问题2:xcorr函数求相关要求两序列长度相同,否则,就不能求有偏估计,是这样么?如果是,应该怎么求啊? |
|