Question

matlab高手帮我看看啊

我有一段程序，怎么画出的图出来总是坐标都是1呀好奇怪，帮忙看看改改啦
Sigma_1 =1.2030;
I_1 =1.1200;
I_2 =1.0863;
I_5 =0.0830;
x=0.4;
T =1.1837e-007;
rho=0.5;
Q=6;
e=1.60218*10^(-19);
z =1.5972e-032;
syms d;
Sigma1=Sigma_1+d*I_2;
I5=Sigma1-I_1;
r_1=-(Sigma1+I5)/(2*Sigma1*I5)+[(Sigma1+I5)^2+16*(1+x)*Sigma1*I5/x^2]^(1/2)/(2*Sigma1*I5);
r_2=(1/r_1+Sigma1)^(1/2)+(1/r_1+I5)^1/2;
P_r=1/(2*T*rho)*(e*Q^(2+x)*r_1^(x/2)*r_2^(2+x)*z^(x/2))^(1/(1+x));
Id=10^(-9);
E=1-T/(d+T-Id/(2*rho*P_r));
d=0:0.001:0.5;
plot(subs(E),subs(P_r))

谢谢啦

Answer 1

这样表达就可以了。 

Sigma_1 =1.2030; 

I_1 =1.1200; 

I_2 =1.0863; 

I_5 =0.0830; 

x=0.4; 

T =1.1837e-007; 

rho=0.5; 

Q=6; 

e=1.60218*10^(-19); 

z =1.5972e-032; 

syms d; 

Sigma1=Sigma_1+d*I_2; 

I5=Sigma1-I_1; 

r_1=-(Sigma1+I5)/(2*Sigma1*I5)+[(Sigma1+I5)^2+16*(1+x)*Sigma1*I5/x^2]^(1/2)/(2*Sigma1*I5); 

r_2=(1/r_1+Sigma1)^(1/2)+(1/r_1+I5)^1/2; 

P_r=1/(2*T*rho)*(e*Q^(2+x)*r_1^(x/2)*r_2^(2+x)*z^(x/2))^(1/(1+x)); 

Id=10^(-9); 

E=1-T/(d+T-Id/(2*rho*P_r)); 

d=0:0.001:0.5; 

format long 

E=subs(E) 

plot(E-1,subs(P_r)) 

结果： 

E = 

Columns 1 through 5 

1.00000004663732 1.00000004671297 1.00000004678889 1.00000004686507 1.00000004694152 

Columns 6 through 10 

1.00000004701824 1.00000004709521 1.00000004717244 1.00000004724993 1.00000004732766 

......

Answer 2

这是因为你算得的横坐E，基本上都是1－10^－8左右，而坐标上显然不能显示出这么精确的数值来，所以只有都显示1了，
有两种解决方法，一种是直接把E的数值变一下，像下面改动E
E=1-T/(d+T-Id/(2*rho*P_r))*10^-7;
就可以了，还有一种就是在坐标显示上作文章了，改动一下E：
E=-T/(d+T-Id/(2*rho*P_r));
然后再在坐标里insert->title,输入"1-",然后再将这个标题拖到横坐标前注明一下，横坐标的值都要加上这个"1-"