信号实验报告二.docxVIP

四?实验结果与分析 (-)实验所需验证性程序 NT N-1 _母賊 X伙)=工班沙呼=》＞(讥“ % (1)离散傅立叶变换函数 心0 心0 % dft.m function [Xk]=dft(xn,N) n=[0:l:N-l]; k=[0: N-1J; WN=exp(-j*2*pi/N); nk=n*k; WNnk=WN.A(nk); Xk=xn*WNnk; % 采用矩阵相乘的方法 i N-\ % (2)离散傅立叶反变换函数 “ “0 function [xn]=idft(Xk,N) n=[0:l:N-l]; k=[0:l:N-l]; WN=exp(-j*2*pi/N); nk=n*k; WNnk=WN.A(-nk); xn=(Xk*WNnk)/N; % (3)实序列分解函数 %实序列可分解为共扼对称分量 xec = (l/2)*[x(n) + x((N-n))N] % 和共扼反对称分量 xoc = (1/2)* [x(n) - X((N -n))N] function [xec,xoc]=circevod(x) N=length(x); n=0:(N-l); xec=().5*(x+x(mod(-n,N)+l)); % mod取余，例如 mod(4,5)=4, mod(-4,5)= mod(5-4,5)=l xoc=0.5*(x-x(mod(-n,N)+1)); % (4)序列的循环移位函数y(n)=双⑺-加))N function y=cirshftt(x,m,N) if length(x)N error(N mustbe = the length of x) %要求移位周期大于信号长度 end x=fx zeros( 1 ,N-length(x))]; n=[0:l:N-l]; n=mod(n-m?N); y=x(n+l); % (5)计算两序列的循环卷积函数 %循环卷积的长度N应大于max(length(x 1 ),length(x2)) function y=circonv(xl,x2,N) if length(xl)N error(*N must be = the length of x T) end if length(x2)N error(N must be = the length of x2) end xl=[xl zeros( 1 ,N-length(x 1))]; %将序列补0成为N长序列 x2=[x2 zeros( 1 ,N-length(x2))]; m=[0:l:N-l]; x3二x2(mod(-m,N)+l); %该语句的功能相当于序列翻褶，延拓，取主值序列 H=zeros(N,N); for n=l:l:N %得到只3序列的循环移位矩阵 H(n,:)=cirshftt(x3,n-1 ,N); end y=xl*H* %用矩阵相乘的方法得到结果 %(6)序列的傅里叶变换及离散傅里叶变换计算 N=5; n=0:4; x=[ones(l,5)]; %产生矩形序列 k=():999;w=(pi/5()())*k; X=x*(exp(-j*pi/500)).A(n*k); %计算序列的傅立叶变换 Xe=abs(X); subplot(3,2,1 );stem(n,x);ylabel(x(n)‘)； subplot(3,2,2);plot(w/pi,Xe);ylabel(*IX(ejw)r); %画出序列的傅立叶变换 X=dft(x,N); %进行5 点 DFT magX=abs(X); k=(O:length(magX)1-1 )*N/length(magX); subplot(3,2,3);stem(n,x);ylabel(x(n)); subplot(3,2,4);stem(k,magX); axis(lO,5,O,5]);ylabelCIX(k)r); % (7)序列产生及其DFT N=20; n=O:N-l; x=2*cos(0.25*pi*n)+cos(0.65*pi*n); subplot(2,1,1 );stem(n,x);title(N=20 时的信号)； Y=dft(x,N); kl=0:N-l;wl=2*pi/N*kl; subplot(2,l ,2);stem(w l/pi,abs(Y));title(信号的频谱)； (二)实验程序代码及图形结果分析 1、序列的离散傅里叶变换实现 分别对复序列，实序列，实偶序列，实奇序列，虚奇序列进行离散傅里叶变换，得到实 验结果并对艾特点进行分析。实验所需序列自选。 驚复序列傅里叶变换 xl= [li,lj z 2i,2j z 3i] n=5 XKl=dft(xl.n) % XK1=0.0000 + 9?0000i 9021 - l.OOOOi