数学实验答案2

%根据以下公式计算y，计算结果作为函数值返回；n通过形参传入
function y=fun1(n)
clc
%n=11;
%ans=1.8333……
y=0;
for k=1:n
fm=0;
for g=1:k
fm=fm+g;
end
y=y+1/fm;
end

%计算并输出给定数组s(长度为9)中每相邻两个元素之平均值的平方根之和。
function y=fun2(s)
clc
y=0;
%s=[12 34 4 23 34 45 18 3 11];
%ans =35.9510
for k=1:length(s)-1
y=y+sqrt((s(k)+s(k+1))/2);
end

%求小于形参n同时能被3与7整除的所有自然数之和的平方根，并作为函数值返回。
function y=fun3(n)
clc
%n=1000;
%ans=153.9091
N=1:n;
f=find(rem(N,3)==0 & rem(N,7)==0);
y=0;
for k=1:length(f)
y=y+f(k);
end
y=sqrt(y);

%根据以下公式求P的值
function y=fun4(m,n)
clc
%m=12;n=8;
%ans=495
fz=1;fm1=1;fm2=1;
for k=1:m
fz=fz*k;
if k<=n
fm1=fm1*k;
end
if k<=m-n
fm2=fm2*k;
end
end
y=fz/fm1/fm2;

%将字符串s中的前导*号全部删除，中间和尾部的*号不删除。
function y=fun5(s)
clc
%s='*******A*BC*DEF*G****';
%ans=A*BC*DEF*G****
for k=1:length(s)
if s(k)~='*'
t=k-1;
break
end
end
s(1:t)=[]

%求出向量s的最大元素在数组中的下标并存放在k所指的存储单元中。
function y=fun6(s)
clc
%s=[876 675 896 101 301 401 980 431 451 777];
%k=7 m=980
m=s(1);
for g=1:length(s)
if s(g)>m
m=s(g);
k=g;
end
end
k
m

%将两个两位数的正整数a、b合并形成一个整数放在c中。
%合并的方式是：将a数的十位和个位数依次放在c数的百位和个位上
% b数的十位和个位数依次放在c数的千位和十位上。
function y=fun7(a,b)
clc
%a=45;b=12;
%ans=1425
qw=fix(b/10)*1000;
bw=fix(a/10)*100;
sw=rem(b,10)*10;
gw=rem(a,10);
y=qw+bw+sw+gw;

%计算并输出给定整数n的所有因子(不包括1与n自身)之和。
%规定n的值不大于1000。
function y=fun8(n)
clc
%n=856;
%ans=763
s=0;
for k=2:n-1
if rem(n,k)==0
s=s+k;
end
end
s

%将大于形参 m 且紧靠 m 的 k 个素数存入x所指的数组中。
function y=fun9(m,k)
clc
%m=17;k=5;
%ans=19 23 29 31 37
n=0;
while n m=m+1;
if rem(m,2:m-1)~=0
m
n=n+1;
end
end

%使数组左下半三角元素中的值全部置成0
function y=fun10(a)
clc
%a=[1 9 7;2 3 8;4 5 6];
d=size(a);
for k=1:d(1)
for g=1:d(2)
if k>=g
a(k,g)=0;
end
end
end
a

%函数的功能是删除字符串s中的所有空格
function y=fun11(s)
clc
%s='asd af aa z67';
%ans=asdafaaz67
for k=length(s):-1:1
if s(k)==' '
s(k)=[];
end
end
s

%对长度为7个字符的字符串s，
%除首、尾字符外，将其余5个字符按ASCII码降序排列
function y=fun12(s)
clc
%s='CEAedca';
%ans =CedcEAa
sz=abs(s);
szz=sz(2:end-1);
for k=1:length(szz)
[a,b]=max(szz);
sz(k+1)=szz(b);
szz(b)=[];
end
s=setstr(sz)

%将s所指字符串中下标为偶数的字符删除
%串中剩余字符形成的新串放在y中
function y=fun13(n)
clc
%s='ABCDEFGHIJK'
%ans =ACEGIK
for k=length(s):-1:1
if rem(k,2)==0
s(k)=[];
end
end
y=s;

%比较两个字符串的长度，函数返回较长的字符串
%若两个字符串长度相同，则返回第一个字符串。
function y=fun14(a,b)
clc
%a='beijing';b='shanghai';
%ans=shanghai
if length(a)>=length(b)
y=a;
else
y=b;
end

%删除字符串s中所有的*号
function y=fun15(s)
clc
%s='****A*BC*DEF*G*******';
%ans=ABCDEFG
for k=length(s):-1:1
if s(k)=='*'
s(k)=[];
end
end
s

%按分数的高低排列学生的记录，高分在前。
function y=fun16(s)
clc
s=['GA005',85;'GA003',76;'GA002',69;'GA004',85;'GA001',96;'GA007',72;'GA008',64;'GA006',87; 'GA015',85;'GA013',94;'GA012',64;'GA014',91;'GA011',90;'GA017',64;'GA018',64;'GA016',72];
sl=s(:,6);
y1_5=s(:,1:5);
d=size(s);
y6=sl;
c=sl;
for k=1:d(1)
[a,b]=max(sl);
y1_5(k,1:5)=s(b,1:5);
c(k)=a;
y6=int2str(c);
sl(b)=0;
end
y=strcat(y1_5,32,32,32,y6);

%求出能整除x且不是偶数的各整数,并按从小到大的顺序放在y中。
function y=fun17(x)
clc
x=35;
a=find(rem(x,1:x)==0)

%将低于平均分的人数n和低于平均分的分数放在y中返回。
function [y,n]=fun18(s)
clc
s=[10 20 30 40 50 60 70 80 90];
p=0;
for k=1:length(s)
p=p+s(k);
end
p=p/length(s);
[a,b]=find(sn=length(a)
y=s(b);

%计算并输出3到n之间(含3和n)所有素数的平方根之和。
function y=fun19(n)
clc
n=100;
y=0;
for k=3:n
if rem(k,2:k-1)~=0
y=y+sqrt(k);
end
end

%求出该学生的平均分放在y中
function y=fun20(s)
clc
s={'GA005',[85.5,76,69.5,85,91,72,64.5,87.5]};
y=0;
t=0;
for k=length(s{2}):-1:1
y=y+s{2}(k);
t=t+1;
if t==8
y=y/8;
break
end
end

1
%只删除字符串前导和尾部的*号，串中字母之间的*号都不删除
function y=fun1(s)
clc
% s='****A*BC*DEF*G*******'
% ans= A*BC*DEF*G
n=length(s) % 形参n给出了字符串的长度
h=0; % 形参h给出了字符串中前导*号的个数
e=0; % 形参e给出了字符串中最后*号的个数
for k=n:-1:1
if s(k)=='*'
e=e+1;
s(k)=[];
else
break
end
end
n1=length(s);
for k=1:n1
if s(k)~='*'
t=k;
break
end
h=h+1;
end
h
e
s(1:k-1)=[]

2
%把指定分数范围内的学生数据放在b所指的数组中
%分数范围内的学生人数由函数值返回
function y=fun2(n1,n2)
clc
% n1=60;n2=69;
% b=GA002 GA008 GA012 GA017 GA018
% y=5
s=['GA005',85;'GA003',76;'GA002',69;'GA004',85;'GA001',96;'GA007',72;'GA008',64;'GA006',87; 'GA015',85;'GA013',94;'GA012',64;'GA014',91;'GA011',90;'GA017',64;'GA018',64;'GA016',72];
d=size(s);
for k=d(1):-1:1
if abs(s(k,6))n2
s(k,:)=[];
end
end
b=s;
b(:,6)=[]
d=size(b);
y=d(1);

3
%将两个两位数的正整数a、b合并形成一个整数放在c中
%合并的方式是: 将a数的十位和个位数依次放在c数的个位和百位上
% b数的十位和个位数依次放在c数的十位和千位上。
function y=fun3(a,b)
clc
%a=45;b=12;
%c = 2514
qw=rem(b,10)*1000;
bw=rem(a,10)*100;
sw=b-rem(b,10);
gw=fix(a/10);
c=qw+bw+sw+gw
4
%将两个两位数的正整数a、b合并形成一个整数放在c中
%合并的方式是: 将a数的十位和个位数依次放在c数的个位和百位上
% b数的十位和个位数依次放在c数的十位和千位上。
function y=fun3(a,b)
clc
%a=45;b=12;
%c = 2514
qw=rem(b,10)*1000;
bw=rem(a,10)*100;
sw=b-rem(b,10);
gw=fix(a/10);
c=qw+bw+sw+gw
5
%除了尾部的*之外，将字符串中其他*号全部删除。
function y=fun5(s)
clc
s='****A*BC*DEF*G*******';
%ans=ABCDEFG*******
for k=length(s):-1:1
if s(k)~='*'
t=k;
break
end
end
for g=t:-1:1
if s(g)=='*'
s(g)=[];
end
end
s
6
%其功能是计算并输出当x<0.97时下列多项式的值，直到abs(s(n)-s(n-1))<0.000001为止
%sn=1+0.5*x+0.5*(0.5-1)*x^2/1/2+0.5*(0.5-1)*(0.5-2)*x^3/1/2/3+…
% +0.5*（0.5-1）…（0.5-n+1）*x^n/(1*2*3*…*n)
function y=fun6(n)
clc
n=100;
%s =1.1000
x=0.21;
s=1;
for k=1:n
jc=1;
fz=1;
for g=1:k
jc=jc*g;
fz=fz*(0.5-g+1);
end
s=s+fz/jc*x^k;
end
s

7
%把高于等于平均分k的学生数据放在y中，输出y和k。
function [y,k]=fun7(s)
clc
s=['GA005',85;'GA003',76;'GA002',69;'GA004',85;'GA001',96;'GA007',72;'GA008',64;'GA006',87; 'GA015',85;'GA013',94;'GA012',64;'GA014',91;'GA011',90;'GA017',64;'GA018',64;'GA016',72];
d=size(s);
k=0;
for g=1:d(1)
k=k+abs(s(g,6));
end
k=k/d(1);
for gg=d(1):-1:1
if abs(s(gg,6)) s(gg,:)=[];
end
end
y1=s(:,1:5);
fs=int2str(s(:,6));
y=strcat(y1,32,32,32,fs);

8
%将两个两位数的正整数a、b合并形成一个整数放在c中
%合并的方式是：将a数的十位和个位数依次放在c数的十位和千位上
% b数的十位和个位数依次放在c数的百位和个位上。
function y=fun8(a,b)
clc
a=45;b=12;
%ans=5142
qw=mod(a,10)*1000;
bw=floor(b/10)*100;
sw=fix(a/10)*10;
gw=rem(b,10)*1;
c=qw+bw+sw+gw
9
%将字符串中的前面的*号全部移到字符串的尾部。
function y=fun9(s)
clc
s='*******A*BC*DEF*G****';
%ans=A*BC*DEF*G***********
for k=1:length(s)
if s(k)~='*'
t=k-1;
break
end
end
if t~=0
for g=1:length(s)-t
s(g)=s(g+7);
end
for gg=length(s)-t+1:length(s)
s(gg)='*';
end
end
s
10
%使矩阵右上半三角元素中的值乘以m
function y=fun10(a,m)
clc
a=[1 9 7;2 3 8;4 5 6];
m=2;
d=size(a);
for k=1:d(1)
for g=1:d(2)
if k<=g
a(k,g)=a(k,g)*m;
end
end
end
a
11
%找出数组元素x中最大的值和它所在的下标
% 最大的值和它所在的下标通过形参传回。
function [y,k]=fun11(x)
clc
x=[1 22 3 6 554 4 752 8 9 10 11 52 78 96 45 12 45 36 79 45 610 45 500 16 14];
y=x(1);
k=1;
for g=2:length(x)
if x(g)>y
y=x(g);
k=g;
end
end
k
12
%把低于平均分的学生数据放在y中
%低于平均分的学生人数k，返回y和k。
function [y,k]=fun12(s)
clc
s=['GA005',85;'GA003',76;'GA002',69;'GA004',85;'GA001',96;'GA007',72;'GA008',64;'GA006',87; 'GA015',85;'GA013',94;'GA012',64;'GA014',91;'GA011',90;'GA017',64;'GA018',64;'GA016',72];
d=size(s);
k=0;
for g=1:d(1)
k=k+abs(s(g,6));
end
k=k/d(1);
for gg=d(1):-1:1
if abs(s(gg,6))>=k
s(gg,:)=[];
end
end
y1=s(:,1:5);
fs=int2str(s(:,6));
y=strcat(y1,32,32,32,fs);
13
%使字符串的前导*号不得多于n个;若多于n个，则删除多余的*号
%若少于或等于n个，则什么也不做
%字符串中间和尾部的*号不删除。
function y=fun13(s,n)
clc
n=4;
%n=8;
s='*******A*BC*DEF*G****';
for k=1:length(s)
if s(k)~='*'
t=k-1;
break
end
end
if n s(1:t-n)=[];
end
s
14
%除了字符串前导的*号之外,将串中其它*号全部删除。
function y=fun14(s)
clc
s='****A*BC*DEF*G*******';
for k=1:length(s)
if s(k)~='*'
t=k;
break
end
end
for g=length(s):-1:t
if s(g)=='*'
s(g)=[];
end
end
s

15
%将所有大于1小于整数m的非素数输出
function y=fun15(n)
clc
n=17;
N=zeros(1,n);
for k=4:n
for g=2:k-1
if rem(k,g)==0;
N(k)=k;
break
end
end
end
y=find(N>0);

16
%将所有大于1小于整数m的非素数输出
function y=fun15(n)
clc
n=17;
N=zeros(1,n);
for k=4:n
for g=2:k-1
if rem(k,g)==0;
N(k)=k;
break
end
end
end
y=find(N>0);

17
%计算级数和s=∑x^n/n! (由0—n)
function y=fun17(n)
clc
n=10;
x=0.3;
y=1;
jc=1;
for k=1:n
for g=1:k
jc=jc*g;
end
y=y+x^k/jc;
jc=1;
end

18
%将两个两位数的正整数a、b合并形成一个整数放在c中
%合并的方式是：将a数的十位和个位数依次放在c数的个位和百位上
% b数的十位和个位数依次放在c数的千位和十位上。
function y=fun18(a,b)
clc
a=45;b=12;
qw=fix(b/10)*1000;
bw=rem(a,10)*100;
sw=rem(b,10)*10;
gw=floor(a/10);
c=qw+bw+sw+gw

19
%计算并输出下列级数和s=1/(1*2)+1/(2*3)+…+1/(n/(n+1))
%不就等于 1-1/(n+1) 嘛
function y=fun19(n)
clc
n=10;
%s=1-1/(n+1)
y=0;
for k=1:n
y=y+1/k/(k+1);
end

20
%将字符串s中所有下标为奇数位置上的字母转换为大写（若该位置上不是字母, 则不转换）。
function y=fun20(s)
clc
s='ab4cefg';
%ans =Ab4cEfG
for k=1:2:length(s)
if s(k)>='a' & s(k)<='z'
s(k)=char(s(k)-32);
end
end
s

%1. 函数fun1的功能是：将s所指字符串中ASCII值为偶数的字符删除，串中剩余字符形成一个新串放在t所指的数组中。
function y=fun1(s)
clc
%s='ABCDEFG12345';
%ans =ACEG135
sz=abs(s);
sc=find(mod(sz,2)==0);
sz(sc)=[];
t=setstr(sz);
y=t;
%2. 假定输入的字符串中只含有字母和*号。请编写函数fun2，它的功能是：将字符串尾部的*号全部删除，前面和中间的*号不删除。
function y=fun2(s)
clc
%s='****A*BC*DEF*G*******'
%ans =****A*BC*DEF*G
sz=abs(s);
for k=0:length(sz)
if sz(end)==abs('*')
sz(end)=[];
else
break
end
end
y=setstr(sz);
%3. 请编写函数fun3, 函数的功能是:统计一行字符串中单词的个数,作为函数值返回。 规定所有单词由小写字母组成, 单词之间由若干个空格隔开。
function y=fun3(s)
clc
%s='good me I hi hello';
%ans =5
sz=abs(s);
m=1;
for k=1:length(sz)
if sz(k)==' '
m=m+1;
end
end
y=m;
%4. 编写函数fun4,函数的功能是求出小于或等于n的所有素数并放在y中, 函数返回所求出的素数y和其个数k。
function [y,k]=fun4(n)
clc
N=1:n;
for g=n:-1:3
for k=g-1:-1:2
if rem(g,k)==0
N(g)=[];
break
end
end
end
N(1)=[];
k=length(N)
y=N;
%5. 编写函数fun5，它的功能是计算： y=(ln(1)+ln(2)+……+ln(m))^0.5
function y=fun5(m)
clc
y=0;
for k=1:m
y=y+log(k);
end
y=y^0.5;
%函数的功能是: 删去一维数组中所有相同的数, 使之只剩一个。数组中的数已按由小到大的顺序排列,函数返回删除后数组中数据的个数。
function y=fun6(s)
clc
%s=[2 2 2 3 4 4 5 6 6 6 6 7 7 8 9 9 10 10 10];
f=s(end);
for k=length(s)-1:-1:1
if s(k)==f
f=s(k);
s(k+1)=[];
else
f=s(k);
end
end
y=length(s);
%7 统计各年龄段的人数
function y=fun7(s)
clc;
s=round(rand(1,100)*120);
d=zeros(1,11);
for k=1:length(s)
switch floor(s(k)/10)
case 0,
d(11)=d(11)+1;
case 1
d(1)=d(1)+1;
case 2
d(2)=d(2)+1;
case 3
d(3)=d(3)+1;
case 4
d(4)=d(4)+1;
case 5
d(5)=d(5)+1;
case 6
d(6)=d(6)+1;
case 7
d(7)=d(7)+1;
case 8
d(8)=d(8)+1;
case 9
d(9)=d(9)+1;
otherwise
d(10)=d(10)+1;
end
end
disp(['d[0]=',num2str(d(11))])
r='d[ ]';
for g=1:9
m=num2str(g);
r(3)=m;
disp([r,num2str(d(g))])
end
disp(['d[10]=',num2str(d(10))])
8.将s所指字符串中所有下标为奇数位置上的字母转换为大写（该位置上不是字母，则不转换）。
function y=fun8(s)
clc
s='abc4EFg'
for k=1:2:length(s)
if s(k)>='a' & s(k)<='z'
s(k)=char(s(k)-32);
end
end
%9.求出矩阵所有元素的平均值并作为函数值在函数中输出。
function y=fun9(a)
clc
%a=[0 1 2 7 9;1 9 7 4 5;2 3 8 3 1;4 5 6 8 2;5 9 1 4 1];
%ans = 4.2800
d=size(a);
y=0;
for k=1:d(1)
for g=1:d(2)
y=y+a(k,g);
end
end
y=y/(d(1)*d(2));
%10.将M行N列的字符矩阵中的字符数据, 按行的顺序依次放到一个字符串中，然后将不是字符的删除。
function y=fun10(s)
clc;
%syms a b c d e f g h o;
%s=[8 a b;c 5 d;3 e 4;f 0 g;5 h o];
%ans =abcdefgho
d=size(s);
s1=reshape(s',1,d(1)*d(2));
for k=length(s1):-1:1
for g=0:9
if s1(k)==g
s1(k)=[];
end
end
end
t=s1';
y=t(1);
for k=2:length(t)
y=strcat(y,t(k));
end
%移动字符串中的内容, 移动的规则如下: 把第1个字符移到第m个字符, 把第2到第m+1个字符
%若第k到字符已移到字符串的最后,则第k+1到最后的字符移到字符串的前部。
function y=fun11(s)
clc
%s='ABCDEFGHIJK'
m=input('请输入m的值：')
ss=s;
if m>=1
for k=1:m-1
s(m-k)=ss(end-k+1);
end
for g=m:length(ss)
s(g)=ss(g-m+1);
end
s
else
disp('m值有误！')
end
%12.计算并输出给定10个数的方差 y=(0.1∑（xk-x'）^2)^0.5 x'= 0.1∑xk
function y=fun12(s)
clc
s=[95.0 89.0 76.0 65.0 88.0 72.0 85.0 81.0 90.0 56.0];
x_=0;
y=0;
for k=1:10
x_=x_+s(k);
end
x_=0.1*x_;
for g=1:10
y=y+(s(k)-x_)^2;
end
y=(0.1*y)^0.5
%13.求出矩阵a所有元素之和，作为函数值返回。
function y=fun13(a)
clc
%a=[1 3 5 7 9;2 9 9 9 4;6 9 9 9 8;1 3 5 7 0];
%sum =115
d=size(a);
sum=0;
for k=1:d(1)
for g=1:d(2)
sum=sum+a(k,g);
end
end
sum
%14.利用以下所示的简单迭代方法求方程：cos(x)-x=0 的一个实根。
function y=fun14(a)
%迭代步骤如下：
%（1）取 x1 初值为a；
%（2）x0 = x1，把x1的值赋给 x0；
%（3）x1 = cos(x0)，求出一个新的x1；
%（4）若x0 - x1的绝对值小于 0.000001，执行步骤（5），否则执行步骤（2）；
%（5）所求x1就是方程 cos(x)-x=0 的一个实根，作为函数值返回。
%程序将输出结果Root=0.739085。
clc
x1=0;
x0=x1;
x1=cos(x0);
while abs(x0-x1)>0.000001
x0=x1;
x1=cos(x0);
end
y=x1;
%15.将两个两位数的正整数a、b合并形成一个整数放在c中。
function y=fun15(a,b)
%合并的方式是：将a数的十位和个位数依次放在c数的千位和十位上， b数的十位和个位数依次放在c数的百位和个位上。
clc
%a=45;b=12;
%ans=4152
qw=floor(a/10)*1000;
bw=floor(b/10)*100;
sw=rem(a,10)*10;
gw=rem(b,10);
y=qw+bw+sw+gw;
%16.统计在s字符串中 'a' 到 'z' 26个字母各自出现的次数，并依次放在pp所指数组中。
function y=fun16(s)
clc
%s='abcdefgabcdeabc';
%pp=3 3 3 2 2 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
p=abs('a':'z');
pp=zeros(1,26);
for k=1:length(s)
for g=1:16
if abs(s(k))==p(g)
pp(g)=pp(g)+1;
end
end
end
pp
%17.使矩阵左下三角元素中的值乘以n 。
function y=fun17(a)
clc
%a=[1 9 7;2 3 8;4 5 6]
%ans =
% 3 9 7
% 6 9 8
% 12 15 18
d=size(a);
for k=1:d(1)
for g=1:d(2)
if k>=g
a(k,g)=a(k,g)*3;
end
end
end
a
%18.计算并输出n(包括n)以内能被5或9整除的所有自然数的倒数之和。
function y=fun18(n)
clc
n=20;
sum=0;
for k=1:n
if rem(k,9)==0
sum=sum+1/k;
else
switch rem(k,5)
case 0
sum=sum+1/k;
otherwise
sum=sum+0;
end

end
end
sum
%19.判断字符串是否为回文
function y=fun19(s)
clc
%s='level';
s='123312';
for k=1:length(s)
if s(k)~=s(end+1-k)
disp('NO')
break
end
if k==length(s)
disp('YES');
end
end
%20.求出 1 到 1000 之内能被 7 或11整除、但不能同时被 7 和 11 整除的所有整数并将它们放在 a所指的数组中，
%通过 n返回这些数的个数。
function y=fun20(s)
clc
a=[1:1000];
for k=1000:-1:1
if rem(k,7)~=0 & rem(k,11)~=0 %不能被7或11整除
a(k)=[];
end
if rem(k,7)==0 & rem(k,11)==0 %同时被7和11整除
a(k)=[];
end
end
a
n=length(a)