clear;
syms?sita;
alfa=3*pi/180;
y?=?0.5*cos（sita）;
%?pu=0.043;qu=0.081;au=0.61;bu=1.6;cu=0.48;du=1.1;
pu=0.01;qu=0.05;au=1;bu=1;cu=1;du=0.6;
%?zu?=?pu*（1?+?cos（sita））^au?*（1?-?cos（sita））^bu?+qu*（1?+?cos（sita））^cu*（1?-?cos（sita））^du;
zu?=?pu*（sin（sita））^au?*（sin（sita））^bu?+qu*（sin（sita））^cu*（1?-?cos（sita））^du;
dzu=diff（zu）;
%?pl=0.022;ql=0.130;al=2.6;bl=1;cl=0.7;dl=1.7;
pl=0.022;ql=0.130;al=2.6;bl=1;cl=0.7;dl=1.7;
zl?=?pl*（1?+?cos（sita））^al*?（1?-cos（sita））^bl-ql*（1?+cos（sita））^cl*（1?-cos（sita））^dl;
dzl=diff（zl）;
Vu=（sin（sita）+2*zu）^2*（sin（sita-alfa）+sin（alfa））^2/（4*dzu^2*sin（sita）^2+sin（sita）^4）;
Vl=（sin（sita）+2*zl）^2*（sin（sita-alfa）+sin（alfa））^2/（4*dzl^2*sin（sita）^2+sin（sita）^4）;
Cpu=1-Vu;
Cpl=1-Vl;
dCLu=Cpu*（dzu*sin（alfa）-sin（sita）*cos（alfa）/2）;
dCLl=Cpl*（dzl*sin（alfa）-sin（sita）*cos（alfa）/2）;
%?CLu=integral（@（sita）dCLu0pi）;CLl=integral（@（sita）dCLlpi2*pi）;
%?CL=CLu+CLl;
CLu=0;CLl=0;
%?for?sita=0:0.001:pi
%?????CLu=CLu+0.001*subs（dCLu）;
%?????CLl=CLl+0.001*subs（dCLl）;
%?end
%?CL=CLu+CLl;
ezplot（yzu）hold?on;
grid?on;
%?ezplot（sitaCpl）hold?on;
%?ezplot（sitadCLu）hold?on;
%?ezplot（sitadCLl）hold?on;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1154??2017-12-14?15:27??m文件\airfoil_752.m

?????文件???????5981??2017-12-22?11:12??m文件\blade_752_20171219.m

?????文件???????2832??2017-12-15?10:37??m文件\blade_752_performance.m

?????目錄??????????0??2018-05-22?14:54??m文件

-----------?---------??----------?-----??----

?????????????????9967????????????????????4