matlab难题
有木有哪位了解matlab的大侠，可以帮小妹看一下这段程序错在哪里吗？
clc;
clear all;
%A站位置
x_a=10000;
y_a=10000;
z_a=0;
%B站位置
x_b=-10000;
y_b=-11000;
z_b=0;
%测角误差
deta_ele=0.001;
deta_ath=0.001;
deta_sta=5;
deta_deta=6;
syms x y;
%测角高度
z=10000;
theta_a =0.7834;
phi_a =0.6145;
Ra=sqrt((x-x_a)^2+(y-y_a)^2+(z-z_a)^2);
Rb=sqrt((x-x_b)^2+(y-y_b)^2+(z-z_b)^2);
Cax=(x-x_a)/Ra;
Cay=(y-y_a)/Ra;
Caz=(z-z_a)/Ra;
Cbx=(x-x_b)/Rb;
Cby=(y-y_b)/Rb;
Cbz=(z-z_b)/Rb;
A=[-(sin(theta_a))^2/(y-y_a) [cos(theta_a)]^2/(x-x_a) 0;-(z-z_a)*cos(theta_a)/Ra^2 -(z-z_a)*sin(theta_a)/Ra^2 cos(phi_a)/Ra;Cbx-Cax Cby-Cay Cbz-Caz];
A=vpa(A,4)
Pdx=inv(A)*[deta_ath^2+1/((x-x_a)^2+(y-y_a)^2) 0 0;0 deta_ele^2+1/Ra^2 0;0 0 (deta_deta)^2+2]*inv(A');
Pdx=vpa(Pdx,4)
[x,y]=meshgrid(-50:1:50);
GDOP=sqrt(Pdx(1,1)+Pdx(2,2)+Pdx(3,3))
[C,h]=contour(x,y,GDOP)
clabel(C,h)
程序运行到GDOP这个结果时就死机了，一直运算不出GDOP的值，而且也没有出来图形。
希望各位大侠帮帮小妹哦，呵呵
[解决办法]
死机还是出错?如果出国就把error message贴出来
因为就我看Pdx是个矩阵
GDOP=sqrt(Pdx(1,1)+Pdx(2,2)+Pdx(3,3))这一句不应该死机一直busy,
就是个简单的求根函数,怎么会死机?
matlab可以求负数与复数的sqrt,那更没出错可能了
你的现象到底是什么?
[解决办法]
去ｍａｔｌａｂ论坛吧 这里会的人少
