matlab画三维伯德图，bode图

做了一个简单的小例子，有什么疑惑可以留言。

clc
clear
close all
%% 初始化
R_Sta=1;    %R的起始，结束范围以及间隔
R_end=10;
Nap_R=1;
% 数组初始化
num=[R_Sta 1];
den=[5e-4 1];
w={1e-2,1e4};   %确定频率范围，单位：rad/s
[mag,phase,wout]=bode(num,den,w);   % 先要运行一次bode函数以确定输出量的长度

R=R_Sta:Nap_R:R_end;
Num_R=length(R);
Num_w=length(wout);
Z_mag = ones(Num_R,Num_w);
Z_phase = ones(Num_R,Num_w);
[X,Y]=meshgrid(wout,R);

R=R_Sta;
for i=1:Num_R
    num=[R 1]
    den=[5e-3 1]
    [mag,phase,wout]=bode(num,den,w);
    Z_mag(i,:)=20*log10(mag);
    Z_phase(i,:)=phase;
    R=R+Nap_R;
end

% figure(1)
% bode(num,den);
figure(2)
C_mag=Z_mag;
surf(X,Y,Z_mag,C_mag)
set(gca,'XScale','log')
colorbar    
shading interp   % 去掉网格
xlabel('rad/s');ylabel('R');zlabel('mag')
title('Magnitude')
figure(3)
C_phase=Z_phase;
surf(X,Y,Z_phase,C_phase)
set(gca,'XScale','log')
colorbar    
shading interp   % 去掉网格
xlabel('rad/s');ylabel('R');zlabel('Phase')
title('Phase')

输出图形如下：

本文链接：https://blog.csdn.net/qq_36697422/article/details/103882377

智能推荐

web安全简易规范123

web安全，大公司往往有专门的安全开发流程去保证，有专门的安全团队去维护，而对于中小网络公司，本身体量小，开发同时兼带运维工作，时间精力有限，但是，同样需要做一些力所能及的必要的事情。有时候，安全威胁并不是因为你的防盗窗被人撬开了，而是你晚上睡觉的时候忘了关门，而关上门对开发来说也许只是举手之劳。 1、不要用root，确定使用的中间件和框架是否默认打开了后门我们总会在线上使用部署一些中间件、开源...

css弹性盒模型详解----justify-content

本篇文章详细介绍justify-content 效果演示如下：效果演示如下：效果演示如下: 效果演示如下：效果演示如下...

html5拖放--15行js代码实现两个div内容互换

本文首发于我的个人博客:http://cherryblog.site/ ，欢迎大家前去参观本文项目地址，sortable插件地址：https://github.com/sunshine940326/sortable demo地址：https://github.com/sunshine940326/drag 在写我们后台的管理程序中需要有一个拖放的功能，然后我们有一个这样的功能，实现11个固定且大...

git切换分支报错，不管什么标题名字，都报非法字符，所以就不起名字了。

切换分支的时候，报了标题这么个错误，error: ”xxx did not match any file(s) known to git. 看见不能切换分支，我首先 git status 查看了一下当前状态，如下图然后，就会发现，其实我的这个错误非常明显，就是在我的 beat 分支下有文件修改，所以切换不了。ok，解决方法： 1. 如果修改的这些文件没什么用，完全可以删除。（我这儿的...

Oracle分析函数之LEAD和LAG实际应用

Oracle分析函数之LEAD和LAG实际应用在前几天的工作中按照客户的需求，需要对客户信息进行数据分析，即某人存在多个状态的账号，将客户信息账号状态分析出结果，和客户确认汇报，根据保留规则，保留唯一账号，以保证程序可用性。起初，根据聚合函数进行查询分析，需要写一大串的SQL，即不美观又复杂，很容易产生错误。后续想到Oracle分析函数中的lead和lag，SQL简洁了很多且容易产生报告数据。 ...

matlab画三维伯德图，bode图

智能推荐

猜你喜欢