php中文网 | cnphp.com

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 319|回复: 0

Matlab绘制系统的根轨迹

[复制链接]

2708

主题

2715

帖子

9656

积分

管理员

Rank: 9Rank: 9Rank: 9

UID
1
威望
0
积分
6826
贡献
0
注册时间
2021-4-14
最后登录
2024-6-16
在线时间
680 小时
QQ
发表于 2022-12-2 08:56:34 | 显示全部楼层 |阅读模式
[mw_shl_code=applescript,true]%% g摄动时的参数根轨迹
clear;clc;close all;
%% 基本参数设定
Vpcc=311;  Rv=0.2180;  Xv=0.005;m=1e-4;Ei=306;
f = 50;  wn = 2*pi*f;  fh=0.2;wp=2*pi*fh;theta=90;
      
%% 求解零点极点和开环增益,并绘图
% for g=0.17e-3:0.01e-3:3.0e-3
  for g=1.408e-3:0.01e-3:6.438e-3
%     theta = 0: 2*pi/point_num:2*pi; %角度[0,2*pi]
    A=(Rv^2*Ei*Vpcc*sin(theta)-2*Rv*Xv*Ei*Vpcc*cos(theta)+2*Rv*Xv*Vpcc^2-Xv^2*Ei*Vpcc*sin(theta))/(Xv^2+Rv^2)^2;
    B=((Xv*cos(theta)-Rv*sin(theta))*Ei*Vpcc)/(Xv^2+Rv^2);
    C=((Rv^2-Xv^2)*Ei*Vpcc*cos(theta)+(Xv^2-Rv^2)*Vpcc+2*Rv*Xv*Ei*Vpcc*sin(theta))/(Xv^2+Rv^2)^2;
    D=-(Ei*Vpcc*Xv*sin(theta)+Ei*Vpcc*Xv*cos(theta))/(Xv^2+Rv^2);
    a=2*wp-C*g^2*wp;
    b=B*2*pi*m*wp-wp^2*C*g^2+wp^2;
    c=(A*D*2*pi*m*g^2+B*2*pi*m-B*2*pi*m*C*g^2)*wp^2;
%     d=2*pi*m*wp;
    p=[1 a b c];
    r=roots(p);
    p1=r(1);p2=r(2);p3=r(3);
%     i=1:1:100;
%     a(i,1)=z(1);a(i,2)=z(2);
%     a(i,3)=p(1);a(i,4)=p(2);a(i,5)=p(3);a(i,6)=p(4);
%     scatter(real(z1),imag(z1),'ko');%画出散点图
    scatter(real(p1),imag(p1),'mx');
    hold on;
    scatter(real(p2),imag(p2),'bx');
    hold on;
    scatter(real(p3),imag(p3),'kx');
    hold on;
    grid on;
end
% legend('z1','z2','p1','p2','p3','p4','FontSize',12); % 图形注释
% title('kp变化时参数根轨迹');
% xlabel('实轴 /s^{-1})','FontName','宋体');
% ylabel('虚轴 /s^{-1})','FontName','宋体');
set(gca,'fontsize',7);%图中所有字体设为10榜
set(gca,'box','on','linewidth',1)%图窗四边边框全显示,并设置1榜线型
% axis([-1000,50,-700,700])[/mw_shl_code]





上一篇:每隔2周休息1天工厂日历生成
下一篇:卡尔曼滤波
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|php中文网 | cnphp.com ( 赣ICP备2021002321号-2 )

GMT+8, 2024-6-18 13:37 , Processed in 0.237359 second(s), 38 queries , Gzip On.

Powered by Discuz! X3.4 Licensed

Copyright © 2001-2020, Tencent Cloud.

申明:本站所有资源皆搜集自网络,相关版权归版权持有人所有,如有侵权,请电邮(fiorkn@foxmail.com)告之,本站会尽快删除。

快速回复 返回顶部 返回列表