mvdr01.m

来自「窄带信号mvdr波束形成,这几天的学习结果」· M 代码 · 共 17 行

close all;
clear all;
clc
N=1024;
dt=0.0001;
df=1/(N*dt);
f1=2000;
t=dt*(-N/2:N/2-1);
f=df*(-N/2:N/2-1);
snr=0;
am=sqrt(2)*10^(snr/20);
s=am*exp(j*f1*2*pi*t);
noise=randn(1,N)+j*randn(1,N);
S=s+noise;
x=fftshift(fft(fftshift(S)));
% x=fft(S);
figure; plot(f,abs(x));