第1页共14页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第1页共14页C语言程序设计实验报告学院：资源加工与生物工程专业：无机非金属材料工程班级：0501姓名：金光龙学号：0302050107指导老师：金胜明制作日期：05年7月9日第2页共14页第1页共14页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第2页共14页C语言程序设计说明书目录一程序设计目的及任务；二程序设计的内容；三源程序；四程序说明；五实践心得；一程序设计目的及任务第3页共14页第2页共14页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第3页共14页1、学会数据文件的建立在实验室试验和现场生产过程中，经常要存储、处理许多数据（试验数据和生产数据），在没有在线检测的条件下，这些数据都要进行人工的输入，建立数据文件。以前，我们都是编写建立数据文件的程序来建立数据文件，当某个数据输入错误时就要重新输入或编写数据文件修改程序，比较麻烦。直接在编辑状态下（TC,TP等）编写数据文件简单、直观、实用。2、理解并且学会运用文件操作函数3、学会用线性回归法建立数学模型二程序设计的内容一、数据文件1．文件操作函数（1）文件的打开fopen#includeFILE*fopen(char*filename,char*type);filename—文件名，字符串变量；type—文件打开的方式(r,w,a)函数功能：打开由文件filename指定的文件，并把一个stream与它联系起来。返回一个指针，该指针用来识别随后操作的文件stream。若成功地完成打开，则返回最近打开的文件stream，若出现错误，则返回NULL。(2)文件的关闭fclose#includeintfclose(FILE*stream);filename—文件名，字符串变量；type—文件打开的方式(r,w,a)函数功能：关闭指定的流式文件。（3）fscanf—从流中执行格式化输入#includeintfscanf(FILE*stream,char*format[,argument,…]);函数功能：从指定的输入流中按规定的格式把数据读入由地址参数&arg1,…&argn给定的单元中。（4）fprintf—发送格式化输出到流#includeintfprintf(FILE*stream,char*format[,argument,…]);函数功能：按规定的格式把数据输出到指定的流式文件stream中。例：#include“stdio.h”main(){FILEfp,*fopen();floatd,date[10][2];inti,j;第4页共14页第3页共14页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第4页共14页fp=fopen(“试验数据.dat”，“r”);for(i=0;i<10;i++)or(j=0;j<2;j++){fscanf(fp,“%f”,&d);date[i][j]=d;}fclose(fp);}例如，建立存储下列数据的数据文件，并对它进行读、写操作。序号Na2OK2OCaOSiO2弹性模量17.268.1612.68727.263626.767.6611.68727.258236.267.1610.68727.245647.267.6610.68747.200156.767.1612.68747.313066.268.1611.68747.256077.267.1611.68767.255886.768.1610.68767.198896.267.6612.68767.3086二、线性回归法建立数学模型1建立数学模型的方法（1）机理分析法（2）数据分析法回归分析法---静态数模时间序列分析法—动态数模（系统辨识）（3）综合分析法2建立线性回归模型的步骤（1）第一步y=f(x1,x2,⋯,xn)y=β0+β1x1+β2x2+⋯βnxn第5页共14页第4页共14页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第5页共14页一组试验数据（m组，m>n）y1=β0+β1x11+β2x12+⋯βnx1n⋮yi=β0+β1xi1+β2xi2+⋯βnxin⋮ym=β0+β1xm1+β2xm2+⋯βnxmn目标：求βj的估计值bj。使yi−(b0+b1xi1+⋯+bnxin)=εi(预报误差）最小。用最小二乘法，使Q=∑i=1m(yi−y¿i)2=极小值∂Q∂b0=2∑i=1m[yi−(b0+b1xi1+¿⋅¿+bnxin)]=0⋮∂Q∂bk=2∑i=1m[yi−(b0+b1xi1+¿⋅¿+bnxin)]xik=0⋮mb0+b1∑i=1mxi1+⋯+bn∑i=1mxin=∑i=1myib0∑i=1mxi1+b1∑i=1mxi12+⋯+bn∑i=1mxinxi1=∑i=1myixi1b0∑i=1mxik+b1∑i=1mxi1xik+⋯+bn∑i=1mxinxik=∑i=1myixikS00=m,S01=∑i=1mxi1,⋯,S0n=∑i=1mxin,d0=∑i=1myiS10=∑i=1mxi1,S11=∑i=1mxi12,⋯,S1n=∑i=1mxinxi1,d1=∑i=1myixi1⋮Sk0=∑i=1mxik,Sk1=∑i=1mxi1xik,⋯,Skn=∑i=1mxinxik,dk=∑i=1myixik⋮第6页共14页第5页共14页编...
