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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » matlab 旋轮线轨迹模拟 程序错误
北京石油化工学院2026年研究生招生接收调剂公告
31/1返回列表
查看: 813  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

guofucheng07

金虫 (正式写手)

[交流] matlab 旋轮线轨迹模拟 程序错误

从网上下载一个程序，想看看，可惜有错误，大神给修改一下呗，多谢啦

% 旋轮线轨迹模拟
% 圆锥齿轮传动参数
m=2;z7=11;z8=36;
dt7=10.5;dt8=35;
bt=65.25;            % 行星轮轴线与XOY平面夹角
% 搅拌杆外点旋转半径
l=65;
hd=pi/180;
% 圆锥齿轮几何尺寸
r7=0.5*m*z7;
% r7v=r7/cos(dt7*hd);
r8=0.5*m*z8;
% r8v=r8/cos(dt8*hd);
k1=r7+r8;
k2=(r7+r8)/r7;
% 1----二维旋轮线参数计算
for i=1:1:360
    x=k1*cos(i*hd)+l*cos(k2*i*hd);
    y=k1*sin(i*hd)+l*sin(k2*i*hd);
    xl1(i,=[i x y];
end
for i=1:1:720
    x=k1*cos(i*hd)+l*cos(k2*i*hd);
    y=k1*sin(i*hd)+l*sin(k2*i*hd);
    xl2(i,=[i x y];
end
for i=1:1:1080
    x=k1*cos(i*hd)+l*cos(k2*i*hd);
    y=k1*sin(i*hd)+l*sin(k2*i*hd);
    xl3(i,=[i x y];
end
for i=1:1:1440
    x=k1*cos(i*hd)+l*cos(k2*i*hd);
    y=k1*sin(i*hd)+l*sin(k2*i*hd);
    xl4(i,=[i x y];
end
% 绘制二维旋轮线
figure(1);                % 生成第1个图形窗口
subplot(2,2,1);           % 选择第1个子窗口
plot(xl1(:,2),xl1(:,3))
grid                      % 绘制网格线
title('二维旋轮线(\theta =360度)')
subplot(2,2,2);           % 选择第2个子窗口
plot(xl2(:,2),xl2(:,3))
grid                      % 绘制网格线
title('二维旋轮线(\theta =720度)')
subplot(2,2,3);           % 选择第3个子窗口
plot(xl3(:,2),xl3(:,3))
grid                      % 绘制网格线
title('二维旋轮线(\theta =1080度)')
subplot(2,2,4);           % 选择第4个子窗口
plot(xl4(:,2),xl4(:,3))
grid                      % 绘制网格线
title('二维旋轮线(\theta =1440度)')
% 2----三维旋轮线参数计算
for i=1:1:360
    x=(k1*cos(i*hd)+l*cos(k2*i*hd))*cos(bt*hd);
    y=(k1*sin(i*hd)+l*sin(k2*i*hd))*cos(bt*hd);
    z=sqrt(x^2+y^2)*tan(bt*hd);
    xlx1(i,=[i x y z];
end
for i=1:1:720
    x=(k1*cos(i*hd)+l*cos(k2*i*hd))*cos(bt*hd);
    y=(k1*sin(i*hd)+l*sin(k2*i*hd))*cos(bt*hd);
    z=sqrt(x^2+y^2)*tan(bt*hd);
    xlx2(i,=[i x y z];
end
for i=1:1:1080
    x=(k1*cos(i*hd)+l*cos(k2*i*hd))*cos(bt*hd);
    y=(k1*sin(i*hd)+l*sin(k2*i*hd))*cos(bt*hd);
    z=sqrt(x^2+y^2)*tan(bt*hd);
    xlx3(i,=[i x y z];
end
for i=1:1:1440
    x=(k1*cos(i*hd)+l*cos(k2*i*hd))*cos(bt*hd);
    y=(k1*sin(i*hd)+l*sin(k2*i*hd))*cos(bt*hd);
    z=sqrt(x^2+y^2)*tan(bt*hd);
    xlx4(i,=[i x y z];
end
% 矩阵cs(:,j)表示第j列的各行元素
% 绘制三维旋轮线
figure(2);                % 生成第2个图形窗口
subplot(2,2,1);           % 选择第1个子窗口
plot3(xlx1(:,2),xlx1(:,3),xlx1(:,4))
grid
title('三维旋轮线(\theta =360度)')
subplot(2,2,2);           % 选择第2个子窗口
plot3(xlx2(:,2),xlx2(:,3),xlx2(:,4))
grid
title('三维旋轮线(\theta =720度)')
subplot(2,2,3);           % 选择第3个子窗口
plot3(xlx3(:,2),xlx3(:,3),xlx3(:,4))
grid
title('三维旋轮线(\theta =1080度)')
subplot(2,2,4);           % 选择第4个子窗口
plot3(xlx4(:,2),xlx4(:,3),xlx4(:,4))
grid
title('三维旋轮线(\theta =1440度)')
% 输出数据
disp ' '
disp '          ==========   圆锥齿轮传动参数   ==========';
fprintf(1,'                   行星轮齿数    z7 = %3.0f \n',z7);
fprintf(1,'                   分度圆半径    r7 = %3.3f 度 \n',r7);
fprintf(1,'                   分度圆锥角   dt7 = %3.3f 度 \n',dt7);
fprintf(1,'                   中心轮齿数    z8 = %3.0f \n',z8);
fprintf(1,'                   分度圆半径    r8 = %3.3f 度 \n',r8);
fprintf(1,'                   分度圆锥角   dt8 = %3.3f 度 \n',dt8);
fprintf(1,'   中心轮与行星轮当量半径之比    rb = %3.0f \n',r8/r7);
disp ' '
['  行星轮公转',' 动点x坐标',' 动点y坐标',' 动点z坐标']
[xlx1(:,1),xlx1(:,2),xlx1(:,3),xlx1(:,4)]
回复此楼

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2014-01-03 17:25:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guofucheng07

金虫 (正式写手)
程序m代码
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 07-搅拌机工作头旋轮轨迹分析与仿真-1.m
    • 2014-01-03 17:28:25, 3.35 K
2楼2014-01-03 17:28:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
guofucheng073楼
2014-01-25 10:04   回复  
引用回帖:
2楼: Originally posted by guofucheng07 at 2014-01-03 17:28:44 程序m代码

谢谢分享
相关版块跳转 我要订阅楼主 guofucheng07 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 278求调剂 +4 范婷娜 2026-04-04 5/250 2026-04-04 08:45 by 来看流星雨10
[考研] 081700化学工程与技术 一志愿中海洋 323 求调剂学校 +11 披星河 2026-04-03 11/550 2026-04-04 08:30 by leqianduo
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 调剂 +5 asdasdassda 2026-04-03 6/300 2026-04-03 20:27 by 岸上的一条鱼
[考研] 321求调剂 +10 认真求上学 2026-04-02 10/500 2026-04-03 16:17 by lijunpoly
[考研] 320求调剂 +5 振—TZ 2026-04-02 5/250 2026-04-03 14:42 by fxue1114
[考研] 0854求调剂 +3 assdll 2026-04-03 3/150 2026-04-03 14:32 by fxue1114
[考研] 085600专硕材料与化工348分求调剂 +10 上学啦！ 2026-04-01 11/550 2026-04-03 14:13 by 百灵童888
[考研] 建环，能源，土木老师路过看一看！！！ +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[考研] 320求调剂 +3 农业工程与信息� 2026-04-03 3/150 2026-04-03 11:40 by 土木硕士招生
[考研] 326求调剂 +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 0703求调剂 +4 zizimo 2026-03-31 4/200 2026-04-01 16:04 by yanflower7133
[考研] 283求调剂 +9 A child 2026-03-28 9/450 2026-04-01 14:20 by Jaylen.
[考研] 085600，321分求调剂 +13 大馋小子 2026-03-31 13/650 2026-04-01 12:35 by chemdavid
[考研] 一志愿西电085401数一英一299求调剂 六级521 +4 爱吃大鸭梨 2026-03-31 4/200 2026-03-31 11:51 by 搏击518
[考研] 08工科求调剂286 +5 tgs_001 2026-03-28 5/250 2026-03-31 08:18 by 一只好果子?
[考研] 福建理工大学材料学院先进合金团队招收考研调剂学生 +3 大华金商都 2026-03-30 4/200 2026-03-31 01:04 by 方英俊602
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3 kaobao456 2026-03-29 4/200 2026-03-30 23:16 by nBu锂
[考研] 086000生物与医药调剂 +5 Feisty。 2026-03-28 9/450 2026-03-29 12:02 by longlotian
信息提示
关闭
请填处理意见
关闭