s_function.m

simulink仿真噪声模型

function [sys,x0,str,ts] = s_function(t,x,u,flag)
global counter in_put w step out_put hs_size hs L r;
switch flag,
case 0,
    [sys,x0,str,ts]=mdlInitializeSizes;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
case 1,
     sys=0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 case 3,
     step=u(2);
     counter=counter+1;
     in_put=[u(1);in_put(1:(hs_size-1))];
     x=in_put'*hs;
     r=[x;r(1:(L-1))];
     w=w-step*u(3)*r;
     sys=w'*r;
    
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
case { 2, 4, 9 },
    sys = [];
otherwise
    error(['Unhandled flag = ',num2str(flag)]);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [sys,x0,str,ts]=mdlInitializeSizes
global counter in_put w step out_put hs_size hs L r;
L=300;           %表示滤波器的长度
in_put=zeros(hs_size,1);
out_put=zeros(L,1);
w=zeros(L,1);         %L表示滤波器的权系数
r=zeros(L,1);


hs_size=8;       %次级声道的传递函数长度
hs=[0.9321 0.6500 -0.2000 0.3900 0.6277 0.4140 0.1040 -0.1000]' %次级通道数学模型
%hs=zeros((hs_size-1),1);hs=[hs;1];  %次级通道数学模型


sizes = simsizes;
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 1;
sizes.NumInputs      = 3;
sizes.DirFeedthrough = 1;
sizes.NumSampleTimes = 1;   % at least one sample time is needed
sys = simsizes(sizes);
x0  = [];
str = [];
ts  = [0 0];