目前尚不清楚我们的目标是什么。你从A, B, C, D,并试图找出的分区是基于知道维度w, u, z,和y ?或者你从分区和试图形成一个党卫军对象代表分区的系统?很难说因为sys示例的输出,所以实在是没有任何分区。
在不确定系统矩阵分区问题:如何找到矩阵A、B1、B2, C1, C2,这里,D12, D21, D22摊位吗?
15的观点(30天)
显示旧的评论
大家好!
我得到不确定系统使用这段代码:
a1 =尿素的(“a1”p_p2,“范围”,p_p2 - 0.0070 p_p2 + 0.0070);
a2 =尿素的(“a2”p_p3,“范围”,p_p3 - 0.0128 p_p3 + 0.0128);
b0 =尿素的(“b0”p_z1,“范围”,p_z1 - 0.0019 p_z1 + 0.0019);
b1 =尿素的(“b1”p_z2,“范围”,p_z2 - 0.0145 p_z2 + 0.0145);
(0 = 1;a2 a1];
B = [0;1);
C = (b1 b0);
sys = ss (A, B, C, 0);
sys.InputGroup。ActualIn = 1;
sys.OutputGroup。ActualOut = 1;
[Msysδ]= lftdata(系统);
(当然p_p2, p_p3 p_z1 p_z2有价值,但是我为了简单起见省略方程)
所以Msys表单:
但是我搜索一种分区,B, C, D, B1, B2, C1, C2,这里,D12, D21, D22摊位,这种形式:
我所知道的是,应该是一样的,D22摊位应该是0,但是如何实现这个分区?
由于任何人都可以帮助我。
弗朗西斯科
接受的答案
保罗
2021年2月2日
假设您有一个模型
xdot = * x + B1 * w + B2 * u
z = C1 * x +这里* w + D12 * u
y = x + C2 * D21 * w + D22摊位* u
这个系统可以写成:
xdot = * x + [B1 B2] * [w; u] = x + B * * [w; u]
[z; y] = (C1, C2) * x +[这里D12; D21 D22摊位]* [w; u] = x + D C * * (w; u);
假设你有n个状态,w是西北的一个向量输入,和uν输入向量,y和z是新西兰输出向量,向量的纽约输出。然后我们必须
B1 = B(: 1:西北);
B2 = B (:, nw + 1:西北+ν);%或B (nw + 1::结束)
C1 = C(1:新西兰:);
C2 = C(新西兰+ 1:新西兰+纽约);%或C(新西兰+ 1:最终,);
这里= D(1:新西兰,1:西北);
D12 = D(1:新西兰,nw + 1:西北+ν);%或D(1:新西兰,nw + 1:结束)
D21 = D(新西兰+ 1:新西兰+纽约,1:西北);%或D(新西兰+ 1:结束,1:西北)
D22摊位= D(新西兰+ 1:新西兰+纽约,nw + 1:西北+ν);%或D(新西兰+ 1:最终,nw + 1:结束)
w的维度,u, z,和y必须已知确定的分区,B, C, D兼容的维度。
7评论
穆罕默德Yakoob说道
2022年2月17日
你好亲爱的再次
上面是我的代码但不是最后一个,我处理uncertian系统,这也添加一些defctes我的代码,我想detrmine他H_inf来设计一个控制器
答案(1)
保罗
2022年2月17日
编辑:保罗
2022年2月19日
% %植物Paramiters
t = 0.0005;
f = 60;
我们= 2 *π* f;
Vdc = 300;
C_st =尿素的(“C_st”15 * 10 ^ 6,“比例”,-40年,40);
L_t =尿素的(“L_t”2 * 10 ^ 3,“比例”,-30,30);
R_l = 40;
% %实现单相微型智能电网" contious时间形式
基于代码的编写,输出方程Y是Y =残雪+杜+ 0 w
% x_dot = AX +不+ dW;Y =残雪+ DU。
% X = (il Vg) ';U =(甚短波);d = (ig);Y = (Vg)。
Aop = [0 1 / L_t;1 / C_st 0);
防喷器= [1 / L_t;0);
d = [0;1 / C_st];
警察= [0 1];
夹住= [0 0];
% B = (Bop d);
H-inf控制扰动输入,w,通常是“前”的控制输入,u。
% sys1 = ss (Aop,(防喷器d),警察,夹住);
sys1 = ss (Aop, [d防喷器],警察,夹住);
此时sys1只有两个输入(w, u)和一个输出,y。
大小(sys1)
我不明白下一行。如果sys2应该sys1离散近似,那么它应该使用汇集()。因为sys2并不相关的实际定义分区的问题,我要把sys2 sys1,这样没有任何问题后续sys2引用。
% sys2 = ss (Aop,(防喷器d),警察,夹住,Ts);
sys2 = sys1;
广告= sys2.A;
Bd = sys2.B;
Cd = sys2.C;
Dd = sys2.D;
%确定尺寸
ns =大小(Aop, 1);%的状态数
ν=大小(防喷器,2);%的致动器输入
纽约=大小(警察,1);%的数量测量输出
目前为止一切都很顺利。现在以下行表明,应该有尽可能多的控制输出,z,因为有测量输出y(为什么使用nq、而不是新西兰?)。但有一个定义的输出方程方程z的x, w, u, z是不包括在sys1作为输出。
nq =大小(警察,1);%的数量调节输出
西北=ν;%的干扰
因为我不知道z应该是什么,我只是假设z = 2 * y为目的的说明。在这种假设下,输出方程z z = 2 *警察* x + 2 *夹住* [w; u]。重建与z作为第一个输出和y sys1第二输出,和重复的所有后续代码。
sys1 = ss (Aop, [d防喷器],[2 *警察;警察],[2 *夹住;夹住);
大小(sys1)
sys2 = sys1;
广告= sys2.A;
Bd = sys2.B;
Cd = sys2.C;
Dd = sys2.D;
%确定尺寸
ns =大小(Aop, 1);%的状态数
ν=大小(防喷器,2);%的致动器输入
纽约=大小(警察,1);%的数量测量输出
nq =纽约;%的数量调节输出
西北=ν;%的干扰
现在的分区
% % 9矩阵表示
一个=广告;
B1 = Bd(:, 1:西北);
B2 = Bd (:, nw + 1:西北+ν);
C1 = Cd (1: nq,:);
C2 = Cd (nq + 1: nq +纽约);
这里= Dd (1: nq, 1:西北);
D12 = Dd (1: nq、nw + 1:西北+ν);
D21 = Dd (nq + 1:结束,1:西北);
D22摊位= Dd (nq + 1: nq + ny, nw + 1:西北+ν);
disp (“完成”)
