数学软件与实验课后题.doc

数学软件与数学实验作业 3．计算，要求： （1）用循环结构并使用上体的M函数文件 ①s=0; for k=1:10 f=f1(k); s=s+f; end s 结果：s = 4500244 ②function f=f4(n) f=0; for i=1:n f=f+f1(i); end （2）用help命令查询sum和prod的用法，并用这两个指令完成计算。 n=input(enter a number:); sum=0; for i=1:n sum=sum+prod([1:i]); end sum 4．用尽可能多的方法编写Matlab函数mysum1.m来计算和式f(n)=.其中，n为函数的输入变量，函数输出为f（n）的值，分析指出他们执行效率的高低，给出理由。 ①n=input(enter a number:); s=0; for i=1:n t=log(1+1/(i*i)); s=s+t; end s ②function f=mysum1(n) f=0; k=0; while kn k=k+1; f=f+log(1+1/(k*k)); end ③function f=mysum11(n) if n==1 f=log(2); elseif n1 f=mysum11(n-1)+log(1+1/(n*n)); else error(wrong input parameter!); end ④function f=mysum12(n) k=1:n; s=log(1+1./(k.*k)); f=sum(s); 5．、用尽可能多的方法编写Matlab函数mysum1.m来计算如下和式： f(n)=1/(1*4)+1/(4*7)+…+1/((3n-2)*（3n+1）) 其中，n为函数的输入变量，函数输出为f(n)的值。 ①function f=mysum2(n) f=0; for k=1:n t=1/((3*k-2)*(3*k+1)); f=f+t; end ②function f=mysum21(n) f=0; k=0; while kn k=k+1; t=(3*k-2)*(3*k+1); f=f+1/t; end ③function f=mysum22(n) if n==1 f=1/4; elseif n1 t=(3*n-2)*(3*n+1); f=mysum22(n-1)+1/t; else error(wrong input parameter!); end ④function f=mysum23(n) k=1:n; t=(3.*k-2).*(3.*k+1); f=sum(1./t); 6．写一个程序用于得到阶乘不超过10^100的最小整数。 for k=1:5000000000000 f=1; for i=1:k f=f*i; end if(f=10^100) continue; end break; end k 结果：k = 70 7．写一个函数rs=f(s),对输入放入字符串变量s,删除其中的小写字母，然后将原来的大写字母变为小写字母，得到rs的返回 function rs=fy(s) k=length(s); p=[]; for i=1:k if s(i)=as(i)=z p=[p,i]; elseif s(i)=As(i)=Z s(i)=char(s(i)+32); end end s(p)=[]; rs=s; 8. 设计函数min_element.m,其功能是在一个二维矩阵中找出其最小元素，函数定义如下：[minEle,row,column]=min_element(matrix) 其中，matrix是一个二维矩阵，minEle是返回的最小元素的值，row是其所在行，coulmn是其所在列。要求： 在调用函数时，要写出调用的语句； 不适用Matlab函数min编写； (3)使用min函数和find函数编写。 function [min,row,column]=min_element(x) n=size(x); min=x(1); for i=1:n(1) for j=1:n(2) if x(i,j)min min=x(i,j); [row,column]=find(x=min); end