数据大概是这样，x=1 2 3 4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29

y=96.552
96.552
96.552
96.552
96.552
96.552
96.552
89.655
82.759
82.759
79.31
79.31
75.862
75.862
75.862
72.414
72.414
72.414
65.517
62.069
62.069
58.621
55.172
48.276
48.276
20.69
10.345
6.897
3.448
然后我想做个双对数坐标轴的双线性拟合，然后可以进行双段幂律分布分析。
结果这样

数据我看了，你可以做成log，瞬间就能看到断点在25了。

那意思是自己划分？大神知道这个编程吗？类似这样
>> x=[ ];  / 输入修正的书信间隔时间
>> y=[ ];  / 输入累积概率
>> loglog(x,y,'ko')  / 画出双对数下坐标图（k表示黑色，o表示圆圈）
>> hold on  / 保留刚才所画图表，以便继续在此图画出拟合直线
>>  a=polyfit(log(x),log(y),1)  / 求拟合直线的参数（一次项和常数项）

a =

   -1.0700    5.9525  

>> b=2.71828^5.9525  /  常数项转化

b =

  384.7124

>> x=100:100000;  /  根据上面所作的图指定x的取值范围
>> y=b*x.^ -1.0700 ;  / 根据刚才的计算输入x和y的关系式
>> plot(x,y,'k')  / 在双对数坐标下画出拟合直线
我怎么也画不出来两条直线在图上
求大神给个具体过程

你可以将breakpoint test的程序敲出来，相当复杂。这些基本的命令和运算没办法。你是求数字的a的值，这个我不大会，我们都是用矩阵运算出来的b=（x‘x）^-1 *(x'y).然后其他的程序也是这么做出来的。你用的方法，对我来说，感觉太老旧了。

你对计量的理解有些问题。所有问题都是有专业计量方法的，不是你看着像两条就是两条的。所以，一般过程是，1，回归，2，系数稳定性检测，3，发现过不了检测说明不稳定，4怀疑有断点，5专业断点检测，6能得到几个断点并且知道断点的位置，7分段再回归或者加入二元变量再回归，8再检测回归的系数稳定性。最后稳定了，再画拟合线。

那这个断点应该是25，（看出来的）专业断点检测是画回归目测还是有一定的软件结果输出（R2?)，假设断点就是25，你所说的分段再回归或者加入二元变量再回归具体步骤怎么做？还有所得到的R2是在哪里体现的。谢谢大神指导。

大神，你所说的每一步具体怎么做可以教教我吗

大哥，你找本实证计量的书，看看学学就搞定了。

{:3_60:}老师催得紧，菜鸟一枚，才开始接触这个东西，老师要求国庆前必须做出来平时还要上课。大神，有时间帮忙做一下吗，可以给报酬，我真的没有办法了[cry]，唐突了。

你去项目交易区里面找人把

你好，我遇到了跟你一样的问题，请问你这个问题你搞定了么，能否分享一下代码呢，感谢大神