matlab文件数据格式化读写小程序设计

文件数据格式化读写

一、实验目的

通过编辑和运行简单的matlab数据读写程序，熟悉和掌握数据格式化输入输出方法。

二、实验内容：

1、文件的打开与关闭

1)创建标准的对话框并通过交互式操作取得文件名

[FileName,PathName,FilterIndex]=uigetfile(FilterSpec,DialogTitle,Def aultName)

FileName：返回的文件名

PathName：返回的文件的路径名

FilterIndex：选择的文件类型

FilterSpec：文件类型设置

DialogTitle：打开对话框的标题

DefaultName：默认指向的文件名

例：

a.最简单的形式，只设置了一种可选择的文件类型

[FileName,PathName] = uigetfile('*.m','Select the M-file');

b.设置多种文件类型

[filename, pathname] =

uigetfile({'*.txt';'*.xls';'*.mat';'*.*'},'File Selector');

c.所有文件类型

[FileName,PathName] = uigetfile('*.*');

d.路径文件名合并

file = fullfile(pathname, filename)

或者file = [pathname, filename]

2）打开文件

在读写文件之前，必须先用fopen函数打开或创建文件，并指定对该文件进行的操作方式。fopen函数的调用格式为：fid=fopen（文件名，‘打开方式’）

说明：其中fid用于存储文件句柄值，如果返回的句柄值大于0，则说明文件打开成功。文件名用字符串形式，表示待打开的数据文件。常见的打开方式如下：‘r’：只读方式打开文件（默认的方式），该文件必须已存在。

‘r+’：读写方式打开文件，打开后先读后写。该文件必须已存在。

‘w’：打开后写入数据。该文件已存在则更新；不存在则创建。

‘w+’：读写方式打开文件。先读后写。该文件已存在则更新；不存在则创建。

‘a’：在打开的文件末端添加数据。文件不存在则创建。

‘a+’：打开文件后，先读入数据再添加数据。文件不存在则创建。

另外，在这些字符串后添加一个“t”，如‘rt’或‘wt+’，则将该文件以文本方式打开；如果添加的是“b”，则以二进制格式打开，这也是fopen函数默认的打开方式。

例：fid = fopen(file,'rt');

2）关闭文件

文件在进行完读、写等操作后，应及时关闭，以免数据丢失。关闭文件用fclose 函数，调用格式为：

sta = fclose(fid)

说明：该函数关闭fid所表示的文件。sta表示关闭文件操作的返回代码，若关闭成功，返回0，否则返回-1。如果要关闭所有已打开的文件用fclose(‘all’)。

2、文本文件的读写操作

1）读文本文件

fscanf函数可以读取文本文件的内容，并按指定格式存入矩阵。其调用格式为：

[A，COUNT]=fscanf(fid，format，size)

说明：其中A用来存放读取的数据，COUNT返回所读取的数据元素个数，fid 为文件句柄，format用来控制读取的数据格式，由%加上格式符组成，常见的格式符有：d（整型）、f（浮点型）、s（字符串型）、c（字符型）等，在%与格式符之间还可以插入附加格式说明符，如数据宽度说明等。size为可选项，决定矩阵A中数据的排列形式，它可以取下列值：N（读取N个元素到一个列向量）、inf（读取整个文件）、[M，N]（读数据到M×N的矩阵中，数据按列存放）。2）写文本文件

fprintf函数可以将数据按指定格式写入到文本文件中。其调用格式为：

fprintf（fid，format，A）

说明：fid为文件句柄，指定要写入数据的文件，format是用来控制所写数据格式的格式符，与fscanf函数相同，A是用来存放数据的矩阵。

例文件数据导入

test1.txt 的文件中以下为文件内容：

你好欢迎来到花果山

1 11 111 1111

2 22 222 2222

3 33 333 3333

4 44 444 4444

5 55 555 5555

方法：

[filename, pathname]= uigetfile('*.*');

file=fullfile(pathname, filename);

A=importdata(file);

A.data=

1 11 111 1111

2 22 222 2222

3 33 333 3333

4 44 444 4444

5 55 555 5555

A.textdata=

'你好' '欢迎来到' '防灾科技学院' ' https://www.360docs.net/doc/aa19130403.html,/ '

test2.txt 的文件中以下为文件内容：

你好

1 11 111 1111

欢迎来到

2 22 222 2222

防灾科技学院

3 33 333 3333

https://www.360docs.net/doc/aa19130403.html,/

4 44 444 4444

5 55 555 5555

方法：

[filename, pathname]= uigetfile('*.*'); %文件查找窗口

file=fullfile(pathname, filename); %合并路径文件名

fidin=fopen(file,'rt'); % 打开文件

fidout=fopen('dataout.txt','w'); % 创建输出文件

while ~feof(fidin) % 判断是否为文件末尾

tline=fgetl(fidin); % 从文件读行

if ~isempty(tline) % 判断是否空行

[m,n]=size(tline);

flag=1;

for i=1:n %判断一行中有没有字符（- +.和空格键除外）

if ~(tline(i)==' '| tline(i)=='.'| tline(i)=='-'| tline(i)=='+'| (double(tline(i))>=48&&double(tline(i))<=57))

flag=0;

break;

end

end

if flag==1 % 如果是数字行，把此行数据写入文件fprintf(fidout,'%s\n',tline); % 如果是数字行，把此行数据写入文件dataout.txt

end

end

end

fclose('all');

A=importdata('dataout.txt'); % 将生成的dataout.txt文件导入工作--

% 空间，变量名为A，

3、角度格式转换

通常角度存贮格式：dd.mmss

例如：22.3511183900为22°35′11.1839″

转换：度分秒－>度

function degree = dms2degree(jiaodu)

%度分秒(dd.mmss)－>度

degree = fix(jiaodu);

mimute = fix((jiaodu-degree)*100);

second = (jiaodu-degree-mimute/100)*10000;

degree = degree+mimute/60+second/3600;

三、实验任务

通过对实验内容的练习，完成以下任务：

1）读取文件data.txt中的数据，并将其中经纬度和大地高数据按照格式存贮到BLH.txt文件中，要求格式为：序号经度纬度大地高，且经纬度格式为度

读取data文件程序如下：

clear;

clc;

[filename, pathname]= uigetfile('*.*'); % 文件查找窗口

file=fullfile(pathname, filename); % 合并路径文件名

fidin=fopen(file,'rt'); % 打开文件

fidout=fopen('dataout.txt','w'); % 创建输出文件

while ~feof(fidin) % 判断是否为文件末尾

tline=fgetl(fidin); % 从文件读行

if ~isempty(tline) % 判断是否空行

[m,n]=size(tline);

flag=1;

for i=1:n % 判断一行中有没有字符（- +.和空格键除外）if ~(tline(i)==' '| tline(i)=='.'| tline(i)=='-'| tline(i)=='+'| (double(tline(i))>=48&&double(tline(i))<=57))

flag=0;

break;

end

end

if flag==1 % 如果是数字行，把此行数据写入文件fprintf(fidout,'%s\n',tline); % 如果是数字行，写入文件dataout.txt end

end

end

fclose('all');

A=importdata('dataout.txt'); % 将生成的dataout.txt文件导入工作空间，变量名为A

jiaodu=A(:,2:3); %提取A中所有行和2、3列

degree=dms2degree(jiaodu); %度分秒转化成度

B=[A(:,1),degree(:,2),degree(:,1),A(:,4)]; %组成矩阵B

fid=fopen('BLH.txt','wt'); %创建文件BLH

fprintf(fid,'序号WGS84 经度WGS84 纬度大地高\n' ); %将文本写入到指定的文件中

[m,n]=size(B); %输入矩阵B, m,n取行列号

for i=1:m

for j=1:n

if j==n

fprintf(fid,'%f\n',B(i,j)); %当j==n时，将数据写入到B矩阵else

fprintf(fid,'%f\t',B(i,j));

end

end

end

fclose(fid);

2) 编写函数：实现角度格式转换，度－>度分秒(dd.mmss)

函数程序如下：

function dms=degree3dms(du)%度－>度分秒(dd.mmss)

degree=fix(du);

min=fix((du-degree)*60);

second=(((du-degree)*60-min)*60);

dms=degree+min/100+second/10000;