matlab中的freqs和freqz函数

1.freqs

模拟滤波器的频率响应

语法:h = freqs(b,a,w)

[h,w] = freqs(b,a)

[h,w] = freqs(b,a,f)

freqs(b,a)

1.1描述:

freqs 返回一个模拟滤波器的H(jw)的复频域响应(拉普拉斯格式)

h = freqs(b, a, w) 根据系数向量计算返回模拟滤波器的复频域响应

freqs 计算在复平面虚轴上的频率响应h,角频率w确定了输入的实向量,因此必须包含至少一个频率点。

[h, w] = freqs(b, a) 自动挑选200个频率点来计算频率响应h

[h, w] = freqs(b, a, f) 挑选f个频率点来计算频率响应h

1.2例子:

找到并画出下面传递函数的频率响应

Matlab代码:a = [1 0.4 1];

b = [0.2 0.3 1];

w = logspace(-1, 1);

freqs(b, a, w);

logspace 功能:生成从10的a次方到10的b次方之间按对数等分的n个元素的行向量

n如果省略,则默认值为50。h=freqs(b,a,w);

mag = abs(h);phase = angle(h);

subplot(2,1,1), loglog(w,mag);

subplot(2,1,2), semilogx(w,phase);

f = w/(2*pi);mag = 20*log10(mag);phase = phase*180/pi;

2.freqz

MATLAB提供了专门用于求离散系统频响特性的函数freqz()

调用freqz()的格式有以下两种:

2.1[H,w]=freqz(B,A,N)

B和A分别为离散系统的系统函数分子、分母多项式的系数向量,N为正整数,返回量H则包含了离散系统频响 在 0——pi范围内N个频率等分点的值,向量w则包含范围内N个频率等分点。调用中若N默认,默认值为512。

2.2[H,w]=freqz(B,A,N,’whole’)

该调用格式将计算离散系统在0—pi范内的N个频率等分店的频率响应的值。因此,可以先调用freqz()函数计算系统的频率响应,然后利用abs()和angle()函数及plot()函数,即可绘制出系统在 或 范围内的频响曲线。

例:绘制如下系统的频响曲线

H(z)=(z-0.5)/z

MATLAB命令如下:B=[1 -0.5];

A =[1 0];

[H,w]=freqz(B,A,400,'whole');

H是频率响应的幅度,w是0—pi内的400个点Hf=abs(H);

Hx=angle(H);

clf

figure(1)

plot(w,Hf)

title('离散系统幅频特性曲线')

figure(2)

plot(w,Hx)

title('离散系统相频特性曲线')

幅频特性:

相频特性:

这样画出来的是线性的,要想获得db格式的幅度,需要转换 20*log10(Hf)

之后再画就是db格式的

也可以直接用freqz(b,a,w)这样就会画出幅频响应和相频响应,幅频响应直接是db格式的幅度。

matlab freqz-m,Matlab函数freqs和freqz相关推荐

  1. Matlab的freqz()函数和窗函数法fir1()一起使用

    窗函数法设计FIRDF使用Matlab的fir1()函数十分方便,DSP课堂上已经讲过,此处不再赘述.关心的是: 使用         hn = fir1(M, wc):%M为滤波器阶数M,长度N, ...

  2. MATLAB 函数之 fir1,freqz

    fir1 Window-based finite impulse response filter design 基于窗口的有限脉冲响应滤波器设计 b = fir1(n,Wn) b = fir1(n,W ...

  3. matlab空间重采样,matlab重采样函数

    然后在 Matlab 软件平台下,利用函数 wavread 对语音信号进行采样,记住采样频率 和采样点数.首先画出语音信号的时域波形;然后对语音号进行快速傅里叶变换 (FFT 变换...... col ...

  4. matlab里butter,matlab中butter函数

    利用切比雪夫最佳一致逼近法设计一多阻带陷波器 qiebijinduobs 在 MATLAB 中,可以用函数 y=filter(p,d,x)实现差分方程的仿真,也可以用函 数 y=conv(x,h).. ...

  5. Matlab中的lsqcurvefit函数的使用

    Matlab中的lsqcurvefit函数的使用 lsqcurvefit函数 调用示例 lsqcurvefit函数 非线性曲线拟合是已知输入向量xdata和输出向量ydata,并且知道输入与输出的函数 ...

  6. matlab语言unique,Matlab的unique函数的C++实现

    Matlab中的unique函数,实现的是去除重复元素,只保留一个,且剩下的非重元素按大小排列: C++中stl::unique函数与其不同之处在于:是去除相邻的重复元素,且不改变向量大小,把重复元素 ...

  7. 【 MATLAB 】使用 filter 函数产生给定线性时不变系统的单位阶跃响应

    先给出filter的大致介绍:[ MATLAB ]filter 函数介绍(一维数字滤波器) 上篇博文写了产生一个系统的脉冲响应的博文,有兴趣可以看看:[ MATLAB ]使用 impz 函数计算并画出 ...

  8. matlab拟合未定义函数或变量,关于matlab的未定义函数或变量的问题

    1.首先需要知道matlab中查看数据类型函数是class,可以先在命令行窗口中help class.2.在命令行窗口中输入"a=123 class(a)".3.按回车键可以看到变 ...

  9. matlab 数学库,matlab数学函数库

    (n) 求 n 的阶乘 如何用 matlab 配方 没有发现 matlab 有这一命令,不过我们可以调用 maple 的命令,调用方法如下: 首先加载 maple 中的 student 函数库,加载. ...

  10. 【Matlab 控制】求函数极限

    Matlab 求极限 求函数的0极限: >> syms x; >> y1=(4*x^3-2*x^2+x)/(3*x^2+2*x); >> limit(y1,x,0) ...

最新文章

  1. LuoguP2617 Dynamic Rankings (动态主席树学习理解)
  2. [转]关于安卓与ios的推送系统,我说说自己的看法。
  3. SQL语句中order_by_、group_by_、having的用法区别
  4. Angular jasmine如何从detectChange触发refreshView进而执行到Component的hook实现
  5. [js] 使用ajax请求真的不安全吗?为什么?
  6. 安装TPCC-MySQL报错
  7. java 1000到9999_java – 可以0.99999999999乘以四舍五入到1.0?
  8. curl_multi实现并发
  9. 【CSS】【7】CSS的盒子模型
  10. OpenPose人体姿态识别项目
  11. ant design 时间控件清空值
  12. 贴片电阻电容封装与尺寸,功率,耐压,温度等的关系
  13. 拆机专用磁力桌垫:保证一颗螺丝也不漏网
  14. PredRNN++:网络结构和代码解读
  15. 怎样彻底删除 mac 上的 Adobe 相关文件?
  16. #第四章 星光灿烂4.1苍穹
  17. 在react框架中实现打印过程中踩的坑
  18. MYSQL字段属性之int() 和 tinyint()的区别
  19. LeetCode | 575. Distribute Candies
  20. Dessert(dfs)

热门文章

  1. Java多线程编程 深入详解
  2. 思科交换机配置试题_思科交换机配置命令大全 思科交换机配置手册
  3. 【Pytorch-手写字体识别】手写字体识别项目
  4. PHP实现查找msn联系人
  5. win10电脑性能测试软件,win10系统测试计算机性能的图文方法
  6. 微信公众号服务器管理员权限,公众号管理员和运营者的区别是什么?公众号管理员需要承担责任吗?...
  7. JavaWeb项目--【在线音乐播放器】onlineMusicPlayer
  8. docker版 Transmission以及qbittorrent 下载路径没有权限问题
  9. XPS Silverlight Reader
  10. android 垂直方向布局,Android详解4种基本布局