| 查看: 1551 | 回复: 6 | ||
[求助]
向量与矩阵的乘积 已有1人参与
|
|
a=[x,y,z]; b=[cos(500)*t,-sin(500)*t,0;sin(300)*t,cos(300)*t,0;0,0,1]; c=[3cos(1000)*t;3sin(1000)*t;0]; d=[9cos(500)*t;9sin(500)*t;20*t]; a=b*c+d 目的是想画出a中x,y,z中的三维图？求助 |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 试验设计与分析及参数优化-（美）C.F.Jeff Wu, Michael Hamada
已经有505人回复
- 非线性分析中的ansys跟踪显示
已经有83人回复
- Fortran90数组学习
已经有1人回复
- 一种简单已用的热加工图画法
已经有22人回复
- 【分享】2011考研数学二大纲
已经有9人回复
- 【分享】10年数2大纲！！
已经有13人回复
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2016-01-27 00:05:34
3楼2016-01-27 09:43:40
wang_chuan
银虫 (正式写手)
- 应助: 64 (初中生)
- 金币: 1786
- 散金: 1900
- 红花: 9
- 帖子: 512
- 在线: 90.7小时
- 虫号: 953632
- 注册: 2010-02-07
- 性别: GG
- 专业: 人工智能与知识工程
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
topliu85: 金币+10, ★★★★★最佳答案, 很好很好。如果我还想把同一时刻对应的x,y做一个筛选，比方说x^2+y^2<=16,怎么添加到程序中？ 2016-01-28 10:01:35
感谢参与，应助指数 +1
topliu85: 金币+10, ★★★★★最佳答案, 很好很好。如果我还想把同一时刻对应的x,y做一个筛选，比方说x^2+y^2<=16,怎么添加到程序中？ 2016-01-28 10:01:35
|
t=0:0.1:3; b=cell(31,1); c=cell(31,1); d=cell(31,1); a=cell(31,1); i=1; for t=0:0.1:3 b{i,1}=[cos(500)*t,-sin(500)*t,0;sin(300)*t,cos(300)*t,0;0,0,1]; c{i,1}=[3*cos(1000)*t;3*sin(1000)*t;0]; d{i,1}=[9*cos(500)*t;9*sin(500)*t;20*t]; a{i,1}=b{i,1}*c{i,1}+d{i,1}; i=i+1; end x=zeros(31,1); y=zeros(31,1); z=zeros(31,1); for i=1:31 x(i,1)=a{i,1}(1,1); y(i,1)=a{i,1}(2,1); z(i,1)=a{i,1}(3,1); end plot3(x,y,z); |
4楼2016-01-27 22:26:14
wang_chuan
银虫 (正式写手)
- 应助: 64 (初中生)
- 金币: 1786
- 散金: 1900
- 红花: 9
- 帖子: 512
- 在线: 90.7小时
- 虫号: 953632
- 注册: 2010-02-07
- 性别: GG
- 专业: 人工智能与知识工程
|
b=cell(31,1); c=cell(31,1); d=cell(31,1); a=cell(31,1); i=1; for t=0:0.1:3 b{i,1}=[cos(500)*t,-sin(500)*t,0;sin(300)*t,cos(300)*t,0;0,0,1]; c{i,1}=[3*cos(1000)*t;3*sin(1000)*t;0]; d{i,1}=[9*cos(500)*t;9*sin(500)*t;20*t]; a{i,1}=b{i,1}*c{i,1}+d{i,1}; i=i+1; end x=zeros(31,1); y=zeros(31,1); z=zeros(31,1); for i=1:31 x(i,1)=a{i,1}(1,1); y(i,1)=a{i,1}(2,1); z(i,1)=a{i,1}(3,1); end x_new=pi; y_new=pi; z_new=pi; for j=1:31 if x(j,1)^2+y(j,1)^2<=16 if x_new(1,1)==pi && y_new(1,1)==pi && z_new(1,1)==pi x_new=x(j,1); y_new=y(j,1); z_new=z(j,1); else x_new=[x_new;x(j,1)]; y_new=[y_new;y(j,1)]; z_new=[z_new;z(j,1)]; end end end plot3(x_new,y_new,z_new); |
5楼2016-01-28 10:55:19
6楼2016-01-28 21:50:19
wang_chuan
银虫 (正式写手)
- 应助: 64 (初中生)
- 金币: 1786
- 散金: 1900
- 红花: 9
- 帖子: 512
- 在线: 90.7小时
- 虫号: 953632
- 注册: 2010-02-07
- 性别: GG
- 专业: 人工智能与知识工程
7楼2016-01-29 08:06:14
20