| 查看: 866 | 回复: 2 | |||
[交流]
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)] |
回复此楼» 猜你喜欢
十年后又回来了，论文投稿求助
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
fileCode有新解读？
已经有17人回复
-
这样的filecode谁见过
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在matlab中做一点的空间轨迹的三维图，程序报错
已经有18人回复
- MATLAB关于lsqcurvefit程序几个错误，求大神指导一下
已经有6人回复
- matlab程序错误
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
山东征女友，坐标济南
+1/160
-
一边期待，一边忐忑：发榜前的科研人正在想什么？
+5/120
-
别研究filecode了，祝都有一个好结果：）
+2/112
-
三(三甲基硅基)甘油醚（CAS：6787-10-6 ）
+1/88
-
3-丁烯-1-醇（CAS：627-27-0）
+1/87
-
三(三甲基硅基)硼酸酯（CAS：4325-85-3）
+1/87
-
1,3,2-Dioxaphospholane, 2-(2,2,3,3,3-pentafluoropropoxy)-
+1/86
-
丁烯二醇（CAS: 497-06-3）
+1/86
-
哈工大（深圳）何自开团队诚招2027级化学工程与技术专业博士生
+1/28
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/26
-
中国石油大学（华东）计算机科学与技术学院/智能化新能源系统设计和开发方向 特聘副教
+2/14
-
双一流高校招收博士（现代技术管理/环境技术经济管理）及硕士研究生（长期有效）
+1/13
-
瑞典林雪平大学 高峰课题组 诚聘 氮化物基无机材料方向 博士后1名
+1/7
-
福建师范大学海峡柔性电子学院-2027年9月入学博士-考核制-电催、能源
+1/7
-
希望大家都高中！！
+1/7
-
北理工集成电路杰青团队 | 诚招科助理
+1/4
-
氨氮空白偏高！求指教！
+1/4
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/3
-
北理工集成电路杰青团队 | 诚招科助理
+1/3
-
27届申博求助
+1/1
| 程序m代码 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 07-搅拌机工作头旋轮轨迹分析与仿真-1.m
2014-01-03 17:28:25, 3.35 K
2楼2014-01-03 17:28:44
简单回复
2014-01-25 10:04
回复
谢谢分享
