如何从极坐标图绘制3 d图吗

8视图(30天)
大卫餐桌
大卫餐桌 2022年8月9日
评论道: 威廉•罗斯 2022年8月17日
我看着策划Titinium缓慢的概要文件的速度(逆)解决christoffell方程。
我已经设法得到一个极坐标图的但是想想象这个3 d。我为我极坐标图完成。
c11 = 162000000000;
c12 = 92000000000;
c44 = 47000000000;
c13 = 69000000000;
c33 = 181000000000;
c66 = 35000000000
ρ= 4430;
(c11、c12 c13 C = 0, 0, 0;
c12 c11、c13, 0, 0, 0;
c13、c13 c33, 0, 0, 0;
c44 0, 0, 0, 0, 0;
c44 0, 0, 0, 0, 0;
0,0,0,0,0,c66];
phase_vel3 = 0 (1181);
r = 0;
θ= -ππ/ 180:π
r = r + 1;
n = [cos(θ),0,罪(θ)];
L_11 = (C (1,1) * n (1) ^ 2 + C (6,6) * n (2) ^ 2 + C (5,5) * n (3) ^ 2 + 2 * C (1,6) * n * n (2) (1) + 2 * C (1、5) * (1) n * n (3) + 2 * C (5、6) * n * n (2) (3));
L_12 = (C (1,6) * n (1) ^ 2 + C (2,6) * n (2) ^ 2 + C (4、5) * n (3) ^ 2 + (C (1、2) + C (6,6)) * (1) n * n (2) + (C (1、4) + C (5、6)) * (1) n * n (3) + (C (4、6) + C (2、5)) * n (2) * n (3));
L_13 = (C (1、5) * n (1) ^ 2 + C (4、6) * n (2) ^ 2 + C (3、5) * n (3) ^ 2 + (C (1、4) + C (5、6)) * (1) n * n (2) + (C (1、3) + C (5,5)) * (1) n * n (3) + (C (3、6) + C (4、5)) * n (2) * n (3));
L_22 = (C n (6,6) * (1) ^ 2 + C (2, 2) * n (2) ^ 2 + C (4, 4) * n (3) ^ 2 + 2 * C (2,6) * n * n (2) (1) + 2 * C (4、6) * n * n (3) (1) + 2 * C (2、4) * n * n (2) (3));
L_23 = (C (5、6) * n (1) ^ 2 + C (2、4) * n (2) ^ 2 + C (3、4) * n (3) ^ 2 + (C (4、6) + C (2、5)) * (1) n * n (2) + (C (3、6) + C (4、5)) * (1) n * n (3) + (C (2、3) + C (4, 4)) * n (2) * n (3));
L_33 = (C (5,5) * n (1) ^ 2 + C (4, 4) * n (2) ^ 2 + C (3,3) * n (3) ^ 2 + 2 * C (4、5) * n * n (2) (1) + 2 * C (3、5) * n * n (3) (1) + 2 * C (3、4) * n * n (2) (3));
Christoffel_mat = [L_11 L_12 L_13;
L_12、L_22 L_23;
L_13、L_23 L_33];
(ev, d) = eig (Christoffel_mat);
% eigvecs =极化向量
pl = ev (:, 3);
% eigvals ~缓慢(阶段)
phase_vel3 (r) =√rho. / d (3、3));% quasi-longitudinal
r = 0;
θ= -ππ/ 180:π
r = r + 1;
n = [cos(θ),0,罪(θ)];
L_11 = (C (1,1) * n (1) ^ 2 + C (6,6) * n (2) ^ 2 + C (5,5) * n (3) ^ 2 + 2 * C (1,6) * n * n (2) (1) + 2 * C (1、5) * (1) n * n (3) + 2 * C (5、6) * n * n (2) (3));
L_12 = (C (1,6) * n (1) ^ 2 + C (2,6) * n (2) ^ 2 + C (4、5) * n (3) ^ 2 + (C (1、2) + C (6,6)) * (1) n * n (2) + (C (1、4) + C (5、6)) * (1) n * n (3) + (C (4、6) + C (2、5)) * n (2) * n (3));
L_13 = (C (1、5) * n (1) ^ 2 + C (4、6) * n (2) ^ 2 + C (3、5) * n (3) ^ 2 + (C (1、4) + C (5、6)) * (1) n * n (2) + (C (1、3) + C (5,5)) * (1) n * n (3) + (C (3、6) + C (4、5)) * n (2) * n (3));
L_22 = (C n (6,6) * (1) ^ 2 + C (2, 2) * n (2) ^ 2 + C (4, 4) * n (3) ^ 2 + 2 * C (2,6) * n * n (2) (1) + 2 * C (4、6) * n * n (3) (1) + 2 * C (2、4) * n * n (2) (3));
L_23 = (C (5、6) * n (1) ^ 2 + C (2、4) * n (2) ^ 2 + C (3、4) * n (3) ^ 2 + (C (4、6) + C (2、5)) * (1) n * n (2) + (C (3、6) + C (4、5)) * (1) n * n (3) + (C (2、3) + C (4, 4)) * n (2) * n (3));
L_33 = (C (5,5) * n (1) ^ 2 + C (4, 4) * n (2) ^ 2 + C (3,3) * n (3) ^ 2 + 2 * C (4、5) * n * n (2) (1) + 2 * C (3、5) * n * n (3) (1) + 2 * C (3、4) * n * n (2) (3));
Christoffel_mat = [L_11 L_12 L_13;
L_12、L_22 L_23;
L_13、L_23 L_33];
(ev, d) = eig (Christoffel_mat);
% eigvecs =极化向量
pl = ev (:, 3);
% eigvals ~缓慢(阶段)
phase_vel3 (r) =√rho. / d (3、3));% quasi-longitudinal
极地(0:π/ 180:2 *π,1. / phase_vel3)
任何3 d援助如何想象这将不胜感激。
谢谢
戴夫:)

接受的答案

威廉•罗斯
威廉•罗斯 2022年8月9日
(移动从评论]回答我原本的计划
请尽可能简单的代码片段,清楚地说明了你想做什么。这总是一个好主意,为了提高ods获得有用的答案。
我无法告诉你想阴谋在第三维度。
Matlab不能直接画在圆柱形或球坐标,据我所知。但它有一个内置的函数从圆柱形转换到笛卡尔,下面我用。
z = 0:200;θ= z *π/ 25;r = 1 - z / 400;
[x, y, z] = pol2cart(θ,r, z);%圆柱形转换为笛卡尔坐标
图;
次要情节(121);极地(θ,r);%极坐标图
次要情节(122);plot3 (x, y, z);% 3 d图
我希望有帮助。
4评论

登录置评。

更多的答案(0)

类别

找到更多的在极地的情节帮助中心文件交换

标签

s manbetx 845


释放

R2021b

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!