ch2example17sfun.m

通信程序源代码和模型通信程序源代码和模型通信程序源代码和模型

function [sys,x0,str,ts] =ch2example17Sfun(t,x,u,flag,Amp,Freq,Phase) 
% 正弦波信号源
switch flag, 
    case 0       % flag=0 初始化 
       [sys,x0,str,ts]=mdlInitializeSizes; 
    case 3       % flag=3 计算输出 
       sys=mdlOutputs(t,Amp,Freq,Phase); 
    case {1, 2, 4, 9 }    % 其他作不处理的flag 
       sys=[];            % 无用的flag时返回sys为空矩阵
otherwise % 异常处理 
       error(['Unhandled flag = ',num2str(flag)]); 
end 
% 主函数结束 
% 子函数实现（1）初始化函数----------------------------------
function [sys,x0,str,ts] = mdlInitializeSizes  % 
sizes = simsizes;           % 获取SIMULINK仿真变量结构
sizes.NumContStates  = 0;   % 连续系统的状态数为0
sizes.NumDiscStates  = 0;   % 离散系统的状态数为0
sizes.NumOutputs     = 1;   % 输出信号数目是1 
sizes.NumInputs      = 0;   % 输入信号数目是0
sizes.DirFeedthrough = 0;   % 该系统不是直通的 
sizes.NumSampleTimes = 1;   % 这里必须为1 
sys = simsizes(sizes); 
str = [];           % 通常为空矩阵 
x0  = [];           % 初始状态矩阵x0 （零状态情况） 
ts  = [0 0];        % 表示连续取样时间的仿真
% 初始化函数结束 

% 子函数实现（2）系统输出方程函数-----------------------------
function sys = mdlOutputs(t,Amp,Freq,Phase) 
sys = Amp*sin(2*pi*Freq*t+Phase);% 这里写入系统的输出方程矩阵形式即可
% 修改这个函数可以得到任意的波形输出
% 系统输出方程函数结束