实验七-常微分方程解法 幻灯片.pptVIP

数值实验六 数值积分与数值微分 实验报告要求 实验报告七必须在12月22日前上交 算法使用编程语言不限； 电子版试验报告与程序代码发送至邮箱： zhangkun@ahau.edu.cn 文件名称：实验七张三 Matlab基础（实验一） 基本命令 基本数学运算 符号解 函数文件格式 函数调用 文件保存位置 基本命令 1. help 帮助命令 help format %查看format命令的帮助； 2. clear 清除变量命令 clear x %清除变量x(无参数时清理工作空间） 3. 清理命令窗口 clc %清屏 4. format 输出格式 format long %15位有效数字格式输出 format short %4-7位有效数字格式输出 format compact %紧凑格式(输出不加空行） 基本数学运算 乘法：* 除法：/ 乘方：^ 根号：sqrt(x) 正弦：sin(x) 余弦：cos(x) 自然对数：log(x) 自然指数：exp(x) 符号解 1. vpa 控制运算精度 vpa(表达式,运算精度) vpa(pi,20) %显示?至100位有效数字 2. solve 方程的符号解 solve( 表达式 , x) 符号解 例 解方程 解 在MATLAB工作窗口输入命令: y=solve(x^3-sin(x)-12*x+1=0,x) vpa(y,5） y=roots([2,0,-1,-1]) 函数文件格式 函数文件由function语句引导，其格式为： function [输出形参表]=函数名(输入形参表) % 注释说明部分 函数体： 例： function [k,xk,yk,p]=jhnewtonqx(x0,ddmax) % 牛顿切线法求非线性方程的根 … 注：其中函数名的命名规则与变量名相同。输入形参为函数的输入参数，输出形参为函数的输出参数。当输出形参多于1个时，则应该用方括号括起来。 函数调用 函数调用的一般格式是： [输出实参1,实参2,...]=函数名(输入实参1,实参2,...） 例： [k,xk,yk,p]=jhnewtonqx(1.5,20) 文件保存位置 设置当前目录位置 cd C:\Matlab7\work %设置Matlab当前目录为： C:\Matlab7\work Matlab基础（实验二） 1. 矩阵输入 2. 特殊矩阵 3. 误差分析 4. 内置函数 矩阵输入 a=[] %建立空矩阵a ,可在workspace中编辑 a=[1,2,3;4 5 6] % [分行“,”或空格 分列 “;” ] a=1:5 % 建立1×5矩阵(向量) a=(1,2,3,4,5)’ a=0:pi/5:2*pi % (pi= ?) 建立1×11矩阵（向量）a ：初始元0，步长?/5 , 终止元2 ? 【注意 “ ：”的用法】 常用矩阵 1、零矩阵 z=zeros(3,4) %产生 3×4零矩阵z； zeros(5) %产生 5阶零矩阵； 2、单位矩阵 E=eye(6) %产生6阶单位矩阵E； 2、幺矩阵 ones(3,2); %产生3行2列元素为1的矩阵； 3、随机矩阵 r=rand(7) %产生0~1间分布的随机矩阵r s=round(rand(7)*30) %产生7阶0~30间均匀分部的随机矩阵s 特殊矩阵 1、Hilbert矩阵 h=hilb(5) %产生 5阶Hilbert矩阵h h=sym(hilb(5)) %产生 5阶Hilbert符号矩阵h 2、Vandermonde矩阵 v=vander([1,2,3,4,5]) %产生5阶Vandermonde矩阵v 3、魔方矩阵 m=magic(3) %产生3阶魔方矩阵 3、Toeplitz矩阵 t=toeplitz(0:-1:-5,0:5) 矩阵运算 a=magic(3);b=round(10*rand(3));r=(1:3); a+b %矩阵加法 a-b %矩阵减法 a*b %矩阵乘法 n=inv(a) % inv(a) =a^(-1)= a的逆 a/b % a/b=a*inv(b)=a*b^(-1) x=a\r,a*x % x=a\r=inv(a)*r=a^(-1)*r c=a % a= a的（共轭）转置