仿真平臺:matlab繪制曲面圖像_第1頁
仿真平臺:matlab繪制曲面圖像_第2頁
仿真平臺:matlab繪制曲面圖像_第3頁
仿真平臺:matlab繪制曲面圖像_第4頁
仿真平臺:matlab繪制曲面圖像_第5頁
已閱讀5頁,還剩9頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權,請進行舉報或認領

文檔簡介

1、一、設計目的熟悉曲面圖形的編制,通過常量參數(shù)的變化觀察圖形變化的規(guī)律。二、設計內(nèi)容已知二次曲面方程為x2a2+y2b2+z2c2=d,要求討論參數(shù)a,b,c,d對其形狀的影響并畫出圖形。要求:1.為便于觀察,繪圖時要求加入網(wǎng)格2.實現(xiàn)當a,b,c,d獨立增大或同時增大的圖形變化3.分析觀察的結果4.編制實訓報告三、設計思路x2a2+y2b2+z2c2=dz=±cd-x 2a2-y2b2利用matlab畫出二次曲面的圖形,編程如下:a=input('a=');b=input('b='); c=input('c=');d=input(

2、9;d='); %輸入?yún)?shù), xgrid=linspace(-abs(a*sqrt(d),abs(a*sqrt(d),50); %建立x網(wǎng)格坐標50為網(wǎng)格線數(shù)目ygrid=linspace(-abs(b*sqrt(d),abs(b*sqrt(d),50); %建立y網(wǎng)格坐標x,y=meshgrid(xgrid,ygrid); %確定50×50個點的x,y網(wǎng)格坐標z=c*sqrt(d-y.*y/b2-x.*x/a2);u=1; %u=1,表示z要取正值z1=real(z); %取z的實部z1for k=2:49; %以下7行程序的作用是取消z中含虛數(shù)的點for j=2:49if

3、 imag(z(k,j)=0;%imag函數(shù)取z的虛部z1(k,j)=0;endif all(imag(z(k-1:k+1,j-1:j+1)=0z1(k,j)=NaN;endendendsurf(x,y,z1); %畫空間曲面hold on%使當前軸及圖形保持而不被刷新,在此基礎上畫圖if u=1z2=-z1;surf (x,y,z2);%u=1時加畫負半面 axis(-abs(a*sqrt(d),abs(a*sqrt(d),-abs(b*sqrt(d),abs(b*sqrt(d),-abs(c*sqrt(d),abs(c*sqrt(d);%分別表示xyz軸坐標的刻度范圍endaxis equ

4、al;%將坐標軸改成相等比例xlabel('x'),ylabel('y'),zlabel('z')hold offfigure四、結果與分析1、 結果a=3 b=4 c=5 d=1只有一個量增大: a增大a=3.5 b=4 c=5 d=1a=4 b=4 c=5 d=1a=5 b=4 c=5 d=1a=6 b=4 c=5 d=1 d增大a=3 b=4 c=5 d=2兩個量同時增大: a、b同時增大a=6 b=7 c=5 d=1a=6 b=8 c=5 d=1 a、d同時增大a=3.5 b=4 c=5 d=2三個量同時增大: a、b、c同時增大a=4

5、b=5 c=6 d=1a=6 b=8 c=10 d=1 a、b、d同時增大a=4 b=5 c=5 d=2四個量同時增大:a=5 b=6 c=7 d=3a=6 b=8 c=10 d=2特殊情況:a=3 b=3 c=3 d=1a=4 b=4 c=4 d=1a=3 b=3 c=3 d=22、 分析當d>0,a,b,c都取實數(shù)時,若a=b=c,x2a2+y2b2+z2c2=d是半徑為ad的球面,當d增大或a,b,c同時增大相同倍數(shù)時,x2a2+y2b2+z2c2=d仍是球面,但是其半徑增大,當a,b,c改變不同的倍數(shù)時,x2a2+y2b2+z2c2=d變成橢球面;若a,b,c不相等,x2a2+y

6、2b2+z2c2=d是在x,y,z軸上截距分別為ad, bd, cd的橢球面,當a,b,c,d獨立增大或同時增大時,都會引起橢球面在x,y,z軸上含有相應參量的截距增大,從而使橢球的形狀發(fā)生改變。五、體會記得以前參加數(shù)學建模大賽時,因為不會使用matlab而無法建模,碰到很多困難,最終半途而廢。一直都希望能夠有機會學習一下matlab,幸運的是,仿真平臺與工具應用實踐這一課程讓我對matlab有了初步的a了解。所以,首先,我要感謝學校能夠給我提供此次學習機會,感謝老師悉心的教育和指導。在學習matlab的過程中,由于以前有一定的的C語言和VFP基礎,進展得還是比較順利,但是也還是遇到了很多的困

7、難,主要原因是自己知識面太窄,喜歡根據(jù)以前程序言語的規(guī)則來推測matlab語言,不夠細心。修修改改不知道有多少遍,程序終于可以正常運行了,很激動很有成就感,卻也有很多無奈。我知道,自己做得還不夠好,有的語言運用不夠恰當,運用的很多的程序算法都是知其然不知其所以然,對于整道題解法的邏輯思路也不是很完美,但是我以后一定會更加努力學好matlab的。另外,我很喜歡matlab,覺得它用起來很方便,和其他語言的兼容性比較強,圖像處理功能很強大。我覺得,這次實訓就像是一把鑰匙,為我打開了一扇思維和能力的大門,我一定會從中看到更多更多的風景?!凹埳系脕斫K覺淺,絕知此事要躬行 ”,以此自勉。附件(源程序)a

8、=input('a='); b=input('b='); c=input('c=');d=input('d='); xgrid=linspace(-abs(a*sqrt(d),abs(a*sqrt(d),50); ygrid=linspace(-abs(b*sqrt(d),abs(b*sqrt(d),50); x,y=meshgrid(xgrid,ygrid); z=c*sqrt(d-y.*y/b2-x.*x/a2);u=1; z1=real(z); for k=2:49; for j=2:49if imag(z(k,j)=0z1(k,j)=0;endif all(imag(z(k-1:k+1,j-1:j+1)=0z1(k,j)=NaN;endendendsurf(x,y,z1),hold onif u=1z2=-z1;surf (x,y,z2); axis(-abs(a*sqrt(d),abs(a*sqrt(d),-abs(b*sqrt(d),. abs

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論