function y=mgr 号16.61-MGR with radial MsD onboard 8 Prof. How 9 Call using lines at the bottom after the "return f nargin 1; w=2; end 0=1: g radial offset K=8: M=l, random model parameters g state space model of the dynamics g state is [R A=[01;-(2*K/M-w^2)0]; 号 weird matlab form Sys=SS(A, B, C, D) T=[0:.01:5] 8 Linear model SIMulation RR=lsim(sys, ones (size(T)),T) g basic offset is Ro, so y actual=RO+R figure(1);号c1f p1ot(T,R0+RR(:1),[min(T)max()],R0+R0/(2*K/M/w^2-1)*[11]," Linewidth',2); title(['Freg =' num2str(w)]) axis([min(T)max(r)0 1. 5*max(Rr(:,1))+RO]) xlabel('time) ylabel('R') return subplot (311); mgrl(1); subplot(312); mgr1(2): subplot(313)i mgr1(3.8)