dsp上机参考.docVIP

1. DFT定义，矩阵形式运算，循环运算 1 2. fir1 , kaiserord , Kaiser , freqz 2 3. cheb1ord , cheb1 2 4. 系统函数的二阶分解，tf2sos,tf2zp,sort,unique 3 5. 零极点形式到直接型的转换 3 6. 阶数估计（表10.2）,hann,fir1，根据定义加窗 4 7. freqz , real , imag , abs , angle , unwrap , phase 4 8. fir1 , hamming窗 5 9. 直接型转零极点型，tf2zp 6 10. butter,buttord 6 11. 部分分式展开（逆），从传输函数计算冲击响应，residuez,impz 7 12. freqz , real , imag , abs , phase , angle , unwrap 7 13. ellipord , ellip , freqz 8 14. fir1 , hann , 分贝转换 9 15. 信号能量计算 9 16. 各IIR阶数估计，buttord , cheb1ord , cheb2ord , ellipord 9 17. square , freqz 10 18. 线性卷积，此处用fft实现 11 19. 圆周卷积，要求用基于DFT的方法实现 11 20. 部分分式展开，逆z变换，residuez 12 21. freqz 13 22. 信号通过系统的输出，filter 13 23. 信号通过系统的输出，filter 13 24. 重要提示 14 DFT定义，矩阵形式运算，循环运算 % Writing a MATLAB program to compute 128-point DFT of % the following sequence, you must firstly use % DFT definition (directly computing DFT) to compute % and use MATLAB function to test the result. % Plot the two results in one figure. % g[n] = sin(pi * n / 4) , 0 = n = 31 N = 32 ; M = 128 ; n = 0:(N - 1) ; m = 0:(M - 1) ; g = sin(pi * n / 4) ; %DFT definition %matrix expression of DFT G1 = g * exp(-2 * j * pi / M).^(n * m) ; %following is the direct conversion %from definition to code % for ii = 0:M-1 % G1(ii + 1) = 0 ; % for jj = 0:N-1 % G1(ii + 1) = G1(ii + 1) + g(jj + 1) * exp(-2 * j * pi / M * ii * jj) ; % end % end subplot(2 , 1 , 1) ; plot(m , abs(G1)) ; %fft realized by MATLAB G2 = fft(g , M) ; subplot(2 , 1 , 2) ; plot(m , abs(G2)) ; fir1 , kaiserord , Kaiser , freqz % Using the function fir1 and window of Kaiser, % design a linear-phase FIR lowpass filter % meeting the following specifications: % passband edge frequency=2kHz, % stopband edge frequency=2.5kHz, % passband ripple δp=0.005, % stopband rippleδs=0.005, % and sampling rate of 10kHz. % Plot its gain and phase responses % and check if it meets the specifications? fsamp = 10000; fcuts = [2000 2500]; mags = [1 0]; devs = [0.005 0.005]; [n,Wn,beta,ftype] = kaiserord(fcu