清華數(shù)學實驗第一章 MATLAB 使用入門 III

1、數(shù)據(jù)顯示格式字符串變量操作圖形文件的輸入/輸出寶石切割實驗MATLAB 使用入門 III計算機浮點數(shù)表示法6.378137 103 6.378137 e +003 2.99792458 105 2.99792458e +005 例1.24近似計算地球表面積,并以不同格式輸出數(shù)據(jù) R=6378.137; S=4*pi*R2 format long,S format bank,Sformat short 短格式 5.1121e+008format long 長格式 5.112078933958109e+008format bank 銀行格式 511207893.40數(shù)據(jù)顯示格式MATLAB內(nèi)部常數(shù)

2、常 數(shù)返 回 值ans默認變量名，保存最近的結(jié)果eps浮點數(shù)相對精度(2.2204e-016)realmax最大浮點數(shù)(1.7977e+308)realmin最小浮點數(shù)(2.2251e-308)pi圓周率（3.1416）i,j虛數(shù)單位inf無限大NaN不合法的數(shù)值，非數(shù)值數(shù)據(jù)顯示格式例1.25立方倍積問題:求作一立方體,使其體積為已知立方體的二倍.如果已知立方體體積為V0=1,要作的立方體體積為V=2,則所求立方體高度為 21/3在命令窗口中直接使用命令 h=vpa(2(1/3),40)數(shù)據(jù)結(jié)果為 h = 1.259921049894873190666544360283296555281MAT

3、LAB具有一種稱為VPA類型的數(shù)據(jù)顯示格式.VPA的全稱是可變精度浮點算法（Variable precision arithmetic）.顯示S的D位數(shù)方法如下R=vpa(S，D)字符串變量：將字符串用單引號括起來賦值給變量拼接字符串命令: strcat(T1,T2,T3)將整數(shù)轉(zhuǎn)換成字符命令: int2str( )例1.26 三個名人Euler，Elizabeth，Plato職業(yè)分別是mathematician，movie star，philosopher，編寫程序正確聯(lián)接他們的名字和職業(yè)并輸出n1=Euler;n2=Elizabeth;n3=Plato；p1=mathematician;p

4、2=movie star;p3=philosopher；s1=strcat(n1,-,p1), s2=strcat(n2,-,p2),s3=strcat(n3,-,p3)例1.27 生肖問題08年是鼠年還是狗年？n=input(input n:=);S=鼠牛虎兔龍蛇馬羊猴雞狗豬;k=rem(n-4,12)+1;s=S(k);s=strcat(int2str(n), 年是, s,年)input n:=2006 s= 2006年是狗年Input n:=2008 s=2008年是鼠年復變函數(shù)圖形設有復平面上單位圓域內(nèi)變化的變量 以 u = z 的實部函數(shù)繪圖并輸出圖形文件r=linspace(0,1

5、,20);theta=linspace(-pi,pi,25);z=r*exp(i*theta);x=real(z);y=imag(z);mesh(x,y,x),hold oncolormap(0 0 1)mesh(x,y,-ones(size(x)axis off輸出圖形文件1.在圖形窗口用鼠標點擊窗口左上方的菜單欄“file”，選擇下拉菜單中的“Export”. (在高版本中選擇“Save As”)2.在對話框中選擇文件類型“bmp”格式，將圖形文件命名為riman. 寶石加工廠需要將石料內(nèi)已知位置上的一塊精品切割出來.石料尺寸為191412(cm3) ,需要切割的精品尺寸為542(cm3)

6、,且二者的左側(cè)面、前面、底面相互平行,距離分別為6cm,7cm,9cm。寶石加工問題：切割加工費用:垂直切割3元/cm2 ,水平切割5元/cm2。試確定切割的工作流程,使切割的總費用最少。 (為了減少旋轉(zhuǎn)刀具的次數(shù)，要求同向切割連續(xù)兩次后再旋轉(zhuǎn)刀具)19 14 125 4 2S1=2(1412)5 14 125 4 25 4 125 4 25 4 25 4 2S2=2(512)S3=2(54)全部方案比較石料尺寸(cm)石料:長寬高=a1a2a3；精品:長寬高= b1b2b3；六種方案的切割面積計算a2a3 b1a3 b1b2a2a3 b1a2 b1b3a1a3 b2a3 b1b2a1a3 a

7、1b2 b2b3a1a2 a2b3 b1b3a1a2 a1b3 b2b3(1,2,3)(1,3,2)(2,1,3)(2,3,1)(3,1,2)(3,2,1)1 左、右2 前、后3 上、下a1=19;a2=14;a3=12;b1=5;b2=4;b3=2;p=1 2 3;1 3 2;2 1 3;2 3 1;3 1 2;3 2 1;f1=3;3;5;f2=3;5;3;f3=5;3;3;A=a2*a3,b1*a3,b1*b2; a2*a3,b1*a2,b1*b3; a1*a3,b2*a3,b1*b2; a1*a3,a1*b2,b1*b3; a1*a2,a2*b3,b1*b3; a1*a2,a2*b3,

8、b2*b3;q(1)=A(1,:)*f1;q(2)=A(2,:)*f2;q(3)=A(3,:)*f1;q(4)=A(4,:)*f2;q(5)=A(5,:)*f3;q(6)=A(6,:)*f3;qmin,index=min(q); operate=p(index,:)expense=2*qminMATLAB程序(mlab0.m)練習與思考題1.設地球半徑為6400(km),試分析“坐地日行八萬里”這句話是否有道理2.如何顯示出圓周率的100位數(shù)字？3.分析復變函數(shù)計算中z=r*exp(i*theta);x=real(z);y=imag(z);的功能作用4. 2007年農(nóng)歷是丁亥年,排列天干(甲乙丙丁戊己庚辛壬癸)/地支(子丑寅卯辰巳午未申酉戍亥)，用字符串數(shù)組編寫程序,要求輸入年份,輸出農(nóng)歷紀年5. 下面程序是繪拋射曲線族,分析每行命令功能alpha=pi*(1: