matlab复数向极坐标转换_matlab中如何将极坐标函数转化为直角坐标函数并绘图

alpha=pi/6;sigmaH=48.6%最大水平地应力;sigmav=53.68%垂向地应力;K=93500%切变模量;v=0.218%泊松比;k=2.128;fori=1:100r=i;forj=1:100n=pi/50;theta=n*j;sigma=(sigmaH+sigmav)/2+(s...

alpha=pi/6;

sigmaH=48.6%最大水平地应力;

sigmav=53.68%垂向地应力;

K=93500%切变模量;

v=0.218%泊松比;

k=2.128;

for i=1:100

r=i;

for j=1:100

n=pi/50;

theta=n*j;

sigma=(sigmaH+sigmav)/2+(sigmaH-sigmav)*cos(2*alpha)/2;

tau=(sigmaH-sigmav)*sin(2*alpha)/2;

K1=sigma*(2*pi)^0.5;

K2=tau*(2*pi)^0.5;

sigmax(i,j)=K1*(3*cos(0.5*theta)+cos(2.5*theta))/(4*(2*pi*r)^0.5)- K2*(7*sin(0.5*theta)+sin(2.5*theta))/(4*(2*pi*r)^0.5);

sigmay(i,j)=K1*(5*cos(0.5*theta)-cos(2.5*theta))/(4*(2*pi*r)^0.5)- K2*(sin(0.5*theta)-sin(2.5*theta))/(4*(2*pi*r)^0.5);