把機翼圖形mapping成一個圓 Joukowski Transformation Airfoil

guangwei27

請問大大 Joukowski_Transformation .Joukowski Airfoil
下面程式碼可以做出一個機翼,但我想把它mapping成figure(2) 一個圓型,請問大大要怎麼做?謝謝!

>> l=30;
>> h=2;
>> t=3;
>> alpha=5;
>> alpha = alpha / 180 * pi;%攻角;換算成徑度
>> c = l/4;
>> m = (4/(3*sqrt(3)))*t*c/l;
>> n = h/2;
>> u0 = (m^2 + n^2)^0.5;
>> a = ((m+c)^2+n^2)^0.5;
>> x = linspace(-l/2,l/2,500);
>> figure(1)
>> cla;
>> hold on;
>> yy = ((l^2)/4*(1+(l^2)/(16*(h^2)))-x.^2).^0.5 - (l^2)/(8*h);
>> dely = 0.385*t*(1-2*x./l).*((1-(2*x./l).^2).^0.5);
>> y = [yy+dely ; yy-dely];
>> plot(x,y(1,:),x,y(2,:));
>> axis equal;

ChaChing

Ref:
http://www.mathworks.com/matlabc ... foil-transformation  ???

guangwei27

假如我若要延續我的程式碼的話,要怎麼做改寫 謝謝