MATLAB

第二题写一个计算n!的M 函数文件
function m=myprod(n)
m=1;
for k=1:n
m=m*k;
end
第三题计算∑=10
1!k k ，要求：
（1）用循环结构并使用上题的M 函数文件；
（2）用help 命令查询sum 和prod 的用法，并用这两个指令完成计算；
①s=0
for k=1:10
s=s+myprod(k);
end
s
②for k=1:10
h=1:k;
p(k)=prod(h);
end
s=sum(p)
第四题：用尽可能多的方法编写Matlab 函数mysum1.m 来计算和式f(n)=∑=n i 1ln(1+2
1i )..其中，n 为函数的输入变量，函数的输出值为f(n)的值。

分别指出他们执行效率的高低，给出理由
① function f=mysum1(n)
m=0;
k=0;
while k<n
k=k+1;
m=m+log(1+1/(k*k));
end
m
②function f=mysum1(n)
f=0;
for k=1:n
h=1:k;
p(k)=prod(h);
f=sum(p);
end
f
③ function f=mysum1(n)
if n==1
f(1)=log(2);
else if n>1
f(n)=mysum1(n-1)+log(1+1/n^2);
数学软件与数学实验作业
专业：10级信息与计算科学 班级 一 班 姓名李亚伟 学号 10140201036
else
error('wrong input parsmeter')
end
第五题用尽可能多的方法编写Matlab函数mysum2.m来计算和式：其中，n为函数的输入变量，函数的输出值为f(n)的值。

①function m=mysum2(n)
n=input('请输入一个整数：')
f=0;
for k=1:n
f=f+1/((3*k-2)*(3*k+1));
end
f
②function m=mysum2(n)
n=input('请输入一个整数')
m=0;
k=0;
while k<n
k=k+1;
m=m+1/((3*k-2)*(3*k+1));
end
m
③ function m=mysum2(n)
n=input('请输入一个整数')
if n==1
m=1/4;
else if n>1
m=mysum2(n-1)+1/((3*k-2)*(3*k+1));
else
error('wrong input parameter!')
end
m
④ function m=mysum2(n)
n=input('请输入一个整数')
k=1:n;
s=1./((3.*k-2).*(3.*k+1));
m=sum(s);
⑤ function m=mysum2(n)
n=input('请输入一个整数')
syms k;
s=symsum(1./((3.*k-2).*(3.*k+1)),1,n);
m=double(s);
10的最大整数。

第六题.写一个程序用于得到阶乘不超过100
k=1;
m=1;
while m<=10
k=k+1;
m=m*k;
end
k-1
第七题.写一个函数rs=f(s),对输入的字符串s,删除其中的小写字母，然后将原来的大写字母变为小写字母，得到rs的返回。

例如，rs=”aBcdE,Fg?”,则rs=”be,f?”.
function rs=f(s)
s=input('请输入一串字符:',’s’);
k=length(s);
p=[];
for i=1:k
if s(i)>='a'&s(i)<='z'
p=[p,i];
elseif s(i)>='A'&s(i)<='Z'
s(i)=char(s(i)+32);
end
end
s(p)=[];
rs=s;
第八题.设计函数min_element.m其功能是在一个二维矩阵中找到其最小元素，函数定义如下：
[minEle,row,column]= min_element(matrix)
其中，matrix是一个二维矩阵，minEle是返回的最小元素的值，row是其所在行，column是其所在列。

要求：
（1）在调用函数时，要写出调用的语句；
（2）不使用Matlab函数mia编写；
（3）使用min函数和find函数编写。

function [minEle,row,column]=min_element(matrix)
n=size(matrix);
minEle=matrix(1,1);
row=1;
column=1;
for i=1:n(1)
for j=1:n(2)
if matrix(i,j)<minEle
row=i;
column=j;
minEle=matrix(i,j);
end
end
end。