| 查看: 1168 | 回复: 1 | ||
[求助]
一个椭球旋转问题
|
|
想让椭球长轴指向特定一个方向。 比如：[x,y,z]=ellipsoid(0,0,0,10,1,1)吧，设定一个向量，P=[0.1, 0.5, 0.1]，或者P=[1 1 0]; 想让长轴和P的方向相同。 搞了一下午，转懵了  |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有12人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有19人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有24人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 已知两个坐标系之间x与x',y与y'轴的夹角，怎么求z与z'轴之间的夹角
已经有7人回复
- 关于热源多载荷步施加问题
已经有5人回复
- 电机内磁场问题
已经有3人回复
- 求助大神画晶体的热椭球图，急！
已经有10人回复
- 用diamond画图时出现以下两个问题怎么解决
已经有11人回复
- 求助ansys软件球体周向约束的问题
已经有6人回复
- 求助crystalmaker画晶体图
已经有3人回复
- 求一个截断椭圆的面积和绕轴旋转形成截断椭球的体积。
已经有5人回复
- 请教单晶出图的问题
已经有3人回复
- 这个案例能否采用轴对称进行简化？？
已经有7人回复
- 求助一个关于计算磁阻的定积分
已经有4人回复
- fluent模拟焊接温度场，热源怎么是这样的，像是没出来呢？？？？？
已经有22人回复
- 一个单晶结构的判断，很有意思的
已经有9人回复
- 关于热椭球参数问题
已经有4人回复
- 如何用gambit画椭圆体
已经有8人回复
- 【求助】四分之一波片问题
已经有6人回复
- 求助双椭球热源模型的加载问题
已经有11人回复
- 用diamond沿这个方向作图是什么轴，用英文怎么说？
已经有5人回复
- 晶体的物理性与对称性
已经有232人回复
- Diamond能画出晶体在旋转的GIF图吗？
已经有11人回复
- 结构中的游离水，位于二重轴附近，旋转后二个O原子相距离太近，不知如何处理？
已经有3人回复
- 焊接热源的施加问题？
已经有5人回复
|
根据管理员大人的建议，附上代码，可用于椭球的旋转指向。 [x,y,z] = ellipsoid(0,0,0,10,1,1); %% major axis L and new direction P DIM = 3; L = [1 0 0]; P = [1 1 0]; P = P/norm(P); %% Find rotation matrix t = find(L, 1, 'first'); beta = eye(DIM-1); beta = [beta(1:t-1,  ; zeros(1,DIM-1); beta(t:end, ];[orthBasis, R] = qr([L' beta]); orthBasis = orthBasis * R(1); t = find(P, 1, 'first'); beta = eye(DIM-1); beta = [beta(1:t-1, ; zeros(1,DIM-1); beta(t:end, ];[orthBasis2, R] = qr([P' beta]); orthBasis2 = orthBasis2 * R(1); R = orthBasis/orthBasis2; %% for i = 1:size(x,1) for j = 1:size(x,2) xyz = [x(i,j) y(i,j) z(i,j)] * R; x(i,j) = xyz(1); y(i,j) = xyz(2); z(i,j) = xyz(3); end end figure; hold on mesh(x, y, z) L = [10 0 0] * R; quiver3(0, 0, 0, 1.2*L(1), 1.2*L(2), 1.2*L(3), 0, 'r', 'LineWidth', 3); axis equal grid on view(3) |
2楼2015-06-02 10:03:51
50