| 查看: 1469 | 回复: 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邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问三维金相立体图如何画出来
已经有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人回复
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
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 |
4楼2014-12-29 12:46:15
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
wzp0098: 金币+10, ★有帮助, 多谢，能帮忙改改吗？ 2014-12-30 16:56:21
wzp0098: 金币+10, ★有帮助, 多谢，能帮忙改改吗？ 2014-12-30 16:56:21
|
我查了一下，好像surf是平面图，你这个是三维的，所以不能表示 http://bbs.pinggu.org/thread-1171696-1-1.html |
5楼2014-12-30 15:18:45
20