matlab中调用simulink,M文件中如何调用SIMULINK

%% Fuzzy Tunning PID Control

clear all;

close all;

a=newfis('fuzzpid');

f1=800;

a=addvar(a,'input','e',[-6*f1,6*f1]); % Parameter e

a=addmf(a,'input',1,'NB','trimf',[-6*f1,-6*f1,-4*f1]);

a=addmf(a,'input',1,'NM','trimf',[-6*f1,-4*f1,-2*f1]);

a=addmf(a,'input',1,'NS','trimf',[-4*f1,-2*f1,0]);

a=addmf(a,'input',1,'Z','trimf',[-2*f1,0,2*f1]);

a=addmf(a,'input',1,'PS','trimf',[0*f1,2*f1,4*f1]);

a=addmf(a,'input',1,'PM','trimf',[2*f1,4*f1,6*f1]);

a=addmf(a,'input',1,'PB','trimf',[4*f1,6*f1,6*f1]);

f2=800;

a=addvar(a,'input','ec',[-6*f2,6*f2]);% Parameter ec

a=addmf(a,'input',2,'NB','trimf',[-6*f2,-6*f2,-4*f2]);

a=addmf(a,'input',2,'NM','trimf',[-6*f2,-4*f2,-2*f2]);

a=addmf(a,'input',2,'NS','trimf',[-4*f2,-2*f2,0]);

a=addmf(a,'input',2,'Z','trimf',[-2*f2,0,2*f2]);

a=addmf(a,'input',2,'PS','trimf',[0,2*f2,4*f2]);

a=addmf(a,'input',2,'PM','trimf',[2*f2,4*f2,6*f2]);

a=addmf(a,'input',2,'PB','trimf',[4*f2,6*f2,6*f2]);

f3=1;

a=addvar(a,'output','kp',[-