matlab freqz half,[转载]Matlab freqz函数使用

freqz函数计算线性系统的频率响应,包括幅频响应和相频响应,基本输入为线性系统的AMMA模型系数向量,一个典型的AMMA模型为

% jw -jw -jmw

% jw B(e) b(1) +

b(2)e + .... + b(m+1)e

% H(e) = ---- = ------------------------------------

% jw -jw -jnw

% A(e) a(1) +

a(2)e + .... + a(n+1)e

其中,向量A为IIR系统分母向量,对FIR系统为1;向量B为分子向量

函数使用零散傅立叶变换(fft)计算系统响应,计算范围为0~pi,即单位圆的上半部分,因为默认情况下对于实数系数系统响应关于DC对称.计算公式为上式的频域离散化形式.计算完成后在对目标解析度使用算术插值方法进行向量填充,填充形式取决于附加输入参数信息.

freqz的基本调用格式为freqz(b,a),此时函数内部默认进行512点算术插值,所以输出是512点的幅度值和对应的0~pi自变量向量,如果需要改变向量长度便于绘图,则可以附加维度信息,例如freqz(b,a,N)则进行N点插值

如果系统不是实系数的而是复数系统,则系统响应不再对称,0~pi的范围不足以反映整体特性,此时有两种选择:

1. 绘制整个主值区(-pi~pi),调用格式是增加字串型参数'whole',例如

freqz(b,a,n,'whole');

2. 绘制定制区域,则直接给定范围,该范围为给定自变量向量,例如 freqz(b,a,[-pi:2*pi/n:pi])

使用上述调用格式,输出为-pi~pi范围内的离散傅立叶变换计算结果,如果为绘图方便要标识频率信息,则调用时需要附加采样频率信息,例如freqz(b,a,[-fs/2:(fs/n):fs/2],fs);

freqz的基本使用方法是利用函数返回值进行后续分析处理,此时函数不会自动绘制曲线,如果直接调用而不指定返回接收变量,则函数自动按照通常标准绘制幅频/相频曲线,幅频曲线为对数格式

可以使用type freqz分析源代码

matlab freqz half,[转载]Matlab freqz函数使用相关推荐

lyapunov函数 matlab,科学网-[转载]Matlab的Lyapunov、Sylvester和Riccati方程的Matlab求解-吴雄君的博文...
一.连续Lyapunov方程连续Lyapunov方程可以表示为 Lyapunov方程来源与微分方程稳定性理论,其中要求C为对称正定的n×n方阵,从而可以证明解X亦为n×n对称矩阵,这类方程直接求解比较 ...
MATLAB阶段性方程组,[转载]matlab 解方程组
1.解方程最近有多人问如何用matlab解方程组的问题,其实在matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法: (1)x=i ...
matlab水汽计算公式,[转载]matlab 解方程组
1.解方程最近有多人问如何用matlab解方程组的问题,其实在matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法: (1)x=i ...
matlab灰度归一化,[转载]matlab图像处理为什么要归一化和如何归一化？
matlab图像处理为什么要归一化和如何归一化?一.为什么归一化 1. 基本上归一化思想是利用图像的不变矩寻找一组参数使其能够消除其他变换函数对图像变换的影响.也就是转换成唯一的标准形式以抵抗仿射变换 ...
matlab 比较日期,[转载][Matlab]关于时间的函数的不完全总结
一.日期和时间三种表示格式 (1)日期字符串: NowDate=date NowDate = 08-Aug-2008 year(NowDate) %提取日期的年份 ans = 2008 (2)连续的 ...
matlab中select,[转载]MATLAB阈值获取函数ddencmp、thselect、wbmpen和w
crit(t) wdcbm的调用格式有以下两种: (1)[THR,NKEEP]=wdcbm(C,L,ALPHA); (2)[THR,NKEEP]=wdcbm(C,L,ALPHA,M); 函数wdcbm ...
matlab定步长ode,[转载]matlab ode45 函数传自定义参数用法及定步长ode
要用的时候总是忘记,这回给把它写在这里! %%程序1 arg1 = 2; arg2 = 1; [T,Y] = ode45('vdp1000',[0 10],[2 0], [], arg1, arg2) ...
matlab msgbox 换行,[转载]Matlab/GUI笔记
转自 http://www.kylen314.com/archives/412 不显示坐标刻度: set(gca,'xtick',[]) =============================== ...
matlab画图模糊,[转载]matlab中模糊工具箱的使用
用 Matlab 中的 Fuzzy 工具箱做一个简单的模糊控制,流程如下: 1.创建一个 FIS (Fuzzy Inference System ) 对象, a = newfis(fisName,fi ...
matlab xcorr lags,[转载]matlab中xcorr的用法
Matlab中用于计算自相关函数的指令是xcorr.比如矩阵A=[1 2 3]; xcorr(A)=3.0000 8.0000 14.0000 8.0000 3.0000 自相关函数是信号间隔的函数, ...

matlab freqz half,[转载]Matlab freqz函数使用

matlab freqz half,[转载]Matlab freqz函数使用相关推荐

最新文章

热门文章