function v1=svenson_gz_joint(x);
load Data.mat;
 %svenson 模型
bd=zeros(size(p));
d=zeros(size(cf));
j=1;
num=63;
b0=x(1);
b1=x(2);
b2=x(3);
b3=x(4);
t1=x(5);
t2=x(6);

for i=1:num
for s=t(i,:)
if s==0
break;
else
m1=-s/t1;
m2=-s/t2;
r=b0+b1*((exp(m1)-1)/m1)+b2*((exp(m1)-1)/m1-exp(m1))+b3*((exp(m2)-1)/m2-exp(m2));
b=exp(-r*s);
end
d(i,j)=b;
j=j+1;
end
j=1;
end
m=d.*cf;
for i=1:num
for h=m(i,:)
if h==0
break;
else bd(i)=bd(i)+h;
end
end
end


v1=zeros(size(du));
dd=sum(1./du);
w=du;
for i=1:num
v1(i)=1/w(i)/dd*(p(i)-bd(i))/(num*1000)^(0.5);%為防止函數值過小，在優化時過早收斂乘上1000 不影響結果
end                                          


save jjh.mat;