MATLAB的答案

c2d函数“zoh”方法的公式

271(30天)
Tamas Toth-K
Tamas Toth-K 2020年10月4日
回答: 拉斐尔安图内斯2022年10月14日
这就跟你问声好！
我的问题,说我有一个连续时间PI控制器:
我应用函数汇集( , ,“zoh”),然后似乎matlab使用他替换:
。
这是完全的好,但是如果我有一个传递函数:
然后应用汇集功能汇集时( , zoh),看来以前替换在这种情况下更改(s ! = (z 1) / Ts)取决于的价值 。我想知道为什么Matlab这,和它的代数公式如果可能的话是什么?
谢谢你！

在回答这个问题。

接受的答案

保罗
保罗 2020年10月4日
编辑:保罗 2022年4月11日
跑:
我怀疑在你的第一个案例C替换你引用的形式只适用,因为C .一般来说,ZOH近似不使用替换。尽管可能没有实现这种方式,genaral形式ZOH近似可以实现如下所示,相比和Matlab产生
P = 5, Ti = 6; t = 1;%示例数据
C =特遣部队(P * [Ti 1], [Ti 0]);
Cz =汇集(C, Ts,“zoh”);
Cznew = minreal(汇集(C *特遣部队(1 0 [1]),Ts,“冲动”)*特遣部队([1],[1 0],Ts) / Ts);
(Cz Cznew)
ans =从输入1输出:5 z - 4.917 - - - - - - - - - - - - z - 1输入2输出:5 z - 4.917 - - - - - - - - - - - - z - 1样品时间:0.1秒的离散传递函数。
R = 10, Te = 5;W =特遣部队(1,[R * Te 1]);
Wz =汇集(W, Ts,“zoh”);
Wznew = minreal(汇集(W *特遣部队(1 0 [1]),Ts,“冲动”)*特遣部队([1],[1 0],Ts) / Ts);
(Wz Wznew)
ans =从输入1输出:0.001998 - - - - - - - - - - - -从输入2输出z - 0.998: 0.001998 - - - - - - - - - - - - z - 0.998样品时间:0.1秒的离散传递函数。
你引用的替换,s = (z 1) / Ts,远期矩形规则,这似乎是一个容许,虽然非法,汇集的输入方法。
编辑:矩形规则是没有实现。当指定方法输入字符串,只有第一个字符的问题,所以“向前”实际上是“呸”
H =特遣部队(1,(1 1 1));
t = 0.1;
H1 =汇集(H, Ts,“呸”);
H2 =汇集(H, Ts,“前进”);
H3 =汇集(H, Ts,“福”);
H1, H2
ans = 0静态增益。
H1 - H3
ans = 0静态增益。
4评论
显示3年长的评论
Tamas Toth-K
Tamas Toth-K 2020年10月7日
非常感谢你,!后说,他同我又看,意识到许多替代变量是不正确的,在修复它是相同的:)

登录置评。

更多的答案(4)

拉斐尔安图内斯
拉斐尔安图内斯 2022年10月3日
n = 1.501877347 e + 06;
125.1564456 d = 4.255319149 [1 e + 03 e0 + 03];
g =特遣部队(n, d)
广州=汇集(g、0.000369“zoh”)
极(广州)
零(广州)
1评论
保罗
保罗 2022年10月3日
如果这个答案应该是一个新问题,我建议删除这个答案,发布一个新的问题(更多细节的问题到底是什么,除了代码)。

登录置评。

拉斐尔安图内斯
拉斐尔安图内斯 2022年10月4日
n = 1.502 e + 06;
d = [1 4.255 125.1 e + e + 03 03];
g =特遣部队(n, d)
广州=汇集(g、0.0000125“zoh”)
特遣部队(z =“z”);
c = (z - 0.9998) / (z 1)
步骤(反馈(广州* c / (1 + gz * c)))
拉斐尔安图内斯
拉斐尔安图内斯 2022年10月4日
跑:
% | - - - - - - - - - - - - - - - - - Funcao de transferencia em节奏数字低音- - - - - - - - - - - |
n = 1.502 e + 06;
d = [1 4.255 125.1 e + e + 03 03];
g =特遣部队(n, d)
g = 1.502 e06 - - - - - - - - - - - - - - - - - - - - - - - - ^ 2 + 4255年代+ 125100连续时间传递函数。
rlocus (g)
% | - - - - - - - - - - - - - - - - - - -Controlador em节奏数字低音- - - - - - - - - - - - - - - - - - - |
s =特遣部队(“年代”);
c = 1.1 * ((s + 32.5) / (1.2 s +))
c = 1.1 + 35.75 - - - - - - - - - - - - - s + 1.2连续时间传递函数。
% | - - - - - - - - - - - - - - - - - - - - - - - -注意de Malha fechada连续- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
女朋友= g * c / (1 + g * c)
女朋友= 1.652 e06 s ^ 4 + 7.086 e09 ^ 3 + 4.437 e11年代^ 2 + 7.24 e12汽油(s + 8.061 e12汽油- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ^ 6 + 8512年代^ 5 + 2.003 e07 s ^ 4 + 8.194 e09 ^ 3 + 4.619 e11年代^ 2 + 7.279 e12汽油(s + 8.083 e12汽油连续时间传递函数。
极(gf)
ans = 6×1
1.0 e + 03 * -4.2254 -3.7896 -0.4338 -0.0328 -0.0296 -0.0012
零(gf)
ans = 4×1
1.0 e + 03 * -4.2254 -0.0325 -0.0296 -0.0012
步骤(gf, 0.5)
% | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ZOH - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
广州=汇集(g、0.0000125“zoh”)
广州z = 0.0001153 + 0.0001133 - - - - - - - - - - - - - - - - - - - - - - - - z z ^ 2 - 1.948 + 0.9482样品时间:1.25 e-05秒离散传递函数。
极(广州)
ans = 2×1
0.9996 - 0.9486
零(广州)
ans = -0.9824
% | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Discretizacao - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
cz =汇集(c, 0.0000125,“tustin”)
cz z = 1.1 - 1.1 - - - - - - - - - - - - z - 1样品时间:1.25 e-05秒离散传递函数。
cprewarping =汇集(0.0000125 c, c2dOptions (“方法”,“tustin”,“PrewarpFrequency”1))
cprewarping z = 1.1 - 1.1 - - - - - - - - - - - - z - 1样品时间:1.25 e-05秒离散传递函数。
cmpz =汇集(c, 0.0000125,“匹配”)
cmpz z = 1.1 - 1.1 - - - - - - - - - - - - z - 1样品时间:1.25 e-05秒离散传递函数。
% | - - - - - - - - - - - - - - - - - - - - - - - -注意de Malha fechada discreta - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
地球科学=广州* cz / (1 + gz * cz)
地球科学= 0.0001268 z ^ 5 - 0.0003761 z ^ 4 + 0.0002492 z ^ 3 + 0.0002407 z z ^ 2 - 0.0003587 + 0.0001181 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - z ^ 6 - 5.896 z ^ 5 + 14.48 z ^ 4 - 18.97 z ^ 3 + 13.98 z z ^ 2 - 5.493 + 0.8992样品时间:1.25 e-05秒离散传递函数。
极(德国)
ans =
1.0002 + 0.0000我0.9995 + 0.0004 0.9995 - 0.0004 0.9946 0.9486我0.9539 + 0.0000 + 0.0000 + 0.0000
零(德国)
ans = 5×1
-0.9824 1.0000 0.9996 0.9996 0.9486
% | - - - - - - - - - - - - - - - - - - - - - - - - -Resposta ao degrau (aproximacoes discretas) - - - - - - - - - - - - - - - - - - |
步骤(反馈(cmpz *广州,1),“r”,0.5)
持有在
步骤(反馈(cprewarping *广州,1),“b”,0.5)
步骤(反馈(cmpz *广州,1),‘G’,0.5)
传奇(“Tustin”,“Tustin prewarping”,“合成”)
% | - - - - - - - - - - - - - - - - - - - - - - - -数字低音x discreto - - - - - - - - - - - - - - - - - - - - - - - - - - - - - |
步骤(反馈(cmpz *广州,1),“r”,0.5)
持有在
步骤(火焰杯,“b”,0.5)
步骤(反馈(cmpz *广州,1),“r”,0.05)
持有在
步骤(火焰杯,“b”,0.05)
拉斐尔安图内斯
拉斐尔安图内斯 2022年10月14日
n = -20.27;
d = (1 0 0。-66.93);
g =特遣部队(n, d)
s =特遣部队(“年代”);
pi1 = 1.3 * s + 6;
g1f =反馈(-13.05 * g软件)
辅助= (-0.1628 * s ^ 2 + 9.81) / s ^ 2
皮= 0.8 * s + 0.7;
g2f =反馈(皮* g1f *辅助,1)
步骤(g2f)

类别

找到更多的在z变换在帮助中心和文件交换

标签

下载188bet金宝搏


释放

R2018b

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!