function [sys,x0,str,ts]=ch2example23Sfunshow(t,x,u,flag,v0,theta0,L,m,g,k) 
% 單擺的動(dòng)畫顯示
switch flag, 
    case 0                 % flag=0 初始化 
       sizes = simsizes;           % 獲取SIMULINK仿真變量結(jié)構(gòu)
       sizes.NumContStates  = 0;   % 連續(xù)系統(tǒng)的狀態(tài)數(shù)是0
       sizes.NumDiscStates  = 0;   % 離散系統(tǒng)的狀態(tài)數(shù)是0
       sizes.NumOutputs     = 0;   % 輸出信號(hào)數(shù)目是0
       sizes.NumInputs      = 1;   % 輸入信號(hào)數(shù)目是0
       sizes.DirFeedthrough = 1;   % 該系統(tǒng)是直通的 
       sizes.NumSampleTimes = 1;   % 這里必須為1 
       sys = simsizes(sizes); 
       str = [];           % 通常為空矩陣 
       x0  = [];           % 初始狀態(tài)矩陣x0
       ts  = [-1 0];       
   case 3                     % flag=3 計(jì)算輸出 
       sita=u;
       x=L*sin(sita);         % 計(jì)算擺錘的位置
       y=-L*cos(sita);
       plot([0,x],[0,y],'-o');% 動(dòng)畫作圖   
       text(.5,-.5,['t=',num2str(t)]);% 顯示參數(shù)
       text(.5,-.6,['L=',num2str(L)]); 
       text(.5,-.7,['m=',num2str(m)]); 
       text(.5,-.8,['k=',num2str(k)]); 
       axis equal;       % 坐標(biāo)縱橫比例相同
       axis([-L L -L L]);% 坐標(biāo)范圍固定
       set(gcf,'DoubleBuffer','on');% 雙緩沖避免作圖閃爍
       drawnow;                     % 立即執(zhí)行作圖命令
   case {1,2,4,9}    % 其他作不處理的flag 
       sys=[];       % 無用的flag時(shí)返回sys為空矩陣
otherwise % 異常處理 
       error(['Unhandled flag = ',num2str(flag)]); 
end