MATLAB经典试题研讨.pptVIP

matlab 基本语句 1.循环语句forfor??i=s1:s3：s2????循环语句组end首先给i赋值s1；然后，判断i是否介于s1与s2之间；如果是，则执行循环语句组,i=i+s3（否则，退出循环.）；执行完毕后，继续下一次循环。 例：求1到100的和 sum=0???for?i=1:1:100?????sum=sum+i???end?注：for循环可以通过break语句结束整个for循环. .循环语句while 例：sum=0;i=1;??????while(i=100)???????sum=sum+i;i=i+1;??????end?????? Break Rem（）求除法后的剩余数 对于 x0，求 。（ 提示：理论结果为 ） 〖目的〗 符号变量的限定性定义的作用。 〖解答〗 syms k; x=sym(x,positive); f_k=2/(2*k+1)*((x-1)/(x+1))^(2*k+1); s=simple(symsum(f_k,k,0,inf)) %结果与理论值lnx相符！ s = log(x) 2.3.3 符号积分 【例2.3-9】求 。 intf=int(f,v) 给出f 对指定变量v的不定积分 intf=int(f,v,a,b) 给出f对指定变量v的定积分 返回从a到b的定积分f，积分变量是为v syms x f= x*log(x) s=int(f,x) s=simple(s) f = x*log(x) s = (x^2*(log(x) - 1/2))/2 s = x^2*(log(x)/2 - 1/4) 例2.3-11求积分 syms x y z; f=x^2+y^2+z^2; F2=int(int(int(f, z, sqrt(x*y), x^2*y), y, sqrt(x), x^2), x, 1, 2) 第一重积分 第二重积分 Warning: Explicit integral could not be found. F2 =(14912*2^(1/4))/4641 - (6072064*2^(1/2))/348075 + (64*2^(3/4))/225 + 1610027357/6563700 VF2=vpa(F2) %默认32位有效数字 VF2 = 224.92153573331143159790710032805 教材中48位有效数字 (1) 冒号生成法： array=a: inc: b 向量 a---数组的第一个元素　 inc---采样点之间的间隔, 即步长. 最后一个元素不一定等于b,其大小为b’=a +inc*fix((b-a)/inc); 步长可以省略, 默认为1; inc可以取正数或负数, 但要注意当取正时,要保证ba, 数组最后一个元素不超过b, 取负时ba, 最后一个元素不小于b. 特点: 等差数列 a b inc0 a b Inc0 3.2.1 一维数组的创建和寻访 1. 递增/递减一维数组的创建: b’ b’ (2) 线性(或对数)定点采样法 在设定的总点数下, 均匀采样生成一维行数组．格式为： 　　　　　x=linspace (a, b, n) 特点：方便对数组的元素的总个数进行控制，最后一个元素肯定是设定之值． a---数组初始值 b---数组最后值 n---数组总个数 增量的计算---inc=(b-a)/(n-1) 等效于: x=a: (b-a)/(n-1) :b 1. 递增/递减一维数组的创建: 　　　　　x=logspace (a, b, n) 要求在闭区间 上产生具有10个等距采样点的一维数组。试用两种不同的指令实现。 目的〗 数值计算中产生自变量采样点的两个常用指令的异同 %方法一 t1=linspace(0,2*pi,10) %方法二 t2=0:2*pi/9:2*pi %要注意采样间距的选择，如这里的2*pi/9. t1 = Columns 1 through 7 0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 Columns 8 through 10 4.8869 5.5851 6.2832 t2 = Columns 1 through 7 0 0.6981 1.3963 2.0