【正文】 數(shù)合并collect(S)對(duì)符號(hào)表達(dá)式S中默認(rèn)變量（findsym(S,1)）同次冪項(xiàng)的系數(shù)合并例10 將按自變量和展開。解：syms x t f fx ft。f=x*(x*(xt)^2+12)。ft=collect(f,t)。fx=collect(f)。輸出結(jié)果ft=x^2*t^22*x^3*t+x*(x^3+12)fx=x^42*x^3*t+x^2*t^2+12*x符號(hào)表達(dá)式的化簡(jiǎn)符號(hào)表達(dá)式的化簡(jiǎn)可由函數(shù)simple 和simplify實(shí)現(xiàn)，調(diào)用格式為simple(S) simplify(S)例11化簡(jiǎn) 解：syms x f f1 f2。f=cos(x)^2+sin(x)^2。f1=simplify(f)。 f2=simple(f)。輸出結(jié)果f1 =1f2 =1符號(hào)表達(dá)式的通分符號(hào)表達(dá)式的通分函數(shù)為numden，其調(diào)用格式為[n d]=numden(S) 將符號(hào)表達(dá)式S轉(zhuǎn)換為分子和分母都是整系數(shù)的最佳多項(xiàng)式。例12 對(duì)表達(dá)式進(jìn)行通分。解：syms x y f。 f=x/y+y/x。 [n d]=numden(f)輸出結(jié)果n =x^2+y^2d =x*y符號(hào)表達(dá)式的替換函數(shù)subs用特定符號(hào)替換表達(dá)式中的某一特定符號(hào)，其調(diào)用格式為subs(S,old,new) 用新的符號(hào)變量new代替符號(hào)表達(dá)式S中的變量old。例13 syms x y z。y=sin(x)。z=subs(y,x,pi)輸出結(jié)果z = 符號(hào)表達(dá)式的轉(zhuǎn)換符號(hào)函數(shù)sym()可將一個(gè)數(shù)字轉(zhuǎn)換為符號(hào)表達(dá)式，而函數(shù)single()恰好相反，它能把一個(gè)常數(shù)轉(zhuǎn)換為一個(gè)數(shù)值，格式如為single(f)。例14 f=sym(39。1sqrt()+39。) f = 1sqrt()+f = single(f)f = 3 符號(hào)函數(shù)的運(yùn)算這里只簡(jiǎn)單介紹符號(hào)函數(shù)的復(fù)合運(yùn)算和反函數(shù)運(yùn)算，MATLAB函數(shù)分別為pose 和finverse。pose函數(shù)的調(diào)用格式為pose(f,g) 返回當(dāng)f=f(x)和g=g(y)時(shí)的復(fù)合函數(shù)f(g(y))。pose(f,g,z) 返回的復(fù)合函數(shù)以z為自變量。finverse的調(diào)用格式為g=finverse(f) 返回符號(hào)函數(shù)f的反函數(shù)。g=finverse(f,v) 返回的符號(hào)函數(shù)表達(dá)式的自變量為v。例15 syms x y z t u。 f=1/(1+x^2)。 g=sin(y)。 則 fg=pose(f,g) fg = 1/(1+sin(y)^2) gf=pose(g,f) gf = sin(1/(1+x^2)) fgz=pose(f,g,t) fgz = 1/(1+sin(t)^2)例16 syms x y。 f=2*x1。則g=finverse(f) g = 1/2+1/2*xh=x^2+y。g=finverse(h,y)g =x^2+y多項(xiàng)式在MATLAB中，使用行向量來表示多項(xiàng)式的系數(shù)，行向量中各元素按多項(xiàng)式系數(shù)由高到低排列。即的系數(shù)向量為。因此，在MATLAB中，將多項(xiàng)式問題轉(zhuǎn)化為矢量問題。生成多項(xiàng)式的基本函數(shù)為poly，poly2sym將多項(xiàng)式轉(zhuǎn)換為矢量形式。多項(xiàng)式的創(chuàng)建主要有以下三種方法。1 系數(shù)矢量的直接輸入法例1 創(chuàng)建多項(xiàng)式。解：p=[1 4 0 2]。p=poly2sym(p)p =x^34*x^2+22 特征多項(xiàng)式輸入法多項(xiàng)式也可由通過矩陣的特征多項(xiàng)式取得，為階方陣，得特征多項(xiàng)式為。的特征多項(xiàng)式創(chuàng)建的命令為：p=poly(A)，p=poly2sym(poly(A))。例2 求的特征多項(xiàng)式。解：A=[1 2 3。4 4 6。7 8 9]。p=poly(A)p = %特征多項(xiàng)式的系數(shù)p=poly2sym(poly(A))p =x^314*x^228*x12 %特征多項(xiàng)式。3 由根創(chuàng)建多項(xiàng)式將多項(xiàng)式的根用向量表示，即，則p=poly(r)表示以、….、為根的多項(xiàng)式。例3 特征多項(xiàng)式的根為1，2，3，求該多項(xiàng)式。解： r=[1 2 3]。 p=poly(r)。 p=poly2sym(p)p =x^36*x^2+11*x6 多項(xiàng)式的運(yùn)算這里只介紹多項(xiàng)式的運(yùn)算包含求值、求根、多項(xiàng)式的四則運(yùn)算，至于多項(xiàng)式的其它形式運(yùn)算不作介紹。多項(xiàng)式求值 求多項(xiàng)式值有兩種方算法，一種按數(shù)值規(guī)則運(yùn)算，對(duì)應(yīng)的函數(shù)為 polyval；另一種按矩陣的運(yùn)算規(guī)則運(yùn)算，對(duì)應(yīng)的函數(shù)為polyvalm。函數(shù)polyval的調(diào)用格式為polyval(p,x) 求多項(xiàng)式p在x點(diǎn)的值，這里x可以是一個(gè)向量。polyvalm(p,x) 求多項(xiàng)式p在矩陣x的值。例1 求多項(xiàng)式在1，2，3點(diǎn)的值。解：p=[1 4 0 2]。x=[1 2 3]。 f=polyval(p,x)f = 1 6 7例2 求多項(xiàng)式對(duì)于矩陣的值。解：p=[1 4 0 2]。A=[1 2。7 8]。f=polyvalm(p,A)f = 83 102 357 440多項(xiàng)式求根MATLAB有兩種進(jìn)行多項(xiàng)式求根的方法，一種是直接調(diào)用求根函數(shù)roots，另一種是先把矩陣轉(zhuǎn)化為伴隨矩陣，然后求其特征值。在MATLAB中，多項(xiàng)式的根用列向量表示。例3 求多項(xiàng)式的根。解：p=[1 1 1 1]。r=roots(p)r = + 本題也可以這樣求解：p=[1 1 1 1]。A=pan(p)。r=eig(A)r = + 多項(xiàng)式的四則運(yùn)算多項(xiàng)式的加減運(yùn)算，同向量的加減運(yùn)算。值得注意的地方是，應(yīng)以高次看齊，低次要看成高次的。例4 求多項(xiàng)式和多項(xiàng)式的和與差。解：p=[1 1 1 1]。q=[0 1 1 0]。s=p+q。 s=poly2sym(s)s =x^3+2*x^2+1d=pq。d=poly2sym(d)d =x^3+2*x+1多項(xiàng)式的乘法由conv函數(shù)實(shí)現(xiàn)，多項(xiàng)式的除法由deconv實(shí)現(xiàn)。conv函數(shù)的調(diào)用格式為：c=conv(p,q) 求多項(xiàng)式p和q的乘積c。deconv函數(shù)的調(diào)用格式為：d=deconv(p,q)求多項(xiàng)式p和q的商c。例5 求多項(xiàng)式和多項(xiàng)式的積與商。解：p=[1 1 1 1]。q=[1 1 0]。c=conv(p,q)。 c=poly2sym(c)c =x^5xd=deconv(p,q)。 d=poly2sym(d)d =x+2參考文獻(xiàn)王家文 王皓 編程基礎(chǔ) 北京，機(jī)械工業(yè)出版社，2006。李海濤 鄧櫻 MATLAB程序設(shè)計(jì)教程 北京，高等教育出版社，2002。聞新 周露 張鴻 MATLAB科學(xué)圖形構(gòu)建基礎(chǔ)與應(yīng)用（）, 北京，科學(xué)出版社，2002。Shoichiro Nakamura (梁恒 劉曉艷 等譯) 科學(xué)計(jì)算引論（第二版）, 北京，電子工業(yè)學(xué)出版社，2002。焦光虹 主編，數(shù)學(xué)實(shí)驗(yàn)，哈爾濱，哈爾濱工業(yè)大學(xué)出版社，2000。肖海軍，數(shù)學(xué)實(shí)驗(yàn)初步，北京，科學(xué)出版社，2007。王沫然 ，北京，電子工業(yè)學(xué)出版社，2001。王向東 戎海武 文翰 數(shù)學(xué)實(shí)驗(yàn)，北京，高等教育出版社，2004。李繼成 主編 數(shù)學(xué)實(shí)驗(yàn) ，西安，西安交通大學(xué)出版社，2003。宋世德 郭滿才 數(shù)學(xué)實(shí)驗(yàn)，北京，高等教育出版社，2002。任玉杰 孫文惠 高等數(shù)學(xué)（一元微積分及其數(shù)學(xué)實(shí)驗(yàn)），北京，機(jī)械工業(yè)出版社，2004。任玉杰 高等數(shù)學(xué)（多元微積分及其數(shù)學(xué)實(shí)驗(yàn)），北京，機(jī)械工業(yè)出版社，2005。薛定宇 陳陽泉 高等應(yīng)用數(shù)學(xué)問題的MATLAB求解，第二版，北京，清華大學(xué)出版社，2008。薛毅 數(shù)值分析與實(shí)驗(yàn)，北京，北京工業(yè)大學(xué)出版社，2004。飛思科技產(chǎn)品研發(fā)中心，北京，電子工業(yè)學(xué)出版社，2003。吳贛昌 高等數(shù)學(xué)（理工類，第三版），北京，中國人民大學(xué)出版社，2009。習(xí)題11 利用help命名查看syms、diff、int、taylor等函數(shù)的用法。2 計(jì)算。3 利用符號(hào)運(yùn)算表達(dá)多項(xiàng)式，并將用替代。4 編寫函數(shù)文件，并畫出其圖形。5 繪制基本初等函數(shù)及二次曲面圖形。解1 略。解2 a=sin(5/12*pi)/(sqrt(2)+exp(sqrt(3)))a =解3 syms x p a。p=[1 3 2 1 1]。p=poly2sym(p)p=x^43*x^32*x^2+x1p1=subs(p,x,a+1)p1=(a+1)^43*(a+1)^32*(a+1)^2+a解4 function y=f(x)if x=0 y=0。elseif x=1 y=x。else y=1。endfplot(39。f(x)39。,[10 10])5 略。47