频率响应分析--伯德图（matlab中bode函数的使用）

1 伯德图是什么

伯德图
bode根据系统动力学原理能画出频率响应（幅度响应（dB）和相位响应（deg：度数）），如下图所示

% 画一个系统函数的频域响应图
H = tf([-2 3.6 -32 1.7],[1 3.3 16.8 25.4 40 13 ]);
bode(H)

2 bode函数的使用
这是MATLAB文档中关于bode函数的所有使用语句，下面对每一条语句进行学习，每个语句都分为
1 语句
2 解释
3 程序例子
4 图片结果

bode(sys)%画一个伯德图
bode(sys1,sys2,...,sysN)% 在同一图上绘制多个动力系统的频率响应
bode(sys1,LineSpec1,...,sysN,LineSpecN)%为图中的每个系统指定颜色、线条样式和标记。
bode(___,w)%绘制系统响应以w指定的频率。
[mag,phase,wout] = bode(sys)
[mag,phase,wout] = bode(sys,w)
[mag,phase,wout,sdmag,sdphase] = bode(sys,w)

- bode(sys1,sys2,…,sysN)

在同一图上绘制多个动力系统的频率响应。所有系统必须有相同数量的输入和输出。

sys1=tf([-2 3.6 -32 1.7],[1 3.3 16.8 25.4 40 13 ]);
sys2=tf([-2 3.6],[1 3.3 16.8 25.4]);
sys3=tf([-2],[1 3.3 16.8]);
bode(sys1,sys2,sys3)%画多个伯德图
grid%画格子

- bode(___,w)

绘制系统响应以w指定的频率。如果w是形式为{wmin,wmax}的单元格数组，则bode绘制在wmin和wmax之间的频率范围内的响应。如果w是一个频率向量，则bode绘制每个指定频率的响应曲线。

sys1=tf([-2 3.6 -32 1.7],[1 3.3 16.8 25.4 40 13 ]);
sys2=tf([-2 3.6],[1 3.3 16.8 25.4]);
sys3=tf([-2],[1 3.3 16.8]);
w=[1e-2,1e+5];
w=[1e+2,1e+5];
w={0.01,100}
bode(sys1,'k.',sys2,'k-',sys3,'k*',w)%画多个伯德图
grid%画格子

参考一下常用的画图符号

这张图来自这篇文章：点击进入

[mag,phase,wout] = bode(sys,w)
返回幅度、相位、频率向量，我们获得了数据，就可以对数据进行处理，也可以用自己的画图方法和单位。
- 下面用这个例子作为总结，也是对前面所有的复习：

%% ---------------------------------
%2017 zhao  3/5 jie gaosi wavelet filter
clear, close all
sys1=tf([-2 3.6 -32 1.7],[1 3.3 16.8 25.4 40 13 ]);
sys2=tf([-2 3.6],[1 3.3 16.8 25.4]);interval=100;%w的另一种定义：频率限制在1e-2-1e+3之间，整个图像用interval个点来画
w = logspace(-2, 3, interval);
% ********************************************************************
figure(1)
bode(sys1,'k.',sys2,'k-',w)%画多个伯德图
% ********************************************************************
w={0.01,100}%w的一种定义：频率限制在0.01-100
figure(2)
num =[-2 3.6 -32 1.7];
den = [1 3.3 16.8 25.4 40 13 ];
[mag, phase, w] = bode(num, den, w);%用bode(sys1），会将mag转化为数组，导致无法画图
magdB = 20*log10(mag);%输出的为hz，转化为rad/s
% dBmax = 5*ones(1,interval); %画一下横线，可以实现控制上下限的效果
% dBmin = -50*ones(1,interval); %点数必须与 w 中频率点的数目相等
subplot(2, 1, 1)
% semilogx(w, magdB, '.', w, dBmax, 'r-', w, dBmin, 'r:')
semilogx(w, magdB, '-')
grid
title('BodeDiagram of G(s)')
xlabel('Frequency (rad/sec)')
ylabel('Gain dB')
subplot(2, 1, 2)
% pmax = 150*ones(1,interval);
% pmin = -150*ones(1,interval);
% semilogx( w, phase, '.', w, pmax, '--', w, pmin, ':')
semilogx( w, phase, '-')
grid%画格子

频率响应分析--伯德图（matlab中bode函数的使用）相关推荐

matlab里butter,matlab中butter函数
利用切比雪夫最佳一致逼近法设计一多阻带陷波器 qiebijinduobs 在 MATLAB 中,可以用函数 y=filter(p,d,x)实现差分方程的仿真,也可以用函数 y=conv(x,h).. ...
MATLAB中round函数的使用
MATLAB中round函数的使用调用格式: Y = round(x) 在matlab中round函数是一个四舍五入的函数,下面举个例子来说: 这个就看起来很明了,我们调用round(a)函数,对a ...
matlab中numden函数,numden函数 MATLAB里面numden函数有什么用？
MATLAB里面numden函数有什么用? Examples:[n,d]=numden(sym(4/5))returnsn=4andd=5.[n,d]=numden(x/y y/x)returnsn= ...
matlab中频域信号IFFT,MATLAB中ifft函数用法、性质、特性-以及与fft的组合应用全面深入解析(含程序)...
MATLAB中ifft函数用法.性质.特性,以及与fft的组合应用全面深入解析(含程序) 前言在我之前的<MATLAB中fft函数用法.性质.特性.缺陷全面深入解析(含程序)>中,我已经 ...
MATLAB中freqz函数的用法
本文引用自乱七八糟<MATLAB中freqz函数的用法> FREQZ 是计算数字滤波器的频率响应的函数 [H,W] = FREQZ(B,A,N)returns the N-point c ...
ssbdemod在matlab中的意义,matlab中modulate函数
(7)利用 MATLAB 中的 demodulate 和 modem.qamdemod 函数生成解调器对 16QAM 信号的解调,并将十六进制信号转化成二进制比特流信息. (8)用得到比特流信息除以. ...
matlab函数乘阶跃,MATLAB中FFT函数和点乘、乘、单位阶跃函数
注意问题:1, 对于单位阶跃函数u(t-b)和Dirac函数delta(t-a),可以借用Maple函数库中的定义,他们分别是Heaviside(t-b),Dirca(t-a). 2,点乘和点除代表矩 ...
Matlab中bwmorph函数的使用
Matlab中bwmorph函数的使用 Matlab中提供了一个基于形态学的处理函数,即以膨胀.腐蚀等操作为基础,其语法格式如下: bw2=bwmorph(bw1,operation,n); 其中bw ...
9.matlab中repmat函数
来源: matlab中repmat函数的用法 - CSDN博客 https://blog.csdn.net/anqier1009/article/details/5214978 B = repmat( ...

频率响应分析--伯德图（matlab中bode函数的使用）

频率响应分析--伯德图（matlab中bode函数的使用）相关推荐

最新文章

热门文章