Question

matlab 编程出现了小错误，求高手指点！！！（高斯光束的matlab 模拟）

高斯光束的三维模拟
clear all
N=200;
w0=200e-9;

rho=linspace(0,3*w0,N);
theta=linspace(0,2*pi,N);
x=meshgrid(rho,theta);
y=pol2cart(theta,rho);
z=exp(-2*rho.^2./w0.^2);
surf(x,y,z);
colormap hot;

shading interp;
xlabel('x /m');
ylabel('y /m');
zlabel('intensity /a.u.');
axis()
title('焦点附近三维光场分布')
colorbar;

【跳出来的错误】
Z must be a matrix, not a scalar or vector.

我知道这话是说，Z是矩阵/向量，不是标量。但是我不知道哪里错了，求指导，O(∩_∩)O谢谢!!

Answer 1

你的函数是，z 与theta无关，实际上与rho的关系用平面曲线表示就够了，硬要画三维图，就改成下面的了
close all;
clear;
clc;
N=200;
w0=200e-9;

theta=linspace(0,2*pi,N);
rho=linspace(0,3*w0,N);
dz=exp(-2*rho.^2/w0^2);
[x,y]=meshgrid(theta,rho);
z=zeros(N);
for i=1:N
z(i,:)=dz;
end
surfc(x,y,z)

colormap hot;

shading interp;
xlabel('x /m');
ylabel('y /m');
zlabel('intensity /a.u.');
axis()
title('焦点附近三维光场分布')
colorbar;

追问

为什么是这样的？？

追答

你的函数关系确定的啊，看一下是不是漏项了

Answer 2

surf(x,y,z)里x,y可以是向量，但是第三项z必须是一个矩阵；
z里的每一点表示坐标为(x,y)上曲面的高度。