2matlab的數(shù)值計算-免費閱讀
【正文】 ,x0,y0,39。,5,5) 例 2： f(x)=100(x2x12)2+(ax1)2在 x1=a, x2=a2處有最小值 function f=xun(x,a) f=100*(x(2)x(1).^2).^2+(ax(1)).^2。 [t,x]=ode23(39。 2 3 4]。b=[1。 方程解 (a 39。) ans = x^4 + 2 x^3 + 3 x^2 + 4 x + 5 b=polyder(a) b = 4 6 6 4 poly2str(b,39。b=[4 5 6]。4 5 6。,pi) 3e+000 g 定長輸出符 sprintf(39。) f 浮點輸出符 sprintf(39。 for i=1:14,s(i)=setstr(s(i)(‘a’‘A’)),end,s 請考慮輸出 如果改為： s=‘I am a teacher39。i am a teacher39。 c1=a.\b。 a.*b ans = 2 8 18 4 15 30 49 72 90 a=[1 2 3。tril: 抽取主下三角 。4,5,6。b=[1。 z=[x。 （ 4）拼接 行拼接：兩個矩陣的行數(shù)必須相同 x=[1 2 3。4 5 6。7 8 9]。? 指令修改 可以用 A(?,?)= ? 來修改。 變量名盡可能不要重復，否則會覆蓋 。sqrt(3) 3+5i] 矩陣元素 符號的作用 逗號和分號的作用 ?逗號和分號可作為指令間的分隔符， matlab允許多條語句在同一行出現(xiàn)。 注意： matlab嚴格區(qū)分大小寫字母 ，因此 a與 A是兩個不同的變量。7 8 9]。4 5 6。4 5 6。4 5 6。 a=[1 2 3。2]。c(:)=a(:) 矩陣的變向 rot90:旋轉 。 b=[2 4 6。7 9 10]。I am a teacher39。))。) fprintf(obj,‘format’,‘cmd’) %obj可以省略 sscanf調用格式： A = sscanf(s,format) A = sscanf(s,format,size) sprintf調用格式： [s,errmsg] = sprintf(format,A,...) 控制符 符號 說 明 示 例 \b 退格符 sprintf(‘Who am I?\b’) （ Who am I） \n 回車換行符 sprintf(‘Who am I\n?’) \t 橫向走紙符 sprintf(39。%.2e39。 This is a good example.39。p=poly(a) p = r=roots(p) r = —— 顯然 r是矩陣 a的特征值 當然我們可用 poly令其返回多項式形式 p2=poly(r) p2 = matlab規(guī)定多項式系數(shù)向量用行向量表示，一組根用列向量表示 ,convs多項式乘運算 例 :a(x)=x2+2x+3。 命令格式： polyder(p): 求 p的微分 polyder(a,b): 求多項式 a,b乘積的微分 例： a=[1 2 3 4 5]。b=[8。 例 : x1+2x2=1 2x1+3x2=2 3x1+4x2=3 a=[1 2。 ? b x = x = 0 21xx321 = 433221 ax = b 當方程數(shù)少于未知量個數(shù)時 ,即不定 情況 ,有無窮多個解存在。f (xn,yn) n=0,1,2… y(x0)=y0 hyyxxyydtdy nnnnnn ????? ??? 111Runge Kutta法 龍格 庫塔法：實際上取兩點斜率 的平均 斜率來計算的，其精度高 于歐拉算法 。求解微分方程 t0=0。 plot(t,x)。yy=polyval(p,xx)。 當不能很快地求出所需中間點的函數(shù)時，插值是一個非常有價值的工具。,[0,0],[ ],[ ],sqrt(2)) 七、數(shù)據(jù)分析與插值函數(shù) max —— 各列最大值 mean —— 各列平均值 sum —— 各列求和 std —— 各列標準差 var —— 各列方差 sort —— 各列遞增排序 八、擬合與插值 1. 多項式擬合 x0=0::1； y0=[ ]。x(2)*(1x(1)^
環(huán)評公示相關推薦
鄂ICP備17016276號-1