關(guān)于matlab中sfunction參數(shù)設(shè)置的兩個例子講解詳盡受益匪淺

1、例10.1設(shè)系統(tǒng)的微分方程為:試建立系統(tǒng)仿真模型。例10.2利用Simuliiik構(gòu)建函數(shù)|II|線尸5P+16。例10.3利用Simuliiik仿真求1=xln(l+x)dx。例10.4有初始狀態(tài)為0的二階微分方程xF.2x40.4滬0.2u(t),其中u(t)是單位階躍函數(shù)，試建立系統(tǒng)模型并仿真。例10.5PID控制器是在H動控制中經(jīng)常使用的模塊，在工程應(yīng)用中其標(biāo)準(zhǔn)的數(shù)學(xué)模型為gs/N)EG)其中采用了一階環(huán)節(jié)來近似純微分動作，為保證有良好的微分近似的效果，一般選NN10。試建立PLD控制器的模型并建立子系統(tǒng)。例10.6利用使能子系統(tǒng)構(gòu)成一個1E弦半波整流器。例10.7利用觸發(fā)子系統(tǒng)將一鋸

2、齒波轉(zhuǎn)換成方波。例108采用S函數(shù)實現(xiàn)模塊尸nx,即模塊的功能是把一個輸入信號n倍以后再輸出。(1)利用MATLAB語言編寫S函數(shù)。程序如下：咎犬大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大咎S函數(shù)timesn.m,其輸出是輸入的n倍咎犬大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大functionsys,x0fstruts=timesn(tufflagfn)switchflagfcase0$初始化sysfxOzstr,ts=mdllnitializeSizes;case3$計算輸出量sys=

3、mdlOutputs(tfxzuzn);case2,4,9$未使用的flmg值sys=;otherwise冬出錯處理error(1Unhandieflag=1,num2str(flag);end咎犬大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大%mdlInitializeSizes:當(dāng)flag為0時進(jìn)行整個系統(tǒng)的初始化咎犬大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大大 functionsyszxOfstruts=mdllnitializeSizes(T)E調(diào)用函數(shù)simsizes以創(chuàng)建結(jié)構(gòu)體si

4、zessizes=simsizes;$用初始化信息填充結(jié)構(gòu)體sizessizes.NumContStates=0;sizes.NumDiscStates=0;sizes.NumOutputs=l;sizes.NumInputs=l;sizes.DirFeedthrough=l;sizes.NumSampleTimes=l;$根據(jù)上面的設(shè)置設(shè)定系統(tǒng)初始化參數(shù)sys=simsizes(sizes);$給其他返回參數(shù)賦值。x0=;str=;ts=-lz0;$初始化子程序結(jié)束紅無連續(xù)狀態(tài)$無離散狀態(tài)%有一個輸出量E有一個輸入信號%輸出量中含有輸入量E設(shè)置初始狀態(tài)為零狀態(tài)$將str變量設(shè)置為空字符串$假

5、定繼承輸入信號的采樣周期咎犬犬大大大大大大大犬大大大大大大大犬犬犬大大大大大犬犬犬大大大大大大犬犬犬犬大大大大大犬犬大大大大%mdlOutputs：當(dāng)flag值為3時,計算輸出量咎犬犬大大大大大大大犬大大大大大大大犬犬犬大大大大大犬犬犬大大大大大大犬犬犬犬大大大大大犬犬大大大大functionsys=mcilOutputs(tfxfufn)sys=n*u;$輸出量計算子程序結(jié)束。(2)模塊的封裝與測試。例109采用S函數(shù)來構(gòu)造非線性分段函數(shù)。y=s33_伍_3尸22-(x-5尸1xllx33x44x55x6(1)利用MATLAB語言編寫S函數(shù)。程序如下：functionsys,xOfsts=s

6、function(xrurflag)switchflagfcase0,sysfxOfstrfts=mdllnitializeSizes;case3Zsys=mcilOutputs(tfxfu);case2,4,9 sys=;otherwiseerror(1Unhandiedflag=1rnum2str(flag);endfunctionsys,x0,ts=mdllnitializeSizessizes=simsizes;sizes.NumContStates=0;sizes.NumDiscStates=0;sizes.NumOutputs=1;sizes.NumInputs=1;sizes.DirFeedth匚ough=1;sizes.NumSampleTimes=1;sys=simsizes(sizes);x0=