二阶广义积分离散方法遇到了难题

Question

目前正在研究，正弦波经过二阶广义积分后得到90的移相，然后离散化时出现问题，离散后波形冲上了云霄之上，采样时间为60khz，matlab自带的tustin等方法都试了，全部超调，瞅瞅我的问题好么？
wt=2pi50

wt =

314.1593

prs=tf([wtwt],[1,wt,wtwt])

prs =

      9.87e04

---

s^2 + 314.2 s + 9.87e04

Continuous-time transfer function.

prz = c2d(prs,1/50000,'tustin')

prz =

9.839e-06 z^2 + 1.968e-05 z + 9.839e-06

      z^2 - 1.994 z + 0.9937

s函数为

但是当我把采样频率调到6khz时，波形正常，不知道问题在哪里

Answer 1

matlab 2020b仿真结果正常

wt=2*pi*50;
prs = tf([wt*wt],[1,wt,wt*wt]);
ts = 1/60e3;   % 采样时间
prz = c2d(prs,ts,'tustin');
t = 0:ts:0.2;
u = sin(100*t);
y = lsim(prz,u,t);

figure
plot(t,u,t,y)