| 查看: 1204 | 回复: 18 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
wzp0098木虫 (著名写手)
|
[求助]
用matlab 画三维图，请高手指教. 已有2人参与
|
||
|
程序如下，可是没有办法运行，请高手指教，重谢。 W=-3:0.1:3;K=-3:0.1:3i;V=-3:0.1:3; P = []; for n=1:length(W); Ky=W(n); for m=1:length(K); Kx=K(m); for v=1:length(V); Kz=V(v); w = 1; K2 = K1+2*w; x = 2*(sin(Kx)+sin(Ky)); r1 = 8; r2 = 0.001*r1; T = 3*sin(Kz); A=[ -r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B = [-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y = A\B; P(n, m, v) = imag(Y(3)); % <---------------------------P向量中的元素赋值 end end end [X,Y,Z]=meshgrid(W,K,V); surf(X,Y,Z,P); shading interp |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有9人回复
-
最失望的一年
已经有17人回复
-
为什么nbs上溴 没有产物点出现呢
已经有6人回复
-
求推荐博导
已经有4人回复
-
存款400万可以在学校里躺平吗
已经有34人回复
-
求助一下有机合成大神
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
疑惑？
已经有5人回复
wzp0098
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 2524.4
- 散金: 39
- 红花: 13
- 帖子: 1218
- 在线: 226.6小时
- 虫号: 853822
- 注册: 2009-09-22
- 性别: GG
- 专业: 光学
送红花一朵 |
请教高人，我想画出P(n, m, v)在整个四维空间的分布，将它所有的值都在图形里面体现出来，怎么修改？ function slice_4D W=-3:0.05:3;K=-3:0.05:3;V=-3:0.05:3; P = []; for n=1:length(W); Ky=W(n); for m=1:length(K); Kx=K(m); for v=1:length(V); Kz=V(v); w = 1;K1=0; K2 = K1+2*w; x = 2*(sin(Kx)+sin(Ky)); r1 = 8; r2 = 0.001*r1; T = 3*sin(Kz); A=[ -r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B = [-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y = A\B; P(n, m, v) = imag(Y(3)); % end end end [X,Y,Z]=meshgrid(W,K,V); figure;slice(X,Y,Z,P,[-2.5 0 2.5],0,[-2 2]) result = HistRate(P(  );[Pm,I]=max(result(:,2)); figure; result(I,1) Pv=isosurface(X,Y,Z,P,result(I,1)); p=patch(Pv) |
19楼2015-01-11 21:39:56
2楼2014-12-29 09:09:14
wzp0098
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 2524.4
- 散金: 39
- 红花: 13
- 帖子: 1218
- 在线: 226.6小时
- 虫号: 853822
- 注册: 2009-09-22
- 性别: GG
- 专业: 光学
|
我改了一下，程序还是运行不起来，请大虾指教！ W=-3:0.1:3;K=-3:0.1:3;V=-3:0.1:3; P = []; for n=1:length(W); Ky=W(n); for m=1:length(K); Kx=K(m); for v=1:length(V); Kz=V(v); w = 1;K1=0; K2 = K1+2*w; x = 2*(sin(Kx)+sin(Ky)); r1 = 8; r2 = 0.001*r1; T = 3*sin(Kz); A=[ -r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B = [-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y = A\B; P(n, m, v) = imag(Y(3)); % <---------------------------P向量中的元素赋值 end end end [X,Y,Z]=meshgrid(W,K,V); surf(X,Y,Z,P); shading interp |
3楼2014-12-29 12:45:43
wzp0098
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 2524.4
- 散金: 39
- 红花: 13
- 帖子: 1218
- 在线: 226.6小时
- 虫号: 853822
- 注册: 2009-09-22
- 性别: GG
- 专业: 光学
|
我改了一下，程序还是运行不起来，请大虾指教！ W=-3:0.1:3;K=-3:0.1:3;V=-3:0.1:3; P = []; for n=1:length(W); Ky=W(n); for m=1:length(K); Kx=K(m); for v=1:length(V); Kz=V(v); w = 1;K1=0; K2 = K1+2*w; x = 2*(sin(Kx)+sin(Ky)); r1 = 8; r2 = 0.001*r1; T = 3*sin(Kz); A=[ -r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B = [-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y = A\B; P(n, m, v) = imag(Y(3)); % <---------------------------P向量中的元素赋值 end end end [X,Y,Z]=meshgrid(W,K,V); surf(X,Y,Z,P); shading interp |
4楼2014-12-29 12:46:15