Embedded Files
Post date: Aug 22, 2013 9:14:55 PM
Thank Dr. Qian HUANG from Technical University of Denmark for the help
Ref:
Nielsen JK, Rasmussen HK, Denberg M, Almdal K, and Hassager O. Macromolecules 2006;39(25):8844-8853.
Rasmussen HK, Christensen JH, and Gøttsche S. Journal of Non-Newtonian Fluid Mechanics 2000;93(2–3):245-263.
Huang Q, Mednova O, Rasmussen HK, Alvarez NJ, Skov AL, Almdal K, and Hassager O. Macromolecules 2013;46(12):5026-5035.
Huang Q, Alvarez NJ, Matsumiya Y, Rasmussen HK, Watanabe H, and Hassager O. ACS Macro Letters 2013;2(8):741-744.
-----------------------------------------------------------
function resBSW2=G1_G2_BSW2res(para)
n1e=para(1);
n2g=para(2);
G0N1=para(3);
lamda1max=para(4);
lamda2c=para(5);
lamda2max=para(6);
G0N2=para(7);
global w_G1w0_G2w0;
Y=w_G1w0_G2w0;
n=length(Y);
n1=n/3;
n2=n1*2;
w=Y(1:n1);
Y1=Y(n1+1:n2);
Y2=Y(n2+1:n);
resBSW2=zeros(n2,1);
resBSW2(1:n1)=(arrayfun(@(wvecsca) quadv(@(x) x^(1+n1e)/(1+x^2),0,wvecsca*lamda1max)*(wvecsca*lamda1max)^-n1e*n1e*G0N1+quadv(@(x) x^(1-n2g)/(1+x^2),0,wvecsca*lamda1max)*(wvecsca*lamda2c)^n2g*n1e*G0N1 + quadv(@(x) x^(1+n1e)/(1+x^2),0,wvecsca*lamda2max)*(wvecsca*lamda2max)^-n1e*n1e*G0N2+quadv(@(x) x^(1-n2g)/(1+x^2),0,wvecsca*lamda2max)*(wvecsca*lamda2c)^n2g*n1e*G0N2 , w)./Y1)-1;
resBSW2(n1+1:n2)= (arrayfun(@(wvecsca) quadv(@(x) x^(n1e)/(1+x^2),0,wvecsca*lamda1max)*(wvecsca*lamda1max)^-n1e*n1e*G0N1+quadv(@(x) x^(-n2g)/(1+x^2),0,wvecsca*lamda1max)*(wvecsca*lamda2c)^n2g*n1e*G0N1 + quadv(@(x) x^(n1e)/(1+x^2),0,wvecsca*lamda2max)*(wvecsca*lamda2max)^-n1e*n1e*G0N2+quadv(@(x) x^(-n2g)/(1+x^2),0,wvecsca*lamda2max)*(wvecsca*lamda2c)^n2g*n1e*G0N2 , w)./Y2)-1;
-------------------------------
w=LCBPS1M(:,1);
N=size(w,1);
G1w0=LCBPS1M(:,2);
G2w0=LCBPS1M(:,3);
loglog(w,G1w0,w,G2w0)
w(1:25,:)=[];
G1w0(1:25,:)=[];
G2w0(1:25,:)=[];
N=size(w,1);
loglog(w,G1w0,w,G2w0)
lamda1max=30;
lamda2c=0.006;
lamda2max=1600;
G0N0=7600;
G0N1=1392984;
G0N2=907508;
n1e=0.01;
n2g=0.73;
n1e=0.02;
para(1)=n1e;
para(2)=n2g;
para(3)=G0N1;
para(4)=lamda1max;
para(5)=lamda2c;
para(6)=lamda2max;
para(7)=G0N2;
G1w0_G2w0=[G1w0;G2w0];
global w_G1w0_G2w0;
w_G1w0_G2w0=[w;G1w0_G2w0];
fitParams = lsqnonlin (@G1_G2_BSW2res, para);
n1e_fit=fitParams(1);
n2g_fit=fitParams(2);
G0N1_fit=fitParams(3);
lamda1max_fit=fitParams(4);
lamda2c_fit=fitParams(5);
lamda2max_fit=fitParams(6);
G0N2_fit=fitParams(7);
G1w_fit = arrayfun(@(wvecsca) quadv(@(x) x^(1+n1e_fit)/(1+x^2),0,wvecsca*lamda1max_fit)*(wvecsca*lamda1max_fit)^-n1e_fit*n1e_fit*G0N1_fit+quadv(@(x) x^(1-n2g_fit)/(1+x^2),0,wvecsca*lamda1max_fit)*(wvecsca*lamda2c_fit)^n2g_fit*n1e_fit*G0N1_fit +quadv(@(x) x^(1+n1e_fit)/(1+x^2),0,wvecsca*lamda2max_fit)*(wvecsca*lamda2max_fit)^-n1e_fit*n1e_fit*G0N2_fit+quadv(@(x) x^(1-n2g_fit)/(1+x^2),0,wvecsca*lamda2max_fit)*(wvecsca*lamda2c_fit)^n2g_fit*n1e_fit*G0N2_fit , w);
G2w_fit = arrayfun(@(wvecsca) quadv(@(x) x^(n1e_fit)/(1+x^2),0,wvecsca*lamda1max_fit)*(wvecsca*lamda1max_fit)^-n1e_fit*n1e_fit*G0N1_fit+quadv(@(x) x^(-n2g_fit)/(1+x^2),0,wvecsca*lamda1max_fit)*(wvecsca*lamda2c_fit)^n2g_fit*n1e_fit*G0N1_fit + quadv(@(x) x^(n1e_fit)/(1+x^2),0,wvecsca*lamda2max_fit)*(wvecsca*lamda2max_fit)^-n1e_fit*n1e_fit*G0N2_fit+quadv(@(x) x^(-n2g_fit)/(1+x^2),0,wvecsca*lamda2max_fit)*(wvecsca*lamda2c_fit)^n2g_fit*n1e_fit*G0N2_fit , w );
w=logspace(-6,5,N);
n1=length(w_G1w0_G2w0);
w1=w_G1w0_G2w0(1:n1/3);
loglog(w1,G1w0,'gs',w1,G2w0,'bs',w,G1w_fit,'y.',w,G2w_fit,'k.')
Google Sites
Report abuse