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];
?????????