數學軟件

江蘇科技大學數理學制 作業:課件內的練習,共26題,做好后 同一個word文件內.文件名:信息專業:X_學號后三位例如

在規定時間之前發送到作業要求:●標注練習題目抄寫,用Mathtype輸入公用WORD2003做練習,WORD文件要求：無網1江蘇科技大函數運本章介 的函數運算M文件inline講eval講函2江蘇科技大首先介 符號運符號計算基符號方程求微分方程3江蘇科技大比數值運符號運運算對象數運算結果數例結果為結果為a=是符號常4江蘇科技大第1號運符號對的一種數據類型 非數值字符常量,變量,函數或表達常量,變量,函數或表達

建立單個符號量,調用格式符號量名=sym('符號字符串1次可定義多個符號變量建立符號表達5江蘇科技大格式方法如格式定義變量指令sym('變量名空例如 >>x=sym('x')%定義x為符號空定義多個變量指令

syms

x2功能 定義x1,x2,…,xn為符號變量例如:>>symsabxt定義a,b,x,t為符號變6江蘇科技大格式格式

定義方式先定義符號變量AA注:這里b,x并未定義為符號變注意不

>>>>symsx>>f='2*x'

這是字符串不是符號表達式可用eval計算不可進行表達式計>>clear,>>clear,>>eval('x+4')ans=8江蘇科技大符號表達式的運>>g=sym('3*x^2+5*x-ansansans=ans=ans=(3*x^2+5*x-8江蘇科技大例計算A= 例計算A= 1 需要符號計算symsksymsk %求行列y=factor(D);%因式分D=3*k^2+k^3y=k^2*(3+k)9江蘇科技大多項式插

問題問題:已知(xi,yi),求n-1次多項式pn-1(x)滿足pn-1(xi)=基本方法

n

nl

x- x-這里pn-1(x)的計算很復雜

下面用符號計算出pn-1(x)的表達式江蘇科技大(程序

x=[12345問題的問題的syms

0.14-0.75-0.95- x向量維

x-for

x-

%下面構造

fork=1:j-for

L=L*(t-x(k))/(x(j)-L=L*(t-x(k))/(x(j)-

p

yl

)

江蘇科技大計算結果pp計算x=1.5處的插值>>>>subs(p,1.5)ans=畫出插值多項式圖形

>>>>>>>>江蘇科技大>>symsx

函數顯示方(x*exp(x)-換個顯示方式>>

xexp(x)------------江蘇科技大符號與數值間的轉對表達式進行計算格式:eval(f

還有>>>>syms>>>>>>f=

符號表達>>clear,>>clear,syms>>>>subs(f,x,pi/4)ans=y=

數值結說明:eval(s),wheresisastringtothestringasan江蘇科技大重要指令 R= %替換S中所有的符號R=subs(S,new) Rsubs(S,old,new)%用new替換S中old符號變量>>>>symsa>>subs(a+b,a,4)ans=4+b>>clear;syms???Undefinedfunctionorvariable>>syms>>subs(exp(a*t),a,-4)ans=exp(-4*t)a要定義江蘇科技大例>>symsthgA>>t1=solve('v*sin(A)*t-1/2*g*t^2=-t1t1=符號表達[1/2/g*(2*v*sin(A)-2*(v^2*sin(A)^2+2*g*h)^(1/2))]>>>>tdtd數值結江蘇科技大第2聯函數格式:y=inline('函數','變量1','變量例如f(t,a,bsin2t)e-atb|t>>y=inline('sin(t)^2*exp(-a*t)-yyInliney(t,a,b)=sin(t)^2*exp(-a*t)-下面計算y(1,2,3)的值>>ansans-江蘇科技大例如>>G=inline('GInlineG(a,x,y)=a*exp(x)>>ans-江蘇科技大第3節中,還有兩個演算函數eval:對字符串計算feval:對 江蘇科技大eval格式:y= str必須是字符串

最常用:求字符串的值表達式的任何合法指令 語句例>>a=eval(例>>a=eval('sqrt(2)')%a=>>eval('a=sqrt(2)')%a=江蘇科技大例題:計算字符串>>clear;>>結結果y=>>clear;該用法同該用法同符號表達式>>syms>>>>>>符號表y0.7071數值結結果:ans=

定義字符串f;下面計算江蘇科技大feval 求fun(x)feval('fun',x):求'fun'給定的函數值, 求fun(x)>>a=feval('sqrt'>>a=feval('sqrt',2)a= fun(xy,z)江蘇科技大函 格式:函數名=@(變量) 例如:>>f=@(x,y)x^2+y^2;%屬性是函數句柄>>f(1,2)%ans=%還可以輸入數組>>>>>>

不是這個表達式ans= 江蘇科技大的例子c=0.1:0.1:0.5;%c作為參數傳 %工作區中的變量a=-.3;b=-f=@(x)1./((x+a).^2+0.1)+1./((x+b).^2+0.4)-f([-101])%用 a=-.3;b=-f=@(x)1./((x+a).^2+0.1)+1./((x+b).^2+0.4)-江蘇科技大 計算

f(x)=cosx0.5+(1+x2),f(x)=ln(x)ea=學號后