2matlab的數值計算(留存版)

　　

【正文】 230 5 10 15 20 25 303210123六、函數優化 尋優函數： fmin —— 單變量函數 fmins —— 多變量函數 constr —— 有約束條件 無約束條件 例 1： f(x)=‘x2+3x+2’在 [5 5]區間的最小值 f=fmin(39。) 插值的定義 —— 是對某些集合給定的數據點之間函數的估值方法。, t0, tf, x0) plot(t,x), figure(2),plot(x(:,1),x(:,2)) 命令格式 : [T,Y] = ODE23(ODEFUN,TSPAN,Y0) 建立 m文件 function dxdt=wf(t,x) dxdt=[x(2)。2]。 b x=(a39。x39。,pi) _ _ _ _ _ _ _ sprintf(39。end,s 請考慮輸出 （ 2）字符串轉換 函數 意義 abs 字符串到 ASCII轉換 dec2hex 十進制到十六進制轉換 fprintf 把格式化的文本寫到文件中或顯示屏上 fscanf 從文件中讀入有格式數據 hex2dec 十六進制字符串轉換成十進制數 hex2num 十六進制字符串轉換成 IEEE浮點數 int2str 整數轉換成字符串 lower 字符串轉換成小寫 num2str 數字轉換成字符串 setstr ASCII轉換成字符串 sprintf 把數據格式化，寫給一個字符串 sscanf 用格式控制，把字符串轉換成數字 str2mat 字符串轉換成一個文本矩陣 str2num 字符串轉換成數字 upper 字符串轉換成大寫 rad=。b=[4 5 6]。 a=[1 2 3。3]。7 8 9]。4 5 6。 ?用于選出矩陣指定行、列及元素。 ?循環語句 matlab函數創建矩陣 空陣 [ ] — matlab允許輸入空陣，當一項操作無結果時，返回空陣。7 8 9]。 y=ones(3)。c=a*b c =14 32 23 d=[1。4 5 6。 z=a.^2 z = z=a.^b z = （ 5）向量函數 一般標量函數都可以用于數組，此時函數作用于數組的每個元素： x=1:5。 t=[‘A circle of radius ’,rad] %是否正確 fscanf調用格式： A = fscanf(fid,format) [A,count] = fscanf(fid,format,size) fprintf調用格式： fprintf(obj,39。%39。) p = 4 x^4 + 13 x^3 + 28 x^2 + 27 x + 18 a=[1 2 3]。 a)1 a 39。 x=a\b x=pinv(a)?b x = x = 0 0 432321321xxx21= ax = b 五、微分方程求解 微分方程求解的仿真算法有多種，常用的有 Euler（歐拉法）、 Runge Kutta(龍 格 庫塔法。x(2)*(1x(1)^2)x(1)]。 當不能很快地求出所需中間點的函數時，插值是一個非常有價值的工具。 plot(t,x)。f (xn,yn) n=0,1,2… y(x0)=y0 hyyxxyydtdy nnnnnn ????? ??? 111Runge Kutta法 龍格 庫塔法：實際上取兩點斜率 的平均 斜率來計算的，其精度高 于歐拉算法 。 例 : x1+2x2=1 2x1+3x2=2 3x1+4x2=3 a=[1 2。 命令格式： polyder(p): 求 p的微分 polyder(a,b): 求多項式 a,b乘積的微分 例： a=[1 2 3 4 5]。 This is a good example.39。) 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。I am a teacher39。 b=[2 4 6。2]。4 5 6。4 5 6。 注意： matl