中南大学matlab课程设计

《中南大学matlab课程设计》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、一、《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分布

2、情况(需要用到subplot函数)；(3)用plot函数给出沿α取向线(φ1=0~90，φ＝45，φ2＝0)的f分布情况。开始流程图读取Data.txtSlice函数绘制图形Pcolor函数绘制图形Contour函数绘制图形Plot函数绘制图形结束解：（1）将文件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=z

3、eros(19,19,19);while~feof(fid)tline=fgetl(fid);data=str2num(tline);line=line+1;ifmod(line,20)==1phi2=(data/5)+1;phi=1;elseforphi1=1:19f(phi1,phi,phi2)=data(phi1);endphi=phi+1;endendfclose(fid);将以上代码保存为readtext.m文件并在MATLAB中运行，运行结果如下图所示：将以下代码保存为code1_1.m文件：fopen('readtext.m

4、');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])%运用slice函数绘制图形运行结果如右图所示（2）将以下代码保存为code1_2_1.m文件：fopen('readtext.m');readtext;fori=1:19subplot(5,4,i)pcolor(f(:,:,i))%运用pcolor函数绘制图形end运行结果如右图所示将以下代码保存为code1_2_2.m文件：fopen('readtext.m')

5、;%运用contour函数绘制图形readtext;fori=1:19subplot(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')%运用plot函数绘制图形text(60,6,'phi=45phi2=0')运行结果如下图所示：1.编程实现以下科学计算算法，并举一

6、例应用之。（参考书籍《精通MATLAB科学算法》，王正林等著，电子工业出版社，2009年）“多项式拟合”。思考：多项式拟合是用多项式拟合曲线的一种方式，低次数下运用此方法符合较好，但较高次数下波动太大，失去真实性。1.1多项式曲线拟合概述对给定数据点（xi,yi）(i=1,2,...N），构造m次多项式，P（x）=+（m

7、：A=multifit(X,Y,m)其中：X为实验数据点的x坐标向量；Y为实验数据点的y坐标向量；m为拟合多项式的次数；A为拟合多项式的系数向量。1.2多项式曲线拟合编程流程图输入向量X，Y，多项式次数m.M=N得出X的项数为N，Y的项数为M建立长度为2m+1的零向量C和长度为2m的零向量b，变量j＝0,k＝0。j<2m+1YN显示输入不正确k

8、1.3多项式曲线拟合的MATLAB程序代码functionA=multifit(X,Y,m)%离散实验数据点的多项式曲线拟合%实验数据点的x坐标向量：X%实验数据点的y坐标向量：Y%拟合多项式