matlab程序設計(2)(編輯修改稿)
2025-02-08 18:53 本頁面
【文章內容簡介】 直至不滿足為止。這時將結束 for語句的執(zhí)行，而繼續(xù)執(zhí)行 for語句后面的語句。 例 已知求 y的表達式，當 n=100時 ,求 y的值。 211niyi?? ?程序如下： y=0。 n=100。 for i=1:n y=y+1/i/i。 end y 在實際 MATLAB編程中，為提高程序的執(zhí)行速度，常用向量運算來代替循環(huán)操作： n=100。 i=1:n。 f=1./i.^2。 y=sum(f) 例 求定積分 3 sin( ) d6xy e x x? ?????程序如下： a=0。b=3*pi。 n=1000。 h=(ba)/n。 x=a:h:b。 f=exp(*x).*sin(x+pi/6)。 for i=1:n s(i)= (f(i)+f(i+1))*h/2 end s=sum(s) 事實上，MATLAB提供了有關數(shù)值積分的標準函數(shù)，實際應用中可直接調用這些函數(shù)求數(shù)值積分。 注意： for語句更一般的格式： for 循環(huán)變量 =矩陣表達式 循環(huán)體語句 end 說明： 執(zhí)行過程是依次將矩陣的各列元素賦給循環(huán)變量，然后執(zhí)行循環(huán)體語句，直至各列元素處理完畢。實際上， “ 表達式 1:表達式 2:表達式 5” 是一個僅為一行的矩陣 (行向量 )，因而列向量是單個數(shù)據(jù)。 例 已知 5個學生 4門功課的成績，求每名學生的總成績。 s=0。 a=[65,76,56,78。 98,85,74,85。 76,67,78,79。 98,58,42,75。 67,89,76,87]。 for k=a s=s+k。 End disp(s39。)。 while語句 1)格式： while (條件 ) 循環(huán)體語句 end 2)說明： (1)執(zhí)行過程：若條件成立 ,則執(zhí)行循環(huán)體語句，執(zhí)行后再判斷條件是否成立 ,如果不成立則跳出循環(huán) 。 (2)與循環(huán)結構相關的還有一個 break語句，當在循環(huán)體內執(zhí)行到該語句時，程序將跳出循環(huán)。該語句一般與 if語句配合使用。 例 已知求 y的表達式 ,當 n=100時 ,求 y的值。 211niyi?? ?y=0。 i=1。 while 1 f=1/i/i。 y=y+f。 if i==100 break。 end i=i+1。 end y 循環(huán)嵌套 1)如果一個循環(huán)結構的循環(huán)體又包括一個循環(huán)結構，就稱為循環(huán)的嵌套，或稱為多重循環(huán)結構。 2)要特別注意內、外循環(huán)之間的關系，以及各語句放置的位置，不要搞錯。 3)可跳出循環(huán)的指令有 break和 continue 在執(zhí)行 keyboard指令后需要返回正在執(zhí)行程序的指令是 return 。 例 用篩選法求某自然數(shù)范圍內的全部素數(shù)。 程序如下： m=input(39。m=39。)。 p=2:m。 for i=2:sqrt(m) n=find(rem(p,i)==0amp。p~=i)。 p(n)=[]。 end p 函數(shù)文件 函數(shù)文件的基本結構 1)函數(shù)文件由 function語句引導 。 2)基本結構： function 輸出形參表 =函數(shù)名 (輸入形參表 ) 注釋說明部分 函數(shù)體語句 (1)函數(shù)文件名 : 函數(shù)文件名與函數(shù)名可以不相同 。當兩者不同時， MATLAB將忽略函數(shù)名而確認函數(shù)文件名，故調用時使用函數(shù)文件名。 (2) 注釋說明部分： ①緊隨函數(shù)文件引導行之后以 %開頭的第一 注釋行； ②第一注釋行及之后連續(xù)的注釋行； ③與在線幫助文本相隔一空行的注釋行。 3)說明： 例 編寫函數(shù)文件求半徑為 r的圓的面積和周長。 function [s,p]=fcircle(r) %CIRCLE calculate the area and perimeter of a circle of radil r %r 圓半徑 %s 圓面積 %p 圓周長 s=pi*r*r。 p=2*pi*r。 將以上函數(shù)文件以文件名 c:\matlab\work下，
畢業(yè)設計相關推薦
鄂ICP備17016276號-1