如何画拟合直线

我用股票市场价格数据研究它的无标度特性，做出了节点度分布图，想在双对数坐标下，给节点度分布图画一条拟合直线，想如图1这样。 可是我画出的度分布的拟合线是图2这样。
我的代码如下：
p=polyfit([0:M],P_DeD,1)   
z=polyval(p,[0:M]);
loglog([0:M],P_DeD,'.');
hold on
plot(log([0:M]),z,'r');
    
其中 p =  -0.0004    0.0283
  
  不知道为什么拟合出来是曲线不是直线啊？请各位指教

最后一个代码是plot(log([0:M]),log(z),'r'); ; 不是plot(log([0:M]),z,'r'); 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。

fantuanxiaot 发表于 2014-11-20 12:42
最后一个代码是plot(log([0:M]),log(z),'r');  ; 不是plot(log([0:M]),z,'r'); 。。。。。。。。。。。。。 ...

我把它改成log(z)后，画出的图中，就没有拟合线了，只有原来的点。

p=polyfit([0:M],P_DeD,1)   
z=polyval(p,[0:M]);
这两句代码，我尝试了各种log组合，比如 ：p=polyfit(log([0:M]),log(P_DeD),1) ,或者z=polyval(p,log([0:M]))，拟合结果都不对。

plot(log([0:M]),log(z),'r-')

fantuanxiaot 发表于 2014-11-20 20:23
plot(log([0:M]),log(z),'r-')

改为这句，仍然没有拟合线画出来啊。还是只有原来的点。

plot(log([0:M]),log(P_DeD),'.');
p=polyfit([0:M],P_DeD,1)   
z=polyval(p,[0:M]);
hold on
plot(log([0:M]),log(z),'r-');
修改后拟合线是这条红线，仍然不是直线，哪位大侠帮忙看下啊

这个问题解决了，因为我的数据中，有不少0值，取对数产生NAN，无法拟合。去掉数据中的0值，就可以画出拟合直线了。感谢大家的关注。