matlab作业答案.docVIP

MATLAB数据处理与应用 班级：07-机械12班 学号 姓名：邹常玖 第一次作业 clear; fid1= fopen (nodes.txt,rt); fid2=fopen(nodepairs_num.dat,wt+); num=0; while feof(fid1)==0 line=fgetl(fid1); num=num+1; if isempty(line)~=1 node(num,1:11)=str2num(line(1:11)); x_co(num)=str2num(line(14:25)); y_co(num)=str2num(line(28:41)); z_co(num)=str2num(line(44:57)); x_coor=x_co(num); if x_coor=1.0e-6 x_coor=1; end y_coor=y_co(num); if y_coor=1.0e-6 y_coor=1; end z_coor=z_co(num); if z_coor=1.0e-6 z_coor=1; end xyz(num)=x_coor*y_coor*z_coor; Lnode(num,1:57)=line; end end fclose(fid1); m=0; for i=1:num cor_x=x_co(i); cor_y=y_co(i); cor_z=z_co(i); if cor_x=1.0e-6 cor_x=1; end if cor_y=1.0e-6 cor_y=1; end if cor_z=1.0e-6 cor_z=1; end XYZ=cor_x*cor_y*cor_z; L=find(XYZ==xyz); if length(L)==2 if x_co(L(1))==x_co(L(2))y_co(L(1))==y_co(L(2))........ z_co(L(1))==z_co(L(2)) m=m+1; fprintf(fid2,%s\n,Lnode(L(1),1:11)); fprintf(fid2,%s\n,Lnode(L(2),1:11)); xyz(L(1))=rand(1)+rand(1); xyz(L(2))=rand(1)+rand(1); end end if length(L)=3 %disp(There is three or more nodes whose coordinates are identical); %disp(i); if x_co(L(1))==x_co(L(2))y_co(L(1))==y_co(L(2))........ z_co(L(1))==z_co(L(2)) fprintf(fid2,%s\n,Lnode(L(1),1:11)); fprintf(fid2,%s\n,Lnode(L(2),1:11)); xyz(L(1))=rand(1)+rand(1); xyz(L(2))=rand(1)+rand(1); m=m+1; end if x_co(L(1))==x_co(L(3))y_co(L(1))==y_co(L(3))........ z_co(L(1))==z_co(L(3)) fprintf(fid2,%s\n,Lnode(L(1),1:11)); fprintf(fid2,%s\n,Lnode(L(3),1:11)); xyz(L(1))=rand(1); xyz(L(3))=rand(1); m=m+1; end if x_co(L(2))==x_co(L(3))y_co(L(2))==y_co(L(3))........