很急 帮帮忙啊！

lambda=808e-7; %LD泵浦光的波长(cm)
wp0=500e-4; %LD泵浦光的束腰(cm)
I0=500; %沿激光棒轴向单位长度内的泵浦功率(W/cm^2)
r0=0.15; %激光棒的横截面半径(cm)
alpha=3; %激光棒对泵浦光的吸收系数
x=linspace(-r0,r0,100); %x坐标取值100个点(cm)
y=linspace(-r0,r0,100); %y坐标取值100个点(cm)
[xx1,yy1]=meshgrid(x,y); %第一束LD泵浦光的坐标(cm)
xx2=xx1*cos(2*pi/3)+yy1*sin(2*pi/3); %第二束LD泵浦光的x坐标(cm)
yy2=-xx1*sin(2*pi/3)+yy1*cos(2*pi/3); %第二束LD泵浦光的y坐标(cm)
xx3=xx1*cos(-2*pi/3)+yy1*sin(-2*pi/3); %第三束LD泵浦光LD的x坐标(cm)
yy3=-xx1*sin(-2*pi/3)+yy1*cos(-2*pi/3); %第三束LD泵浦光LD的y坐标(cm)
d1=sqrt(r0^2-xx1.^2)+yy1; %第一束LD泵浦光通过晶体的距离(cm)
d2=sqrt(r0^2-xx2.^2)+yy2; %第二束LD泵浦光通过晶体的距离(cm)
d3=sqrt(r0^2-xx3.^2)+yy3; %第三束LD泵浦光通过晶体的距离(cm)
wp1=wp0*sqrt(1+(lambda.*yy1/pi/(wp0^2))); %第一束LD泵浦光束宽的分布(cm)
wp2=wp0*sqrt(1+(lambda.*yy2/pi/(wp0^2))); %第二束LD泵浦光束宽的分布(cm)
wp3=wp0*sqrt(1+(lambda.*yy3/pi/(wp0^2))); %第三束LD泵浦光束宽的分布(cm)
I1=sqrt(2/pi)*I0./wp1.*exp(-2.*xx1.^2./(wp1).^2-alpha.*d1); %第一束LD泵浦光的光强分布
I2=sqrt(2/pi)*I0./wp2.*exp(-2.*xx2.^2./(wp2).^2-alpha.*d2); %第二束LD泵浦光的光强分布
I3=sqrt(2/pi)*I0./wp3.*exp(-2.*xx3.^2./(wp3).^2-alpha.*d3); %第三束LD泵浦光的光强分布
II=I1+I2+I3; %LD三向泵浦时的光强分布
I=II/max(max(II));
surf(xx1,yy1,I);
shading interp;grid off;box on;
%subplot(2,1,2);
figure;
contour(xx1,yy1,I,20);
colorbar('vert');
figure;
plot(xx1(23,:),I(

23,:));




??? Error using ==> surf at 78
X, Y, Z, and C cannot be complex.  怎么改。

:'(  帮个忙啦  做毕业设计呢。。。。

surf(xx1,yy1,I);
这句话错了，里头的I不能是复数的。

surf(xx1,yy1,I);
这句话错了，里头的I不能是复数的。

加个绝对值吧

lambda=808e-7; %LD泵浦光的波长(cm)
wp0=500e-4; %LD泵浦光的束腰(cm)
I0=500; %沿激光棒轴向单位长度内的泵浦功率(W/cm^2)
r0=0.15; %激光棒的横截面半径(cm)
alpha=3; %激光棒对泵浦光的吸收系数
x=linspace(-r0,r0,100); %x坐标取值100个点(cm)
y=linspace(-r0,r0,100); %y坐标取值100个点(cm)
[xx1,yy1]=meshgrid(x,y); %第一束LD泵浦光的坐标(cm)
xx2=xx1*cos(2*pi/3)+yy1*sin(2*pi/3); %第二束LD泵浦光的x坐标(cm)
yy2=-xx1*sin(2*pi/3)+yy1*cos(2*pi/3); %第二束LD泵浦光的y坐标(cm)
xx3=xx1*cos(-2*pi/3)+yy1*sin(-2*pi/3); %第三束LD泵浦光LD的x坐标(cm)
yy3=-xx1*sin(-2*pi/3)+yy1*cos(-2*pi/3); %第三束LD泵浦光LD的y坐标(cm)
d1=sqrt(r0^2-xx1.^2)+yy1; %第一束LD泵浦光通过晶体的距离(cm)
d2=sqrt(r0^2-xx2.^2)+yy2; %第二束LD泵浦光通过晶体的距离(cm)
d3=sqrt(r0^2-xx3.^2)+yy3; %第三束LD泵浦光通过晶体的距离(cm)
wp1=wp0*sqrt(1+(lambda.*yy1/pi/(wp0^2))); %第一束LD泵浦光束宽的分布(cm)
wp2=wp0*sqrt(1+(lambda.*yy2/pi/(wp0^2))); %第二束LD泵浦光束宽的分布(cm)
wp3=wp0*sqrt(1+(lambda.*yy3/pi/(wp0^2))); %第三束LD泵浦光束宽的分布(cm)
I1=sqrt(2/pi)*I0./wp1.*exp(-2.*xx1.^2./(wp1).^2-alpha.*d1); %第一束LD泵浦光的光强分布
I2=sqrt(2/pi)*I0./wp2.*exp(-2.*xx2.^2./(wp2).^2-alpha.*d2); %第二束LD泵浦光的光强分布
I3=sqrt(2/pi)*I0./wp3.*exp(-2.*xx3.^2./(wp3).^2-alpha.*d3); %第三束LD泵浦光的光强分布
II=I1+I2+I3; %LD三向泵浦时的光强分布
I=II/max(max(II));
surf(xx1,yy1,abs(I));           
shading interp;grid off;box on;
%subplot(2,1,2);
figure;
contour(xx1,yy1,abs(I),20);
colorbar('vert');
figure;
plot(xx1(23,:),I(23,:));

可以运行，不知道是不是你想要的结果。

zhanglu 发表于 2012-5-17 18:56
加个绝对值吧

lambda=808e-7; %LD泵浦光的波长(cm)

太感谢了  不管行不行！