How do I call "frd" for discrete time?

Question

Mark 2014-11-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/163598-how-do-i-call-frd-for-discrete-time

评论： Mark 2014-11-21

在 MATLAB Online 中打开

Following along this example for mu analysis: Robust Performance Analysis

I have M, and blk from

[M,NDelta,blk] = lftdata(ClosedLoop);

Making a frequency vector and defining my discrete time step as:

freq=logspace(-1,2,100);
dt = 0.05;

I now want to call frd as:

M_g = frd(M,freq,dt);

but get the error: "The frequency units must be specified as one of the strings 'rad/TimeUnit', 'cycles/TimeUnit', 'rad/s', 'Hz', 'kHz', 'MHz', 'GHz', or 'rpm'."

I have no idea how to specify these units though... Help?!?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Sebastian Castro 2014-11-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/163598-how-do-i-call-frd-for-discrete-time#answer_159693

编辑：Sebastian Castro 2014-11-20

在 MATLAB Online 中打开

Hi Mark,

I think you may need to discretize the system M_g first using c2d before using the frd command.

The following commands worked for me:

M_disc = c2d(M,0.05);
M_g = frd(M_disc,freq,'rad/s');

After displaying M_g, at the very end I get the following:

Sample time: 0.05 seconds
Discrete-time frequency response.

Best,

Sebastian

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Mark 2014-11-21

在 MATLAB Online 中打开

Well, my system is already in discrete time with a sample time of 0.05 before I call lftdata, so the use of c2d is unnecessary. Unfortunately, calling

M_g = frd(M_disc,freq,dt)

yields the same error. I believe the third argument is required for discrete time. Calling

M_g = frd(M_disc,freq,dt,'rad/s')

also returns an invalid syntax error.

请先登录，再进行评论。