CODE:
function zbzbzbzb
data=[ 0.04556 0.04352 0.099 0.03592 % c1,q1,c2,q2
0.10253 0.08542 0.12602 0.08207
0.2056 0.17075 0.27154 0.16137
0.6889 0.30201 0.87702 0.27525
1.51941 0.38387 1.81968 0.34115
2.40038 0.45855 2.72298 0.41266
3.37753 0.51955 3.82469 0.45595
6.32386 0.60045 6.74548 0.54047
9.51436 0.64661 9.89625 0.59228
12.75747 0.68528 13.15723 0.62842
16.07689 0.7131 16.59229 0.63979];
xdata=[data(:,1),data(:,3)];
ydata=[data(:,2),data(:,4)];
x0=[0.2,0.2,0.2,9314,106,0.6]; %n11,n12,k12,n22,n21,k21
[x,resnorm] = lsqcurvefit(@myfun,x0,xdata,ydata)
y=myfun(x,xdata)
end
function y=myfun(x,xdata)
c1=xdata(:,1);
c2=xdata(:,2);
n11=x(1);
n12=x(2);
k12=x(3);
n22=x(4);
n21=x(5);
k21=x(6);
k1=0.45967;
n1=1/0.33824;
k2=0.40308;
n2=1/0.38477;
q1=(k1*c1.^(1/n11 +1/n1 ))./(c1.^(1/n11 )+k12*c2.^(1/n12 ) );
q2=(k2*c2.^(1/n22 +1/n2 ))./(c2.^(1/n22 )+k21*c1.^(1/n21 ) );
y=[q1,q2];
end