s函数的问题.谁能帮我把这段程序变成s函数的形式啊,以t为输入,y为输出,clcclearRad90=deg2rad(90);%度数到弧度的转换.1弧度=180/pi度,即1rad=57.29°根据文献51制作出跟定的桨矩角曲线Rad60=deg2ra

Question

s函数的问题.
谁能帮我把这段程序变成s函数的形式啊,以t为输入,y为输出,
clc
clear
Rad_90=deg2rad(90); %度数到弧度的转换.1弧度=180/pi 度,即1rad=57.29° 根据文献【51】制作出跟定的桨矩角曲线
Rad_60=deg2rad(60);
Rad_50=deg2rad(50);
Rad_40=deg2rad(40);
Rad_30=deg2rad(30);
Rad_20=deg2rad(20);
Rad_10=deg2rad(10);
y=[];
for t=0:36
if 0t
k=((Rad_90-0)/(8-0));
x0=8;y0=Rad_90;
b=y0-k*x0;
y(t)=k*t+b;
elseif 8=t
y(t)=Rad_90;
elseif 10t
k=((Rad_60-Rad_90)/(13-10));
x0=13;y0=Rad_60;
b=y0-k*x0;
y(t)=k*t+b;
elseif 13=t
y(t)=Rad_60;
elseif 15t
k=((Rad_40-Rad_60)/(17-15));
x0=17;y0=Rad_40;
b=y0-k*x0;
y(t)=k*t+b;
elseif 17=t
y(t)=Rad_40;
elseif 19t
k=((Rad_50-Rad_40)/(21-19));
x0=21;y0=Rad_50;
b=y0-k*x0;
y(t)=k*t+b;
elseif 21=t
y(t)=Rad_50;
elseif 23t
k=((Rad_30-Rad_50)/(25-23));
x0=25;y0=Rad_30;
b=y0-k*x0;
y(t)=k*t+b;
elseif 25=t
y(t)=Rad_30;
elseif 27t
k=((Rad_10-Rad_30)/(29-27));
x0=29;y0=Rad_10;
b=y0-k*x0;
y(t)=k*t+b;
elseif 29=t
y(t)=Rad_10;
elseif 31t
k=((Rad_90-Rad_10)/(34-31));
x0=34;y0=Rad_90;
b=y0-k*x0;
y(t)=k*t+b;
elseif 34

Answer 1

这个很简单啊,为什么自己不试一试呢
找个S 函数模板
修改输出函数的代码,输入维数0,输出维数1,状态维数0
function sys=mdlOutput(t,x,u)
Rad_90=deg2rad(90); %度数到弧度的转换.1弧度=180/pi 度,即1rad=57.29° 根据文献【51】制作出跟定的桨矩角曲线
Rad_60=deg2rad(60);
Rad_50=deg2rad(50);
Rad_40=deg2rad(40);
Rad_30=deg2rad(30);
Rad_20=deg2rad(20);
Rad_10=deg2rad(10);
if 0t
k=((Rad_90-0)/(8-0));
x0=8;y0=Rad_90;
b=y0-k*x0;
sys=k*t+b;
elseif 8=t
sys=Rad_90;
elseif 10t
k=((Rad_60-Rad_90)/(13-10));
x0=13;y0=Rad_60;
b=y0-k*x0;
sys=k*t+b;
elseif 13=t
sys=Rad_60;
elseif 15t
k=((Rad_40-Rad_60)/(17-15));
x0=17;y0=Rad_40;
b=y0-k*x0;
sys=k*t+b;
elseif 17=t
sys=Rad_40;
elseif 19t
k=((Rad_50-Rad_40)/(21-19));
x0=21;y0=Rad_50;
b=y0-k*x0;
sys=k*t+b;
elseif 21=t
sys=Rad_50;
elseif 23t
k=((Rad_30-Rad_50)/(25-23));
x0=25;y0=Rad_30;
b=y0-k*x0;
sys=k*t+b;
elseif 25=t
sys=Rad_30;
elseif 27t
k=((Rad_10-Rad_30)/(29-27));
x0=29;y0=Rad_10;
b=y0-k*x0;
sys=k*t+b;
elseif 29=t
sys=Rad_10;
elseif 31t
k=((Rad_90-Rad_10)/(34-31));
x0=34;y0=Rad_90;
b=y0-k*x0;
sys=k*t+b;
elseif 34