矩阵运算处理,包括FFT,IFFT,FFTSHIFT等

计算全息close all clc clear A=zeros(64) A(15:20,20:40)=1 A(15:50,20:25)=1 A(45:50,20:40)=1 A(30:34,20:35)=1 ％ ppp=exp(rand(64)*pi*2*i) A=A.*ppp ％ Author s email: zjliu2001@163.com figure imshow(abs(A),[]) Fa=fft2(fftshift(A)) Fs=fftshi

：滤波 fs=400 采样频率,必须大于两倍基带信号最高频率 ts=1/fs 采样时间间隔 T=5 时间窗大小 t=0:ts:T -hf1=zeros(1,N) hf1(f<=B1&f>=-B1)=1 y1=hf1.*fftshift(y) figure subplot(211) plot(f,abs(fftshift(y))/N,f,hf1, r-- ) axis([-fs/2,fs/2,0,1.2]) title( 原信号

例3：分析占空比为ta的矩形脉冲频谱 T=1 时间窗大小 ta=0.5 占空比 N=256 t=linspace(-T,T,N) ts=2*T/N fs=1/ts -x=zeros(1,N) x(N/2-round(N/2*ta):N/2+round(N/2*ta))=ones(1,2*round(N/2*ta)+1) y=fft(x,N) figure plot(t,x) title([ 时域信号图 占空比为 ,num