杂烩收集.docVIP

例2 用双重循环编写程序生成4阶Hilbert矩阵. 《Mobile Fading Channels: Modelling, Analysis, Simulation》 for H end end while循环语句适用于由条件控制的循环算法的实现.反复执行某些操作，直到条件不满足为止.其具体格式为 while〈条件〉， 循环 end 例3 产生裴波拉奇数列：1，1，2，5，8，13，21，….要求数列之中最后一个数不超过1000.数列产生的递推公式为： [1 1]； I 1； while , end 下面是将for语句和if语句结合起来一达到特殊的目的的一个例子. 在实现循环算法时常用到转移语句.Matlab中的转换语句是以if〈条件〉break，end的格式来定义的. 例4 用if…break语句跳出while循环.计算无穷级数的近似值.要求误差不超过. while if end s 这段程序执行的结果，，误差界限不超过.如果没有转移语句，程序见按照公式求和并一直进行下去永不休止，形成死循环. 3、switch——case——end分支结构 switch语句是多分支选择语句，虽然在某些场合switch的功能可以由if语句的多层嵌套来完成，但是会使程序变得复杂和难于修改维护，而利用switch语句构造多分支选择时显得更加简单明了、容易理解。 switch的基本用法为： switch 表达式 case 常量表达式1 语句块1 case 常量表达式2 语句块2 case 常量表达式n，常量表达式n+1, 常量表达式n+2… 语句块n otherwise 语句块n+1 end switch 后面的表达式可以为任何类型，如字符串、矩阵等 当表达式的值与case后面的某个表达式的值相等时，就执行这个case后面的语句块，如果所有的常量表达式的值都与此表达式的值不匹配时，就执行otherwise后面的语句块。 与if语句不同的是，各个case和otherwise语句出现的先后顺序并不会影响程序的执行结果，将otherwise语句放到case语句的前面也是合法的。 每个case后面的常量表达式可以有多个，而且可以是不同的类型，每个case中常量表达式的值是可以重复的，这在语法上没有错误，只是在执行时后面那些符合条件的case语句将被忽略，不会起作用。 每次只执行一个语块，执行完一个语句块后就退出switch语句。 例5 我国新税法规定：个体工商户的生产、经营所得和对事业单位的承包经营、承租经营所得应交纳的个人所得税为：全年收入中应纳税所得税部分 税率（%） 1、不超过5000元的 5 2、超过5000元至10000元的部分 10 3、超过10000元至30000元的部分 20 4、超过30000元至50000元的部分 30 5、超过50000元的部分 35 function y shui x x input input x: n fix x/1000 ; switch n case 0,1,2,3,4 y x*0.05; case 5,6,7,8,9 y x*0.1; case 10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29 y x*0.2; case 30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49 y x*0.3; otherwise y x*0.35; end 38 Matlab中如何求解广义积分？即积分限到有无穷的或者有歧 异点的积分 瑕积分 ？ Matlab的quad/quad8只能作定积分，广义积分需要自己来写程序逼近， 流程大概如下： 1.设定收敛限epsi 2.把为inf/-inf或者歧义点的积分限设置为一个初始值,k 1 3.计算定积分Q k 4.然后朝着inf/-inf或者歧义点移动一个步长，然后计算定积分 Q k+1 5.判断 abs Q k+1 -Q k QuasiMonteCarlo或者Oscillatory]。广义积分数值求解是个较困难的问题，通常数值积分方法都不适用。若先编一个M函数 function y fun6 x y exp sin x -x.^2/50 ; 再用命令 quad8 fun6,-inf,inf 结果为 ans NaN 是无穷大，错误。为能正确求解，考虑到 但问题是去多大？比如取，用命令 quad8 fun6,-1e+10,1e+10 结果为 ans 6.8135e+005 是一个明显错误的大数。正确的方法是，先取一个适当大的（比如10），计算一个值；然后将以一适当的倍