DOC

matlab

By Calvin Hill,2014-05-25 14:43
10 views 0
matlab2

22fXxx()??12

gXx()10???1

??

null
null
null

    ?

    z0.3?hg

    z?0.0188??r0.74

    z

    ??z

    ?zz?zwhile fr

    Eflag~=min

    z0.10.85(0.2)????1

     1.

    function [G,hg,a,L,b,Gm,hgm,am,bm]=fknow x(2)=r%x=[0.6718;0.46542;1]; and(1,

    1); G=4080*9.8; %N

     hg=0.845; %m

    phan_a=2.1; %m

    f=x(2)L=3.95; %m

    b=L-a; *x(1)*

    L/(b+xGm=9290*9.8; %N

    (1)*hghgm=1.17; %m

    ); am=2.95; %m

    bm=L-am; %2

    function w=fw

    w=0.75;

    3.

    function x=fint(z) [G,hg,a,L,b,Gm,hgm,am,bm]=fknow; c=[1]; phan_x(1)=z; r=(1-xflag=0; (2))*x(

     1)*L/(a-x(1)*hg); %

    c(14,1) phan_f_m=x(2)*x(1)*L/(bm+x(1)*hgm); %

    =0.85* phan_r_m=(1-x(2))*x(1)*L/(am-x(1)*hgm); %

    phan_ Ef=b/L/(x(2)-phan_f*hg/L) ; %

    r-x(1)- Er=a/L/(1-x(2)+phan_r*hg/L) ; %

    0.07; Ef_m=bm/L/(x(2)-phan_f_m*hgm/L) ; %

     Er_m=am/L/(1-x(2)+phan_r_m*hgm/L) ; %

     %c=[]; c(15,1)

     =0.85*

    phan_ %0.15

     if x(1)>=0.15 && x(1)<=0.3 f_m-x(

     c(1,1)=x(1)-phan_f-0.08; 1)-0.07

     c(2,1)=phan_f-x(1)-0.08; ;

     c(3,1)=x(1)-phan_r-0.08;

     c(4,1)=phan_r-x(1)-0.08; c(16,1)

     c(5,1)=x(1)-phan_f_m-0.08; =0.85*

     c(6,1)=phan_f_m-x(1)-0.08; phan_

     c(7,1)=x(1)-phan_r_m-0.08; r_m-x

     c(8,1)=phan_r_m-x(1)-0.08; (1)-0.0

     c(9,1)=0.7-Ef; 7;

     c(10,1)=0.7-Er;

     c(11,1)=0.7-Ef_m;

     c(12,1)=0.7-Er_m; %

     c(13,1)=0.85*phan_f-x(1)-0.07;

     0.30

     elseif x(1)>=0.30 && x(1)<=0.61 for

     c(1,1)=0.74*phan_r-x(1)+0.018; i=1:le

     c(2,1)=0.74*phan_r_m-x(1)+0.018; ngth(c

     c(3,1)=0.7-Ef; )

     c(4,1)=0.7-Er;

     c(5,1)=0.7-Ef_m; flag1=

     c(6,1)=0.7-Er_m; 1;

     c(7,1)=0.85*phan_f-x(1)-0.07;

     c(8,1)=0.85*phan_r-x(1)-0.07; if

     c(9,1)=0.85*phan_f_m-x(1)-0.07; c(i,1)>

     c(10,1)=0.85*phan_r_m-x(1)-0.07; 0

     flag1=

    0; %0.61

     elseif x(1)>=0.61 && x(1)<=0.8

     c(1,1)=0.74*phan_r-x(1)+0.018; end

     c(2,1)=0.74*phan_r_m-x(1)+0.018;

     c(3,1)=0.7-Ef; flag=fl

     c(4,1)=0.7-Er; ag*fla

     c(5,1)=0.7-Ef_m; g1;

     c(6,1)=0.7-Er_m;

     end end

    end flag=1; %

    x=x'; Er=a/

    L/(1-x4.

    function f=myfun2(x) (2)+ph[G,hg,a,L,b,Gm,hgm,am,bm]=fknow; an_r*

    hg/L) ;w=fw; %

    z=0.15:0.01:0.30; f=0; %for i=1:length(z)

     phan_f=x(2,1)*z(i)*L/(b+z(i)*hg); %

     phan_r=(1-x(2,1))*z(i)*L/(a-z(i)*hg); %

     phan_f_m=x(2,1)*z(i)*L/(bm+z(i)*hgm); %

    Ef_m= phan_r_m=(1-x(2,1))*z(i)*L/(am-z(i)*hgm); %

     bm/L/f=f+w*((phan_f_m-z(i))^2+(phan_r_m-z(i))^2)+(1-w)*((phan_f-z(i))^2+(phan(x(2)-p

    han_f_r-z(i))^2);%

    end _m*hg

    m/L) ;5.

    function [c,ceq]=mycon2(x) [G,hg,a,L,b,Gm,hgm,am,bm]=fknow; %phan_f=x(2)*x(1)*L/(b+x(1)*hg); %

    phan_r=(1-x(2))*x(1)*L/(a-x(1)*hg); %

    phan_f_m=x(2)*x(1)*L/(bm+x(1)*hgm); % phan_r_m=(1-x(2))*x(1)*L/(am-x(1)*hgm); %

    Er_mEf=b/L/(x(2)-phan_f*hg/L) ; %

    z=0.3:=am/L/(1-x(2)+phan_r_m*hgm/L) ; % c=[]; 0.01:0.

    61; %0.15

     c(1,1)=x(1)-phan_f-0.08; f=0;

     c(2,1)=phan_f-x(1)-0.08; for

     c(3,1)=x(1)-phan_r-0.08; i=1:le

     c(4,1)=phan_r-x(1)-0.08; ngth(z

     c(5,1)=x(1)-phan_f_m-0.08; )

     c(6,1)=phan_f_m-x(1)-0.08;

     c(7,1)=x(1)-phan_r_m-0.08; phan_

     c(8,1)=phan_r_m-x(1)-0.08; f=x(2,

     c(9,1)=0.7-Ef; 1)*z(i)

     c(10,1)=0.7-Er; *L/(b+

     c(11,1)=0.7-Ef_m; z(i)*h

     c(12,1)=0.7-Er_m; g);

     c(13,1)=0.85*phan_f-x(1)-0.07; %

     c(14,1)=0.85*phan_r-x(1)-0.07;

     c(15,1)=0.85*phan_f_m-x(1)-0.07;

     c(16,1)=0.85*phan_r_m-x(1)-0.07;

    ceq=[];

    6. function f=myfun3(x) [G,hg,a,L,b,Gm,hgm,am,bm]=fknow; phan_

    r=(1-xw=fw; %

Report this document

For any questions or suggestions please email
cust-service@docsford.com