【正文】 輸入的是素數(shù)!39。while i if rem(m,i)==0 disp(39。)endif m=0|mround(m)~=0 warning(39。2。您所求解的方程組的解是：39。x=A\b。)。endnn=m。while hheps s2=s2+hh。endh1=h。h=H。)else t=[]。 otherwise money=p*q*。 case q=3000amp。 case q=1000amp。linear39。Method is nearest39。Method is cubic39。Method is linear39。else f=exp(2*x)。elseif 0=xamp。)。您描述的可能是雞雞雞！39。 ji_ya(s)（b）格式2if 邏輯表達式 執(zhí)行語句1else 執(zhí)行語句1end例4：根據(jù)嘴型判斷是否是雞。您描述的是鴨鴨鴨!39。)。function ji_ya(s)if s==39。%功能：求一個矩陣所有元素的平均值function average=junzhi(A)s=sum(sum(A))。函數(shù)M文件需要接受參數(shù)的輸入和輸出，只處理輸入?yún)?shù)傳遞的數(shù)據(jù)，并把處理結(jié)果作為函數(shù)輸出參數(shù)返回給指定接收變量。subplot(2,2,3)[c,h]=contour(z,6)。腳本M文件matlab代碼按順序組成的命令序列集合，不接受參數(shù)的輸入和輸出，一般用于實現(xiàn)一個相對獨立的功能，如某種分、繪圖或求解等。sin39。39。)（b）使用insert操作添加圖例：（a）使用函數(shù)legendx = pi:pi/20:pi。plot(years,value)xlabel(39。添加直線和箭頭：例：指出高斯曲面的最高點和最低點%先畫出圖形[X,Y,Z]=peaks(30)。正弦曲線39。ko39。x=2*pi::2*pi。string39。正弦曲線39。sin(x)39。subplot(2,2,4)[c,h]=contourf(z,8)。subplot(2,2,1)contour(z)subplot(2,2,2)[c,h]=contour(z,[,])。Z = X.*exp(X.^2Y.^2)。clabel(c,h)。})figure(3)pie3(x)%三維餅狀圖（4）等值線圖命令：contour，clabel，contourf例10：繪制z=peaks的等值線圖z=peaks。,39。figure(1)pie(x)figure(2)pie(x,{39。data=[，，]。%建立三角形網(wǎng)格figure(1)trimesh(tri,x,y,z)figure(2)trisurf(tri,x,y,z)三、特殊圖形（1）直方圖命令：hist例7:%均勻分布直方圖x=rand(500,1)。Z = X.*exp(X.^2Y.^2)。z=x.*y。z=x.^2+y.^2。x=exp(*t).*cos(t)。 x=2*cos(t)。（2）數(shù)據(jù)曲線圖命令：plot(y),plot(x,y)例2：，，請畫出銷售收入增長趨勢圖。sqrt(x)139。sin(x)39。 mean(R39。 x=[93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55]。 %產(chǎn)生均值為10的1行10000列的泊松分布隨機數(shù)； mean(y1) %檢驗均值 var(y1) %檢驗方差 y2=random(39。,0,1,1,10000)。 %在0到1之間產(chǎn)生200個均勻分的隨機數(shù) mean(a1) %計算其期望 var(a1) %計算其方差 std(a1) %計算其標準差（4）協(xié)方差、相關(guān)系數(shù)命令：cov()%協(xié)方差corrcoef()%相關(guān)系數(shù)例9： cov(a1) %計算a1的協(xié)方差（方差） cov(a1,a2) %計算a1和a2的協(xié)方差 corrcoef(a1,a2) %計算a1和a2的相關(guān)系數(shù) b1=[1 2 3 4 5 6 7]。plot(cs2)。 sum(x)例6：求每個學生的總成績；例7：cs1=[3 4 2 5 6 3]。english。四、數(shù)理統(tǒng)計（1）數(shù)據(jù)排序、最值、均值、查找①排序命令：sort(x)%向量x按曾序排列；[y,ind]=sort(x)%向量x增序排列為y向量，并輸出下表向量；例1：已知某5個同學的數(shù)學、外語、計算機成績分別為math=[86 78 82 91 85]。,39。x39。x*Dy+2*yexp(x)=039。,39。x39。) ans = x^2*t+C1由于系統(tǒng)默認的自變量是t，故在求解時將x看作了常數(shù)。（4）dsolve(‘equation’,’condition’,’v’) %給出微分方程初值問題的解，表示為v的函數(shù)。具體算法實現(xiàn)將在第5講中講述。1,1]。 [r,val]=fsolve(myfun,x0)結(jié)果為：r = val = * 也可以用fsolve函數(shù)求解矩陣方程；例8：求矩陣方程的一個根，初始解為x0=[1,1。（2）非線性方程組①直接解法直接采用fsolve函數(shù)進行求解。 B=[2。 X1=pinv(A)*B %利用偽逆法 X2=lsqnonneg(A,B) %利用最小二乘法⑤欠定方程AX=B的通解matlab在求解欠定方程組時，只給出其中一組解，不會給出方程組的通解，但是可以用rref命令化為行階梯（查看增廣矩陣的秩和系數(shù)矩陣的秩）來分析方程是否有通解。 B=[1。④超定方程的近似解法偽逆法：X=pinv(A)*B最小二乘法：X=lsqnonneg(A,B)例5：利用偽逆法和最小二乘法求解超定方程組 A=[1 2 3。1]。36 192 180。R X=B，則X= R\(R39。2。 。2]。14 4 5]。7]。1 3 7。 b=[1。,[0 2])（2）其它迭代解法如二分法、牛頓切線法、不動點迭代法、弦截法等，這些算法的實現(xiàn)將在第五講應(yīng)用舉例中講述。,[1 2]) x=fzero(39。,) x=fzero(39。x*exp(2*x)839。x*3^x139。z_y=int(z,y,x^2,x)。int(x^3*atan(x),0,1)例4:計算廣義積分syms x。dx=diff(cos(a*x^21),x)。diff(cos(a*x^21),x)例3 求關(guān)于a的導數(shù)syms a x。limit((x^21)/(x1),x,1)例2:求syms x。)。x^22*x339。)syms x yv=x^22*y^2+x*y2*uv+2（2）符號表達式的四則運算syms x y。1 1 1 0] %建立4行4列的矩陣d=[1 2 3 4。6。第二講：Matlab基本運算一、數(shù)值類型（1）變量①命名規(guī)則(a):變量名以英文字母開始，由字母、數(shù)字、下劃線組成，不能有空格和標點符號；(b):區(qū)分大小寫；(c):變量名長度不能超過31位，超過部分將被忽略；(d):某些常量也可作變量名，如虛數(shù)單位i。