机翼翼型Matlab编程.docVIP

编程界面类似下图，再加一个图形显示框和一个生成三维坐标数据的按钮。 该模块的功能是根据用户指定的中线形式以及相关的几何参数，选择相应计算公式来计算叶片截面上、下轮廓的坐标数据点，并将计算结果分别存入数组中。 b——叶片弦长 ???叶片转折角 叶片安装角和最大厚度不需要。 2．原始翼型数据见txt文件，两种即可，文件名NACA4309和NACA4412，文件中存放的是翼型上缘和下缘的二维坐标数据。Matlab将从此文件中读取二维坐标数据。 注意：其中x坐标为弦长b的百分数。 3．翼型中线形式，圆弧形和抛物线型两种。 翼型中线通常采用圆弧形或抛物线形。圆弧形中线可由两段或一段圆弧构成。 对于图1所示的双圆弧中线， (0xa) （1） 同理，R2圆弧在xoy坐标中的方程为 (axb) （2） 设叶型的弯曲角为θ，则设计中一般取，，。 图1 双圆弧中线 （2）对于由抛物线构成的翼型中线，其方程为 （3） 图2 叶片截面的上、下轮廓曲线上任一点坐标值都确定 如图2所示，沿弯曲的翼型中线有关坐标点的法线方向叠加原始翼型数据便可以获得所需的叶片截面形状。叶片截面的上、下轮廓曲线上任一点坐标值可以由下式确定： 叶片截面上轮廓： （4） 叶片截面下轮廓： （4） 式（4）中x，ye的值可由所选原始翼型提供型值点数据代入计算；yz值根据不同形式的中线分别由式（1）、（2）、（3）确定；中线上任意点的切线方向与x轴之间的夹角φ可由下式计算： （5） 4.我想要的程序是能够读入原始坐标参数，然后通过计算得到最终的坐标参数，然后存储为txt或者dat格式的文件。 如果你感觉数据缺失也没有关系，你只需要帮我把整体程序做一下，然后我自己再修改。 clc clear %x，ye的值可由所选原始翼型提供型值点数据代入计算,这里假设二者均已知 %先计算yz的值，以式(1)为例，其他类似 %假设θ、a、b(若未知，再讨论) sita=（θ具体的值）; a=（a的具体值）; b=（b的具体值）; %计算yz的值 syms x yz dyz fai yz=((a/sin(0.6*sita)^2-(x-a)^2))^0.5-a/(tan(0.6*sita)); yzvalue=subs(yz,x);%注意这里的x就指的是一个具体的x的值，由于不知道x的具体形式，这里只讨论了x为一个具体数值的情况，若x为数组编程稍加修改即可 %下面计算φ dyz=diff(yz,x); fai=subs(atan(dyz),x);%这里的x也跟上述一样是具体的数值，不是参数 %下面叶片截面上轮廓xu，yu xu=x-ye*sin(fai); yu=yzvalue+ye*cos(fai); %叶片截面下轮廓与此类似 数据的导出比较简单，若能调通，再编写 建议取一个x及ye的值代入计算，看运行结果是否满足要求 另我对题目的理解可能有问题，再联系