西北工业大学《MATLAB基础与应用》课件-第7章数值计算基本概念.pptxVIP

Matlab与数值分析西北工业大学《MATLAB基础与应用》第7章 数值计算基本概念数的表示方法定点形式（定点数)R进制定点数形式： (αnαn-1…α1α0.β1β2…βm) R=αnRn + αn-1Rn-1 + … + α0R0 + β1R-1 + β2R-2 + … + βmR-m 其中左边为整数部分，右边为小数部分其中每个数字αi , βj 是介于0与 R-1之间的正整数；上边是数的书写形式，下边是该数所表达的值，它等于每位数字与其单位乘积之和，该数的总位数称为字长。数的表示方法浮点形式（浮点数)R进制浮点数形式： Rp ?(0.d1d2…dm) R=Rp ?(d1R-1 + d2R-2 + … + dmR-m) 其中 P 称为阶码，(0.d1d2…dm) R称为尾数其中0 ≤di ≤ R-1, i=0,1,…,m 当d1≠0时，称为规格化数当d1=0时，称为不规格化数尾数的位数称为字长数的表示方法例：十进制数的表示(37.21829)10=102(0.3721829)10=103(010 定点数 浮点规格化数 浮点非规格化数采用浮点形式表达数，可以拓宽数值的表达范围，化简特大或特小的定点数的写法；同一个数其进位制的基数R越大，数的字长就越短，但运算规则越复杂，反之亦然。数的表示方法二进制记数法（计算机中采用二进制表示法）国际通用标准，双精度二进制64位(8字节)，其中11个指数位，52个数值位和1个符号位符号阶码尾数11位1位52位其中：d=0或1，d1?0，s为阶码数的表示方法浮点数标准：IEEE 754浮点数：IEEE 754（1985）: （可在google上下载该文档）；在Matlab帮助文件搜索栏中输入：“IEEE 754”，可得到IEEE 754的介绍。数的表示方法Matlab所有数值都使用64bit表示，根据IEEE 754标准规定，表示形式如下：第1位是符号指示位(S) ；接下来11位是指数部分(C) ；最后成为尾数的52位小数部分(F) ；指数的基是2；符号位(S)指数位(S)尾数位(S)1位11位52位+/--1023-102410^(-308)-10^(308)数的表示方法例：计算以下64位二进制数对应的十进制0 100000000011 1011100100010000000000000000000000000000000000000000解：符号位S=0,所以为整数；指数位C=1*2^10+0*2^9+…+1*2^0=1024+2+1=1027所以指数部分为2^(1027-1023)=2^4;尾数位F=1*0.5^1+0*0.5^2+1*0.5^3+…+0*0.5^52=2961/4096所以对应的十进制数为数的表示方法例：计算以下64位二进制数对应的十进制0 100000000011 1011100100010000000000000000000000000000000000000000解：clear all;x64=0100000000011110000000000000000000000000000000000000000000000000;s = bin2dec(x64(1))c = bin2dec(x64(2:12))m = bin2dec(x64(13:64))f = sym(0.5).^(52)*mx = (-1)^(s)*2^(c-1023)*(1+f)数的表示方法例：将“7.5”转换为64位浮点数二进制表示解：S=0e=2+1023=1025f = 0.5+0.25+0.1250 10000000001 1011100100010000000000000000000000000000000000000000clear all;close all; clcformat long % 浮点数display(浮点数的二进制表示)x = 7.5fid = fopen(conv.bin, w); fwrite(fid, x, double);fclose(fid);fid = fopen(conv.bin, r);B = fread(fid, 1, ubit64);fclose(fid);Bin = dec2bin(B, 64)数的表示方法浮点数的性质浮点数是什么数？实数？ 有理数？ 有限小数？有多少个不同的浮点数？264（64位，每位有两个状态，“0”和“1”）浮点数是由264个有限小数（包含整数）构成的集合？错。IEEE 定义了一些异常值，inf （无穷）和 NaN（“非数字”）浮点数精度是多少？eps=2-52 = 2.2204E-16最大的浮点数是多少？realmax =（2-esp）×21023