【电路仿真】基于Simulink钟摆自由控制含Matlab源码.zip

【电路仿真】基于Simulink钟摆自由控制含Matlab源码

1 简介

【电路仿真】基于Simulink钟摆自由控制含Matlab源码

2 部分代码

function [sys,x0,str,ts] = simpendzzy(t,x,u,flag,dampzzy,gravzzy,angzzy)

%SFUNTMPL General M-file S-function template

% With M-file S-functions, you can define you own ordinary differential

% equations (ODEs), discrete system equations, and/or just about

% value of the FLAG, the current state vector, X, and the current

% input vector, U.

%

% FLAG RESULT DESCRIPTION

% ----- ------ --------------------------------------------

% 0 [SIZES,X0,STR,TS] Initialization, return system sizes in SYS,

% initial state in X0, state ordering strings

% in STR, and sample times in TS.

% 1 DX Return continuous state derivatives in SYS.

% 2 DS Update discrete states SYS = X(n+1)

%

% The state vectors, X and X0 consists of continuous states followed

% by discrete states.

%

% Optional parameters, P1,...,Pn can be provided to the S-function and

% used during any FLAG operation.

%

% When SFUNC is called with FLAG = 0, the following information

% should be returned:

%

% actual length for all other flags will be equal to the

% length of the input, U.

% SYS(5) = Reserved for root finding. Must be zero.

% SYS(6) = Direct feedthrough flag (1=yes, 0=no). The s-function

% has direct feedthrough if U is used during the FLAG=3

% call. Setting this to 0 is akin to making a promise that

% U will not be used during FLAG=3. If you break the promise

% then unpredictable results will occur.

% SYS(7) = Number of sample times. This is the number of rows in TS.

%

% (period, offset) information. Where m = number of sample

% times. The ordering of the sample times must be:

%

% TS = [0 0, : Continuous sample time.

% 0 1, : Continuous, but fixed in minor step

% sample time.

% PERIOD OFFSET, : Discrete sample time where

% PERIOD > 0 & OFFSET < PERIOD.

% -2 0]; : Variable step discrete sample time

% specified in TS. When specifying than one

% sample time, you must check for sample hits explicitly by

% seeing if

% abs(round((T-OFFSET)/PERIOD) - (T-OFFSET)/PERIOD)

% is within a specified tolerance, generally 1e-8. This

% tolerance is dependent upon your model's sampling times

% and simulation time.

%

% You can also specify that the sample time of the S-function

% is inherited from the driving block. For functions which

% $Revision: 1.18 $



%

% The following outlines the general structure of an S-function.

%

switch flag,

