| 查看: 2585 | 回复: 6 | ||
[求助]
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文件的仿真图片  我用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;   |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 数据如何到matlab的workspace中？
已经有42人回复
- matlab/simulink
已经有5人回复
- SIMULINK模块间连线如何加粗，添加颜色，请各位大侠指教！
已经有3人回复
- simulink的状态空间模块怎么输入变矩阵
已经有12人回复
- matlab中是否可以用simulink实现遗传算法
已经有6人回复
- 【求助】求助simulink累加模块
已经有5人回复
CSU_HUSTer
新虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2705.1
- 红花: 3
- 帖子: 713
- 在线: 169.4小时
- 虫号: 2836351
- 注册: 2013-11-29
- 性别: GG
- 专业: 工程热物理相关交叉领域
2楼2014-06-17 19:31:56
3楼2014-06-19 12:53:01
CSU_HUSTer
新虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2705.1
- 红花: 3
- 帖子: 713
- 在线: 169.4小时
- 虫号: 2836351
- 注册: 2013-11-29
- 性别: GG
- 专业: 工程热物理相关交叉领域
4楼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;   |
5楼2014-06-19 15:05:55
寒域110
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 584.2
- 散金: 20
- 帖子: 128
- 在线: 43.1小时
- 虫号: 3009529
- 注册: 2014-03-02
- 性别: GG
- 专业: 电力电子学
6楼2015-02-15 15:57:20
7楼2015-12-15 09:09:39
10