Bode Plot axes changes with sampling time?

Question

John 2014-4-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/126399-bode-plot-axes-changes-with-sampling-time

评论： Arkadiy Turevskiy 2014-4-28

Why does the Bode Plot 'bode()' frequency axis change in magnitude if the sampling time 'dt' is changed in magnitude? If dt is changed from 1e-3, to 1e-4, the axis labels go from 10^1 --> 10^4, to 10^2 --> 10^5. But sampling time shouldn't change the axis magnitude, just the accuracy.

dt = 1e-4;
t = (0 : dt : 50)';
uc = chirp(t, 0.1, t(end), 100);
yc = circshift(uc, 100); %phase shift the output
data = iddata(yc, uc, dt);
TF= etfe( data );
bode(TF);

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

John 2014-4-25

Nobody else sees this when running the code above? It seems a basic thing to not work...

请先登录，再进行评论。

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

Answer 1

Arkadiy Turevskiy 2014-4-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/126399-bode-plot-axes-changes-with-sampling-time#answer_134063

编辑：Arkadiy Turevskiy 2014-4-21

Take a look at the examples section here .

The FFT frequency is basically an inverse of sampling time, so when you make dt 10 times smaller, frequency goes 10 time higher, as shown in the example referenced above.

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

请先登录，再进行评论。

Answer 2

John 2014-4-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/126399-bode-plot-axes-changes-with-sampling-time#answer_134064

Right, but efte automatically sets the axes, and it also has the sampling time as an input via 'iddata'. It should already compensate for that, since there is not other provision for entering the time vector.

2 个评论
显示无隐藏无

Arkadiy Turevskiy 2014-4-21

If you take a look at a signal processing book like Discrete-Time Signal Processing by Oppenheim and Schafer, and look at frequency-domain representation of sampling, you will see that all the frequency plots are defined with sampling frequency, 1/T Hz , or 2pi/T rad/sec.

The bode plot will show the right content at the right frequencies, but the range will always extend to (1/t)/2 - half of the sampling frequency. This is the right behavior. If you have t=1e-4, you can still capture frequency content happening at just below half of sampling frequency, 5 Khz. With sampling time 1e-3, you can only capture frequency content up to 500 Hz.

Think about it this way: if your transfer function is a very narrow passband at 100 Hz, no matter what dt is (1e-3 or 1e-4), you will see your passband at 100Hz, but the range will extend to either 500Hz or 5,000 Hz.

John 2014-4-22

But that's just it, the passband is not fixed, and its value varies with dt -- of course, it shouldn't change just because of dt.

Hence the confusion of what 'etfe' is doing.

In other words, as i change 'dt', the plot shown on the screen is exactly the same and remains in place (i.e. it is not squeezed or stretched), and the axis magnitude just changes under it.

For example, the passband could be 100 in one case, and 1000 in another case as dt changes by a factor 10.

That shouldn't happen...

请先登录，再进行评论。