中南大学材料专业MATLAB实践4号..doc
文本预览下载声明
一、《MATLAB程序设计实践》Matlab基础
班级：
学号：
姓名：
（2）用pcolor或contour函数分别给出(φ2＝0, 5, 10, 15, 20, 25, 30, 35 … 90)切面上f分布情况(需要用到subplot函数)；
(3) 用plot函数给出沿α取向线(φ1=0~90，φ＝45，φ2＝0)的f分布情况。
备注：data.txt数据格式说明
将文件Data.txt内的数据按照要求读取到矩阵f(phi1,phi,phi2)中，
代码如下：
fid=fopen(data.txt); %读取数据文件Data.txt
for i=1:18
tline=fgetl(fid);
end
phi1=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;
else
for phi1=1:19
f(phi1,phi,phi2)=data(phi1);
end
phi=phi+1;
end
end
fclose(fid);
建立新的脚本，将以上代码保存为m文件，命名为readtext.m，在MATLAB运行后。其变量f结果为：
(1). 用Slice函数给出其整体分布特征
程序代码：
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]) %运用slice函数绘制图形
将程序代码保存为code1_1.m文件，其运行结果：
(2).. 用subplot函数和pcolor函数给出(φ2＝0, 5, 10, 15, 20, 25, 30, 35 … 90)切面上f分布情况；
程序代码：
fopen(readtext.m);
readtext; for i=1:19
subplot(5,4,i)
pcolor(f(:,:,i)) %运用pcolor函数绘制图形
end
将上述代码保存为code1_2_1.m文件，其运行结果：
. 使用subplot函数和contour函数给出(φ2＝0, 5, 10, 15, 20, 25, 30, 35 … 90)切面上f分布情况；
程序代码：
fopen(readtext.m); %运用contour函数绘制图形 readtext;
for i=1:19
subplot(5,4,i)
contour(f(:,:,i))
end
将上述代码保存为code1_2_2.m文件，运行结果：
（3）. 用plot函数给出沿α取向线(φ1=0~90，φ＝45，φ2＝0)的f分布情况。
程序代码：
fopen(readtext.m);
readtext;
plot([0:5:90],f(:,10,1),-bo) %运用plot函数绘制图形
text(60,6,\phi=45 \phi2=0)
将上述代码保存为code1_3.m文件，运行结果：
二 《MATLAB程序设计实践》科学计算（04）
班级
１、编程实现以下科学计算算法，并举一例应用之。（参考书籍《精通ＭＡＬＡＢ科学计算》，王正林等著，电子工业出版社，２００９年）“切比雪夫逼近”
其中为次切比雪夫多项式，具体表达可通过递推得出：
，
它们之间满足如下的正交关系：
在实际应用中，可根据所需的精度来截取有限的项数，切比雪夫级数中的系数由下式决定：
流程图：
程序源代码（m文件）：
function f = Chebyshev(y,k,x
显示全部