汪长城《科学计算与MATLAB语言》第二章 Matlab程序设计.pptVIP

第二章 程序设计 声明和致谢 课件取自网络电子课件，只进行了很少的编辑和加工，对不知名的网络电子课件作者表示感谢！ 前言 现代大地测量领域的大型解算软件，如GPS和InSAR解算软件等，一般按以下方式完成解算： 数据文件编辑 参数文件配置 软件解算并保存结果到文件 读取解算结果文件进行分析和可视化显示（绘图） 文件操作是大型程序设计中很重要的一个内容，绘图是结果表达的重要手段 提纲 1.文件操作 1.1 文件的打开与关闭 1.2 文件的读写操作 1.3 数据文件定位写操作 2. 程序设计 2.1 M文件 2.2 函数文件 2.3 程序设计 2.4 程序调试 Matlab环境下的文件与其它系统一样，也有二类文件组成，一是文件，又称M文件，另一类是数据文件。系统除提供了文件的一般管理功能外，还提供了对数据文件进行操作的特殊功能函数。 数据文件的打开与关闭 Matlab提供了对数据文件建立、打开、读、写以及关闭等一系列函数，数据文件一般存放在磁盘等介质上，用文件名标识，系统对文件名没有特殊要求。文件数据格式有二种形式，一是二进制格式文件，二是ASCII文本文件，系统对这两类文件提供了不同的读写功能函数。 打开文件 在读写文件之前，必须先用fopen命令打开一个文件，并指定允许对该文件进行的操作。文件操作结束后，应及时关闭文件，以免数据的丢失或误修改。fopen函数格式为： Fid= fopen(filename，permission) 打开文件 其中filename为文件名，permission为文件格式，可以是下列格式之一： ‘r’ 打开文件，读数据，文件必须存在。 ‘w’ 打开文件，写数据，若文件不存在，系统会自动建立。 ‘a’ 打开文件，在文件末尾添加数据。 ‘r+’ 打开文件，可以读和写数据，文件必须存在。 ‘w+’ 打开文件，供读与写数据用，会删除以前数据。 ‘a+’ 打开文件，供读与添加数据用。 ‘W’ 打开文件供写数据用，无自动刷新功能。 ‘A’ 打开文件供添加数据用，无自动刷新功能。 打开文件 例如：打开一个名为std.dat的数据文件并进行读操作，其命令格式为： Fid=fopen( ‘std.dat’, ’r’ ) 上述打开格式均为二进制格式，如果想用ASCII文本格式，则必须在格式字符串中加上字符t，例如用’r t’表示以ASCII格式打开供读操作的数据文件。 关闭文件 文件在进行完读、写等操作后，应及时关闭，以保证文件的安全可靠。关闭文件命令格式为： Sta=fclose(Fid) 关闭Fid所表示的文件 Sta表示关闭文件操作的返回代码，若关闭成功，返回0，否则返回–1。 提纲 1.文件操作 1.1 文件的打开与关闭 1.2 文件的读写操作 1.3 数据文件定位写操作 2. 程序设计 2.1 M文件 2.2 函数文件 2.3 程序设计 2.4 程序调试 读二进制数据文件 fread 读二进制数据文件。格式为： [A,COUNT]=fread(Fid,size,precision) 其中A为数据矩阵，COUNT返回所读取的数据元素个数。size为可选项，若不选用则读取整个文件内容，若选用它的值可以是下列值： 读二进制数据文件 N 读取 N个元素到一个列向量。 inf 读取整个文件。 [M,N] 读数据到M×N的矩阵中，数据按列存放。 读二进制数据文件 precision用于控制所读数据的精度格式。 缺省格式为uchar，即无符号字符格式。例如：Fid=fopen(‘std.dat’, ‘r’); A=fread(Fid, 100, ’long’); Sta=fclose(fid); 以读数据方式打开数据文件std.dat，并按长整型数据格式读取文件的前100个数据放入向量A，然后关闭文件。 写二进制数据文件 fwrite 函数以二进制格式向数据文件写数据，其格式为： COUNT=fwrite (Fid, A, precision) 例如：Fid=fopen(‘magic5.bin’, ‘wb’); fwrite(Fid, magic, ’int32’); 上述语句将矩阵magic中的数据写入文件magic5.bin中，数据格式为32位整型二进制格式。 写二进制数据文件 【例】建立一数据文件test.dat，用于存放矩阵A的数据。 已知 A=[-0.6515 -0.2727 -0.4354 -0.3190 -0.9047 -0.7534 -0.4567 -0.3212 -0.4132 -0.3583 -0.9264 -0.8