实验四 程序编写

《实验四 程序编写》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、实验4MATLAB程序编写一、实验目的1．学会编写简单函数m文件，调用及调试函数m文件；2．学会MATLAB的简单程序编写。二、实验设备1．电脑2．MATLAB软件三、实验内容1．从键盘输入一个4位整数，按如下规则加密后输出。加密规则：每位数字都加上7，然后用和除以10的余数取代该数字；然后将第1位数与第3位数交换，第2位数与第4位数交换。程序：a=input('请输入一个四位整数:');a1=fix(a/1000);%千a2=rem(fix(a/100),10);%百a3=rem(fix(a/10),10);%十a4=rem(a,10);%个%加密a1=a1+7;a

2、2=a2+7;a3=a3+7;a4=a4+7;a1=rem(a1,10);a2=rem(a2,10);a3=rem(a3,10);a4=rem(a4,10);m=a3*1000+a4*100+a1*10+a2;ifa3==0disp(['输出：','0',num2str(m)])elsedisp(['输出：',num2str(m)])end结果：请输入一个四位整数:1234输出：01892.硅谷公司员工的工资计算方法如下：(1)工作时数超过120h者，超过部分加发15%；(2)工作时数低于60h者，扣发700元；(3)其余按每小时84元计发。试编程按输入的工号和该员工

3、的工时数，计算应发工资。程序：clearall;clc;a=input('请输入工号:');b=input('请输入工时:');if60<=b&b<=120m=b*84;endifb<60m=b*84-700;endifb>120m=(b-120)*84*(1+0.15)+120*84;enddisp(['工号为',num2str(a),'的员工应发工资',num2str(m),'元'])结果：请输入工号:2453请输入工时:24工号为2453的员工应发工资1316元3．根据，求的近似值。当n分别取100，1000，10000时，结果是多少？（要求：分别用循环结构和向

4、量运算（使用sum函数）来实现）程序(1):clc;;clearall;n=input('输入n:');v=0;fori=1:nv=[v,1/(i^2)];endm=sum(v);disp(sqrt(6*m))%pi的近似值结果：输入n:10003.1406程序(2):clc;clearall;n=input('输入n:');v=0;fori=1:nv=v+1/(i^2);enddisp(sqrt(6*v))%pi的近似值结果：输入n:10003.14064．根据，求(1)y<3时的最大n值；(2)与(1)的n值对应的y值。考虑程序：clc;clearall;y=0;

5、n=0;while(y<=3)n=n+1;y=y+1/(2*n-1);endy=y-1/(2*n-1);n=n-1;disp(['y<3时的最大的n=',num2str(n)]);disp(['y=',num2str(y)]);结果：y<3时的最大的n=56y=2.99445．考虑一下迭代公式，其中，ａ，ｂ为常数。(1)编写程序求迭代结果，迭代的终止条件为，迭代初值，迭代次数不超过500次。(2)如果迭代过程收敛于r，那么r的准确值是，当(a,b)的值取(1,1)、(8,3)、(10,0.1)时，分别对迭代结果和准确值进行比较。程序：clc;clearall;a=in

6、put('请输入正数a:');b=input('请输入正数b:');n=0;%迭代次数x=1.0;%初值xl=1.0;while(n<=499)x=a/(b+x);n=n+1;h=abs(x-xl);ifh<=10^(-5)o=1;break;endxl=x;enddisp(['x=',num2str(x)]);disp(['n=',num2str(n)]);m=(-b+sqrt(b^2+4*a))/2结果：请输入正数a:123请输入正数b:134x=0.91171n=3m=0.9117