MATLAB程序設(shè)計與應(yīng)用(劉衛(wèi)國編)課后實驗答案

1、實驗一MATLAB運算基礎(chǔ)1.先求下列表達式的值，然后顯示MATLAB工作空間的使用情況并保存全部變量。(1)2sin850z=11+e2(2)z=ln(x+.；1+x2)，22其中x=_02451+2i5(3)e0.3ae-0.3a.0.3+asm(a+0.3)+In22a=3.0,2.9,L,2.9,3.0(4)12Z=1214122t+10t11t2，其中t=0:2t=0&t=1&t=2&t=0&t=1&t=2&t3).*(t.八2-2*t+1)z1=z2=2.已知：1234-4-13-A=34787,B=20336573-27求下列表達式的值：A+6*B和A-B+I（其中丨為單位矩陣）

2、(2)A*B和A.*B(3)A八3和A.八3(4)A/B及BA(5)A,B和A(1,3,：）；B八2解：M文件：A二1234-4;34787;3657;B=13-1；203;3-27;A+6.*B運算結(jié)果：ans=1852-10467105215349ans=1231-3328840671ans=684462309-72596154-5241ans=3.設(shè)有矩陣A和Bj2345-301667891017_69A=1112131415,B=023_41617181920970212223242541311(1)求它們的乘積c。將矩陣C的右下角3X2子矩陣賦給D。(3)查看MATLAB工作空間的使

3、用情況。解：.運算結(jié)果：E=(reshape(1:1:25,5,5);F二3016;17一69;023一4;970;4111；C=E*FH=C(3:5,2:3)C=9315077258335237423520397588705557753890717H=5203977055578907174.完成下列操作：求100,999之間能被21整除的數(shù)的個數(shù)。建立一個字符串向量，刪除其中的大寫字母解：(1)結(jié)果：m=100:999;n=find(mod(m,21)=0);length(n)ans=43(2).建立一個字符串向量例如：ch二ABC123d4e56Fg9;則要求結(jié)果是:ch二ABC123d4

4、e56Fg9;k=find(ch二A&ch二Z);ch(k)=ch=123d4e56g9實驗二MATLAB矩陣分析與處理E1.設(shè)有分塊矩陣A=o3x3L2x3R3x2，其中E、R、0、S分別為單位矩陣、隨機矩陣、零矩陣S2x2和對角陣，試通過數(shù)值計算驗證A2=ER+RSOS2解：M文件如下;輸出結(jié)果：2.產(chǎn)生5階希爾伯特矩陣H和5階帕斯卡矩陣P,且求其行列式的值Hh和Hp以及它們的條件數(shù)Th和Tp,判斷哪個矩陣性能更好。為什么解：M文件如下：輸出結(jié)果：3.建立一個5X5矩陣，求它的行列式值、跡、秩和范數(shù)。解：M文件如下:輸出結(jié)果為：4.已知-29618A=20512-885求A的特征值及特征向

5、量，并分析其數(shù)學意義。解：M文件如圖:scriptLa3Cal|Oj/E-；:5.下面是一個線性方程組：113411x1-0.95x20.67x30.52(1)求方程的解。將方程右邊向量元素b3改為再求解，并比較b3的變化和解的相對變化。(3)計算系數(shù)矩陣A的條件數(shù)并分析結(jié)論。解：M文件如下:輸出結(jié)果：6.建立A矩陣，試比較sqrtm(A)和sqrt(A),分析它們的區(qū)別。解：M文件如下:分析結(jié)果知：sqrtm(A)是類似A的數(shù)值平方根(這可由b1*b1二A的結(jié)果看出)，而sqrt(A)則是對A中的每個元素開根號，兩則區(qū)別就在于此。實驗三選擇結(jié)構(gòu)程序設(shè)計一、實驗?zāi)康恼莆战⒑蛨?zhí)行M文件的方法。

6、掌握利用if語句實現(xiàn)選擇結(jié)構(gòu)的方法。掌握利用switch語句實現(xiàn)多分支選擇結(jié)構(gòu)的方法。掌握try語句的使用。二、實驗內(nèi)容求分段函數(shù)的值。x2+x一6xdispToomanyinputarguments.實驗四循環(huán)結(jié)構(gòu)程序設(shè)計一、實驗?zāi)康恼莆绽胒or語句實現(xiàn)循環(huán)結(jié)構(gòu)的方法。掌握利用while語句實現(xiàn)循環(huán)結(jié)構(gòu)的方法。熟悉利用向量運算來代替循環(huán)操作的方法。二、實驗內(nèi)容兀211111.根據(jù)丁=廠+L+，求n的近似值。當n分別取100、1000、100006122232n2時，結(jié)果是多少要求：分別用循環(huán)結(jié)構(gòu)和向量運算（使用sum函數(shù)）來實現(xiàn)。解：M文件如下：e編輯器-Unlil1ed.9*回岡文件騙輯

7、文本（1）跳轉(zhuǎn)單兀心?X楷Ome勻I口V丨駕謂-1.0+|1.1X1國循環(huán)結(jié)構(gòu)計算“值2y=0;3n=inputm二);4piari=l;n5y=y+l/i/i;6endTpi=sqrt(6y)6995向里方法計算Pi值10n=input(n=、11i二1./(1:r0-2;12E=suni(i);13pi=sqrt町腳本13列IM改寫.：運行結(jié)果如下：1112.根據(jù)y=1+-+-+L+，求:352n-1(1)y3求ff100中：以下是運算結(jié)果：5.若兩個連續(xù)自然數(shù)的乘積減1是素數(shù)，則稱這兩個邊疆自然數(shù)是親密數(shù)對，該素數(shù)是親密素數(shù)。例如，2X3-1=5,由于5是素數(shù)，所以2和3是親密數(shù)，5是親

8、密素數(shù)。求2,50區(qū)間內(nèi)：2929(1)親密數(shù)對的對數(shù)。(2)與上述親密數(shù)對對應(yīng)的所有親密素數(shù)之和。解：M文件:運算結(jié)果為j=23615實驗五函數(shù)文件一、實驗?zāi)康睦斫夂瘮?shù)文件的概念。掌握定義和調(diào)用MATLAB函數(shù)的方法。二、實驗內(nèi)容定義一個函數(shù)文件，求給定復(fù)數(shù)的指數(shù)、對數(shù)、正弦和余弦，并在命令文件中調(diào)用該函數(shù)文件。解：M文件如下:函數(shù)文件：functione,l,s,c二fushu(z)%fushu復(fù)數(shù)的指數(shù)，對數(shù)，正弦，余弦的計算%e復(fù)數(shù)的指數(shù)函數(shù)值%l復(fù)數(shù)的對數(shù)函數(shù)值%s復(fù)數(shù)的正弦函數(shù)值%c復(fù)數(shù)的余弦函數(shù)值e=exp(z);l=log(z);s=sin(z)c=cos(z);命令文件M:z

9、二input(請輸入一個復(fù)數(shù)z二)a,b,c,d二fushu(z)一物理系統(tǒng)可用下列方程組來表示：mcos0-m-sin00_a-0_imsin010cos00iamgi0m-sin002N10_020-cos011N2mg2從鍵盤輸入mi.m2和6的值，求aa?、叫和的值。其中g(shù)取，輸入6時以角度為單位。要求：定義一個求解線性方程組AX二B的函數(shù)文件，然后在命令文件中調(diào)用該函數(shù)文件。解：M文件函數(shù)文件：functionX=fc(A,B)%fcfc是求解線性方程的函數(shù)%AA是未知矩陣的系數(shù)矩陣X=AB；命令M文件：clc;m1=input(輸入m1=);m2二input(輸入m2二)；thet

10、a二input(輸入theta二)；x二theta*pi/180;g=;A=m1*cos(x)-m1-sin(x)0m1*sin(x)0cos(x)00m2-sin(x)000-cos(x)1;B二0;m1*g;0;m2*g;X二fc(A,B)一個自然數(shù)是素數(shù)，且它的數(shù)字位置經(jīng)過任意對換后仍為素數(shù)。例如13是絕對素數(shù)。試求所有兩位絕對素數(shù)。要求：定義一個判斷素數(shù)的函數(shù)文件。解：M文件：函數(shù)文件functionp二prime(p)%輸入p的范圍，找出其中的素數(shù)m=p(length(p);fori=2:sqrt(m)n二find(rem(p,i)=0&p=i);p(n)=;%將卩中能被i整除，而卻

11、不等于i的元素，即下標為n的元素剔除，其余的即為素數(shù)endp;命令文件：clc;p=10:99;p二prime(p);%找出10到99內(nèi)的所有素數(shù)p=10*rem(p,10)+(p-rem(p,10)/10;%將p素數(shù)矩陣每個元素個位十位調(diào)換順序p二prime(p)%再對對換后的素數(shù)矩陣找出所有的素數(shù)運算結(jié)果：4.設(shè)f(x)=11+(x-2)2+0.1(x-3)4+0.01，編寫一個MATLAB函數(shù)文件，使得調(diào)用p二113171137317379779f(x)時，x可用矩陣代入，得出的f(x)為同階矩陣。解：函數(shù)文件:function仁fx(x)%fxfx求算x矩陣下的f(x)的函數(shù)值A(chǔ)=+(

12、x2)八2B=+(x3)八4仁1./A+1./B命令文件：x二input(輸入矩陣x=)仁fx(x)5.已知y=f(40)fGO)+f(20)當f(n)二n+10ln(nN5)時，求y的值。當f(n)=1X2+2X3+3X4+.+nX(n+1)時，求y的值。運算結(jié)果如下：(2).運算結(jié)果如下：實驗六高層繪圖操作一、實驗?zāi)康恼莆绽L制二維圖形的常用函數(shù)。掌握繪制三維圖形的常用函數(shù)。掌握繪制圖形的輔助操作。二、實驗內(nèi)容3sinx1設(shè)y二0.5+Gcosx，在x=02n區(qū)間取101點，繪制函數(shù)的曲線。解：M文件如下:x=linspace(0,2*pi,101);y=+3*sin(x)./(1+x.入2

13、);plot(x,y)運行結(jié)果有：已知y1=x&y2二cos(2x),y3=y1Xy2，完成下列操作:在同一坐標系下用不同的顏色和線型繪制三條曲線。以子圖形式繪制三條曲線。分別用條形圖、階梯圖、桿圖和填充圖繪制三條曲線解：(1)M文件:運行結(jié)果：-I口IxlFlI*dittffInserttoolsdesktoplinicvH-tlp口口出丄T：化-w豪斤甩口目口.運行結(jié)果：3)M文件：卜*1口砂m|4rbitix二-pi:pi/100:pi;y1=x.2；y2二cos(2*x);y3二y1.*y2;subplot(2,2,1);plot(x,y1,b-,x,y2,r:,x,y3,k一);su

14、bplot(2,2,2);bar(x,y1,b);title(y1=x2);subplot(2,2,3);bar(x,y2,r);title(y2=cos(2x);subplot(2,2,4);bar(x,y3,k);title(y3=y1*y2);由上面的M文件，只要依次將bar”改為“stairs”、“stem”、“fill”，再適當更改區(qū)間取的點數(shù)，運行程序即可，即有下面的結(jié)果：lieEJl-LQaikLipindoiHLpEilIaiInitrIIs-ililnhlcfilinJcTlid.】)xosi+nb口i.jddu4-2Q2105wy釦仙牽沖丫3=和冷2-口|*1-口1*1Ei

15、lsEAtIjntrLLwls“血乎LitpItiHelplieEjHTitilstrtJoiLs：Ifdrt卯liniwleljE*；I口QIQ10ID5fl=iC1010OOM4-2024-4-2024-4-20242=8妙）cdekiY：3=y12ME4-202J4202i已知e21ln(x+y/T+x2)在-5WxW5區(qū)間繪制函數(shù)曲線。y=(x+sqrt(pi)/(exp(2).*(x0);plot(x,y)Ini運行結(jié)果：x|0.00.4-0.4-41235e0.6-65aEOFileEditViewtToolsDeslztopWindowHelp由圖可看出，函數(shù)在零點不連續(xù)。繪制極

16、坐標曲線p=asin(b+n6),并分析參數(shù)a、b、n對曲線形狀的影響。解：M文件如下：theta=0:pi/100:2*pi;a=input(輸入a二);b=input(輸入b=);n=input(輸入n二);rho二a*sin(b+n*theta);polar(theta,rho,m)采用控制變量法的辦法，固定兩個參數(shù)，變動第三個參數(shù)觀察輸出圖象的變化。133TO.3JWZIOJ-XIEjIv皿Iui/LIwkjErijrJ*l|lllh：JdiL申HlowlIiilJalkLjoindi勺SEScQDU：；帆呃搖、焜刃口/ficwtI3=b=n=22/15VEFb-2jl=4W2心眼12

17、0301B：i2I020f、1-5嚴-in|x|Til*IdiI.圧i+rInvirLlull,laiklpinJnTl.Cd*.4忌cosptF2Bn=3r4f2jf5丸nFi+TriitVi+rIiiarl圖jlailLipIirr3Tilp赳19：i,3018002102TO-|n|x|?il+Idii.Vi+-Tnvirl:T1klaiklf-linJnTalpuJio-./sQaFal*IJitVi-wI-BEtThT*JjrlicitrpFinJTil-pH印1ZIaj1W30I印/210釣応詁T居|口目|回兇njxjl+EdiI:j.arIrvurtI曲生Qoblcf*inJn|

18、pflL-ljdiLlatliJiikL-?pjoinJ*lpEE日由I旳岷口i用i口kNF*r*3=4J3=20C=X-分析結(jié)果：由這8個圖知道，當a,n固定時，圖形的形狀也就固定了，b只影響圖形的旋轉(zhuǎn)的角度；當a,b固定時，n只影響圖形的扇形數(shù)，特別地，當n是奇數(shù)時，扇葉數(shù)就是n,當是偶數(shù)時，扇葉數(shù)則是2n個；當b,n固定時，a影響的是圖形大小，特別地，當a是整數(shù)時，圖形半徑大小就是a。繪制函數(shù)的曲線圖和等高線。x2豐y2xyz=cosxcosye4其中x的21個值均勻分布-5,5范圍，y的31個值均勻分布在0,10，要求使用subplot(2,1,1)和subplot(2,1,2)將產(chǎn)生

19、的曲面圖和等高線圖畫在同一個窗口上。解：M文件:x=linspace(-5,5,21);y=linspace(0,10,31);x,y二meshgrid(x,y);z二cos(x).*cos(y).*exp(-sqrt(x.入2+y.入2)/4);subplot(2,1,1);surf(x,y,z);title(曲面圖)；subplot(2,1,2);surfc(x,y,z);title(等高線圖)；Figore運行結(jié)果：555510Inixl曲面圖等高線圖:%頭卩動匯銘F00TileEditiewInsertT_oo1eDesktopIVindowHelp繪制曲面圖形，并進行插值著色處理。x

20、=cosscost.c兀c3兀y=cosssint0sJ一,0t);%找出大于數(shù)的個數(shù)p二y/30000%大于的所占百分比運行結(jié)果：將100個學生5門功課的成績存入矩陣P中，進行如下處理:分別求每門課的最高分、最低分及相應(yīng)學生序號。分別求每門課的平均分和標準方差。5門課總分的最高分、最低分及相應(yīng)學生序號。將5門課總分按從大到小順序存入zcj中，相應(yīng)學生序號存入xsxh。提示：上機調(diào)試時，為避免輸入學生成績的麻煩，可用取值范圍在45,95之間的隨機矩陣來表示學生成績。解：M文件:t=45+50*rand(100,5);P=fix(t);%生成100個學生5門功課成績x,l=max(P)%x為每門

21、課最高分行向量，丨為相應(yīng)學生序號y,k二min(P)%y為每門課最低分行向列,k為相應(yīng)學生序號mu二mean(P)%每門課的平均值行向量sig=std(P)%每門課的標準差行向量s=sum(P,2)%5門課總分的列向量X,m=max(s)%5門課總分的最高分X與相應(yīng)學生序號mY,n=min(s)%5門課總分的最低分Y與相應(yīng)學生序號nzcj,xsxh=sort(s)%zcj為5門課總分從大到小排序，相應(yīng)學生序號xsxh運行結(jié)果：某氣象觀測得某日6:0018:00之間每隔2h的室內(nèi)外溫度（0C）如實驗表1所示。實驗表1室內(nèi)外溫度觀測結(jié)果（0C）時間h681012141618室內(nèi)溫度t1室外溫度t2

22、試用三次樣條插值分別求出該日室內(nèi)外6:3018:30之間每隔2h各點的近似溫度（0C）解：T1=interp1(h,t1,spline)%室內(nèi)的3次樣條插值溫度T2二interp1(h,t2,spline)%室外的3次樣條插值溫度運行結(jié)果：1011014.已知Igx在1,101區(qū)間10個整數(shù)采樣點的函數(shù)值如實驗表2所示。實驗表2lgx在10個采樣點的函數(shù)值x1112131415161718191Igx0試求lgx的5次擬合多項式p(x)，并繪制出lgx和p(x)在1,101區(qū)間的函數(shù)曲線。解：M文件：x=1:10:101;y=lg10(x);P二polyfit(x,y,5)y1二polyval

23、(P,x);plot(x,y,:o,x,y1,-*)運行結(jié)果：Warning:Polynomialisbadlyconditioned.AddpointswithdistinctXvalues,reducethedegreeofthepolynomial,ortrycenteringandscalingasdescribedinHELPPOLYFIT.Inpolyfitat80P二這里出現(xiàn)警告是提示不必用5價函數(shù)就已經(jīng)可以完美擬合了，是可以降價擬合。)在1,101的區(qū)間函數(shù)圖像5有3個多項式P1(x)=x4+2x3+4x2+5,P2(x)=x+2,P3(x)=x2+2x+3,試進行下列操作:求

24、P(x)=P(x)+P(x)P(x)。TOC o 1-5 h z123求P(x)的根。當x取矩陣A的每一元素時，求P(x)的值。其中：-11.2-1.4_ HYPERLINK l bookmark4 o Current Document A=0.7523.5052.5當以矩陣A為自變量時，求P(x)的值。其中A的值與第題相同。解：M文件：clc;clear;p1二1,2,4,0,5;p2=1,2;p3=1,2,3;p2=0,0,0,p2;p3=0,0,p3;p4二conv(p2,p3);%p4是p2與p3的乘積后的多項式np4=length(p4);np仁length(p1);p二zeros(

25、1,np4-np1)p1+p4%求p(x)二p1(x)+p2(x)x二roots(p)%求p(x)的根A二-1;2;05;y二polyval(p,A)%x取矩陣A的每一元素時的p(x)值運行結(jié)果：實驗九數(shù)值微積分與方程數(shù)值求解一、實驗?zāi)康恼莆涨髷?shù)值導(dǎo)數(shù)和數(shù)值積分的方法掌握代數(shù)方程數(shù)值求解的方法。掌握常微分方程數(shù)值求解的方法。二、實驗內(nèi)容求函數(shù)在指定點的數(shù)值導(dǎo)數(shù)。xx2X3/(x)二12x3X2,x=1,2,3026x解：M文件：clc;clear;x=1;i=1；仁inline(det(xx入2x八3;12*x3*x2;026*x);whilex=g(i)=f(x);i二i+1;x=x+;%以

26、的步長增加，可再縮小步長提高精度endg;t=1:;dx二diff(g)/;%差分法近似求導(dǎo)f仁dx(1)%x=1的數(shù)值倒數(shù)f2=dx(101)%x=2的數(shù)值倒數(shù)f3=dx(length(g)T)%x=3的數(shù)值倒數(shù)運行結(jié)果：用數(shù)值方法求定積分。I=j2兀cos12+4sin(21)2+1dt的近似值。oI=s2nln(lx)dt01+x2解：M文件：clc;clear;仁inline(sqrt(cos(t.入2)+4*sin(2*t).入2+1);11二quad(f,0,2*pi)g二inline(log(1+x)./(1+x.八2);12二quad(g,0,2*pi)分別用3種不同的數(shù)值方法

27、解線性方程組。6x+5y-2z+5u=一49x-y+4z-u=133x+4y+2z-2u=13x-9y+2u=11解：M文件：clc;clear;A二65-25;9-14-1;3422;3-902;b二413111;x二Aby二inv(A)*bL,U=lu(A);z=U(Lb)求非齊次線性方程組的通解。2x+7x+3x+x=6TOC o 1-5 h z12340%非齊次方程組ifrank(A)二二rank(A,b)ifrank(A)二二ndisp(有唯一解x);x=Ab;elsedisp(有無窮個解，特解x，基礎(chǔ)解系y);x=Ab;y=null(A,r);clc;clear;formatrat

28、A二2731;3522;9417;b=642;x,y=line_solution(A,b)運行結(jié)果：有無窮個解，特解x，基礎(chǔ)解系yWarning:Rankdeficient,rank二2,tol二Inline_solutionat11x二所以原方程組的通解是：-1/11-9/1-2/1-5/111/1110/11X=k1+k0+0，其中k,k為任意常數(shù)。12120105.求代數(shù)方程的數(shù)值解。3x+sinx-ex=0在x二附近的根。0(2)在給定的初值xo=1,y0=1,zo=1下，求方程組的數(shù)值解。sinx+y2+Inz一7=03x+2y一z3+1=0 x+y+z一5=0ans=1289/68

29、2(2).M文件：functionF二fun(X)x=X(1);y=X(2)；z=X(3);F(1)=sin(x)+y入2+log(z)-7;F(2)=3*x+2-z八3+1;F(3)=x+y+z-5;X二fsolve(myfun,1,1,1,optimset(Display,off)運行結(jié)果：求函數(shù)在指定區(qū)間的極值。X3+COSX+XlogXf(x)=在(0,1)內(nèi)的最小值。eX/(X,X)=2X3+4XX3一10XX+X2在0,0附近的最小值點和最小值。12112122解：M文件:function仁g(u)x二u(1);y二u;仁2*x八3+4*x.*y八3-10*x.*y+y八2;clc

30、;clear;formatlong仁inline(x八3+cos(x)+x*log(x)/exp(x);x,fmin1=fminbnd(f,0,1)U,fmin2=fminsearch(g,0,0)運行結(jié)果求微分方程的數(shù)值解。xd2y_dyn-5丄+y=0dx2dxy(0)=0y(0)=0解：M文件：functionxdot二sys(x,y)xdot二y(2);(5*y(2)-y(1)/x;clc;clear;x0二;xf=20;x,y=ode45(sys,x0,xf,00);x,y運行結(jié)果：000000000000000000000000000000000000000000000000000

31、000求微分方程組的數(shù)值解，并繪制解的曲線。y=yyTOC o 1-5 h z123y=-yy4213y=0.51yy312y(0)=0,y(0)=1,y(0)=1123解：令y1二x,y2二y,y3二z;這樣方程變?yōu)椋簓=xzz=0.51xy,自變量是tx(0)=0,y(0)=1,z(0)=1M文件：functionxdot二sys(x,y)xdot二y(2)*y(3);-y(1)*y(3);*y(1)*y(2);clc;clear;t0=0;t仁8;x,y=ode23(sys,t0,tf,0,1,1)plot(x,y)運行結(jié)果：圖形：o.e5|=1IEfli.日4-0.2-.4-ID.6-

32、.BFilftS.-111yiftwlnsort22.0ol=XLndowHlpuZ5gJ實驗十符號計算基礎(chǔ)與符號微積分一、實驗?zāi)康恼莆斩x符號對象的方法。掌握符號表達式的運算法則以及符號矩陣運算。掌握求符號函數(shù)極限及導(dǎo)數(shù)的方法。掌握求符號函數(shù)定積分和不定積分的方法。二、實驗內(nèi)容已知x=6,y=5，利用符號表達式求3+x-y提示：定義符號常數(shù)X二sym(6),y二sym(5)。解：M文件：clearall;clc;x=sym(6);y=sym(5);z=(1+x)/(sqrt(3+x)-sqrt(y)運行結(jié)果：分解因式。(1)x4-y4(2)5135運行結(jié)果：5*13*79化簡表達式。(l)s

33、inPcosP-cosPsinP(2)4+8x+312122X+1解：M文件：clearall;clc;symsbetalbeta2x;fl二sin(beta1)*cos(beta2)-cos(beta1)*sin(beta2);simplify(f1)%(1)問f2=(4*x2+8*x+3)/(2*x+1);simplify(f2)%(2)問運行結(jié)果：已知-010-100-abcP=100,P=010,A=def12001101ghk完成下列運算：。B的逆矩陣并驗證結(jié)果。(3)包括B矩陣主對角線元素的下三角陣。B的行列式值。解：M文件:clearall;clc;symsabcdefghk;p

34、1=010;100;001;p2=100;010;101;A=abc;def;ghk;B二p1*p2*AB1二inv(B)%8的逆矩陣B1*B%驗證逆矩陣結(jié)果B2=tril(B)d二det(B)運行結(jié)果：B二TOC o 1-5 h zd,e,fa,b,ca+g,b+h,c+kB1二-(c*h-b*k)/(a*f*h-b*f*g-c*d*h+c*e*g-a*e*k+b*d*k),(b*f-c*e+f*h-e*k)/(a*f*h-b*f*g-c*d*h+c*e*g-a*e*k+b*d*k),-(b*f-c*e)/(a*f*hb*f*g-c*d*h+c*e*g-a*e*k+b*d*k)(c*g-a*

35、k)/(a*f*h-b*f*g-c*d*h+c*e*g-a*e*k+b*d*k),-(a*f-c*d+f*g-d*k)/(a*f*h-b*f*g-c*d*h+c*e*g-a*e*k+b*d*k),(a*f-c*d)/(a*f*hb*f*g-c*d*h+c*e*g-a*e*k+b*d*k)(a*h-b*g)/(a*f*h-b*f*g-c*d*h+c*e*g-a*e*k+b*d*k),(a*e-b*d-5.用符號方法求下列極限或?qū)?shù)。(1)limxtOx(esinx+1)2(etanx1)sin3x(2)limXT1+VX+1丁=匕凹,求y,y已知A=axtcosxt3lnxdAd2Ad2A刀別求

36、dX喬麗已矢f(x,y)=(x22x)e-x2-y2-xy,求y,fx=0,y=1exexey解：M文件:clearall;clcsymsxtayz;f1=(x*(exp(sin(x)+1)-2*(exp(tan(x)T)/sin(x)八3;%(1)limit(f1)f2=(sqrt(pi)-sqrt(acos(x)/sqrt(x+1);limit(f2,x,T,right)y=(1-cos(2*x)/x;y1=diff(y)y2=diff(y,2)A=axt八3;t*cos(x)log(x);Ax仁diff(A,x,1)At2=diff(A,t,2)Axt二diff(Ax1,t)仁(x入2-

37、2*x)*exp(-x入2-z入2-x*z);Zx二-diff(f,x)/diff(f,z)dfxz=diff(diff(f,x),z);x=sym(0);z=sym(1);%符號運算返回數(shù)值eval(dfxz)運行結(jié)果：-(exp(x入2+x*z+z入2)*(2*x-2)/exp(x入2+x*z+z入2)+(2*x-x入2)*(2*x+z)/exp(x入2+x*z+z入2)/(2*x-x入2)*(x+2*z)ans二4/exp(1)6.用符號方法求下列積分。J匚1+x4+x8(arcsinx)2J1-x2(3)J+s三1dx(4)Jln2ex(1+ex)2dx0 x4+10解：M文件：clear;clc;x二sym(x);f1=1/(1+xP+x八8);%(1)f2=1/(asin(x)廠2/sqrt(1-x入2);%(2)f3=(x八2+1)/(x八4+1);%(3)f4=exp(x)*(1+exp(x)2;%(4)F1二int(f1)F2=int(f2)F3=int(f3,0,inf)F4=int(f4,0,log(2)運行結(jié)果：F1二(3八(1/2)*log(x八2+3八(1/2)*x+1)/12-(3八(1/2