中南大学材料专业matlab实践-4号

《中南大学材料专业matlab实践-4号》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、一、《MATLAB程序设计实践》Matlab基础班级：学号：姓名：表示多晶体材料织构的三维取向分布函数（f＝f（φ1，φ，φ2））是一个非常复杂的函数，难以精确的用解析函数表达，通常采用离散空间函数值来表示取向分布函数，Data.txt是三维取向分布函数的一个实例。由于数据量非常大，不便于分析，需要借助图形来分析。请你编写一个matlab程序画出如下的几种图形来分析其取向分布特征：（1）用Slice函数给出其整体分布特征；19（2）用pcolor或contour函数分别给出(φ2＝0,5,10,15,20,25,30,35…90)切面上f分布情况(需要用到subplot函数)；(3)用plo

2、t函数给出沿α取向线(φ1=0~90，φ＝45，φ2＝0)的f分布情况。19备注：表示以下数据为φ2＝0的数据，即f（φ1，φ，0）此方向表示f随着φ从0,5,10,15,20…到90的变化而变化此方向表示f随着φ1从0,5,10,15,20…到90的变化而变化数据说明部分，与作图无关data.txt数据格式说明将文件Data.txt内的数据按照要求读取到矩阵f(phi1,phi,phi2)中，代码如下：fid=fopen('data.txt');%读取数据文件Data.txtfori=1:18tline=fgetl(fid);endphi1=1;phi=1;phi2=1;line=0;f=

3、zeros(19,19,19);while~feof(fid)tline=fgetl(fid);data=str2num(tline);line=line+1;ifmod(line,20)==1phi2=(data/5)+1;phi=1;19elseforphi1=1:19f(phi1,phi,phi2)=data(phi1);endphi=phi+1;endendfclose(fid);建立新的脚本，将以上代码保存为m文件，命名为readtext.m，在MATLAB运行后。其变量f结果为：(1).用Slice函数给出其整体分布特征程序代码：fopen('readtext.m');readt

4、ext;[x,y,z]=meshgrid(0:5:90,0:5:90,0:5:90);slice(x,y,z,f,[45,90],[45,90],[0,45])%运用slice函数绘制图形将程序代码保存为code1_1.m文件，其运行结果：(2)..用subplot函数和pcolor函数给出(φ2＝0,5,10,15,20,25,1930,35…90)切面上f分布情况；程序代码：fopen('readtext.m');readtext;fori=1:19subplot(5,4,i)pcolor(f(:,:,i))%运用pcolor函数绘制图形end将上述代码保存为code1_2_1.m文件，

5、其运行结果：.使用subplot函数和contour函数给出(φ2＝0,5,10,15,20,25,30,35…90)切面上f分布情况；程序代码：fopen('readtext.m');%运用contour函数绘制图形readtext;fori=1:19subplot(5,4,i)contour(f(:,:,i))end将上述代码保存为code1_2_2.m文件，运行结果：19（3）.用plot函数给出沿α取向线(φ1=0~90，φ＝45，φ2＝0)的f分布情况。程序代码：fopen('readtext.m');readtext;plot([0:5:90],f(:,10,1),'-bo')%

6、运用plot函数绘制图形text(60,6,'phi=45phi2=0')将上述代码保存为code1_3.m文件，运行结果：19二《MATLAB程序设计实践》科学计算（04）班级１、编程实现以下科学计算算法，并举一例应用之。（参考书籍《精通ＭＡＬＡＢ科学计算》，王正林等著，电子工业出版社，２００９年）“切比雪夫逼近”算法说明：当一个连续函数定义在区间[-1,1]上时，它可以展开成切比雪夫级数。即：其中为次切比雪夫多项式，具体表达可通过递推得出：，它们之间满足如下的正交关系：在实际应用中，可根据所需的精度来截取有限的项数，切比雪夫级数中的系数由下式决定：19流程图：结束开始定义sym型变量

7、t以及切比雪夫多项式矩阵T(n)，并规定前两项为1，t输出k项切比雪夫多项式求出逼近x0处六位精度下的值根据前两项多项式计算其系数c(1)，c(2)是输入三项?否由第一项和第二项多项式及其系数计算一级切比雪夫逼近值是i>k+1？i从3渐渐增加到k+1否是由T(i-1)和T(i-2)计算T(i)由T(i)和c(i)计算i级切比雪夫逼近值。由T(i)计算其系数c(i)程序源代码（m文件）：functionf=Ch