运用CST方法拟合NACA翼型

利用CST参数化方法在C语言环境下，拟合NACA翼型数据，以实现较少的变量更加精确的拟合翼型

#include “stdio.h“
#include “stdlib.h“
#include “math.h“
#include “string.h“


void Accuracy（char infile_1[20]int n）
{
	FILE *fp1*fp2*fp3;
	char *middle_1*middle_2*middle_3;
	char *middle_4=“_NACA.txt“;
	char *middle_5=“_CST.txt“;
	char *middle_6=“_Accuracy.txt“;
	char infile_2[20];
	char outfile[20];
	int ij;
	double x_o[200]y_o[200]x_c[200]y_c[200];
	double e[200];
	double dD;


	for （i=0;i<20;i++）
	{
		infile_2[i]=infile_1[i];
	}

	for （i=0;i<20;i++）
	{
		outfile[i]=infile_1[i];
	}


	middle_1=infile_1;
	middle_2=infile_2;
	middle_3=outfile;

	middle_1=strcat（middle_1middle_4）;

	middle_2=strcat（middle_2middle_5）;

	middle_3=strcat（middle_3middle_6）;


	if（（fp1=fopen（middle_1“r“））==NULL）
	{
		printf（“cannot open file\n“）;
		exit（0）;
	}


	if（（fp2=fopen（middle_2“r“））==NULL）
	{
		printf（“cannot open file\n“）;
		exit（0）;
	}


	if（（fp3=fopen（middle_3“w“））==NULL）
	{
		printf（“cannot open file\n“）;
		exit（0）;
	}


	for （j=0;j<2;j++）
	{
		d=0.0;
		for （i=0;i<=n;i++）
		{
			fscanf（fp1“%lf%lf“&x_o[i]&y_o[i]）;
		}


		for （i=0;i<=n;i++）
		{
			fscanf（fp2“%lf%lf“&x_c[i]&y_c[i]）;
		}


		for （i=0;i<=n;i++）
		{
			e[i]=y_c[i]-y_o[i];
			fprintf（fp3“%.10lf%.10lf\n“x_o[i]e[i]）;
		}

		for （i=0;i<=n;i++）
		{
			d+=pow（e[i]2）;
		}

		D=sqrt（d/（double）n）;
		printf（“\n拟合标准差为:%.10lf\n“D）;
		fprintf（fp3“\n拟合标准差为:%.10lf\n\n“D）;

	}

	printf（“\n精度计算完成！\n\n“）;
	
	fclose（fp1）;
	fclose（fp2）;
	fclose（fp3）;

	


}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       6726  2016-05-26 15:32  NACA_CST （data）\Debug\4410_Accuracy.txt

     文件       6633  2016-05-26 15:32  NACA_CST （data）\Debug\4410_CST.txt

     文件        336  2016-05-26 15:32  NACA_CST （data）\Debug\4410_CstCoef.txt

     文件        329  2016-05-26 15:32  NACA_CST （data）\Debug\4410_Inter.txt

     文件       6630  2016-05-26 15:32  NACA_CST （data）\Debug\4410_NACA.txt

     文件      48640  2014-09-23 21:07  NACA_CST （data）\Debug\NACA_CST.exe

     文件     360896  2014-09-23 21:07  NACA_CST （data）\Debug\NACA_CST.ilk

     文件     510976  2014-09-23 21:07  NACA_CST （data）\Debug\NACA_CST.pdb

     文件       1568  2013-12-19 21:16  NACA_CST （data）\NACA_CST\Accuracy.cpp

     文件         39  2013-12-19 21:09  NACA_CST （data）\NACA_CST\Accuracy.h

     文件       5250  2013-12-19 20:36  NACA_CST （data）\NACA_CST\CST.cpp

     文件         42  2013-12-18 21:59  NACA_CST （data）\NACA_CST\CST.h

     文件      13150  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\Accuracy.obj

     文件      28594  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\BuildLog.htm

     文件      24831  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\CST.obj

     文件      11138  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\Inter.obj

     文件      25416  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\main.obj

     文件         65  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\mt.dep

     文件      19633  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\NACA.obj

     文件        663  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\NACA_CST.exe.embed.manifest

     文件        728  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\NACA_CST.exe.embed.manifest.res

     文件        621  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\NACA_CST.exe.intermediate.manifest

     文件     371712  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\vc90.idb

     文件     110592  2014-09-23 21:07  NACA_CST （data）\NACA_CST\Debug\vc90.pdb

     文件       1598  2013-12-19 21:12  NACA_CST （data）\NACA_CST\Inter.cpp

     文件         34  2013-12-18 21:59  NACA_CST （data）\NACA_CST\Inter.h

     文件        802  2013-12-19 21:09  NACA_CST （data）\NACA_CST\main.cpp

     文件       3221  2013-12-19 21:09  NACA_CST （data）\NACA_CST\NACA.cpp

     文件         35  2013-12-18 21:59  NACA_CST （data）\NACA_CST\NACA.h

     文件       4398  2014-09-23 21:07  NACA_CST （data）\NACA_CST\NACA_CST.vcproj

............此处省略12个文件信息