DOC

# matlab

By Calvin Hill,2014-05-25 14:43
12 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