《《matlab實驗報告》word版》由會員分享,可在線閱讀,更多相關(guān)《《matlab實驗報告》word版(6頁珍藏版)》請在裝配圖網(wǎng)上搜索。
1、2011——2012學(xué)年第一學(xué)期
合肥學(xué)院數(shù)理系
實驗報告
課程名稱: MATLAB實驗
實驗項目: 特殊函數(shù)與圖形
實驗類別:綜合性□ 設(shè)計性 √ 驗證性□
專業(yè)班級: 09信息與計算科學(xué)
姓 名: 韋志強 學(xué) 號: 0907011039
實驗地點: 7#604
實驗時間: 2011-10-19
指
2、導(dǎo)教師: 李旭 成 績:
一.實驗?zāi)康模?
1學(xué)會繪制matlab中的一些符號函數(shù)和各種簡單圖形;
2初步了解關(guān)于簡單圖形的簡單程序的編寫;
3懂得腳本文件和運用函數(shù)文件的使用方法。
二.實驗內(nèi)容:
1.ezmesh('(6+2*cos(u))*cos(v)','(6+2*cos(u))*sin(v)','2*sin(u)',[0,2*pi,0,2*pi])
axis equal
hold on
ezmesh('(2*cos(u))*cos(v)','(2*cos(u))*sin(v)','2*sin(u)',[0,2
3、*pi,0,2*pi])
axis equal
title('圓環(huán)面和球面')
2.ezmesh('(2*cos(u))*cos(v)','(2*cos(u))*sin(v)','2*sin(u)',[0,2*pi,0,2*pi])
axis equal
title('球面')
grid off
3.ezmesh('(2*cos(u))*cos(v)','(2*cos(u))*sin(v)','6*sin(u)',[0,2*pi,0,2*pi])
axis equal
title('橢球面')
grid off
4.ezsurf('3*ta
4、n(u)*cos(v)','3*tan(u)*sin(v)','5*sec(u)',[-pi/2,3*pi/2,0,2*pi]);
axis auto
title('雙葉雙曲面')
5.ezsurf('3*sec(u)*cos(v)','3*sec(u)*sin(v)','5*tan(u)',[-pi/2,pi/2,0, 2*pi]);
axis auto
title('單葉雙曲面')
6.[x,y]=meshgrid(-25:1:25,-25:1:25);
z=x.^2/9-y.^2/4;
surf(x,y,z)
title('馬鞍面')
5、
grid on
meshz(x,y,z)
三. 實驗要求
1.學(xué)會并了解Matlab中一些符號函數(shù)的繪圖功能和簡單圖形的繪制,學(xué)會熟練運用各種繪圖語句;
2.在實驗中要求領(lǐng)會各條語句的意義,要求將每個所需要畫出的圖形的代碼寫成函數(shù)文件或者腳本文件。
四 實驗步驟及結(jié)果
1.環(huán)面和球面的組合圖:
[u,v]=meshgrid([0:1/100:2*pi]);
x=(6+2.*cos(u)).*cos(v);
y=(6+2.*cos(u)).*sin(v);
z=2.*sin(u);
mesh(x,y,z);
axis equal;
6、
hold on
[a,b,c]=sphere(20);
surf(2*a,2*b,2*c)
shading flat
2.1,球面:
[x,y,z]=sphere(50);
surf(4.*x+1,4.*y+1,4.*z+1)
2.2,橢球面:
ezsurf('(2*cos(u))*cos(v)','(2*cos(u))*sin(v)','6*sin(u)',
[0,2*pi,0,2*pi])
axis equal
title('橢球面')
2.3,雙葉雙曲面:
ezsurf('3*t
7、an(u)*cos(v)','3*tan(u)*sin(v)','5*sec(u)',[-pi/2,
3*pi/ 2, 0,2*pi]);
axis auto
title('雙葉雙曲面')
2.4,單葉雙曲面:
ezsurf('3*sec(u)*cos(v)','3*sec(u)*sin(v)','5*tan(u)',[-pi/2,pi/2, 0, 2*pi]);
axis auto
title('單葉雙曲面')
5,馬鞍面:
[x,y]=meshgrid (-20:1:20,-20:1:20);
z=x.^2/9-y.^2/
8、4;
meshz(x,y,z)
colormap(gray)
title('馬鞍面')
grid off
6,黎曼函數(shù):
function liman(n)
% n=input ('Enter n:');
x=[];
y=[];
k=1;
for i=2:n
for j=1:i-1
if gcd(i,j)==1
x(k)=j/i;
y(k)=1/i;
k=k+1;
end
end
end
plot(x,y,'.b')
axis([0,1,0,1
9、])
五 實驗總結(jié):
通過這次實驗我學(xué)會了如何熟練運用matlab中各種簡單圖形和一些符號函數(shù)的繪制,也懂得了如何運用函數(shù)文件和腳本文件的使用;我也初步了解了簡單程序的編寫方法。
學(xué)生簽名: 韋志強
2011 年 10月 19日
教師評語及成績
教師簽名:
年 月 日
5