matlab_4数值计算(阅读).pptVIP

创建数组和寻访 数组运算 关系与逻辑操作 第四章 MATLAB的数值计算 —— matlab 具有卓越的数值计算能力 基本内容 数值微积分 代数方程求解 概率分布和统计 多项式运算和卷积 创建矩阵(数组) 直接创建 函数创建rand eye zeros ones magic diagpascal peaks … 矩阵的修改 直接修改 寻访并修改(函数；全下标;单下标)矩阵运算 矩阵加、减（＋,－）运算 规则： ? 相加、减的两矩阵必须有相同的行和列两矩阵对应元素相加减。 ? 允许参与运算的两矩阵之一是标量。标量与矩阵的所有元素分别进行加减操作。 矩阵运算 2. 矩阵乘（?）除运算 规则： A矩阵的列数必须等于B矩阵的行数 标量可与任何矩阵相乘。 a=[1 2 3;4 5 6;7 8 0];b=[1;2;3];c=a*b c =143223 d=[-1;0;2];f=pi*d f = -3.141606.2832矩阵除的运算在线性代数中没有，有矩阵逆的运算，在matlab中有两种矩阵除运算\ / 矩阵除的运算在线性代数中没有矩阵除的运算，有矩阵逆的运算，在matlab中有两种矩阵除运算\ / 3. 矩阵乘方—— a^n,a^p,p^aa ^ p —— a 自乘p次幂a=[1,2,3;4,5,6;7,8,9];a^2ans =30 36 4266 81 96102 126 150 a^0.5ans = 0.4498 + 0.7623i 0.5526 + 0.2068i 0.6555 -0.3487i 1.0185 + 0.0842i 1.2515 + 0.0228i 1.4844 - 0.0385i 1.5873 - 0.5940i 1.9503 - 0.1611i 2.3134 + 0.2717i 4. 矩阵的其它运算 inv —— 矩阵求逆 det —— 行列式的值 eig —— 矩阵的特征值 diag —— 对角矩阵 ’ —— 矩阵转置 sqrtm—— 矩阵开方 5.矩阵的一些特殊操作 矩阵的变维 a=[1:12];b=reshape(a,3,4) c=zeros(3,4);c(:)=a(:) 关系运算 函数极值的数值求解 寻优函数： fminbnd — 求单变量函数 极小值 fminsearch— 单纯形法求多变量函数 极小值 fminunc — 拟牛顿法求多变量函数 极小值 fzero —— 变量函数0值 例　：f(x)=‘x2+3x+2’在[-5 5]区间的最小值 f=fminbnd(x^2+3*x+2,-5,5) 常微分方程数值解 微分方程求解的仿真算法有多种，常用的有Euler（欧拉法）、Runge Kutta(龙 格-库塔法。 Euler法称一步法，用于一阶微分方程 当给定仿真步长时： 所以yn+1 = yn + h·f (xn,yn) n=0,1,2…y(x0)=y0 Runge Kutta法龙格-库塔法：实际上取两点斜率 的平均 斜率来计算的，其精度高 于欧拉算法 。 龙格-库塔法：ode23 ode45例：x+(x2-1)x+x=0 为方便令x1=x，x2=x分别对x1,x2求一 阶导数，整理后写成一阶微分方程组 形式 x1=x2 x2=x2(1-x12)-x1 建立m文件 解微分方程 建立m文件 function xdot=wf(t,x) xdot=zeros(2,1) xdot(1)=x(2) xdot(2)=x(2)*(1-x(1)^2)-x(1) 给定区间、初始值;求解微分方程 t0=0; tf=20; x0=[0 0.25]; [t,x]=ode23(wf, 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);x(2)*(1-x(1)^2)-x(1)]; 求解微分方程 [t,x]=ode23(‘wf’,[0 30],[0;0.25]); plot(t,x); figure(2) plot(x(:,1),x(:,2)) matlab中有两种除运算左除和右除。 对于方程ax+b，a 为an×m矩阵，有三种情 况： ? 当n=m时，此方程成为“恰定”方程 ? 当nm时，此方程成为“超定”方程 ? 当nm时，此方程成为“欠定”方程matlab定义的除运算可以很方便地解上 述三种方程 1.恰定方程组的解 方程ax+b(a为非奇