原创力文档- 1、本文档共19页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
中南大学材料专业MATLAB实践4号
一、《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
文档评论（0）