24小时热门版块排行榜    

查看: 866  |  回复: 2

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)]
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注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的回帖
简单回复
2014-01-25 10:04   回复  
引用回帖:
2楼: Originally posted by guofucheng07 at 2014-01-03 17:28:44 程序m代码

谢谢分享
相关版块跳转 我要订阅楼主 guofucheng07 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +5 a089 2026-08-07 5/250 2026-08-08 14:04 by alian_214
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 关于代码变化问题,想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见