24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 528  |  回复: 1

小琪、

金虫 (小有名气)

[求助] S函数运行出错 求指导

这是我的一个式子 我编写的S函数是这样的:
function [sys,x0,str,ts] =Fz(t,x,u,flag)
switch flag,

  case 0,
    [sys,x0,str,ts]=mdlInitializeSizes(R);
    beta=345.1656;
    amp=0.8*10^(-6);
    i=beta*r;
    R=2.537*besselj(9,i)+0.43472*bessely(9,i)+3.3257*10^(-3)*besseli(9,i)+4.4643*besselk(9,i);

  case 1,
    sys=mdlDerivatives(t,x,u);

  case 3,
    sys=mdlOutputs(t,x,u,R,r);

  case 4,
    sys=[];

  case 9,
    sys=[];

  otherwise
    DAStudio.error(\'Simulink:blocks:unhandledFlag\', num2str(flag));

end

function [sys,x0,str,ts]=mdlInitializeSizes(R)

sizes = simsizes;

sizes.NumContStates  = 0;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 1;
sizes.NumInputs      = 1;
sizes.DirFeedthrough = 1;
sizes.NumSampleTimes = 1;   % at least one sample time is needed

sys = simsizes(sizes);

x0  = [];
str = [];
ts  = [0 0];


function sys=mdlOutputs(t,x,u,R,r)
if u>0
    h=sin(arccos((u-amp*R)/amp/R))+(1-u/amp/R)*arccos((u-amp*R)/R/amp)
else
    h=9*pi*(1-u/R/amp)
end
sys =3.68*10^6*int(h*R*r,r,0.0225,0.03);

function sys=mdlTerminate(t,x,u)

sys =9*pi*(1-u/R/amp);
  下面是我的simulink模型  运行的时候老是出错  式子的R表达式中J,N,I,K 都属于matlab中的贝塞尔函数 可以直接在matlab中表示出来  希望大家帮忙给看一下 谢谢大家了

S函数运行出错 求指导
模拟的式子


S函数运行出错 求指导-1
simulink模型


S函数运行出错 求指导-2
错误
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)

你这个思路应该就是有问题的
那个定积分可以考虑用数值方法计算
bessel函数没想明白怎么处理合适

其实完全不需要用s-function的
2楼2015-04-28 15:55:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小琪、 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 310求调剂 +18 争取九点睡 2026-03-30 18/900 2026-04-03 18:35 by ls刘帅
[考研] 334求调剂 +7 曾仰之 2026-04-03 7/350 2026-04-03 18:21 by ls刘帅
[考研] 285求调剂 +5 AZMK 2026-04-03 8/400 2026-04-03 18:17 by AZMK
[考研] 288求调剂 一志愿哈工大 材料与化工 +38 洛神哥哥 2026-03-31 40/2000 2026-04-03 18:00 by Jimmyandyou
[考研] 一志愿0817化学工程与技术,求调剂 +17 我不是只因 2026-04-02 19/950 2026-04-03 17:08 by lijunpoly
[考研] 274求调剂 +9 顺理成张 2026-04-03 10/500 2026-04-03 15:10 by 啊俊!
[考研] 326分求调剂 +3 于是乎呢 2026-04-01 5/250 2026-04-03 14:23 by 于是乎呢
[考研] 数一英一285求调剂 +7 AZMK 2026-04-03 9/450 2026-04-03 13:03 by ms629
[考研] 319求调剂 +18 太容易1018 2026-04-01 18/900 2026-04-03 11:18 by linyelide
[考研] 295求调剂 +7 愿旅途永远坦然 2026-04-02 7/350 2026-04-03 08:22 by fangshan711
[考研] 土木304求调剂 +4 兔突突突, 2026-04-02 5/250 2026-04-02 21:16 by 兔突突突,
[考研] +4 雾与海 2026-04-02 5/250 2026-04-02 19:16 by 土木硕士招生
[考研] 一志愿华南师范大学-22408计算机-292分-求华南师范大学调剂 +4 爱读书的小鳄鱼 2026-04-02 4/200 2026-04-02 18:35 by 求调剂zz
[考研] 362求调剂 +14 西南交材料专硕3 2026-03-31 14/700 2026-04-02 17:50 by yunlongyang
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 求调剂,一志愿南京师范大学计算机专硕,初试373,六级通过, +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
[考研] 290求调剂 +5 dfffsar 2026-03-29 5/250 2026-04-01 19:45 by 6781022
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
[考研] 生物与医药考研调剂 +5 铁憨憨123425 2026-03-31 5/250 2026-04-01 18:01 by syh9288
[考博] 材料专业申博 +5 杜雨婷dyt 2026-03-29 5/250 2026-03-31 11:19 by oooqiao
信息提示
请填处理意见