数学建模_B题 艾滋病疗法的评价及疗效的预测(答案)(6)

X1=data201(D(1:j(1),2),:); x14=X1(:,4);

t=x14; a1=0.85; a2=1.60;

yw1=(t*7-floor(ceil(t*7/30)/2)*30)*a1+floor(ceil(t*7/30)/2)*30*a2; yw2=(t*7-floor(ceil(t*7/30)/2)*30)*a2+floor(ceil(t*7/30)/2)*30*a1; w1=rem(ceil(t*7/30),2); w2=1-w1;

yw=yw1.*w1+yw2.*w2;

x11=X1(:,1);x12=X1(:,2);x13=X1(:,3);x14=X1(:,4);x15=X1(:,5);x16=yw; XX1=[x11,x12,x13,x14,x15,x16];

X2=data201(D(j(1)+1:j(1)+j(2),2),:);

x21=X2(:,1);x22=X2(:,2);x23=X2(:,3);x24=X2(:,4);x26=X2(:,4)*(1.60+1.85)*7; x25=X2(:,5);

XX2=[x21,x22,x23,x24,x25,x26];

X3=data201(D(j(1)+j(2)+1:j(1)+j(2)+j(3),2),:);

x31=X3(:,1);x32=X3(:,2);x33=X3(:,3);x34=X3(:,4);x36=X3(:,4)*(1.60+0.85)*7; x35=X3(:,5);

XX3=[x31,x32,x33,x34,x35,x36];

X4=data201(D(j(1)+j(2)+j(3)+1:j(1)+j(2)+j(3)+j(4),2),:); x41=X4(:,1);x42=X4(:,2);x43=X4(:,3);x44=X4(:,4); x46=X4(:,4)*(1.60+0.85+1.20)*7; x45=X4(:,5);

XX4=[x41,x42,x43,x44,x45,x46];

data1=XX4; n=size(data1); j=1;

ss(j)=data1(1,1); for i=2:n(1)

if data1(i,1)-data1(i-1,1)>=1 j=j+1; ss(j)=data1(i,1); end end M=ss'; n2=size(M);

x1=data1(:,3); x2=data1(:,6); y=data1(:,5);

26

x=[x1,x2];

rstool(x,y,'purequadratic')

27