系统传递函数的相位和幅值绘图-matlab-使用mesh函数
一、实验目的
1.了解熟悉MATLAB中meshgrid、mesh绘图函数的使用;
2.回忆传递函数,利用meshgrid和mesh函数在三维上绘出传递函数图像。
二、实验原理
传递函数是零初始条件下线性系统响应(即输出)量的拉普拉斯变换(或z变换)与激励(即输入)量的拉普拉斯变换之比,记作G(s)=Y(s)/U(s),其中Y(s)、U(s)分别为输出量和输入量的拉普拉斯变换。
三、实验设计及操作
1、选定自变量幅度mag和相位phs的范围,取相同数量的点;
2、通过meshgrid由向量转变为矩阵M和P;
3、写出传递函数表达式;
4、基于矩阵M和P作相应的计算并绘图。
四、实验结果及分析
作图结果与示例相一致,结果正确。
从三维图像中可以看出高通滤波器,与传递函数计算出的结果相一致,印证了理论结果,实验结果非常完美。
五、实验总结和提升
实现了MATLAB中传递函数的三维局限化,加深了传递函数的理解,为日后的其他信号处理打下基础。
六、源代码
clear
clc
close all
%% 示例
z1=0.5*exp(1j*3);
for mag=0.2:0.07:0.8
for phs=-pi+0.01:pi/20:pi
x=mag*cos(phs);
y=mag*sin(phs);
H=1/(x+y*1j-z1);
subplot(221);plot3(x,y,abs(H),'.');
hold on;grid on;
subplot(222);plot3(x,y,unwrap(angle(H)),'.');
hold on;grid on;
end
end
%% mesh作图
mag=0.2:0.07:0.8;
phs=-pi:pi/20:pi;
[M,P]=meshgrid(mag,phs);
H=1./(M.*exp(P.*1j)-z1);
figure(2)
subplot(121);mesh(M.*cos(P),M.*sin(P),abs(H));
hold on;grid on;
subplot(122);mesh(M.*cos(P),M.*sin(P),angle(H));
hold on;grid on;
系统传递函数的相位和幅值绘图-matlab-使用mesh函数相关推荐
- 基于matlab测幅值,基于MATLAB的分析测试系统
基于MATLAB的分析测试系统 来源:华强电子网 作者:华仔 浏览:663 时间:2016-08-10 14:18 标签: 摘要: 摘 要:介绍基于matlab的分析测试系统的结构,详细阐述了系统的工 ...
- matlab采交流电压幅值,基于MATLAB的电力电子技术仿真分析
9 如上图所示,当晶闸管VT处于断态时,电路中电流Id=0,负载上的电压为0,U2 全部加在VT两端,在触发角α处,触发VT使其导通,U2加于负载两端,由于电感L的存在使电流id不能突变,id从0开始 ...
- matlab如何将函数做成幅度谱,好用的画包络谱和幅值谱matlab函数 | 学步园
画包络函数: function hua_baol (y,fs,style,varargin) %画包络函数是hua_baoluo的代码简化版 %输入参数要么是3个,要么是5个 %当三个的时候 y_hh ...
- 画幅值matlab,好用的画包络谱和幅值谱matlab函数
使用范例: x=@(t) (1+0.5*cos(9*pi*t)).*cos(200*pi*t+2*cos(10*pi*t))+sin(pi*t).*sin(30*pi*t); t=0:0.01:9.9 ...
- matlab fft 画出幅值,画包络谱和幅值谱matlab函数示例代码
使用范例: x=@(t) (1+0.5*cos(9*pi*t)).*cos(200*pi*t+2*cos(10*pi*t))+sin(pi*t).*sin(30*pi*t); t=0:0.01:9.9 ...
- 画伯德图及相位、幅值裕度
s域: 无滞后: num=[122.0501] den=[0.8190 1] g=tf(num,den); margin(g); grid on; 有滞后: >> num=[106.938 ...
- matlab如何求传递函数的幅值_MATLAB中求开环传递函数的幅值裕度、相位裕度、截止频率的margin()函数用法...
一.幅值裕度.相位裕度 相位裕度γ就是Bode图中幅值为0dB时对应的相位加上180°,如图中0dB时的相位是-99.1°,所以该系统的相位裕度就是 -99.1°+180°=80.9°.所谓相位裕度的 ...
- MATLAB中求开环传递函数的幅值裕度、相位裕度、截止频率的margin()函数用法
一.幅值裕度.相位裕度 相位裕度γ就是Bode图中幅值为0dB时对应的相位加上180°,如图中0dB时的相位是-99.1°,所以该系统的相位裕度就是 -99.1°+180°=80.9°.所谓相位裕度的 ...
- matlab复数的相位,复数的幅值和相位
[i]); x[i]=x[i]*180/PI-90; printf("第%d次谐波的相位为为%f \\n", i , x[i]) ; } printf("第%d次谐波的幅 ...
最新文章
- 嵌套中断共享资源问题及解决
- 解决:The application could not be installed: INSTALL_FAILED_SHARED_USER_INCOMPATIBLE
- NPOI读取Excel模板并向其中写入数据
- linux apache设置web访问重定向_从零开始学Linux运维|30.Linux的目录结构
- [转]Teach Yourself Programming in Ten Years——用十年教会自己编程
- 【BZOJ2655】calc,dp+拉格朗日插值法
- Theano at a Glance
- RHEL脚本更换YUM源(sohu)
- C#窗体excel与dbf的导入导出
- 稀土铕配合物掺杂聚苯乙烯荧光微球/含铕配合物聚苯乙烯荧光微球/稀土磁性荧光微球制备
- Web项目 引入Editormd
- arp-scan使用
- 阿里云智能编码插件Cosy,提升开发效率杠杠滴!
- Android默认时区、语言设置
- B站纪录片发布会,发布的不只是纪录片
- 外贸人如何快速学好英语
- ktor启动报错:Module function cannot be found for the fully qualified name 'ApplicationKt.module'
- volo拉丁语_拉丁语:Latin Quotes and Phrases (V)
- 【渝粤教育】广东开放大学 古代汉语 形成性考核 (28)
- 电商常用功能模块的数据库设计
热门文章
- PMP续证有新的证书吗?
- 目标跟踪评估绘图(1):Windows10环境将pycharm跑出来的OTB100数据集的.txt结果文件,转换成.mat文件
- ATC计算机会议,实验室三项成果分别被计算机系统、数据库、分布式计算重要会议 USENIX ATC、ICDE、IPDPS收录...
- Hive 窗口函数 实现原理
- Quartus II11.0破解注意
- 如何快速提取文件名称?
- 收费邮箱有哪些?收费邮箱有哪些真正的价值?
- 怎么用oracle查询一个省,知道区县的代码如何一次得到区县所在的省份和城市?各位朋友帮帮忙。。。谢谢...
- 2022-2028全球与中国荧光增白剂 OB-1市场现状及未来发展趋势
- MapGIS IGServer for java