accumulate222.m

来自「相位累加器」· M 代码 · 共 34 行

function [sys,x0,str,ts] =accumulate(t,x,u,flag,sum)
switch flag,
  case 0,
    [sys,x0,str,ts]=mdlInitializeSizes(sum);
  case 2,
    sys=mdlUpdate(t,x,u,sum);
  case 3,
    sys=mdlOutputs(t,x,u,sum);
  case {1,4,9}
    sys=[ ];
    otherwise
    error(['Unhandled flag = ',num2str(flag)]);
end

function [sys,x0,str,ts]=mdlInitializeSizes(sum)
sizes = simsizes;
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 1;
sizes.NumOutputs     = 1;
sizes.NumInputs      =1;%输入U应该为频率控制字
sizes.DirFeedthrough = 1;
sizes.NumSampleTimes = 1;  
sys = simsizes(sizes);
x0  = [0];%设定初始值，根据需要可以随意修改
str = [ ];
ts  = [1/(1000*64) 0];%1/100是采样间隔，根据需要可以随意修改       

function sys=mdlUpdate(t,x,u,sum)
N=32;%累加器位数，用来控制频率分辨率
sum=x+((2*pi)/(2^N))*u;
sys = sum;
function sys=mdlOutputs(t,x,u,sum)
sys=x;