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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » simulink如何将m文件创建成一个模块,求助呀
71/1返回列表
查看: 2557  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

如瑜?运�

银虫 (初入文坛)

[求助] simulink如何将m文件创建成一个模块,求助呀 已有1人参与

我的m文件原文件如下，我要输出的是一个Fm与时间的关系的一个模块到simulink当中，但是一直都弄不出来，有的说用embedded matlab function，我弄了好久，也没弄出来，我只是想把Fm作为一个激励输入到simulink当中，求高手给个具体的方法和过程，或者帮我编一个s-function函数，我对simulink也不熟，s-function也不会
clc;clear;
syms fs fx;
k=0.6069;
zb=43;
b=11;
rp=72.5;
rrp=5;
E=206000;
i=1;
u=0.3;
kc=0;
P=zeros(1,5);
ks=P;
PP=zeros(1,3);
kk=PP;
ph=PP;
w=0.018;
   for pha=3*2*pi/43:4*pi/43:12*2*pi/43
       s=1+k^2-2*k*cos(pha);
       T=k*(zb+1)*cos(pha)-(1+zb*k^2);
       fs=pi*b*E*rp*s^(3/2);
       fx=4*(1-u^2)*(rp*s^(3/2)+2*rrp*T);
       ks(i)=fs/fx;
       ks(i)
       kc=kc+abs(ks(i));
       kc
       i=i+1;
   end
kp=237190;  %单齿刚度
chd=5.375;  %重合度
nb=600/43;   %摆线齿轮转速
wm=2*pi*nb*zb/60; %频率
t=0:0.001:0.3;
for j=1:1:3
    kk(j)=sqrt(2-2*cos(2*pi*j*(chd-5)))*kp/(pi*j);
    ph(j)=atan((1-cos(2*pi*j*(chd-5)))/sin(2*pi*j*(chd-5)));
    kk(j)
    ph(j)
   
end
    km=kc+kk(1)*cos(wm.*t-ph(1))+kk(2)*cos(2*wm.*t-ph(2))+kk(3)*cos(3*wm.*t-ph(3));
   % plot(t,km)
    Fm=km*w*b/1000;
    plot(t,Fm)
下面这个就是m文件的仿真图片
simulink如何将m文件创建成一个模块,求助呀
我用embedded matlab function编了一个文件，但是出错，下面是我编程和错误的图片，还有simulink的模型
function y = fcn(t)
%#eml
eml.extrinsic('clc','clear','syms');
clc;clear;
%syms fs fx;
k=0.6069;
zb=43;
b=11;
rp=72.5;
rrp=5;
E=206000;
i=1;
u=0.3;
kc=0;
P=zeros(1,5);
ks=P;
PP=zeros(1,3);
kk=PP;
ph=PP;
w=0.018;
   for pha=3*2*pi/43:4*pi/43:12*2*pi/43
       s=1+k^2-2*k*cos(pha);
       T=k*(zb+1)*cos(pha)-(1+zb*k^2);
       fs=pi*b*E*rp*s^(3/2);
       fx=4*(1-u^2)*(rp*s^(3/2)+2*rrp*T);
       ks(i)=fs/fx;
       ks(i)
       kc=kc+abs(ks(i));
       kc
       i=i+1;
   end
kp=237190;  %单齿刚度
chd=5.375;  %重合度
nb=600/43;   %摆线齿轮转速
wm=2*pi*nb*zb/60; %频率
%t=0:0.001:0.3;
for j=1:1:3
    kk(j)=sqrt(2-2*cos(2*pi*j*(chd-5)))*kp/(pi*j);
    ph(j)=atan((1-cos(2*pi*j*(chd-5)))/sin(2*pi*j*(chd-5)));
    kk(j)
    ph(j)
   
end
    km=kc+kk(1)*cos(wm*t-ph(1))+kk(2)*cos(2*wm*t-ph(2))+kk(3)*cos(3*wm*t-ph(3));
   % plot(t,km)
    y=km(t)*w*b/1000;
simulink如何将m文件创建成一个模块,求助呀-1
simulink如何将m文件创建成一个模块,求助呀-2
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-06-17 11:00:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CSU_HUSTer

新虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
matlab FCN即可，设置名字与m文件函数文件名一致
赞一下 回复此楼
2楼2014-06-17 19:31:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

如瑜?运�

银虫 (初入文坛)
引用回帖:
2楼: Originally posted by CSU_HUSTer at 2014-06-17 19:31:56
matlab FCN即可，设置名字与m文件函数文件名一致

好像可以，但是为什么出来的都是好几条平行的直线。就是示波器上的图形显示的是Fm的每个值，每个值是一条直线，没办法跟时间关联起来
赞一下 回复此楼
3楼2014-06-19 12:53:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CSU_HUSTer

新虫 (正式写手)
???????:
3?: Originally posted by ?????? at 2014-06-19 12:53:01
????????????????????????ü?????е???????????????????????????Fm????????????????????????????????????...

???????????д???
赞一下 回复此楼
4楼2014-06-19 14:50:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

如瑜?运�

银虫 (初入文坛)
引用回帖:
4楼: Originally posted by CSU_HUSTer at 2014-06-19 14:50:43
你的输入输出没写对吧...

下面是我的程序，您看看对不，示波器的图显示是那样
function [Fm]=gangdu(t)
syms fs fx;
k=0.6069;
zb=43;
b=11;
rp=72.5;
rrp=5;
E=206000;
i=1;
u=0.3;
kc=0;
P=zeros(1,5);
ks=P;
PP=zeros(1,3);
kk=PP;
ph=PP;
w=0.018;
   for pha=3*2*pi/43:4*pi/43:12*2*pi/43
       s=1+k^2-2*k*cos(pha);
       T=k*(zb+1)*cos(pha)-(1+zb*k^2);
       fs=pi*b*E*rp*s^(3/2);
       fx=4*(1-u^2)*(rp*s^(3/2)+2*rrp*T);
       ks(i)=fs/fx;
       ks(i)
       kc=kc+abs(ks(i));
       kc
       i=i+1;
   end
kp=237190;  %单齿刚度
chd=5.375;  %重合度
nb=600/43;   %摆线齿轮转速
wm=2*pi*nb*zb/60; %频率
%t=0:0.001:0.3;
for j=1:1:3
    kk(j)=sqrt(2-2*cos(2*pi*j*(chd-5)))*kp/(pi*j);
    ph(j)=atan((1-cos(2*pi*j*(chd-5)))/sin(2*pi*j*(chd-5)));
    kk(j)
    ph(j)
   
end
    km=kc+kk(1)*cos(wm*t-ph(1))+kk(2)*cos(2*wm*t-ph(2))+kk(3)*cos(3*wm*t-ph(3));
    plot(t,km)
    Fm=km*w*b/1000;
simulink如何将m文件创建成一个模块,求助呀-3
simulink如何将m文件创建成一个模块,求助呀-4
赞一下 回复此楼
5楼2014-06-19 15:05:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

寒域110

金虫 (小有名气)
楼主，你会用这个模块了吗？求助啊
赞一下 回复此楼
6楼2015-02-15 15:57:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guodeqi0816

新虫 (初入文坛)
楼主您的问题解决了吗   我现在也遇到您这个问题   想把m文件做成类似simulink中的阶跃信号输入模块    可以指点下不
赞一下 回复此楼
7楼2015-12-15 09:09:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 如瑜?运� 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +11 kulium 2026-08-21 14/700 2026-08-24 02:12 by 流流伤
[基金申请] 让我中一个面上吧！ +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 什么时候开奖？ +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭