【正文】 0 0 2A = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 a = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0由ans,所以2. 產(chǎn)生5階希爾伯特矩陣H和5階帕斯卡矩陣P，且求其行列式的值Hh和Hp以及它們的條件數(shù)Th和Tp，判斷哪個(gè)矩陣性能更好。4 13 11]。 A+6.*B AB+eye(3) A*B A.*B A^3 A.^3 A/B B\A [A,B] [A([1,3],:)。t1).*(t.^2)+(t=1amp。 5]。3 65 7]。 解：(1) 結(jié)果：m=100:999。解： M文件如下：輸出結(jié)果為：A = 17 24 1 8 15 23 5 7 14 16 4 6 13 20 22 10 12 19 21 3 11 18 25 2 9d = 5070000t = 65c1 = c2 = cinf = 4. 已知求A的特征值及特征向量，并分析其數(shù)學(xué)意義。要求：(1) 分別用if語(yǔ)句和switch語(yǔ)句實(shí)現(xiàn)。11 2 3 9 7 3。pi=sqrt(6*s)n=100pi = n=1000pi = n=10000pi =2. 根據(jù)，求：(1) y3時(shí)的最大n值。解：M文件：運(yùn)算結(jié)果為：j = 29s = 23615實(shí)驗(yàn)五 函數(shù)文件一、實(shí)驗(yàn)?zāi)康?. 理解函數(shù)文件的概念。解： M文件:function X= fc(A,B)%fc fc是求解線性方程的函數(shù)%A A是未知矩陣的系數(shù)矩陣X=A\B； 命令M文件：clc。m2*g]。B=+(x3).^4。)。n2=input(39。plot(x,y)運(yùn)行結(jié)果有：2. 已知y1=x2，y2=cos(2x)，y3=y1y2，完成下列操作：(1) 在同一坐標(biāo)系下用不同的顏色和線型繪制三條曲線。y1=x.^2。)。,x,y2,39。 title(39。解：M文件如下：clc。當(dāng)a,b固定時(shí)，n只影響圖形的扇形數(shù)，特別地，當(dāng)n是奇數(shù)時(shí)，扇葉數(shù)就是n,當(dāng)是偶數(shù)時(shí)，扇葉數(shù)則是2n個(gè)。title(39。subplot(2,2,2)。)。x=2::2。,9)。grid on。flat39。Yaxis39。)。(2) 分別求每門(mén)課的平均分和標(biāo)準(zhǔn)方差。spline39。 解：M文件：clc。二、實(shí)驗(yàn)內(nèi)容1. 求函數(shù)在指定點(diǎn)的數(shù)值導(dǎo)數(shù)。 %差分法近似求導(dǎo)f1=dx(1) %x=1的數(shù)值倒數(shù)f2=dx(101) %x=2的數(shù)值倒數(shù)f3=dx(length(g)1) %x=3的數(shù)值倒數(shù) 運(yùn)行結(jié)果：f1 = f2 = f3 = 2. 用數(shù)值方法求定積分。3 9 0 2]。)。clear。z=X(3)。 f=2*x.^3+4*x.*y^310*x.*y+y.^2。sys39。4. 掌握求符號(hào)函數(shù)定積分和不定積分的方法。)。syms a b c d e f g h k。 %(2)limit(f2,x,1,39。x=sym(39。n=sym(39。x1=solve(39。y(0)=039。)運(yùn)行結(jié)果：x = C1/exp(t) + C2*exp(2*t) y = C1/exp(t) + C2*exp(2*t) + C3/exp(2*t) z = C2*exp(2*t) + C3/exp(2*t)寧可累死在路上，也不能閑死在家里！寧可去碰壁，也不能面壁。,39。clc。f=log(x)。(1) 計(jì)算。)。clc。(3) 包括B矩陣主對(duì)角線元素的下三角陣。clc。sys39。clear。(1) 在(0,1)內(nèi)的最小值。fzero(39。 y=null(A,39。)。clear。 x=x+。0 5 ]。(2) 求P(x)的根。t2=[ ]。解：M文件:clc。,[4,0,0],39。Xlabel(39。)。edgecolor39。text(,^2*exp(2*),39。)39。surf(x,y,z)。mesh(x,y,z)。曲面圖39。polar(theta,rho,39。x=5::5。subplot(2,2,3)。y3=y1.*y2。r:39。,x,y3,39。二、實(shí)驗(yàn)內(nèi)容1. 設(shè)，在x=0~2π區(qū)間取101點(diǎn)，繪制函數(shù)的曲線。命令文件：clc。解：(1):function f=f(x)f=x+10*log(x^2+5)。p=prime(p)。A=[m1*cos(x) m1 sin(x) 0 m1*sin(x) 0 cos(x) 0 0 m2 sin(x) 0 0 0 cos(x) 1]。)。解：M—文件以下是運(yùn)算結(jié)果：max(f)=437763282635min(f)=899412113528sum(f)=742745601951c1=49c2=2c3=495. 若兩個(gè)連續(xù)自然數(shù)的乘積減1是素?cái)?shù)，則稱(chēng)這兩個(gè)邊疆自然數(shù)是親密數(shù)對(duì)，該素?cái)?shù)是親密素?cái)?shù)。n=39。11 2 3 9 7 3。二、實(shí)驗(yàn)內(nèi)容1. 求分段函數(shù)的值。ch(k)=[]ch =123d4e56g9實(shí)驗(yàn)二 MATLAB矩陣分析與處理1. 設(shè)有分塊矩陣，其中E、R、O、S分別為單位矩陣、隨機(jī)矩陣、零矩陣和對(duì)角陣，試通過(guò)數(shù)值計(jì)算驗(yàn)證。9 7 0。3 2 7]。z4=(t=0amp。z2=1/2*log(x+sqrt(1+x^2)) a=::。B=[1 3 1。n=find(mod(m,21)==0)。解：M文件如圖：輸出結(jié)果為：V = D = 0 0 0 0 0 0 數(shù)學(xué)意義：V的3個(gè)列向量是A的特征向量，D的主對(duì)角線上3個(gè)是A的特征值，特別的，V的3個(gè)列向量分別是D的3個(gè)特征值的特征向量。(2) 輸入百分制成績(jī)后要判斷該成績(jī)的合理性，對(duì)不合理的成績(jī)應(yīng)輸出出錯(cuò)信息。2 3 4 5 6 7]輸入一正整數(shù)n=6 2 3 4 5 6 7ans =Error using == dispToo many input arguments.實(shí)驗(yàn)四 循環(huán)結(jié)構(gòu)程序設(shè)計(jì)一、實(shí)驗(yàn)?zāi)康?. 掌握利用for語(yǔ)句實(shí)現(xiàn)循環(huán)結(jié)構(gòu)的方法。(2) 與(1)的n值對(duì)應(yīng)的y值。2. 掌握定義和調(diào)用MATLAB函數(shù)的方法。m1=input(39。X=fc(A,B) 運(yùn)算結(jié)果：輸入m1=1輸入m2=1輸入theta=30X = 3. 一個(gè)自然數(shù)是素?cái)?shù)，且它的數(shù)字位置經(jīng)過(guò)任意對(duì)換后仍為素?cái)?shù)。f=1./A+1./B。n2=input(39。n2=39。(2) 以子圖形式繪制三條曲線。y2=cos(2*x)。subplot(1,3,3)。r:39。y2=cos(2x)39。theta=0:pi/100:2*pi。當(dāng)b,n固定時(shí)，a影響的是圖形大小，特別地，當(dāng)a是整數(shù)時(shí)，圖形半徑大小就是a。等高線圖39。surf(x,y,z)。 運(yùn)行結(jié)果有：實(shí)驗(yàn)七 低層繪圖操作二、實(shí)驗(yàn)內(nèi)容1. 建立一個(gè)圖形窗口，使之背景顏色為紅色，并在窗口上保留原有的菜單項(xiàng)，而且在按下鼠標(biāo)器的左鍵之后顯示出Left Button Pressed字樣。y=x.^2.*exp(2*x)。運(yùn)行結(jié)果:3. 利用曲面對(duì)象繪制曲面v(x,t)=(+π)。xlabel(39。)。)。shading interp。(3) 5門(mén)課總分的最高分、最低分及相應(yīng)學(xué)生序號(hào)。)%室外的3次樣條插值溫度 運(yùn)行結(jié)果：T1 = Columns 1 through 3 Columns 4 through 6 Column 7 T2 = Columns 1 through 3 Columns 4 through 6 Column 7 4. 已知lgx在[1,101]區(qū)間10個(gè)整數(shù)采樣點(diǎn)的函數(shù)值如實(shí)驗(yàn)表2所示。clear。 解：M文件：clc。(1) 的近似值。b=[4 13 1 11]39。 end else disp(39。format ratA=[2 7 3 1。F(1)=sin(x)+y^2+log(z)7。clc。,[x0,xf],[0 0])。二、實(shí)驗(yàn)內(nèi)容1. 已知x=6,y=5，利用符號(hào)表達(dá)式求提示：定義符號(hào)常數(shù)x=sym(‘6’)，y=sym(‘5’)。a=x^4y^4。p1=[0 1 0。right39。x39。n39。log(x+1)5/(1+sin(x))=239。,39。是狼就要練好牙，是羊就要練好腿。Dy=4*x5*y+3*z39。 解：M文件：clear all。)。二、實(shí)驗(yàn)內(nèi)容1. 級(jí)數(shù)符號(hào)求和。139。解：M文件：clear all。 (2) B的逆矩陣并驗(yàn)證結(jié)果。(1) x4y4 (2) 5135解：M文件：clear all。[x,y]=ode23(39。 clc。))運(yùn)行結(jié)果：X = 909/1073 1735/728 1106/625 6. 求函數(shù)在指定區(qū)間的極值。clear。)。有無(wú)窮個(gè)解，特解x，基礎(chǔ)解系y39。解：M文件：clc。 i=i+1。 2 。）在[1,101]的區(qū)間函數(shù)圖像5. 有3個(gè)多項(xiàng)式P1(x)=x4+2x3+4x2+5，P2(x)=x+2，P3(x)=x2+2x+3，試進(jìn)行下列操作：(1) 求P(x)=P1(x)+P2(x)P3(x)。t1=[ ]。(3) 。Position39。)。meshsurf39。,... [,],39。,2)。39。subplot(2,2,4)。subplot(2,2,1)。title(39。rho=a*sin(b+n*theta)。解：M文件：clc。)。y2=cos(2*x)。plot(x,y2,39。r:39。3. 掌握繪制圖形的輔助操作。ends=sum(g)。(2) 當(dāng)f(n)=12+23+34+...+n(n+1)時(shí)，求y的值。p=10:99。g=。請(qǐng)輸入一個(gè)復(fù)數(shù)z=39。(2) 正數(shù)、零、負(fù)數(shù)的