【文章內容簡介】 ell{d,n} =k 。 k=[]。 end endmmm=zeros(50,k2)。for t=1:50 for j=1:k2 k=cell{t,j}。n=length(k)。if n==0 break。elsegz=zeros(1,n)。 % Coefficients: p1 = 。 p2 = 。 p3 = 。 p4 = 。 p5 = 。 p6 = 。 p7 = 。 p8 = 。 p9 = 。 q1 = 。 q2 = 。 q3 = 。 q4 = 。 q5 = 。 q6 = 。 q7 = 。 q8 = 。 q9 = 。 q10 = 。y11=1/9*p1*k.^9+1/8*p2*k.^8+1/7*p3*k.^7+1/6*p4*k.^6+1/5*p5*k.^5+1/4*p6*k.^4+1/3*p7*k.^3+1/2*p8*k.^2+p9*k。 y2 = q1*k.^9 + q2*k.^8 +q3*k.^7 + q4*k.^6 + q5*k.^5 + q6*k.^4 +q7*k.^3 + q8*k.^2 +q9*k.^1 + q10。 zg=0。 for i=1:n if i==1gz(i)=(10*i+3000+200*k(i))*(y2(i)0)200*(y11(i)0)。zg=zg+gz(i)。 else gz(i)=(10*i+3000+200*k(i))*(y2(i)y2(i1))200*(y11(i)y11(i1))。zg=zg+gz(i)。 endend wgz=(1000+10*(n1))*(1y2(n))。 zhp=y11(n)+k(n)*(1y2(n))。 y=(wgz+zg)/(zhp)。 mmm(t,j)=y。 endendendfor i=1:50for j=1:k2 if mmm(i,j)==0 mmm(i,j)=100。 endendendk3=[b,v]=min(min(mmm))。k4=[b,v]=min(min(mmm39。))。weizhi=find(mmm==mmin)cell{k3,k4}問題二的程序（1）函數方程擬合源代碼： %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%a=[459 362 624 542 509 584 433 748 815 505612 452 434 982 640 742 565 706 593 680926 653 164 487 734 608 428 1153 593 844527 552 513 781 474 388 824 538 862 659775 859 755 649 697 515 628 954 771 609402 960 885 610 292 837 473 677 358 638699 634 555 570 84 416 606 1062 484 120447 654 564 339 280 246 687 539 790 581621 724 531 512 577 496 468 499 544 645764 558 378 765 666 763 217 715 310 851]。a=reshape(a,1,100)。junzhi=mean(a)。biaozhuncha=std(a)。fangcha=var(a)。x1=1:1200。%概率密度y0=normpdf(x1,600,)。%由于考慮到刀具損壞故障占95%, 其它故障僅占5%。y1=y0*+*ones(1,length(y0))。plot(x1,y1),%figure 1 積分擬合后=y2。%對機器零件的概率分布的擬合；y11=y1.*x1。 figure, plot(x1,y11),%figure 2 %對y1積分的概率分布函數y2。 syms x p1 p2 p3 p4 p5 p6 p7 p8 p9 y1=p1*x^8 + p2*x^7 +p3*x^6 + p4*x^5 +p5*x^4 + p6*x^3 +p7*x^2 + p8*x^1 +p9 。y2=int(y1,x,0,39。k39。)。%即y2=1/9*p1*x^9+1/8*p2*x^8+1/7*p3*x^7+1/6*p4*x^6+1/5*p5*x^5+1/4*p6*x^4+1/3*p7*x^3+1/2*p8*x^2+p9*x。%對y11的積分y11=y1(i)*x(i)的概率分布=y3；syms x q1 q2 q3 q4 q5 q6 q7 q8 q9 q10。y11=q1*x^9+q2*x^8+q3*x^7+q4*x^6+q5*x^5+q6*x^4+q7*x^3+q8*x^2+q9*x^1+q10。 y3=int(y11,x,0,39。k39。)。 %即y3=1/10*q1*k^10+1/9*q2*k^9+1/8*q3*k^8+1/7*q4*k^7+1/6*q5*k^6+1/5*q6*k^5+1/4*q7*k^4+1/3*q8*k^3+1/2*q9*k^2+q10*k。%Coefficients:Norm of residuals = 。 p1 = 。 p2 = 。 p3 = 。 p4 = 。 p5 = 。 p6 = 。 p7 = 。 p8 = 。 p9 = 。 %Coefficients:Norm of residuals = q1 = 。 q2 = 。 q3 = 。 q4 = 。 q5 = 。 q6 = 。 q7 = 。 q8 = 。 q9 = 。 q10 = 。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%(2)問題二車床管理程序：1)使用fomfunction [c,ceq]=myfun1(k)ceq=[]。c=[]。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%2)M文件源代碼:function y=chechuang(k)clear， clc， %Coefficients:Norm of residuals = 。 p1 = 。 p2 = 。 p3 = 。 p4 = 。 p5 = 。 p6 = 。 p7 = 。 p8 = 。 p9 = 。 %Coefficients:Norm of residuals = q1 = 。 q2 = 。 q3 = 。 q4 = 。 q5 = 。 q6 = 。 q7 = 。 q8 = 。 q9 = 。 q10 = 。%每次調用時令n=5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20.n=10。k(n+1)=1200。R1=。R2=。zhengpin=0。zhp=zeros(1,n+1)。zhp0=0。zhpn=0。zhpn1=0。for t=1:n+1 y21=1/9*p1*k(t)^9+1/8*p2*k(t)^8+1/7*p3*k(t)^7+1/6*p4*k(t)^6+1/5*p5*k(t)^5+1/4*p6*k(t)^4+1/3*p7*k(t)^3+1/2*p8*k(t)^2+p9*k(t)。y31=1/10*q1*k(t)^10+1/9*q2*k(t)^9+1/8*q3*k(t)^8+1/7*q4*k(t)^7+1/6*q5*k(t)^6+1/5*q6*k(t)^5+1/4*q7*k(t)^4+1/3*q8*k(t)^3+1/2*q9*k(t)^2+q10*k(t)。if t==1 for j=1:n1 zhp0=zhp0+R1*k(j)*R2^(j1)。 end zhp(t)=zhp(t)+(1R2)*zhp0+R2^(n1)*(R1*k(1)+R2*(k(n)k(1)))。 elseif t1amp。t=n1 y32=1/10*q1*k(t1)^10+1/9*q2*k(t1)^9+1/8*q3*k(t1)^8+1/7*q4*k(t1)^7+1/6*q5*k(t1)^6+1/5*q6*k(t1)^5+1/4*q7*k(t1)^4+1/3*q8*k(t1)^3+1/2*q9*k(t1)^2+q10*k(t1)。y22=1/9*p1*k(t1)^9+1/8*p2