数字信号处理-第二次上机Word文档格式.docx

1、N1=16;X1=dft（xn,N1）;figure（1）k=0:N1-1;subplot（2,1,1）stem（k,abs（X1）,.）;xlabel（k）,ylabel（|X1（k）|title（16点DFT）,grid onsubplot（2,1,2）stem（k,angle（X1）,grid onylabel（angle（X1（k）%16点FFTX1=fft（xn,N1）;figure（2）16点FFT%32点DFTN2=32;X2=dft（xn,N2）;figure（3）N2-1;stem（k,abs（X2）,|X（2k）|32点DFTstem（k,angle（X2）,angle（X

2、2（k）%32点FFTX2=fft（xn,N2）;figure（4）32点FFT 运行结果：2. 输出序列x（n）=sin（0.5n+0.2）,0nN-1 ,N自定。计算并输出x（n）的N点DFT：X1k，以及x（n）的2N点DFT：X2k 。观察X1k 和X2k ，能得出什么结论？N=32; %第一次取采样点32个，第二次取采样点2*32=64个N-1;xn=sin（0.5*pi*n+0.2*pi）;Xk1=fft（xn,N）Xk2=fft（xn,2*N）stem（k,abs（Xk1）,|Xk1（k）|）,grid on,title（）2*N-1;stem（k,abs（Xk2）,|Xk2（k

3、）|64点DFT运行结果：3. 用快速卷积算法计算下列两序列的线性卷积序列，并输出结果 x1n=0,2,2,1x2n=1.02n,0n150.98n,16n28输出x1（n）、x2（n）及其FFT信号图形，输出卷积结果。（注意FFT点数应满足循环卷积与线性卷积相等条件）与第一次上机作业中时域线性卷积比较（指计算时间比较，卷积结果应相同）。x1=0,2,2,1; %序列1a1=（1.02）.n;n=16:28;a2=（0.98）.n;x2=a1 a2; %序列2N1=length（x1）;X1=fft（x1,N1）;x1的DFTN2=length（x2）;X2=fft（x2,N2）;|X2（k）

4、|x2的DFTL=N1+N2-1;y1=ifft（fft（x1,L）.*fft（x2,L） %卷积subplot（3,1,1）,stem（x1,nx1subplot（3,1,2）,stem（x2,x2subplot（3,1,3）,stem（y1,yy=x1*x2y2=juanji（x1,x2）y2y=x1*x2直接卷积结果4.若x（n）=x1（n）+jx2（n）,其中x1（n）=cos（n/4）,x2（n）=sin（n/8）。根据DFT的对称性，由X（k）求出X1（k）=DFTx1（n）和X2（k）=DFTx2（n）。（X（k）=DFTx（n）clear N=4;x1=cos（pi*n/4）;x2=sin（pi*n/8）;xn=x1+x2*i;Xk=fft（xn,N）;Xk（N+1）=Xk（1）;k=1:N;Xepk=0.5*（Xk（k）+conj（Xk（N+2-k）;Xopk=0.5*（Xk（k）-conj（Xk（N+2-k）;X1k=Xepk %DFT（x1）X2k=Xopk/j %DFT（x2）