1. 对于直接给出频响函数的情况 这里以滑动平均的频响函数作为例子,滑动窗口为[0, 4]. 上式中M2=4. >> w=0:0.001:2*pi; >> h1=1-exp(-j*w*5); >> h2=1-exp(-j*w); >> h=0.2*h1./h2;                % 注意这里是点除 >> plot(w/pi,abs(h));            % 横轴以pi为单位 >> plot(w/pi,angle

对于多窗口谱估计的理解 目录 对于多窗口谱估计的理解 0. 缘起 1. PMTM 含义 2. 与我们常用的周期谱估计的区别 3. 计算过程 5. 多窗/单窗谱估计结果对比 6. 程序如何生成多窗 - 以sin为例 6.1 生成正交窗的公式代码 6.2 计算是否符合正交 7. Reference 0. 缘起 在语音增强通过改进的谱减法进行谱估计时,使用多窗谱估计法要优于我们一般遇到的周期法谱估计,中文相关的资料不多,这里分享一下我学习理解到的内容. 1. PMTM 含义 P: PSD (Power

%相位角与相关系数曲线 close all; clear all; Samp1=200;  %设置信号的采样精度 Samp2=200;  %设置相位角p分割精度 A=10;%信号幅值 w=1;%信号角频率 na=5;%噪声放大系数 g=zeros(Samp2,Samp1);%产生Samp2 x Samp1的二维零矩阵,每一行对应一个p值时的函数 t = linspace(0,6*pi/w,Samp1); %创建函数向量 p= linspace(0,pi/3,Samp2);%创建相位角向量 f=A

(此帖引至网络资源,仅供参考学习)第一:频谱 一.调用方法 X=FFT(x):X=FFT(x,N):x=IFFT(X);x=IFFT(X,N) 用MATLAB进行谱分析时注意: (1)函数FFT返回值的数据结构具有对称性. 例:N=8;n=0:N-1;xn=[4 3 2 6 7 8 9 0];Xk=fft(xn) →Xk =39.0000            -10.7782 + 6.2929i         0 - 5.0000i    4.7782 - 7.7071i    5.0000

用到了三角窗脉冲序列,各小题的DTFT就不写公式了,直接画图(这里只贴长度M=10的情况). 1. 代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <

一个复数序列可以分解为共轭偶对称和共轭奇对称部分. 代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Probl

前几天看了看博客,从16年底到现在,3年了,终于看书到第9章了.都怪自己愚钝不堪,唯有吃苦努力,一点一点一页一页慢慢啃了. 代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n');

下雨了,大风降温,一地树叶,终于进入冬季了 代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Problem 8

代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Problem 8.45.4 \n\n'); banner()

代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Problem 8.43 \n\n'); banner();

代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Problem 8.42 \n\n'); banner();

代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Problem 8.41.2 \n\n'); banner()

代码: function [wpLP, wsLP, alpha] = bs2lpfre(wpbs, wsbs) % Band-edge frequency conversion from bandstop to lowpass digital filter % ------------------------------------------------------------------------- % [wpLP, wsLP, alpha] = bs2lpfre(wpbs, wsbs)

代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Problem 8.39.2 \n\n'); banner()

代码: function [wpLP, wsLP, alpha] = bp2lpfre(wpbp, wsbp) % Band-edge frequency conversion from bandpass to lowpass digital filter % ------------------------------------------------------------------------- % [wpLP, wsLP, alpha] = bp2lpfre(wpbp, wsbp)

代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Problem 8.37 \n\n'); banner();

上代码: function [wpLP, wsLP, alpha] = lp2lpfre(wplp, wslp) % Band-edge frequency conversion from lowpass to lowpass digital filter % ------------------------------------------------------------------------- % [wpLP, wsLP, alpha] = lp2lpfre(wplp, wslp)

代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Problem 8.35 \n\n'); banner();

今天下了小雨,空气中泛起潮湿的味道,阴冷的感觉袭来,心情受到小小影响. 代码: hp2lpfre子函数 function [wpLP, wsLP, alpha] = hp2lpfre(wphp, wshp) % Band-edge frequency conversion from highpass to lowpass digital filter % -----------------------------------------------------------------------

代码: %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Problem 8.33.2 \n\n'); banner()