2011数学建模 汽车修理问题 - 图文(4)

p0=1/p0; for j=0:c

p(j+1)=(c*(1/(0.739*x)))^j/factorial(j)*p0; end

for j=c+1:6;

p(j+1)=c^c/factorial(c)*(1/(0.739*x))^j*p0; end lq=0;

for k=c+1:6;

lq=lq+(k-c)*p(k+1); end

ls=lq+c*(1/(0.739*x))*(1-p(7));

w=(lq*x*0.0778)./(1-p(7))+12.768*c*x+3.75*lq*x+5.074*c*(1-p(7))+87.9 for n=1:3

A(m,n)=subs(w,x,n) end end

z1=subs(ls,x,3) z2=subs(lq,x,3) 第五问的程序： clear syms x for m=1:3 c=m; p0=0;

for i=0:c-1;

p0=p0+(c*(1/(0.739*x)))^i/factorial(i); end

p0=p0+c^c/factorial(c)*((1/(0.739*x))^(c+1)-(1/(0.739*x))^7)/(1-(1/(0.739*x))); p0=1/p0; for j=0:c

p(j+1)=(c*(1/(0.739*x)))^j/factorial(j)*p0; end

for j=c+1:6;

p(j+1)=c^c/factorial(c)*(1/(0.739*x))^j*p0; end lq=0;

for k=c+1:6;

lq=lq+(k-c)*p(k+1); end

ls=lq+c*(1/(0.739*x))*(1-p(7));

w=(lq*x*0.0778)./(1-p(7))+12.768*c*x+3.75*lq*x+5.074*c*(1-p(7))+87.9 //k=1时的w表达式

w=(lq*x*0.0778*12.5)./(1-p(7))+12.768*c*x+3.75*lq*x+5.074*c*(1-p(7))+87.9

//k=12.5时的w表达式

z=42.375*c*x+12.5*ls*x+293 //问题三的目标函数 for n=1:3

A(m,n)=subs(w,x,n) end end

z0=subs(z,x,3) z1=subs(ls,x,3)

z2=subs(lq,x,3) //对应的t是后来用笔算出来的