程序設(shè)計(jì)上機(jī)

利用Matlab提供的randn函數(shù)生成符合正態(tài)利用Matlab提供的randn函數(shù)生成符合正態(tài)分布隨機(jī)矩陣A，并計(jì)算：A各列元素的均值和標(biāo)準(zhǔn)方差A(yù)的最大元素和最小元素（3）A每行元素的和以及全部元素之和（4）A的每列元素按升序、每行元素降序排序3產(chǎn)生均值為0產(chǎn)生均值為0，方差為1的標(biāo)準(zhǔn)正態(tài)分布的隨機(jī)矩陣A4>>A=randn(10,5)A= （1）A各列元素的（1）A各列元素的均值和標(biāo)準(zhǔn)方差求數(shù)據(jù)序列平均值的函數(shù)是返回一個(gè)行向量，其第i個(gè)元素是Ai列的算術(shù)平5>>mean(A,1)ans= - ->>A=randn(10,5)A=- - - - - - - - - - - - - - - - - -- - - ->>mean(A,2)ans=----（1）A各列元素的均值和標(biāo)準(zhǔn)方差x1,x2對于N個(gè)元素的數(shù)據(jù)序，標(biāo)準(zhǔn)方差的（1）A各列元素的均值和標(biāo)準(zhǔn)方差x1,x2對于N個(gè)元素的數(shù)據(jù)序，標(biāo)準(zhǔn)方差的計(jì)算公式dim=1flag=0按式1flag=1按式2 NN或Nix21Ni默認(rèn)：flag=0其中xNNxii在Matlab中，計(jì)算標(biāo)準(zhǔn)方差的函數(shù)是std(A,flag,std函數(shù)的一般調(diào)用格6（1）A（1）A各列元素的均值和標(biāo)準(zhǔn)方差7>>y2 >>y1 （1）A（1）A各列元素的均值和標(biāo)準(zhǔn)方差8>>y4>>y3實(shí)驗(yàn)隨機(jī)產(chǎn)生100個(gè)200～實(shí)驗(yàn)隨機(jī)產(chǎn)生100個(gè)200～500的隨機(jī)整數(shù)，求這100個(gè)數(shù)的和，平均大于平均數(shù)的個(gè)（3）求最小值、最大（3）把這100個(gè)數(shù)從小到大排9隨機(jī)產(chǎn)生100個(gè)200～500的隨機(jī)產(chǎn)生100個(gè)200～500的隨機(jī)整%產(chǎn)生100個(gè)隨機(jī)forx(i)=fix(rand*(500-%x(i)=round(rand*(500-%求100個(gè)隨機(jī)整數(shù)的平%%-%%-fori=1:100%%-%%-fori=1:100if%把這100個(gè)數(shù)排%YYholdYColumns1throughColumns17throughColumns33throughColumns49throughColumns65throughYColumns1throughColumns17throughColumns33throughColumns49throughColumns65throughColumns81throughColumns97through sumaver=aver_Nansans實(shí)驗(yàn)p(x)axna1實(shí)驗(yàn)p(x)axna1n10pn a1alue=polyval(p,xp5(x)2x53x3p8(x)8x85x4x=1，x=-1，x=3p3(x)3x32xp2(x)x224x3x（2）上列多項(xiàng)式23]、x=4%p5=2x^5+3x^3-%p3=3x^3-%p8=8x^8-%p2=x^2-alue5%p5=2x^5+3x^3-%p3=3x^3-%p8=8x^8-%p2=x^2-alue5p5(1)2x53x34x12alue3p3(1)3x32x1p5=[2030-4p3=[30-2p8=[8000-5003p2=[10-alue8p8(1)8x85x43x28alue2p2(1)x22x%p5=2x^5+3x^3-%p3=3x^3-%p8=8x^8-%p2=x^2-alue5x%p5=2x^5+3x^3-%p3=3x^3-%p8=8x^8-%p2=x^2-alue5 2p5(x)2x53x34x1alue31p3(x)3x32xp5=[2030-4p3=[30-2p8=[8000-5003p2=[10-alue88p8(x)8x85x43x2x=[12;3alue272p2(x)x22求方程用roots命令求方程f=0的根roots=roots(f)%求方程用roots命令求方程f=0的根roots=roots(f)%x0為初始估–––用fzero命令z=fzero(y,x=fzero(F,[a,b])%返回一元函數(shù)F在區(qū)間[a,b個(gè)零點(diǎn)，要求F在區(qū)間端點(diǎn)[x,f,h]=fsole(F,–x返回一元函數(shù)或多元函數(shù)F在自變量x0附近的一個(gè)零f返回對應(yīng)函數(shù)值h回值大于零說明結(jié)果可靠，–8x3108x310（3）求多f(x)3x54x35x2f(x)07.2x（4）已x-0.5075+-0.5075-%p4=x^4+8x^3-10p5=[1800-10];%p4=x^4+8x^3-10grid8x310（3）求多8x310（3）求多的f(x)3x54x35x2f(x)0（4）已 p5=3x^5+4x^3-5x^2-p=[304-5-7.2x-0.3046+-0.3046-cl%p5=3x^5+4x^3-5x^2-y=3*x.^5+4*x.^3-5*x.^2-gridon實(shí)驗(yàn)克萊姆法則解線性方程實(shí)驗(yàn)克萊姆法則解線性方程用求逆矩陣的方法解線性方程ABBAIIB但在MATLAB中，求矩陣的逆非常容易。A的逆矩陣命令“\”求解線性方程組x2y3zx2y3z（5）求解線性方x4y9zx8y27z%%%1112483956A=[123;149;18B=[5-2A1=[523;-249;68A2=[153;1-29;16A3=[125;14-2;18A B56A1 A2 A3 x1x2x3x2yx2y3z（5）求解線性方x4y9zx8y27zA=[123;149;18B=[5-2A B56x=x2yx2y3z（5）求解線性方x4y9zx8y27z解3：A=[123;149;18B=[5-2A B56xx2y3zx2y3z（5）求解線性方x4y9zx8y27z解4：使用solve命令進(jìn)行計(jì)算%%%1112483956symsxy[x,y,z]=solve(‘x+2*y+3*z=5’,‘x+4*y+9*z=-;x=yz=ans 23,-29/2,實(shí)驗(yàn)對于非線性方實(shí)驗(yàn)對于非線性方程組F(X)=0，用fsolve函數(shù)求其數(shù)值fsolve函數(shù)的調(diào)用格––––X為返回的option為最優(yōu)化工具箱的選項(xiàng)設(shè)定，可以使用optimset命令將它們示出來。常用的選項(xiàng)‘iter’表示每步都顯示，‘final’只顯示最終結(jié)果。MaxIter選項(xiàng)：表示最大允許迭代次數(shù)，默認(rèn)為400次。??? （6）求下列方程組在（1，1，1（6）求下列方程組在（1，1，1）附近的解sinxyz2exxyzxyz解先建立函數(shù)文function在給定初值x0=1,y0=1,z0=1，調(diào)用函數(shù)求方程組的解x 將求得的解代回原方程，可以將求得的解代回原方程，可以檢驗(yàn)結(jié)果是否正確，命令如下：x>>q1.0e-006可見得到了較高精度的結(jié)果（7）求下列非線（7）求下列非線性方程x0.7sinx0.2cosyy0.7cosx0.2siny解：1先作圖分析gridonholdezplot('y-初值（7）求下列非（7）求下列非線性方程x0.7sinx0.2cosyy0.7cosx0.2siny解：2.建立函數(shù)文functionF(2)=x(2)-（7）求下列非線性方（7）求下列非線性方程x0.7sinx0.2cosyy0.7cosx0.2siny輸入命令，調(diào)用函數(shù)子程序，設(shè)初值[0,0解>>x實(shí)驗(yàn)A56B實(shí)驗(yàn)A56BB424229（8）722934A24256（9），B設(shè)AB的乘積5求矩312A3（10）>>A^(-4ans2A3（10）>>A^(-4ans解：輸入命令>>A=[12;3求逆陣也A1324>>>>A2A1實(shí)驗(yàn)22n31實(shí)驗(yàn)22n31255n31 15lim53n3nsymsans=3x1例：求極限x3x1例：求極限xlimx4x xlimxsymsxlimit((x^2-3*x-4)/(x-2xxx22limx2xxx22limxx22 lim1xsymsx 例：求極限limxsinx 例：求極限limxsinxxxsymsxsymsxezplot('x*sin(1/x)+sin(x)/x',[-實(shí)驗(yàn)x實(shí)驗(yàn)xyx2cosxy10yx2cosxy10yx2cosxx2cosxx2cos2xcosxx2sinxy2xcosxx2sinx2cosx2xsinxx2xsin2cosx4xsinxx2cosxy2cosx4xsinxx2cosx2cosx4xsinxx26sinx6xcosxx2siny10y972sinx18xcosxx2sin90cosx20xsinxx2cos2cos20*x*sin(x)-symsxyxyx102x109的1階到11yx102(x10)910x918(xy10x918(x10)890x8144(xy90x8144(x10)7720x71008(xsymsxforn=1:11yyx102x的1階到11ans=10*x^9+18*(x-10)^8ans=90*x^8+144*(x-10)^7ans=720*x^7+1008*(x-10)^6ans=5040*x^6+6048*(x-10)^5ans=30240*x^5+30240*(x-10)^4ans=151200*x^4+120960*(x-10)^3ans=604800*x^3+362880*(x-10)^2ans=1814400*x^2+725760*x-7257600ans=3628800*x+725760ansans例：求函數(shù)f(x)sinaxcosbx一階導(dǎo)數(shù)，并求f 例：求函數(shù)f(x)sinaxcosbx一階導(dǎo)數(shù)，并求f 。abf(x)sinaxcossinaxcosbxsinaxcosacosaxcosbxsinaxbsinacosaxcosbxbsinaxsinf acosaxcosbxbsinaxsinbx a1acos cosb bsina sinb aaaa例：求函數(shù)f(x)例：求函數(shù)f(x)sinaxcosbx一階導(dǎo)數(shù)，并求f 。abdscos(a*x)*a*cos(b*x)-dszcos(a/(a+b))*a*cos(b/(a+b))-symsxab例：求由方程2x22xyy2x2y例：求由方程2x22xyy2x2y10令：zf(x,y2x22xyx2y1z4x2yz2x2y4x2yyx2x2y4x2y2x2yds(-4*x+2*y-1)/(-symsxz2z23y ds_x=4*x*y^3ds_y=6*x^2*y^2ds2_x=4*y^3ds2_y=12*x^2*yds2_xyds2_yxz 6x22z 42z12x2212xy22z12xy2symsxyMATLAB軟件求函數(shù)積分的命令是int，它既可以用于計(jì)intfMATLAB軟件求函數(shù)積分的命令是int，它既可以用于計(jì)intff關(guān)于symsf關(guān)于變量v1x C13x123xdxx2dxC12>>syms>>int(sqrt(x))ans=1x2111x2111(x21)(x21)1(x21 1 dxx2dxdx11x3xarctanx3g1/3*x^3-31/3x-x+symsxMATLAB軟件求函數(shù)積分的命令是int，它既可以用于計(jì)MATLAB軟件求函數(shù)積分的命令是int，它既可以用于計(jì)intf,a,f關(guān)于syms定義的符號(hào)變量從abintf,v,a,f關(guān)于變量v的從a42xdx241221>>symsans=xx12011131213161解dx 22xx12011131213161解dx 223xxx0>>sym>>int('(x-x^2)',0,1)ans=21x111x01xx1x21x111x01xx1x1x xxdxdx12222121x1x11122212>>syms>>int(abs(1-x),1,2)ans= dx1dx dx1dxarctan1解1arctanarctan22記號(hào)arctanarctanlim同arctanlim>>syms>