版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

小琪、

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 776.9
红花: 2
帖子: 99
在线: 35.6小时
虫号: 3626547
注册: 2015-01-04
性别: MM
专业: 制造系统与自动化

[求助] 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中表示出来  希望大家帮忙给看一下谢谢大家了

模拟的式子

simulink模型

错误

回复此楼

» 猜你喜欢

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

求助关于狄拉克函数（Dirac Delta function）的一个问题已经有13人回复
运行函数文件出错，Input argument "y" is undefined. 已经有5人回复
matlab求函数积分出错已经有12人回复
【已完结】BP神经网络程序问题已经有3人回复
求各位大神帮我看下一个小小的程序，只需要你一点儿时间。就会帮我很大的忙已经有15人回复
in文件运行失败（lost atoms）求助！！已经有11人回复
ANSYS FLUENT 12.1.2 装完运行FLUENT出错已经有12人回复
如何在matlab里写.m函数，老是出错有谁知道具体格式不已经有4人回复
使用matlab 中的fsolve 老提示出错已经有4人回复
matlab 问题求助，一段程序出错，求纠错已经有7人回复
matable ,quad函数的调用出错。。。已经有5人回复
matlab程序运行出错，求指点修改已经有4人回复
利用matlab中的S-函数实现积分分离PID控制器的仿真已经有4人回复
两道有关固液体系相变，要设计过程求G,H,S状态函数变化的题目已经有23人回复
【求助】IRC出错 Convergence failure -- run terminated. 已经有32人回复
【求助】BSSE of O2运行出错已经有17人回复
【求助】matlab如何求解sinc函数的反函数已经有5人回复
【求助】有很多错误的程序，找不出来了已经有18人回复
【求助】LAMMPS中非周期性边界条件和Born-Mayer-Huggins势函数不会设置已经有4人回复

1楼 2015-04-28 15:21:46

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

信彼南山

木虫 (著名写手)

应助: 33 (小学生)
金币: 4142.9
散金: 1221
红花: 16
帖子: 1178
在线: 233.5小时
虫号: 1133529
注册: 2010-10-27
专业: 导航、制导与传感技术

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

其实完全不需要用s-function的

赞一下(2人)