数学实验答案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)clcy=0;%s=[12 34 4 23 34 45 18 3 11]; %ans =35.9510for 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.9091n=1:n;f=find(rem(n,3)==0 & rem(n,7)==0);y=0;for k=1:length(f) y=y+f(k);endy=sqrt(y);%根据以下公式求p的值function y=fun4(m,n)clc%m=12;n=8; %ans=495fz=1;fm1=1;**=1;for k=1:m fz=fz*k; if k<=n fm1=fm1*k; end if k<=m-n **=***k; endendy=fz/fm1/**;%将字符串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 endends(1:t)=[]%求出向量s的最大元素在数组中的下标并存放在k所指的存储单元中。function y=fun6(s)clc%s=[876 675 896 101 301 401 980 431 451 777]; %k=7 m=980m=s(1);for g=1:length(s) if s(g)>m m=s(g); k=g; endendkm%将两个两位数的正整数a、b合并形成一个整数放在c中。%合并的方式是：将a数的十位和个位数依次放在c数的百位和个位上% b数的十位和个位数依次放在c数的千位和十位上。function y=fun7(a,b)clc%a=45;b=12; %ans=1425qw=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=763s=0;for k=2:n-1 if rem(n,k)==0 s=s+k; endends%将大于形参 m 且紧靠 m 的 k 个素数存入x所指的数组中。function y=fun9(m,k)clc%m=17;k=5; %ans=19 23 29 31 37n=0;while n<k m=m+1; if rem(m,2:m-1)~=0 m n=n+1; endend%使数组左下半三角元素中的值全部置成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 endenda%函数的功能是删除字符串s中的所有空格function y=fun11(s)clc%s='asd af aa z67'; %ans=as**az67for k=length(s):-1:1 if s(k)==' ' s(k)=[]; endends%对长度为7个字符的字符串s，%除首、尾字符外，将其余5个字符按ascii码降序排列function y=fun12(s)clc%s='ceaedca'; %ans =cedceaasz=abs(s);szz=sz(2:end-1);for k=1:length(szz) [a,b]=max(szz); sz(k+1)=szz(b); szz(b)=[];ends=setstr(sz)%将s所指字符串中下标为偶数的字符删除%串中剩余字符形成的新串放在y中function y=fun13(n)clc%s='abcdefghijk' %ans =acegikfor k=length(s):-1:1 if rem(k,2)==0 s(k)=[]; endendy=s;%比较两个字符串的长度，函数返回较长的字符串%若两个字符串长度相同，则返回第一个字符串。function y=fun14(a,b)clc%a='beijing';b='shanghai'; %ans=shanghaiif length(a)>=length(b) y=a;else y=b;end%删除字符串s中所有的*号function y=fun15(s)clc%s='****a*bc*def*g*******'; %ans=abcdefgfor k=length(s):-1:1 if s(k)=='*' s(k)=[]; endends%按分数的高低排列学生的记录，高分在前。function y=fun16(s)clcs=['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;endy=strcat(y1_5,32,32,32,y6);%求出能整除x且不是偶数的各整数,并按从小到大的顺序放在y中。function y=fun17(x)clcx=35;a=find(rem(x,1:x)==0)%将低于平均分的人数n和低于平均分的分数放在y中返回。function [y,n]=fun18(s)clcs=[10 20 30 40 50 60 70 80 90];p=0;for k=1:length(s) p=p+s(k);endp=p/length(s);[a,b]=find(s<p);n=length(a)y=s(b);%计算并输出3到n之间(含3和n)所有素数的平方根之和。function y=fun19(n)clcn=100;y=0;for k=3:n if rem(k,2:k-1)~=0 y=y+sqrt(k); endend%求出该学生的平均分放在y中function y=fun20(s)clcs={'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 endend 1%只删除字符串前导和尾部的*号，串中字母之间的*号都不删除function y=fun1(s)clc% s='****a*bc*def*g*******' % ans= a*bc*def*gn=length(s) % 形参n给出了字符串的长度h=0; % 形参h给出了字符串中前导*号的个数e=0; % 形参e给出了字符串中最后*号的个数for k=n:-1:1 if s(k)=='*' e=e+1; s(k)=[]; else break endendn1=length(s);for k=1:n1 if s(k)~='*' t=k; break end h=h+1;endhes(1:k-1)=[] 2 %把指定分数范围内的**放在b所指的数组中%分数范围内的学生人数由函数值返回function y=fun2(n1,n2)clc% n1=60;n2=69; % b=ga002 ga008 ga012 ga017 ga018 % y=5s=['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))<n1 | abs(s(k,6))>n2 s(k,:)=[]; endendb=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 = 2514qw=rem(b,10)*1000;bw=rem(a,10)*100;sw=b-rem(b,10);gw=fix(a/10);c=qw+bw+sw+gw4%将两个两位数的正整数a、b合并形成一个整数放在c中%合并的方式是: 将a数的十位和个位数依次放在c数的个位和百位上% b数的十位和个位数依次放在c数的十位和千位上。function y=fun3(a,b)clc%a=45;b=12; %c = 2514qw=rem(b,10)*1000;bw=rem(a,10)*100;sw=b-rem(b,10);gw=fix(a/10);c=qw+bw+sw+gw5%除了尾部的*之外，将字符串中其他*号全部删除。function y=fun5(s)clcs='****a*bc*def*g*******'; %ans=abcdefg*******for k=length(s):-1:1 if s(k)~='*' t=k; break endendfor g=t:-1:1 if s(g)=='*' s(g)=[]; endends 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)clcn=100; %s =1.1000x=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;ends7%把高于等于平均分k的**放在y中，输出y和k。function [y,k]=fun7(s)clcs=['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));endk=k/d(1);for gg=d(1):-1:1 if abs(s(gg,6))<k s(gg,:)=[]; endendy1=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)clca=45;b=12; %ans=5142qw=mod(a,10)*1000;bw=floor(b/10)*100;sw=fix(a/10)*10;gw=rem(b,10)*1;c=qw+bw+sw+gw9%将字符串中的前面的*号全部移到字符串的尾部。function y=fun9(s)clcs='*******a*bc*def*g****'; %ans=a*bc*def*g***********for k=1:length(s) if s(k)~='*' t=k-1; break endendif t~=0 for g=1:length(s)-t s(g)=s(g+7); end for gg=length(s)-t+1:length(s) s(gg)='*'; endends10%使矩阵右上半三角元素中的值乘以mfunction y=fun10(a,m)clca=[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 endenda11%找出数组元素x中最大的值和它所在的下标% 最大的值和它所在的下标通过形参传回。function [y,k]=fun11(x)clcx=[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; endendk12%把低于平均分的**放在y中%低于平均分的学生人数k，返回y和k。function [y,k]=fun12(s)clcs=['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));endk=k/d(1);for gg=d(1):-1:1 if abs(s(gg,6))>=k s(gg,:)=[]; endendy1=s(:,1:5);fs=int2str(s(:,6));y=strcat(y1,32,32,32,fs);13%使字符串的前导*号不得多于n个;若多于n个，则删除多余的*号%若少于或等于n个，则什么也不做%字符串中间和尾部的*号不删除。function y=fun13(s,n)clcn=4;%n=8;s='*******a*bc*def*g****';for k=1:length(s) if s(k)~='*' t=k-1; break endendif n<t s(1:t-n)=[];ends14%除了字符串前导的*号之外,将串中其它*号全部删除。function y=fun14(s)clcs='****a*bc*def*g*******';for k=1:length(s) if s(k)~='*' t=k; break endendfor g=length(s):-1:t if s(g)=='*' s(g)=[]; endends15%将所有大于1小于整数m的非素数输出function y=fun15(n)clcn=17;n=zeros(1,n);for k=4:n for g=2:k-1 if rem(k,g)==0; n(k)=k; break end endendy=find(n>0); 16 %将所有大于1小于整数m的非素数输出function y=fun15(n)clcn=17;n=zeros(1,n);for k=4:n for g=2:k-1 if rem(k,g)==0; n(k)=k; break end endendy=find(n>0); 17 %计算级数和s=∑x^n/n! (由0—n)function y=fun17(n)clcn=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;end18%将两个两位数的正整数a、b合并形成一个整数放在c中%合并的方式是：将a数的十位和个位数依次放在c数的个位和百位上% b数的十位和个位数依次放在c数的千位和十位上。function y=fun18(a,b)clca=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+gw19 %计算并输出下列级数和s=1/(1*2)+1/(2*3)+…+1/(n/(n+1)) %不就等于 1-1/(n+1) 嘛 function y=fun19(n)clcn=10;%s=1-1/(n+1)y=0;for k=1:n y=y+1/k/(k+1);end20 %将字符串s中所有下标为奇数位置上的字母转换为大写（若该位置上不是字母, 则不转换）。function y=fun20(s)clcs='ab4cefg'; %ans =ab4cefgfor k=1:2:length(s) if s(k)>='a' & s(k)<='z' s(k)=char(s(k)-32); endends %1.函数fun1的功能是：将s所指字符串中ascii值为偶数的字符删除，串中剩余字符形成一个新串放在t所指的数组中。function y=fun1(s)clc%s='abcdefg12345';%ans =aceg135sz=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*gsz=abs(s);for k=0:length(sz) if sz(end)==abs('*') sz(end)=[]; else break endendy=setstr(sz);%3.请编写函数fun3, 函数的功能是:统计一行字符串中单词的个数,作为函数值返回。 规定所有单词由小写字母组成, 单词之间由若干个空格隔开。 function y=fun3(s) clc %s='good me i hi hello'; %ans =5sz=abs(s);m=1;for k=1:length(sz) if sz(k)==' ' m=m+1; endendy=m;%4.编写函数fun4,函数的功能是求出小于或等于n的所有素数并放在y中, 函数返回所求出的素数y和其个数k。function [y,k]=fun4(n)clcn=1:n;for g=n:-1:3 for k=g-1:-1:2 if rem(g,k)==0 n(g)=[]; break end endendn(1)=[];k=length(n)y=n;%5.编写函数fun5，它的功能是计算： y=(ln(1)+ln(2)+……+ln(m))^0.5function y=fun5(m)clcy=0;for k=1:m y=y+log(k);endy=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 endy=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; otherw**e d(10)=d(10)+1; endendd**p(['d[0]=',num2str(d(11))])r='d[ ]';for g=1:9 m=num2str(g); r(3)=m; d**p([r,num2str(d(g))])endd**p(['d[10]=',num2str(d(10))]) 8.将s所指字符串中所有下标为奇数位置上的字母转换为大写（该位置上不是字母，则不转换）。function y=fun8(s)clcs='abc4efg'for k=1:2:length(s) if s(k)>='a' & s(k)<='z' s(k)=char(s(k)-32); endend%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.2800d=size(a);y=0;for k=1:d(1) for g=1:d(2) y=y+a(k,g); endendy=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 =abcdefghod=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 endendt=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 selse d**p('m值有误！')end%12.计算并输出给定10个数的方差 y=(0.1∑（xk-x'）^2)^0.5 x'= 0.1∑xkfunction y=fun12(s)clcs=[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);endx_=0.1*x_;for g=1:10 y=y+(s(k)-x_)^2;endy=(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 =115d=size(a);sum=0;for k=1:d(1) for g=1:d(2) sum=sum+a(k,g); endendsum %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。clcx1=0;x0=x1;x1=cos(x0);while abs(x0-x1)>0.000001 x0=x1; x1=cos(x0);endy=x1;%15.将两个两位数的正整数a、b合并形成一个整数放在c中。function y=fun15(a,b)%合并的方式是：将a数的十位和个位数依次放在c数的千位和十位上， b数的十位和个位数依次放在c数的百位和个位上。clc%a=45;b=12; %ans=4152qw=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 0p=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 endendpp%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 18d=size(a);for k=1:d(1) for g=1:d(2) if k>=g a(k,g)=a(k,g)*3; end endenda%18.计算并输出n(包括n)以内能被5或9整除的所有自然数的倒数之和。function y=fun18(n)clcn=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; otherw**e sum=sum+0; end endendsum %19.判断字符串是否为回文function y=fun19(s)clc%s='level';s='123312';for k=1:length(s) if s(k)~=s(end+1-k) d**p('no') break end if k==length(s) d**p('yes'); endend %20.求出 1 到 1000 之内能被 7 或11整除、但不能同时被 7 和 11 整除的所有整数并将它们放在 a所指的数组中，%通过 n返回这些数的个数。function y=fun20(s)clca=[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 endan=length(a) 20210311