再舉一個簡單的M文件實(shí)例.

1、8再舉一個簡單的 M 文件實(shí)例：sunooy 發(fā)表于 2006-6-22 0:07:00日本數(shù)學(xué)家角谷提出一個著明的猜想：任何一個大于 1 的自然數(shù)經(jīng)過如下的有限步運(yùn)算，最后結(jié)果為 1(1) 如果是偶數(shù)，將其除以 2(2) 如過是奇數(shù)，將其乘以 3,再加上 1例如 n = 1010/2=55*3+1=1616/2=88/2 = 44/2 = 22/2 = 1根據(jù)上面提供的猜想：建立 M 文件，以便方便的驗(yàn)證對不同自然數(shù)，角谷猜想 是否正確。程序可以參考這里： user2/su nooy/archives/2006/1275717.shtml第一步：建立 M 文件% M file, jiaogu

2、.m%驗(yàn)證對不同的自然數(shù) n,角谷猜想是否成立n=input(Type n:);%從鍵盤輸入 n,并賦值給 ni=1; %標(biāo)記運(yùn)算的步數(shù)f=n; %存放每次運(yùn)算后的 n 的值 while n=1 %當(dāng) n 不等于 1 時就循環(huán)，直到 n= 1時為止。if mod(n,2)=0n=n/2;else n=n*3+1;endi=i+1;f(i)=n;endf%輸出每次運(yùn)算后的 n 的值組成的集合I=i-1 %輸出運(yùn)算的步數(shù),f 集合中會有 i 個元素，但實(shí)際上是只經(jīng)過 i-1 步的運(yùn) 算n %輸出最后一步所得 n 的值，顯然，一般 n = 1，除非猜想不成立時，n 會 取到其它的值！8第二步：保存

3、M 文件: 進(jìn)入 File 菜單，打開 New 菜單，選擇 M file，在保存，注意為了直觀最好將 文件名保存為 jiaogu第三步：調(diào)用 M 文件 在 matlab 窗口輸入調(diào)用命令 jiaogu 回車后顯示 Type n:立即給n 賦不同的 值，馬上得到結(jié)果。比如：n= 10 時，經(jīng)過 6 步可將 10 化為 1Type n:10 f =105I =6n = 1 n= 55 時，經(jīng)過 20 步可將 55 化為 1f = 5516683250447142 .42 n= 100時，經(jīng)過 25 步可將 100 化為 1f =10050257638195829884422113417522613

4、402010 5168421I =25n =1還可給 n 取更多的值或更大的值，總之，建立 M 文件夠，驗(yàn)證工作變得十分方 便了！特別說明：如果把 while n=1 改為 while n=1 &n0 程序會更完美。寫成 while n=1 表 面上看是沒有錯的。但因?yàn)榻枪炔孪氘吘故莻€猜想，我們無法驗(yàn)證所有的數(shù)， 如果存在某個特別大的數(shù)使猜想不成立，它會從n 逐漸的變小，直到負(fù)的無窮大，而不能被化為 1。這樣程序就會進(jìn)入死循環(huán)！所以應(yīng)該寫成while n=1 &n0，這就保證如果有某個數(shù)不滿足角谷猜想，一旦發(fā)現(xiàn)它不是正數(shù)了，就停 止循環(huán)，從而避免了死循環(huán)。如何建立和調(diào)用 M 文

5、件(M-File)呢?(原創(chuàng))sunooy 發(fā)表于 2006-6-21 21:17:001. 建立 M 文件的必要性：對于常用的或很長的程序，如果每次都從頭輸入，則非常費(fèi)時間，也容易出錯將 這些程序做成 M 文件就方便多了 -只要輸入自定義的命令，就可以方便的調(diào)用, 編程的效率大大提高！2. 如何建立 M 文件(1)在記事本(.txt)或 maltab 命令窗口中編輯好 M 文件16842112537618891（為節(jié)省空間，中間的數(shù)略去了）8例:編寫一個輸入半徑，計算圓的面積和球的體積的 M 文件.前兩行是文件頭，一般說明文件的 類型(M-file),文件名(tutex2),和文件的功能 (

6、輸入半徑,計算圓的面積和球的體積),后面是文件體，主要是解決特定問題的程 序代碼% M-file, tutex2.m%輸入半徑，計算圓的面積和球的體積r = input(Type radius:); %從鍵盤輸入 Type radius 的值，并將它賦值給 r.area=pi*rA2;volume=(4/3)*pi*rA3;fprintf(The radius is %12.5fn,r)fprintf(The area of a circle is %12.5fn,area)fprintf(The volume of a sphere is %12.5fn,volume)將上面的代碼保存到工作

7、空間里.保存方法是：File/New/F-ile/保存3. 如何調(diào)用 M 文件？在 matlab 工作窗口直接鍵入 M 文件名，即可調(diào)用.以上面的程序?yàn)槔?tutex2 %鍵入文件名，回車Type radius:%提示輸入半徑,輸入半徑，比如 13.5,回車 顯示如下Type radius:13.5The radius is 13.50000The area of a circle is 572.55526The volume of a sphere is 10305.99470閱讀全文(61) |回復(fù)(0) |引用通告(0) |編輯MATLAB 中有趣的演示實(shí)例sunooy 發(fā)表于 200

8、6-5-7 18:04:008來源:http:/ 為提高讀者對 MATLAB和 Simulink 的興趣，特舉部分以供參考。具體如下，運(yùn)行的時候只要將： “”前面的代碼復(fù)制到 MATLAB 中就可以了，隨之會出現(xiàn)各種各樣的演示實(shí)例，對初學(xué)者幫助 不小哦。平面與立體繪圖graf2d : XY 平面繪圖（火柴棒） graf2d2 : XYZ 立體繪圖（切片）hndlgraf :平面顯示線型處理窗口及命令演示 hndlaxis :平面顯示處理窗口及命令演示 graf3d :立體顯示處理窗口及命令演示復(fù)雜函數(shù)的三維繪圖cplxdemo :復(fù)雜的 XYZ 立體圖形等高線繪制quivdemo :等高線箭頭

9、顯示動畫lorenz : Lorenz 吸引子動畫顯示電影vibes : L-形薄膜振動Fourier 變換sshow sun spots :太陽黑點(diǎn)數(shù)據(jù)的傅里葉分析 fftdemo :分析噪聲序列中兩組數(shù)據(jù)的相關(guān)度數(shù)據(jù)擬合sshow fitdemo :顯示非線性數(shù)據(jù)擬合過程census :預(yù)測世界人口 spline2d :樣條擬合稀疏矩陣 sshow sparsity :降階游戲xpbombs :仿 Windows 系統(tǒng)自帶的掃雷游戲 life :生命發(fā)展游戲三維效果圖klei n1 :膚色三維效果圖 tori4 :四個首尾相接的圓環(huán) spharm2 : 球形和聲 cruller :類似油餅

10、的東西 xpklein : Klein 瓶 bottle modes : L-形薄膜的 12 中模態(tài) logo :MATLAB 的 Logo xpquad :不同比例的巴爾體超四方體 truss :二維桁架的 12 個模模態(tài) tr*el :旅行商問題動畫演示 wrldtrv :在地球儀上演示兩地間的飛行線路makevase :通過點(diǎn)擊鼠標(biāo)來制作花瓶 xpsound :聲音樣本分析funfuns :綜合了找零點(diǎn)，最小化和單輸入函數(shù)積分功能 sshow e2pi : eApi 或者 piAequake :地震波可視化 penny :便士可視化 imageext :改變圖像的映射顏色earthmap

11、 :地球儀優(yōu)化工具箱 bandem :香蕉最優(yōu)化展示 expo-style bananaoptimization sshow filtdem :濾波效果演示 filter effect demo sshow filtdem2 :濾波設(shè)計演示 filter design demo cztdemo : FFT 和 CZT （兩種不同類型的 Z-變換算法）phone :演示電話通聲音的時間與頻率的關(guān)系 sigdemol :離散信號的時頻圖，可用鼠標(biāo)設(shè)置 sigdemo2 :連續(xù)信號的時頻圖，可用鼠標(biāo)設(shè)置 filtdemo :低通濾波器的交互式設(shè)計 moddemo :聲音信號的調(diào)制 sosdemo :

12、數(shù)字濾波器的切片圖神經(jīng)網(wǎng)絡(luò)工具箱neural :神經(jīng)網(wǎng)絡(luò)模塊組 firdemo :二維 FIR 濾波器 nlfdemo :非線性濾波器dctdemo : DCT 演示 mlpdml :利用多層感知器神經(jīng)網(wǎng)絡(luò)擬合曲線動畫mlpdm2 :利用多層感知器神經(jīng)網(wǎng)絡(luò)進(jìn)行 XOR 問題運(yùn)算模糊邏輯工具箱invkine :運(yùn)動逆問題8juggler :跳球戲法 femdemo : FCM slcp :類似倒立擺動畫slcp1 :類似倒立擺動畫 cart and a varying pole slcpp1 :類似倒立擺動畫，有兩個擺，一個可以變化 sltbu :卡車支援 slbb :類似于翹翹板simplif

13、y 把數(shù)學(xué)表達(dá)式化為最簡形式sunooy 發(fā)表于 2006-6-29 23:21:00舉一例：y=1/(xA2+1)A(1/2)-xA2/(xA2+1)A(3/2);這個式子看起來比較煩瑣，不妨用用 simplify 命令對其化簡，看它的廬山真面目到底是什么? y=1/(xA2+1)A(1/2)-xA2/(xA2+1F(3/2); simplify(y)ans =1/(xA2+1)A(3/2)原來是這個，看來 simplify 化簡的效率挺高的。MATLAB 的取整函數(shù):fix(x), floor(x) :,eeil(x) , round(x)sunooy 發(fā)表于2006-5-12 0:16:

14、00(1)fix(x):截尾取整. fix( 3.12 -3.12)ans =3-3floor(x):不超過 x 的最大整數(shù).(高斯取整) floor( 3.12 -3.12)ans =3-4(3) ceil(x):大于 x 的最小整數(shù) ceil( 3.12 -3.12)ans =84-3(4) 四舍五入取整 round(3.12 -3.12)ans =0 round(3.12 -3.12)ans =3-3在 matlab 繪圖中加入圖例說明的辦法sunooy 發(fā)表于 2006-8-1 23:06:00例如:在同一個坐標(biāo)系中畫了兩條不同的曲線，一條線形是空心圓O一條線形是 +,我們希望在圖上做出這樣的圖例說明：O曲線 1+ -曲線 2一般的方法是 繪制好曲線后,在圖形窗口選擇選則 In sert,然后選擇 Lege nd, 這時圖形上就會出現(xiàn)：Odata1+ -data