| 查看: 1457 | 回复: 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 |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问三维金相立体图如何画出来
已经有17人回复
- 如何用matlab导入导入三维数据~
已经有5人回复
- 请问采用AFM扫描的文件如何用matlab画出三维图？
已经有7人回复
- 求matlab功率谱的画法，我这样对么？
已经有6人回复
- 【已完结】matlab仿真图形（交流问题即有奖，呵呵）
已经有6人回复
- 菜鸟求助用matlab如何画三维图
已经有3人回复
- 急用！！！！求画三维拓扑图
已经有10人回复
- 求助已知三维坐标位置，如何画出空间图像
已经有5人回复
- 在matlab中用mesh画三维图存在的问题
已经有11人回复
- MATLAB中画三维双对数图形的问题
已经有10人回复
- MATLAB将三维图和其中一个面的等高线画出来的问题，谢谢了。
已经有7人回复
- 如何用matlab画二维网格
已经有5人回复
- 求助MATLAB做三维图问题
已经有4人回复
- 关于Matlab画三维曲面图的问题，求助高手，感激感激
已经有3人回复
- 请问如何用Matlab画工作空间的轮廓图
已经有10人回复
- 最近用HB2151表达遇到以下情况，恳请高手指教
已经有6人回复
- Matlab三维图中如何实现x，y坐标名称与各自的坐标轴平行？
已经有4人回复
- 怎么用matlab画极坐标三维图？各位大侠帮帮忙啊！
已经有5人回复
- 如何让MATLAB画出来的这个三维图更漂亮？
已经有3人回复
- 弹性模量各向异性的三维图。。。求高手解决！！
已经有10人回复
- 【求助】菜鸟求助 matlab 绘图问题 绘制光谱图
已经有3人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
wzp0098
8楼2014-12-31 09:37:08
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
20