2013-10-24

DirFeedthrough 什么意？

s-function也就是system-function的縮寫。說得簡單，s-function就是用MATLAB所提供的模型不能完全滿足用戶，而提供給用戶自己編寫程序來滿足自己要求模型的接口。 要了解 s-function，必須了解以下知識：

(1)direct feedthrough

(2)dynamically sized inputs

(3)setting sample times and offsets

一。 direct feedthrough

direct feedthrough意思是說系統(tǒng)的輸出或可變采樣時間是否受到輸入的控制。有的系統(tǒng)是受到輸入控制的,如：

y=k*u...全部

s-function也就是system-function的縮寫。說得簡單，s-function就是用MATLAB所提供的模型不能完全滿足用戶，而提供給用戶自己編寫程序來滿足自己要求模型的接口。

要了解 s-function，必須了解以下知識：

(1)direct feedthrough

(2)dynamically sized inputs

(3)setting sample times and offsets

一。

direct feedthrough

direct feedthrough意思是說系統(tǒng)的輸出或可變采樣時間是否受到輸入的控制。有的系統(tǒng)是受到輸入控制的,如：

y=k*u (u是輸入，k是放大因子，y是輸出)

而有的系統(tǒng)輸出是不受到輸入影響，如：

輸出：y=x

dx=u

x表示狀態(tài)

二。

dynamically sized inputs

主要是給出：輸入連續(xù)狀態(tài)數目(size。NumContStates),離散狀態(tài)數目(size。NumDiscStates) ,輸出數目(size。

NumOutputs),輸入數目(size。NumInputs),Direct Feedthrough(size。Dir Feedthrough)。

三。setting sample times and offsets

setting smaple times and offsets主要設置采樣時間。

MATLAB為了用戶方便，已經書寫了S-FUNCTIONS模板函數sfuntmpl。m。為了更好的寫S-Function,該函數sfuntmpl。m內容如下面所示,該函數可以通過which sfuntmpl。

m來查看它的位置。一般是在\toolbox\simulink\blocks\文件夾中,然后通過輸入命令edit sfuntmpl。m來查看M代碼,當然也可以使用type函數。

事例一:連續(xù)系統(tǒng)例子：

function [sys,x0,str,ts] = csfunc(t,x,u,flag)

A=[-0。

09 -0。01; 1 0];

B=[ 1 -7; 0 -2];

C=[ 0 2; 1 -5];

D=[-3 0; 1 0];

switch flag,

case 0

[sys,x0,str,ts]=mdlInitializeSizes(A,B,C,D);

case 1

sys=mdlDerivatives(t,x,u,A,B,C,D);

case 3

sys=mdlOutputs(t,x,u,A,B,C,D);

case { 2, 4, 9 }

sys = [];

otherwise

error(['Unhandled flag = ',num2str(flag)]);

end

function [sys,x0,str,ts]=mdlInitializeSizes(A,B,C,D)

sizes = simsizes;

sizes。

NumContStates = 2; sizes。NumDiscStates = 0; sizes。NumOutputs = 2;

sizes。NumInputs = 2; sizes。DirFeedthrough = 1; sizes。

NumSampleTimes = 1;

sys = simsizes(sizes); x0 = zeros(2,1); str = []; ts = [0 0];

function sys=mdlDerivatives(t,x,u,A,B,C,D)

sys = A*x + B*u; % end mdlDerivatives

function sys=mdlOutputs(t,x,u,A,B,C,D)

sys = C*x + D*u; % end mdlOutputs

事例二:離散系統(tǒng)

function [sys,x0,str,ts] = dsfunc(t,x,u,flag)

A=[-1。

3839 -0。5097 1。0000 0];

B=[-2。5559 0 0 4。2382];

C=[ 0 2。0761 0 7。7891];

D=[ -0。8141 -2。9334 1。

2426 0];

switch flag,

case 0

[sys,x0,str,ts] = mdlInitializeSizes(A,B,C,D);

case 2

sys = mdlUpdate(t,x,u,A,B,C,D); %%%%%%%%%%

case 3

sys = mdlOutputs(t,x,u,A,C,D); %%%%%%%%%%%%%

case 9

sys = [];

otherwise

error(['unhandled flag = ',num2str(flag)]);

end %end dsfunc

function [sys,x0,str,ts] = mdlInitializeSizes(A,B,C,D)

sizes = simsizes;

sizes。

NumContStates = 0; sizes。NumDiscStates = size(A,1); sizes。NumOutputs = size(D,1);

sizes。NumInputs = size(D,2); sizes。

DirFeedthrough = 1; sizes。NumSampleTimes = 1;

sys = simsizes(sizes); x0 = ones(sizes。NumDiscStates,1); str = []; ts = [1 0];

function sys = mdlUpdate(t,x,u,A,B,C,D)

sys = A*x+B*u; %end mdlUpdate

function sys = mdlOutputs(t,x,u,A,C,D)

sys = C*x+D*u; %end mdlUpdate

事例三:離散和連續(xù)的混合型

function [sys,x0,str,ts] = mixedm(t,x,u,flag)

dperiod = 1; doffset = 0;

switch flag

case 0

[sys,x0,str,ts]=mdlInitializeSizes(dperiod,doffset);

case 1

sys=mdlDerivatives(t,x,u);

case 2

sys=mdlUpdate(t,x,u,dperiod,doffset);

case 3

sys=mdlOutputs(t,x,u,doffset,dperiod);

case 9

sys = [];

otherwise

error(['unhandled flag = ',num2str(flag)]);

end % end mixedm %

function [sys,x0,str,ts]=mdlInitializeSizes(dperiod,doffset)

sizes = simsizes; sizes。

NumContStates = 1; sizes。NumDiscStates = 1;

sizes。NumOutputs = 1; sizes。NumInputs = 1; sizes。DirFeedthrough = 0;

sizes。

NumSampleTimes = 2; sys = simsizes(sizes); x0 = ones(2,1); str = [];

ts = [0 0;dperiod doffset];

function sys=mdlDerivatives(t,x,u)

sys = u; % end mdlDerivatives

function sys=mdlUpdate(t,x,u,dperiod,doffset)

if abs(round((t - doffset)/dperiod) - (t - doffset)/dperiod) < 1e-8, sys = x(1);

else sys = []; end % end mdlUpdate

function sys=mdlOutputs(t,x,u,doffset,dperiod)

if abs(round((t - doffset)/dperiod) - (t - doffset)/dperiod) < 1e-8, sys = x(2); else sys = []; end

本帖中對是連續(xù)，離散，連續(xù)與離散混合系統(tǒng)的三個例子加以說明，以至于大家在看下面的例子時能更好的理解。

一。函數mdlInitializeSizes

對于描述連續(xù)型的函數csfunc。m

sizes。NumContStates = 2; sizes。NumDiscStates = 0;

表明本函數是描述連續(xù)型的

sizes。

NumOutputs = 2; sizes。NumInputs = 2;

sizes。DirFeedthrough = 1; sizes。NumSampleTimes = 1;

對于描述離散型的函數dsfunc。

m

sizes。NumContStates = 0; sizes。NumDiscStates = size(A,1);

表明本函數是描述離散型的

sizes。NumOutputs = size(D,1); sizes。

NumInputs = size(D,2);

sizes。DirFeedthrough = 1; sizes。NumSampleTimes = 1;

對于描述混合型的函數mixedm。m

sizes。

NumContStates = 1; sizes。NumDiscStates = 1;

表明本函數是描述混合型的

sizes。NumOutputs = 1; sizes。NumInputs = 1;

sizes。

DirFeedthrough = 0; sizes。NumSampleTimes = 2;

關于變量sizes。numoutputs和sizes。numinputs則由描述的系統(tǒng)的輸出，輸入而定。

關于變量sizes。DirFeedthrough,則可以查看矩陣D,如果D是非空，則為1，否則為0

關于變量sizes。NumSampleTimes，即是ts這個矩陣的行數。

二。函數mdlDerivatives

僅僅出現在含有連續(xù)型的系統(tǒng)中；這部分應該根據系統(tǒng)狀態(tài)來決定。

三。函數mdlUpdate 僅僅出現在含有離散型的系統(tǒng)中；該函數應該根據系統(tǒng)狀態(tài)來決定

四。函數mdlOutputs 該函數在各種類型中都有,應該根據描述系統(tǒng)的輸出來決定;

五。關于書寫s-functions函數的注意之點

(1)我想大家都會覺得上述三個函數都不利于維護,因此在寫s-functions函數時不妨把矩陣 A,B,C,D也考慮成變量,作為輸入參數;

(2)關于在csfuncs。

m中的型如case(2,4,9)這種形式,大家最好不要采用,不利于維護,可以分開來寫

(3)關于結構sizes中的各個變量的值,最好從輸入參量中得到,而不是通過人為判斷來輸入

(4)因此,其實上述三種形式,可以變化成一種形式。

(5)可以通過對相應于你所要求的系統(tǒng)對上述三個函數加以相應的改進來到達自己的要求,因此可以把上述函數當成模板函數。

。收起

總結

以上是生活随笔為你收集整理的matlab中(1 )什么意思,matlab中area(1)什么意？mat – 手机爱问的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。