中南大学MATLAB课程实践

精品文档---下载后可任意编辑MATLAB 课程设计实践目录公共题 1 第一题 3 方 程组 3 解 非线性方程组 5 第二题 6 689精品文档---下载后可任意编辑公共题题目表示多晶体材料织构的三维取向分布函数（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 分布情况。程序流程图程序代码common.m%课程实践公共题目file=fopen('data.txt','r');%No_use 存储没有用的数据for i=1:18精品文档---下载后可任意编辑 No_use=fgetl(file);end%读入数据for i=1:19 %phi2 No_use=fscanf(file,'%f',1); for j=1:19%phi1 for k=1:19%phi f(j,k,i)=fscanf(file,'%f',1); end endend% slice 给出分布特征figure(1);[x,y,z]=meshgrid(0:5:90,0:5:90,0:5:90);slice(x,y,z,f,[45,90],[45,90],[0,45]);%pcolor 给出切面 f 情况figure(2);for i=1:19 subplot(5,4,i); [X,Y]=meshgrid(0:5:90); contour(X,Y,f(:,:,i)); axis ij;end%沿 alpha 取向线分布情况figure(3);plot([0:5:90],f(10,:,1),'-bo');text(60,6,'\phi=45');text(60,5.5,'\phi2=0');运行结果020406080100050100020406080100精品文档---下载后可任意编辑第一题题目编程实现以下科学计算算法，并举一例应用之。（参考书籍《精通ＭＡＬＡＢ科学计算》，王正林等著，电子工业出版社，２００９年） “不动点迭代法和牛顿法非线性方程组求解”1.1 不动点迭代法解非线性方程组算法说明设含有 n 个未知数与 n 个方程的非线性方程组记为：F(x)=0，然后把上述方程组改为便于迭代的等价形式：x=φ(x)，由此就可以构造不动点迭代法的迭代公0500500500500500500500500500500500500500500500500500500500500500500500500500500500500500500500500500500500500500500102030405060708090-20246810=452=0精品文档---下载后可任意编辑式：这样就可以求出非线性方程组的解。调...