summaryrefslogtreecommitdiff
path: root/wing_scripts/get_z.m
diff options
context:
space:
mode:
authorMarius Peter <blendoit@gmail.com>2019-07-15 13:29:36 -0700
committerMarius Peter <blendoit@gmail.com>2019-07-15 13:29:36 -0700
commit32b6a6e80d22a6d6e8449821d3117252de3e89ba (patch)
treec658df6c4cb1f9d13ce62764c14cf2f5672878a0 /wing_scripts/get_z.m
parentcab1c2b9470ddc9099c4458daf2388d30ac12ca6 (diff)
commence class-ification of gui
Diffstat (limited to 'wing_scripts/get_z.m')
-rw-r--r--wing_scripts/get_z.m34
1 files changed, 34 insertions, 0 deletions
diff --git a/wing_scripts/get_z.m b/wing_scripts/get_z.m
new file mode 100644
index 0000000..5387b52
--- /dev/null
+++ b/wing_scripts/get_z.m
@@ -0,0 +1,34 @@
+function z = get_z(x,u)
+
+
+
+if (x < 0 )
+ disp('invalid X')
+end
+
+M = 0.02;
+P = 0.4;
+T = 0.12;
+a0 = 0.2969;
+a1 = -0.126;
+a2 = -0.3516;
+a3 = 0.2843;
+a4 = -0.1015;
+
+if x <P
+ z_camber = M/P^2*(2*P*x - x^2);
+else
+ z_camber = (M/(1-P)^2)*(1 - 2*P +2*P*x - x^2);
+end
+
+%z_camber = M/P^2*(2*P*x - x^2);
+z_thickness = (T/0.2)*(a0*x^.5+a1*x+a2*x^2+a3*x^3+a4*x^4);
+
+if u==1
+ z = z_camber + z_thickness;
+else
+ z = z_camber - z_thickness;
+end
+
+
+
Copyright 2019--2024 Marius PETER