MATLAB基础与应用教程(人民邮电-蔡旭辉)第三章作业参考程序

1、编写一个求圆的面积的函数文件。functions=area(r)%AREAForcalculatingtheareaofaround.%area.ms=pi*r^2;2、编写一个求圆的面积的命令文件Radius=input('Pleaseinputtheradius:');Area=pi*Radius^2;上机练习（5）3、4位整数从键盘输入一个，按如下规则加密后输出。加密规则：每位数字都加上7，然后用和除以10的余数取代该数字；然后将第1位数与第3位数交换，第2位数与第4位数交换。参考程序：x=input('输入一个4位整数：');y=[fix(x/1000),rem(fix(x/100),10),rem(fix(x/10),10),rem(x,10)]z=rem((y+7),10)x=z(3)*1000+z(4)*100+z(1)*10+z(2)4、根据，求：⑴y3时的最大n值；⑵与⑴的n值对应的y值。参考程序：y=0;k=0;whiley3k=k+1;y=y+1/(2*k-1);endn=k-1，y=y-1/(2*k-1)1111352-1yn5、已知，⑴当时，y的值是多少？⑵当时，y的值是多少？参考程序：y=fun(40)/(fun(30)+fun(20))（1）函数文件fun.mfunctionf=fun(n)f=n+log(n^2+5);403020fyff210ln5fnnn12+23+34+1fnnn5、已知，⑵当时，y的值是多少？参考程序：（2）函数文件fun.mfunctionf=fun(n)a=1:n;f=sum(a.*(a+1));或functionf=fun(n)ifn==1f=1*2;elsef=fun(n-1)+n*(n+1);end403020fyff12+23+34+1fnnn6、Fibonacci数列定义如下：f1=1，f2=1，fn=fn-1+fn-2(n2)。编写一个函数文件，求Fibonacci数列的第k项。参考程序：functionfk=fibon(k)%fibon:计算Fibonnaci数列的第k项ifk==1|k==2fk=1;elsefk=fibon(k-2)+fibon(k-1);end7、通过M脚本文件，画出下列分段函数所表示的曲面15457.0117575.015457.0),(215.175.375.0216215.175.375.02112122212212122xxexxexxexxpxxxxxxxx参考程序cleara=2;b=2;clf;x=-a:0.2:a;y=-b:0.2:b;fori=1:length(y)forj=1:length(x)ifx(j)+y(i)1z(i,j)=0.5457*exp(-0.75*y(i)^2-3.75*x(j)^2-1.5*x(j));elseifx(j)+y(i)=-1z(i,j)=0.5457*exp(-0.75*y(i)^2-3.75*x(j)^2+1.5*x(j));elsez(i,j)=0.7575*exp(-y(i)^2-6*x(j)^2);endendendaxis([-a,a,-b,b,min(min(z)),max(max(z))]);colormap(flipud(winter));surf(x,y,z);运行结果