首页 > 代码库 > 《DSP using MATLAB》示例Example5.9
《DSP using MATLAB》示例Example5.9
代码:
n = 0:10; x = 10*(0.8) .^ n; y = x(mod_1(-n,11)+1); %% ----------------------------------------------------------------- %% START a %% ----------------------------------------------------------------- figure(‘NumberTitle‘, ‘off‘, ‘Name‘, ‘Exameple5.9a x sequence‘) set(gcf,‘Color‘,‘white‘); subplot(2,1,1); stem(n,x); title(‘Original sequence x(n)‘); %axis([0,10,-2.5,2.5]); xlabel(‘n‘); ylabel(‘x(n)‘); grid on; subplot(2,1,2); stem(n,y); title(‘Circularly folded sequence y(n)‘); %axis([0,10,-2.5,2.5]); xlabel(‘n‘); ylabel(‘y(n)=x(-n mod 10)‘); grid on; %% ----------------------------------------------------------------- %% END a %% ----------------------------------------------------------------- %% ----------------------------------------------------------------- %% START b %% ----------------------------------------------------------------- X_DFT = dft(x,11); Y_DFT = dft(y,11); % DFT of x and y realX_DFT = real(X_DFT); imagX_DFT = imag(X_DFT); realY_DFT = real(Y_DFT); imagY_DFT = imag(Y_DFT); figure(‘NumberTitle‘, ‘off‘, ‘Name‘, ‘Exameple5.9b x sequence‘) set(gcf,‘Color‘,‘white‘); subplot(2,2,1); stem(n,realX_DFT); title(‘Real {DFT[x(n)]}‘); %axis([0,N,-0.5,1.5]); xlabel(‘k‘); grid on; subplot(2,2,2); stem(n,imagX_DFT); title(‘Imag {DFT[x(n)]}‘); %axis([0,N,-0.5,1.5]); xlabel(‘k‘); grid on; subplot(2,2,3); stem(n,realY_DFT); title(‘Real {DFT[x((-n))11]}‘); %axis([0,N,-0.5,1.5]); xlabel(‘k‘); grid on; subplot(2,2,4); stem(n,imagY_DFT); title(‘Imag {DFT[x((-n))11]}‘); %axis([0,N,-0.5,1.5]); xlabel(‘k‘); grid on; %% ----------------------------------------------------------------- %% END b %% -----------------------------------------------------------------
运行结果:
《DSP using MATLAB》示例Example5.9
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。