数学实验13-matlab编程

1、 数学实验理学院数学学科 李换琴 信息01、02、03班上机时间安排10月15日周五晚上18：3022：0010月29日周五上午10：0012：0011月21日周日上午08：3012：0012月01日周三上午08：0010：00信息04、05、06班上机时间安排10月16日周六上午08：3012：0010月29日周五晚上18：0020：0011月21日周日下午14：3018：0012月04日周六上午10：0012：00计算机01、02、03班上机时间安排10月10日周日上午08：3012：0010月29日周五上午8：009：5011月13日周六上午08：3012：0012月01日周三上午10：

2、0012：00计算机04、05、06班上机时间安排10月10日周日下午14：3018：0010月29日周五晚上20：0022：0011月13日周六晚上18：3022：0012月04日周六上午08：0010：00MATLAB语言编程介绍MATLAB中各种命令可以完成许多单一的任务，对于某些较为复杂的问题，仅靠现有的命令或函数来解决，往往是难以达到目的 。为此，要运用MATLAB编程语言编制程序，形成M-文件。程序是使计算机完成各项运算的命令集，运行一个编制好的程序，计算机会从第一条命令行开始，一行接一行地执行相应的命令，直到终止。程序编写调试完成后，需要存盘，形成永久性文件，可以随时对它进行调用

3、或修改。文件名以字母开头，但不能用专用变量名，如pi,ans,eps等。 程序结构程序一般包含：a=input(请输入数组 an= ); m,n=size(a); M=a(1); for i=2:nif M n2.省略格式： for i=n1:n2这里step=1.例1：求n个奇数和：s=1+3+5+(2n-1) n=input(please input n=); s=0; for i=1:n s=s+(2*i-1); fprintf(i=%.0f, s=%.0fn,i,s); end 程序aa1例2：求阶乘：p=12 3 n=n! n=input(请输入 n= ); p=1; for i=1

4、:n p=p*i; fprintf( i=%.0f, p=%.0fn ,i,p); end aa2.m例3：求e：e=1+1+1/2!+1/3!+1/n! n=input(请输入 n= ); p=1;e=1; for i=1:n p=p*i; p1=1/p; e=e+p1; fprintf( i=%.0f, p=%.0f, e=%.8f n ,i,p,e); end aa3.mM-文件中条件循环命令(while命令)格式： while (condition is true) commands; end 作用：当条件成立时，执行命令集commands,直到条件不成立.例4：求和：s=1+2+n3

5、000 s=0;k=1; while s+k3000 s=s+k; fprintf( k=%.0f, s=%.0fn ,k,s); k=k+1; end aa4.m例5：求e：e=1+1+1/2!+1/3!+ (1/n!=1.0e-8p=p*i; r=1/p; e=e+r; fprintf( i=%.0f, r=%.9f, e=%.9f n ,i,r,e); i=i+1; end aa5.mM-文件中选择控制命令(if命令)单项选择控制格式： if (condition is true) commands; end 作用：若条件成立，则执行命令集 commands. 否则，不执行。例6：求n个

6、实数中最大的数M. a=input(请输入数组 an= ); m,n=size(a); M=a(1); for i=2:n if M0 z=1; elseif x=0 z=0; else z=-1; end y=x2*z aa7.m用M-文件定义函数格式 function y1,y2=funname(x,y,a,n)其中： function 为关键字； y1,y2,为输出列表； funname为自定义的函数名； (x,y,a,n,)为函数参数，用来传递相 关数值. 存盘时，要求用funname作为函数名. 例7：建立符号函数sgn(x) function sn=sgn(x) if x0 sn=

7、1; elseif x=0 sn=0; else sn=-1; end以sgn作为文件名存盘，即建立了函数。调用： 在命令区执行 : sn=sgn(10)或sn=sgn(-2)例8：建立计算(a+b)n,(a-b)n的函数 function p,q=fun1(a,b,n) p=(a+b).n; q=(a-b).n; fun1.m 例编制程序，从键盘输入a,b,n,计算(a+b)n,(a-b)n aa8.ma=input(a=);b=input(b=);n=input(n=)；p,q=fun1(a,b,n);fprintf(a+b)n=%.4f,(a-b)n=%.4fn,p,q)例10：用二分法

8、求函数x2-2=0的正实根.输入初值：输出结果：迭代次数k,方程的近似根c.NoNoNoYesYesYesa=1;b=2;er=1.0e-8;f=inline(x.2-2);k=1;c=(a+b)/2;r=c+1;while abs(r-c)err=c;if abs(f(r)1.0e-10fprintf(k=%.0f,r=%.8fn,k,r);break;endif f(a)*f(r)0 b=r;else a=r;endc=(a+b)/2;fprintf(k=%.0f,r=%.8fn,k,r);k=k+1;enderfenfa.m本次实验三人一个小组，合作完成任务，交一份实验报告。任务：1，2

9、任选1题，3，4任选1题， 练习1第1题任选1小题。实验报告内容：1）任务描述； 2）问题分析及程序（包括符号说明）； 3）运行结果 ；4）结果分析与说明。实验报告提交时间：待下次任务完成以后，将两次内容一并打印，封皮上写清楚班级，学号，姓名。课代表收齐于17周，18周周五晚上7点半至9点半交到我办公室（理科楼225）。1. 每门课程考试阅卷完毕,任课教师都要对各班的考试成绩进行统计,统计内容包括:全班人数,总得分,平均得分,不及格的人数及90分(包括90分)以上的人数.请编制程序解决这一问题,并自给一组数据验证程序的正确性.要求:使用者在提示下通过键盘输入学生成绩,计算机自动处理后,显示需要

10、的结果.2. 根据我国个人所得税计算方法,编制程序,要求:使用者在系统提示下通过键盘输入月工资薪金收入总数,计算机则在屏幕上显示个人所得税额,界面友好,方便使用.个人所得税计算方法: 月个人所得税（月工资薪金收入1600）*适用税率速算扣除数附表：个人所得税税率表（工资、薪金所得适用）个人所得税税率表（工资、薪金所得适用）级 数 全月应纳税所得额 税率（）税率 速算扣除数1 不超过500元的 52 超过500元至2000元的部分 10 253 超过2000元至5000元的部分 15 1254 超过5000元至20000元的部分 20 3755 超过20000元至40000元的部分 25 137

11、56 超过40000元至60000元的部分 30 33757 超过60000元至80000元的部分 35 63758 超过80000元至100000元的部分 40 103759 超过100000元的部分 45 15375首先由计算机随机产生一个 1,100 之间的一个整数，然后由用户猜测所产生的这个数。根据用户猜测的情况给出不同的提示，如果猜测的数大于产生的数，则显示 “High” ，小于则显示 “ Low ” ，等于则显示 “You won！”，同时退出游戏。用户最多有 7 次机会。3. 编写程序实现猜数游戏4. 编写程序完成复印资料给出的问题。a=1.;b=2.;r=1.0e-7;ya=fc1(a);yb=fc1(b);k=1;if abs(ya)r root=a; fprintf(k=%.0f, Root x=%.4f, f(x)=%12.