各位大神好。

我编的程序运行后，两个数值相加，结果本来应该都是实数的，当设定某一参数时，部分结果出现了虚数单位，大概占结果的百分之30%，所以就不能画三维图形了，请问这么处理这个问题。

-0.0000                 -0.0000                 -0.0000

-0.0000 - 0.0000i  -0.0000                 -0.0000

-0.0000                 -0.0000                 -0.0000

-0.0000                 -0.0000                 -0.0000

-0.0000                 -0.0000 + 0.0000i  -0.0000

-0.0000                 -0.0000                 -0.0000

-0.0000                 -0.0000                 -0.0000

-0.0000 - 0.0000i  -0.0000 - 0.0000i  -0.0000 - 0.0000i

-0.0000                 -0.0000                 -0.0000

-0.0000 - 0.0000i  -0.0000                 -0.0000 + 0.0000i

-0.0000                 -0.0000                 -0.0000

-0.0000                 -0.0000                 -0.0000 - 0.0000i

-0.0000                 -0.0000                 -0.0000

-0.0000                 -0.0000                 -0.0000 - 0.0000i

-0.0000                 -0.0000                 -0.0000

-0.0000                 -0.0000 - 0.0000i  -0.0000

-0.0000 + 0.0000i  -0.0000                 -0.0000

-0.0000                 -0.0000                 -0.0000 + 0.0000i

附原程序：

clear

clc

a =0.1;

r=0.6;

t = sqrt(1 - abs(r).^2);

aa = abs(a);

at = abs(t);

ar = abs(r);

ca = conj(a);

ct = conj(t);

cr = conj(r);

x=linspace(-3,3,101);

y=linspace(-3,3,101);

[x,y]=meshgrid(x,y);

b=x+1i*y;

W1=( ar.^4.*( (2.*b+a).^2.*conj(2.*b-a).^2-4.*conj(2.*b-a).*(2.*b+a)+2 ) ...

-2.*t.*cr.*ar.^2.*a.*(  (2.*b+a).^2.*conj(2.*b-a)-2.*(2.*b+a))+2.*ct.*r.*ar.^2.*ca.*( conj(2.*b-a).^2.*(2.*b+a)-2.*conj(2.*b-a)  ) ...

-4.*at.^2.*ar.^2.*aa.^2.*( (2.*b+a).*conj(2.*b-a) -1 )+(t.^2.*cr.^2.*a.^2+t.*cr.*ar.^2).*(2.*b+a).^2 ...

+(ct.^2.*r.^2.*ca.^2+ct.*r.*ar.^2 ).*conj(2.*b-a).^2+(2.*t.*cr.*at.^2.*a.^2.*ca+2.*at.^2.*ar.^2.*ca).*(2.*b+a) ...

-(2.*ct.*r.*at.^2.*a.*ca.^2+2.*at.^2.*ar.^2.*a    ).*conj(2.*b-a) ...

+at.^4.*aa.^4+ct.*r.*at.^2.*ca.^2+ t.*cr.*at.^2.*a.^2+at.^2.*ar.^2  ).*exp(- 2*(b+a ).*conj(b-a)-2*aa.^2);

W2=(  ar.^4.*( (2.*b-a).^2.*conj(2.*b+a).^2-4.*conj(2.*b+a).*(2.*b-a)+2 ) ...

+2.*t.*cr.*ar.^2.*a.*(  (2.*b-a).^2.*conj(2.*b+a)-2.*(2.*b-a))      -2.*ct.*r.*ar.^2.*ca.*( conj(2.*b+a).^2.*(2.*b-a)-2.*conj(2.*b+a)  ) ...

-4.*at.^2.*ar.^2.*aa.^2.*( (2.*b-a).*conj(2.*b+a) -1 )+(t.^2.*cr.^2.*a.^2+t.*cr.*ar.^2).*(2.*b-a).^2 ...

+(ct.^2.*r.^2.*ca.^2+ct.*r.*ar.^2 ).*conj(2.*b+a).^2   -(2.*t.*cr.*at.^2.*a.^2.*ca+2.*at.^2.*ar.^2.*ca).*(2.*b-a) ...

+(2.*ct.*r.*at.^2.*a.*ca.^2+2.*at.^2.*ar.^2.*a    ).*conj(2.*b+a) ...

+at.^4.*aa.^4+ct.*r.*at.^2.*ca.^2+ t.*cr.*at.^2.*a.^2+at.^2.*ar.^2  ).*exp(- 2*(b-a ).*conj(b+a)-2*aa.^2);

W=W1+W2;

surf(x,y,W)

xlabel('W');

ylabel('x');

xlabel('y');