简单独立连接模型（IC）matlab代码

clear;clf;clc
%计算矩阵Q（H，G）
dir=‘C:\MATLAB7\work‘;
A=sprintf（‘%s\\Tomo-gravity\\Data\\Matrix_data\\A1.dat‘dir）;
A=importdata（A）;
H=A（1:2:24:）;
G=A（2:2:24:）;
Q=[H;G];
%求出真实的链路数据
X=sprintf（‘%s\\Tomo-gravity\\X01.dat‘dir）;
Xall=importdata（X）‘;
Xreal=Xall（1:5:720:）;
Xsg=Xall（2:5:720:）;
Xstg=Xall（3:5:720:）;
Xgg=Xall（4:5:720:）;
Xgtg=Xall（5:5:720:）;
Y=A*Xreal;
Yi=Y（1:2:24:）;            
Ye=Y（2:2:24:）;
Y=[Yi;Ye];
%计算参数Pi
S=sum（Yi）;
S=repmat（S121）;
P=（sum（（Yi./S）‘））‘/2016;
%设定参数f
f=0.23;
fa=f;
fb=1-fa;
%计算矩阵theta
temp=zeros（112）;
for i=1:12
    for j=1:12
        if i==1&j==1
             theta=temp;
             theta（11）=P（i）;
         else if i==j
                 T=temp;
                 T（1i）=P（i）;
                 theta=[theta;T];