精品文档---下载后可任意编辑MATLAB 程序设计实践班级：学号：姓名：一、《MATLAB 程序设计实践》Matlab 基础表示多晶体材料织构的三维取向分布函数（f＝f（φ1，φ，φ2））是一个非常复杂的函数，难以精确的用解析函数表达，通常采纳离散空间函数值来表示取向分布函数，Data.txt 是三维取向分布函数的一个实例。由于数据量非常大，不便于分析，需要借助图形来分析。请你编写一个 matlab 程序画出如下的几种图形来分析其取向分布特征：（1）用 Slice 函数给出其整体分布特征；（2）用 pcolor 或contour 函数分别给出(φ2＝0, 5, 10, 15, 20, 25, 30, 35 … 90)切面上 f 分布情况(需要用到subplot 函数)；(3) 用 plot 函数给出沿 α 取向线(φ1=0~90，φ＝45，φ2＝0)的 f 分布情况。备注：数据格式说明解：（1）将文件数据说明部分，与作图无关此方向表示 f 随着 φ1从 0,5,10,15, 20 … 到90 的变化而变化此方向表示 f 随着 φ从 0,5,10,15, 20 … 到90 的变化而变化表 示 以 下 数 据 为φ2＝0 的数据，即f（φ1，φ，0）精品文档---下载后可任意编辑Data.txt 内的数据根据要求读取到矩阵 f(phi1,phi,phi2)中，代码如下：fid=fopen('data.txt');for i=1:18 tline=fgetl(fid);endphi1=1;phi=1;phi2=1;line=0;f=zeros(19,19,19);while ~feof(fid) tline=fgetl(fid); data=str2num(tline); line=line+1;if mod(line,20)==1 phi2=(data/5)+1; phi=1;elsefor phi1=1:19 f(phi1,phi,phi2)=data(phi1);end phi=phi+1;endendfclose(fid);将以上代码保存为 readtext.m 文件并在 MATLAB 中运行，运行结果如下图所示：将以下代码保存为 code1_1.m 文件：fopen('readtext.m');readtext;[x,y,z]=meshgrid(0:5:90,0:5:90,0:5:90);slice(x,y,z,f,[45,90],[45,90],[0,45])运行结果如右图所示：（2）将以下代码保存为 code1_2_1.m 文件：fopen('readtext.m');readtext;for i=1:19 subplot(5,4,i) pcolor(f(:,:,i))End运行结果如右图所示：将以下代码保存为 code1_2_2.m 文件：fopen('readtext.m');readtext;for i=1:19 subplot(5,4,i) contour(f(:,:,i))end运行结果如右图所示：（3）φ1=0~90，φ＝45，φ2＝0 所对应的 f(φ1,φ,φ2)即为 f(:,10,1)。将以下代码保存为 code1_3.m 文件：fopen('readtext.m');readtext;plot([0:5:90],f(:,10,1),'-bo')text(60,6,'\phi=4...
