i0=5.83;r=0.367;t=0.85;m=3880;G=38024; f=0.013;CdA=2.77;pg=7.12;ig3=1.71; n=[815 1207 1614 2012 2603 3006 3403 3804]; B0=[1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7]; B1=[-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291] ; B2=[72.379 36.657 14.524 7.0035 4.4763 2.8593 2.9788 0.71113]; B3=[-5.8629 -2.0553 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449 -0.00075215]; B4=[0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230 -0.000038568]; Tq=-19.313+295.27*n/1000-165.44*((n/1000).^2)+40.874*((n/1000).^3)-3.8445*((n/1000).^4); P=(1/t)*(G*f*ua/3600+CdA*ua.^3/76140); b=B0+B1.*P+B2.*P.^2+B3.* P.^3+B4.*P.^4; plot(ua,Qs,'ko',ua4,f,'LineWidth',2); ylabel('Qs/[1/L(100km)]') |