在科學(xué)計(jì)算中的應(yīng)用(完整版)
【正文】 實(shí) MATLAB提供了求取階乘的函數(shù) factorial(),其核心算法為 prod(1:n),從結(jié)構(gòu)上更簡(jiǎn)單、直觀(guān)，速度也更快。 % 求出各個(gè)點(diǎn)上的函數(shù)值 plot(x,y) plot(x,y,??)%當(dāng)用多個(gè)要求的時(shí)候可以用單引號(hào)把這些要求放在一起?？捎?set( )函數(shù)設(shè)置對(duì)象的屬性，用 get( )函數(shù)獲得對(duì)象的某個(gè)屬性。 % 先生成繪圖用數(shù)據(jù) subplot(2,2,1)%將幾幅圖放在同一個(gè)畫(huà)面上 , stairs(t,y) % 分割窗口，在左上角繪制階梯曲線(xiàn) subplot(2,2,2), stem(t,y) % 火柴桿曲線(xiàn)繪制 subplot(2,2,3), bar(t,y) % 條型圖繪制 subplot(2,2,4), semilogx(t,y) % 橫坐標(biāo)為對(duì)數(shù)的曲線(xiàn) 隱函數(shù)繪制及應(yīng)用 ? 隱函數(shù) 例： ezplot(?x^2 *sin(x+y^2) +y^2*exp(x+y)+5*cos(x^2+y)?) ％ x自選 ezplot(39。 hold on%表示保持原來(lái)的曲線(xiàn)，然后在此基礎(chǔ)上疊加其他圖層 %。 contour3(x,y,z,30) ％三維等高線(xiàn)圖， 30等高線(xiàn)條數(shù) ? 例：試?yán)L制出二元函數(shù) [x,y]=meshgrid(2:.1:2)。zlim39。% 俯視圖 subplot(222), surf(x,y,z), axis([0,31,0,31,0,1])。 % 側(cè)視圖 % subplot(223), surf(x,y,z), axis([0,31,0,31,0,1])。 n=2。 Warning: Divide by zero. Warning: Divide by zero. (Type warning off MATLAB: divideByZero to suppress this warning.) surf(x,y,z), shading flat %修飾其顯 示形式 2 2 2 211( , )( 1 ) ( 1 )z f x yx y x y? ? ?? ? ? ? xx=[2:.1:, ::, ::, ::, ::2]。 ? 其他函數(shù)，光照下 surfl( ),等高線(xiàn) surfc( )，瀑布型 waterfall( ) ? 等高線(xiàn)繪制 contour( ), contour3( ) 例 :Butterworth 濾波器 [x,y]=meshgrid(0:31)。,[10 10]) 三維圖形繪制 ? 三維曲線(xiàn)繪制 ? stem3（三維火柴桿型曲線(xiàn)） , fill3（三維填充圖形） , bar3（三維直方圖）等。 其他二維圖形繪制語(yǔ)句 二維條形圖 bar(x,y)、羅盤(pán)圖 pass(x,y)、羽毛狀圖feather(x,y)、直方圖、極坐標(biāo)圖、階梯圖形、 x半對(duì)數(shù)圖、彗星狀軌跡圖、誤差限圖形、二維填充圖、對(duì)數(shù)圖、磁力線(xiàn)圖、火柴桿圖、 y半對(duì)數(shù)圖。 % 以變步距方式構(gòu)造自變量向量， 注意紅色數(shù)字，應(yīng)該比前面的數(shù)字大一點(diǎn)。 for i=1:length(varargin), a=conv(a,varargin{i})。)。 for i=1: n for j=1:m A1(i,j)=1/(i+j1)。 end [m1,s1]=findsum(145323) m1 = 539 s1 = 145530 ? 無(wú)需修改程序 例： ? 若只給出一個(gè)輸入?yún)?shù)，則會(huì)自動(dòng)生成一個(gè)方陣 ? 在函數(shù)中給出合適的幫助信息 ? 檢測(cè)輸入和返回變量的個(gè)數(shù) ? edit myhilb , 1 /( 1 )ijh i j? ? ?function A=myhilb(n, m) %產(chǎn)生 A=MYHILB(N,M)或 A=MYHILB(N)； if nargout1, error(39。 end end i i = 141 開(kāi)關(guān)結(jié)構(gòu) ?和 C 語(yǔ)言的區(qū)別 – 當(dāng)開(kāi)關(guān)表達(dá)式的值等于某表達(dá)式，執(zhí)行該語(yǔ)句后結(jié)束該結(jié)構(gòu)，不用 break – 當(dāng)需要在開(kāi)關(guān)表達(dá)式滿(mǎn)足若干個(gè)表達(dá)式之一時(shí)執(zhí)行某一程序段，則用單元形式 （用大括號(hào)把這些表達(dá)式括起來(lái)，用逗號(hào)分隔） – otherwise 語(yǔ)句，不是 C語(yǔ)言中的 default（但與之等價(jià)） – 程序的執(zhí)行結(jié)果和各個(gè) case順序無(wú)關(guān) – case 語(yǔ)句中條件不能重復(fù)，否則列在后面的條件將不能執(zhí)行 ? 全新結(jié)構(gòu)（首先試探性執(zhí)行語(yǔ)句 1，若執(zhí)行過(guò)程中有錯(cuò)，將錯(cuò)誤信息賦給保留的 lasterr變量，并終止這段語(yǔ)句的執(zhí)行，轉(zhuǎn)而執(zhí)行語(yǔ)句 2。 for i=1:100000, s=s+1/2^i+1/3^i。 i=1。 isprime(A) ％若向量 A中某個(gè)整數(shù)值為質(zhì)數(shù)，則相應(yīng)位置為 1，其他為零。 floor(A) % 向 inf 方向取整 ans = 1 0 1 2 1 2 ceil(A) % 向 +inf 方向取整 ans = 0 1 2 3 2 1 round(A) %取最近的整數(shù) ans = 0 1 1 2 1 2 fix(A) %向 0 的方向取整 ans = 0 0 1 2 1 1 ? 例： 3x3 Hilbert 矩陣，試用 rat() 函數(shù)變換 A=hilb(3)。[i,j] ％顯示行標(biāo)，列標(biāo) ans = A= 3 1 1 2 3 2 2 4 5 6 3 2 7 8 0 2 3 all(A=5) ％某列元素全大于或等于 5時(shí)，相應(yīng)元素為 1，否則為 0。 A = B2 =
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1