南郵MATLAB上機(jī)實(shí)驗(yàn)二

1、實(shí)驗(yàn)二MATLAB的符號(hào)運(yùn)算符號(hào)計(jì)算的特點(diǎn)：一，運(yùn)算以推理解析的方式進(jìn)行，因此不受計(jì) 算誤差積累問(wèn)題困擾；二，符號(hào)計(jì)算，或給出完全正確的封閉解，或 給出任意精度的數(shù)值解當(dāng)封閉解不存在時(shí)；三，符號(hào)計(jì)算指令的 調(diào)用比擬簡(jiǎn)單；四，計(jì)算所需時(shí)間較長(zhǎng)。在 MATLAB，符號(hào)計(jì)算雖以 數(shù)值計(jì)算的補(bǔ)充身份出現(xiàn)，但涉及符號(hào)計(jì)算的指令使用、運(yùn)算符操 作、計(jì)算結(jié)果可視化、程序編制以及在線幫助系統(tǒng)都是十分完整、便 捷的。一、實(shí)驗(yàn)?zāi)康?、熟悉符號(hào)對(duì)象和符號(hào)表達(dá)式；2、掌握符號(hào)表達(dá)式和符號(hào)函數(shù)的操作;3、學(xué)習(xí)利用符號(hào)代數(shù)方程的求解 。二、實(shí)驗(yàn)根本知識(shí)1、根本運(yùn)算指令；2. 掌握符號(hào)運(yùn)算的應(yīng)用 三、上機(jī)練習(xí)1. 1通過(guò)符

2、號(hào)計(jì)算求y(t) sint的導(dǎo)數(shù)dy。 2然后根據(jù)此結(jié)dt果，求獸和等。dt 11 0dt t 2>> clear allsyms ty=abs(si n( t)df=diff(y)df1=limit(df,t,0,'left')df2=subs(df,'t',sym(pi/2) y =abs(si n( t) df =sig n( si n( t)*cos(t)df1 =-1 df2 =02 x2222.計(jì)算二重積分1 1 (x y )dydx>>clear allsyms x yf=xA2+yA2Fin t=(i nt(i nt(f,

3、y,1,xA2),x,1,2) f = xA2 + yA2Fint =1006/1053. f x2y2(1)求今;求11 fdx. dy1>> clear all syms x y f=xA2+yA2; dfdy2=diff(f,y,2) dfdy2 =>> in tf=i nt(f,x,-1,1) intf =2*yA2 + 2/34. 解方程：3x24x 60>> clear allsyms xeq=solve(3*xA2+4*x-6)eq =-2/3+1/3*22八(1/2)-2/3-1/3*22八(1/2)5. 方程：d3y o dyr 3 2y

4、6dx dx(1)求通解;(2)求滿足y (0) 1,y (0) 2,y(0) 3的特解.>>clear all>> y=dsolve('D3y+3*Dy+2*y-6=0','x')y =3+C1*exp(1+2A(1/2)F(1/3)*(-1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)+C2*exp(- 1/2*(1+2A(1/2)A(1/3)*(-1-(1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)*si 門(1/2*3八(1/2)*(1+2A(1

5、/2)A(1/3)*(1-(1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)+C3*exp(- 1/2*(1+2A(1/2)A(1/3)*(-1-(1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)*cos(1/2*3A(1/2)*(1+2A(1/2)A(1/3)*(1-(1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)>> y二dsolve('D3y+3*Dy+2*y-6=0','y(0)=3,D2y(0)=1,Dy(0)=2','

6、;x')3+(-1/12*(1+2八(1/2)八(1/3)*2八(1/2)+1/6*(1+2八(1/2)八(2/3)+1/3*(1+2 八(1/2)八(1/3)-1/4*(1+2A(1/2)A(2/3)*2A(1/2)*exp(1+2A(1/2)A(1/3)*(-1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)+(1/3*3A(1/2)* (1+2A(1/2)A(1/3)-1/12*3A(1/2)*(1+2A(1/2)A(1/3)*2A(1/2)- 1/6*3A(1/2)*(1+2A(1/2)A(2/3)+1/4*3A(1/2)*(1+2A

7、(1/2)A(2/3)*2 A(1/2)*exp(-1/2*(1+2A(1/2)F(1/3)*(-1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)*sin(1/2*3A(1/ 2)*(1+2A(1/2)A(1/3)*(1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)+(1/12*(1+2A(1 /2)A(1/3)*2A(1/2)-1/6*(1+2A(1/2)A(2/3)- 1/3*(1+2A(1/2)A(1/3)+1/4*(1+2A(1/2)A(2/3)*2A(1/2)*exp(- 1/2*(1+2A(

8、1/2)A(1/3)*(-1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)*cos(1/2*3A(1/ 2)*(1+2A(1/2)A(1/3)*(1- (1+2A(1/2)A(1/3)+(1+2A(1/2)A(1/3)*2A(1/2)*x)226. 求方程x2 y2 1, xy 2 的解。>>clear all>> syms x y>> S=solve('xA2+yA2-1=0','x*y-2=0','x','y')S =x: 4x1 symy: 4

9、x1 sym>> disp('S.x'),disp(S.x),disp('S.y'),disp(S.y)S.x-1/2*(1/2*5A(1/2)+1/2*i*3A(1/2)A3+1/4*5A(1/2)+1/4*i*3A(1/2) -1/2*(1/2*5A(1/2)-1/2*i*3A(1/2)A3+1/4*5A(1/2)-1/4*i*3A(1/2) -1/2*(-1/2*5A(1/2)+1/2*i*3A(1/2)A3-1/4*5A(1/2)+1/4*i*3A(1/2) -1/2*(-1/2*5A(1/2)-1/2*i*3A(1/2)A3-1/4*5A(1/2)-1/4*i*3A(1/2)S.y1/2*5八(1/2)+1/2*i*3八(1/2)1/2*5八(1/2)-1/2*i*3八(1/2)-1/2*5A(1/2)+1/2*i*3A(1/2)-1/2*5A(1/2)-1/2*i*3A(1/2)實(shí)驗(yàn)小結(jié)：本次試驗(yàn)中，是為了讓我們更好的學(xué)習(xí)MATLAB符號(hào)運(yùn)算，讓我們熟練掌握符號(hào)計(jì)算的解題理念和模式，經(jīng)過(guò)本次實(shí)驗(yàn)，能夠使學(xué)生 感受到MATLA優(yōu)缺點(diǎn)，對(duì)以后的解題很有幫助。第一題說(shuō)明了極限和導(dǎo)數(shù)的符號(hào)運(yùn)算，學(xué)會(huì)應(yīng)用 diff 的不同情況 和 limit 的求取極限作用。第二題是對(duì)積分