【正文】 p=polyfit(x0,y0,3) p = xx=0::1。 建立 m文件 function xdot=wf(t,x) xdot=zeros(2,1) xdot(1)=x(2) xdot(2)=x(2)*(1x(1)^2)x(1) 給定區(qū)間、初始值 。?a) ? a39。2 3]。7 8 0]。,pi) 3 e 指數輸出符 sprintf(39。A39。1 3 5。b=reshape(a,3,4) c=zeros(3,4)。 二、矩陣運算 2. 矩陣乘（ ?） 運算 規(guī)則： A矩陣的列數必須等于 B矩陣的行數 標量可與任何矩陣相乘。2]) （ 3）刪除 x=[1 2 3。4 5 6。4 5 6] x=[2 pi/2。 例如 a=[1 2 0。7 8 9]。y] z = 1 2 3 4 5 6 7 8 9 1 1 1 6. 矩陣函數 7. 矩陣分解 1. 矩陣加、減（＋ ,－）運算 規(guī)則： ? 相加、減的兩矩陣必須有相同的行和列兩矩陣對應元素相加減。7,8,9]。4 5 6。 for i=1:14,s(i)=setstr(s(i)(39。%.2f39。7 8 0]。x39。2。 x=fmins(39。b39。 x0=[0 ]39。 x1+2x2+3x3=1 2x1+3x2+4x3=2 a=[1 2 3。 ?x=inv(a)*b ? x=a\b x = x = 2 3 322121xx138 = a x = b 例 : x1+2x2=8 2x1+3x2=13 左除、右除的異同 方程 ax=b ,mn時此時不存在唯一解。 c = (x2+2x+3)(4x2+5x+6) a=[1 2 3]。%.0e39。i am a teacher39。b=[4 5 6]。 flipud:下翻 矩陣的抽取 diag:抽取主對角線 。7 8 0]。 x(:,1)=[ ] x = 2 3 5 6 8 9 size(x) ans = 3 2 Matlab中只能刪除矩陣的整行或整列， x(:,1)=[ ]將矩陣 x的第一列刪除，同樣， x(1,:)=[ ]將矩陣 x的第一行刪除。4 5 6。 注意：只要是賦過值的變量，不管是否在屏幕上顯示過，都存儲在工作空間中，以后可隨時顯示或調用。 matlab函數名必須小寫 。7 8 9]。7 8 9]。f=pi*d f = 0 矩陣除的運算在線性代數中沒有，有矩陣逆的運算，在 matlab中有兩種矩陣除運算 a ^ p —— a 自乘 p次冪 方陣 1的整數 3. 矩陣乘方 —— a^p,p^a 對于 p的其它值 ,計算將涉及特征值 和特征向量，如果 p是矩陣， a是標量 a^p使用特征值和特征向量自乘到 p次 冪；如 a,p都是矩陣， a^p則無意義。1 3 5。 s s = I am a teacher u=abs(s) u = 73 32 97 109 32 97 32 116 101 97 99 104 101 114 8. 文本操作 str1=s(5:10) str1 = a tea y=setstr(u) %函數 setstr y = I am a teacher u=s(10:1:5)輸出內容？ u=s’輸出內容？ 字符串 I’m a teacher如何輸入？ 可以像數組一樣連接字符串 可以用 disp函數打印字符串變量表示的字符串，如 disp(x)。Who am I?\tadf39。[token,rem] = strtok(s) token =This rem = is a good example. matlab語言把多項式表達成一個行向量， 該向量中的元素是按多項式降冪排列的。 poly2str(a,39。2 3。 龍格 庫塔法： ode23 ode45 211 2121 kkyynn ????k1=hf(xn,yn) k2=hf(xn+h,yn+k) 例： x+(x21)x+x=0 為方便令 x1=x， x2=x分別對 x1,x2求一 階導數，整理后寫成一階微分方程組 形式 x1=x2 x2=x2(1x12)x1 1. 建立 m文件 2. 解微分方程 figure(2) plot(x(:,1),x(:,2)) 2 1 0 1 2 32101