Matlab上机题库和详细答案

Matlab上机题库及详细答案例1、用一个简单命令求解线性系统3x1+x2-x3=3.6x1+2x2+4x3=2.1-x1+4x2+5x3=-1.4解：A=[31-1;124;-145];b=[3.6;2.1;-1.4];x=A\b%Inv（A）*bx=1.4818-0.46060.3848例2、用简短命令计算并绘制在0x6范围内的sin(2x)、sinx2、sin2x。解：x=linspace(0,6)y1=sin(2*x),y2=sin(x.^2),y3=(sin(x)).^2;plot(x,y1,x,y2,x,y3)例3：画出指数衰减曲线y1=exp(-t/3)*sin(3*t)和它的包络y2=exp(-t/3),t的取值范围是（0，4pi）。解：t=0:pi/100:4*pi;y2=exp(-t/3);y1=y2.*sin(3*t);plot(t,y1,'-r',t,y2,':b',t,-y2,':b')知识点：1.A(:,j)表示取A矩阵的第j列全部元素；A(i,:)表示A矩阵第i行的全部元素；A(i,j)表示取A矩阵第i行、第j列的元素。2.A(i:i+m,:)表示取A矩阵第i～i+m行的全部元素；A(:,k:k+m)表示取A矩阵第k～k+m列的全部元素，A(i:i+m,k:k+m)表示取A矩阵第i～i+m行内，并在第k～k+m列中的所有元素。A=[1,2,3,4,5;11,12,13,14,15;21,22,23,24,25;31,32,33,34,35];A(2:3,4:5)A(2:3,1:2:5)3.利用end运算符，end表示某一维的末尾元素下标。A=[1,2,3,4,5;11,12,13,14,15;21,22,23,24,25;31,32,33,34,35];A(end,:)%取A最后一行元素A([1,4],3:end)%取A第1、4行中第3列到最后一列元素。4.A(:)将矩阵A的每一列元素堆叠起来，成为一个行向量。例4：使用关系运算和元素的逻辑运算找出大于60小于100的数的位置。解：num=round(rand(1,10)*100);%生成<100的整数，%round(A),A中的元素四舍五入到整数，rand(1,10),产生随机数0.0-1.0的1X10矩阵。n=(num>60)&(num<100)n=n.*numResult=find(n)%查找非0位置，按索引方式as=n(result)%显示这些数字例5：产生对角线上全1，其余为0的2行3列矩阵的命令是_C。A.ones(2,3)B.ones(3,2)C.eye(2,3)D.eye(3,2)例6：已知a=0:4,b=1:5,下面的运算表达式出错的为_D__。A.a+bB.a.*bC.a’*bD.a*b例7：有一个2行2列的元胞数组c，则c(2)指的是__D。A.第1行第2列的元素内容B.第2行第1列的元素内容C.第1行第2列的元素D.第2行第1列的元素例8：已知x=0:10，则x有B个元素。A.10B.11C.9D.12例9：已知数组a=则a(:,end)是指CA.所有元素B.第一行元C.第三列元素D.第三行元素例10：计算x从0到20，y=sin(x)中，π<x<4π范围中，y>0的所有值。解：x=0:20;x1=x.*((x>pi)&(x<4*pi));y=sin(x);y1=sin(x1).*(sin(x1)>0)例11:设A=[123;456]为2×3矩阵,分别建立3×3、3×2和与矩阵A同样大小的零矩阵。解：zeros(3)；zeros(3,2)；A=[123;456]；zeros(size(A))例12：建立随机矩阵：(1)在区间[20,50]内均匀分布的5阶随机矩阵。(2)均值为0.6方差为0.1的5阶正态分布随机矩阵解：>>x=20+(50-20)*rand(5)；>>y=0.6+sqrt(0.1)*randn(5)例13：将101~125等25个数填入一个5行5列的表格中，使其每行每列及对角线的和均为565。解：M=100+magic(5)M=117124101108115123105107114116104106113120122110112119121103111118125102109例14:先建立5×5矩阵A，然后将A的第一行元素乘以1，第二行乘以2，…，第五行乘以5。解A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;11,18,25,2,19];D=diag(1:5);D*A%用D左乘A，对A的每行乘以一个指定常数例15：输入x,y的值，并将它们的值互换后输出。程序如下：x=input('Inputxplease.')；y=input('Inputyplease.');z=x;x=y;y=z;disp(x);disp(y);例16：:求一元二次方程ax2+bx+c=0的根。程序如下：a=input('a=?')；b=input('b=?');c=input('c=?');d=b*b-4*a*c;x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)];disp(['x1=',num2str(x(1)),',x2=',num2str(x(2))]);例17：输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出。解：c=input('请输入一个字符','s');ifc>='A'&c<='Z'disp(setstr(abs(c)+abs('a')-abs('A')));elseifc>='a'&c<='z'disp(setstr(abs(c)-abs('a')+abs('A')));elseifc>='0'&c<='9'disp(abs(c)-abs('0'));elsedisp(c);end例18：某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)：price<200没有折扣200≤price<5003%折扣500≤price<10005%折扣1000≤price<25008%折扣2500≤price<500010%折扣5000≤price14%折扣输入所售商品的价格，求其实际销售价格。解：price=input('请输入商品价格');switchfix(price/100)case{0,1}%价格小于200rate=0;case{2,3,4}%价格大于等于200但小于500rate=3/100;casenum2cell(5:9)%价格大于等于500但小于1000rate=5/100;casenum2cell(10:24)%价格大于等于1000但小于2500rate=8/100;casenum2cell(25:49)%价格大于等于2500但小于5000rate=10/100;otherwise%价格大于等于5000rate=14/100;endprice=price*(1-rate)%输出商品实际销售价格例19已知，当n=100时，求的值。解：程序如下：y=0;n=100;fori=1:n；y=y+1/(2*i-1);End例20：一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数解：form=100:999m1=fix(m/100);m2=rem(fix(m/10),10);m3=rem(m,10);ifm==m1*m1*m1+m2*m2*m2+m3*m3*m3；disp(m)；end例21：从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。程序如下：解：sum=0;cnt=0;val=input('Enteranumber(endin0):');while(val~=0)sum=sum+val;cnt=cnt+1;val=input('Enteranumber(endin0):');endif(cnt>0)Sum；mean=sum/cntend例22：求[100，200]之间第一个能被21整除的整数。程序如下：解：forn=100:200ifrem(n,21)~=0Continue;endbreakend例23:若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求[1,500]之间的全部完数。解：form=1:500；s=0;fork=1:m/2；ifrem(m,k)==0；s=s+k;end；end；ifm==sdisp(m);end；end例24：编写函数文件求半径为r的圆的面积和周长。函数文件如下：解：function[s,p]=fcircle(r)%CIRCLEcalculatetheareaandperimeterofacircleofradiir%r圆半径%s圆面积%p圆周长s=pi*r*r;p=2*pi*r;例25：利用函数文件，实现直角坐标(x,y)与极坐标(ρ,θ)之间的转换。函数文件tran.m：解：function[rho,theta]=tran(x,y)；rho=sqrt(x*x+y*y);theta=atan(y/x);调用tran.m的命令文件main1.m：在x=input('Pleaseinputx=:');y=input('Pleaseinputy=:');[rho,the]=tran(x,y);rho例26：在0≤x≤2区间内，绘制曲线y=2e-0.5xcos(4πx)解：程序如下：x=0:pi/100:2*pi;y=2*exp(-0.5*x).*cos(4*pi*x);plot(x,y)例27：用不同标度在同一坐标内绘制曲线y1=0.2e-0.5xcos(4πx)和y2=2e-0.5xcos(πx)。解：x=0:pi/100:2*pi;y1=0.2*exp(-0.5*x).*cos(4*pi*x);y2=2*exp(-0.5*x).*cos(pi*x);Plotyy(x,y1,x,y2);例28：在0≤x≤2区间内，用不同线型和颜色绘制曲线y=2e-0.5xsin(2πx)及其包络线。解：x=(0:pi/100:2*pi)';y1=2*exp(-0.5*x)*[1,-1];y2=2*exp(-0.5*x).*sin(2*pi*x)；x1=(0:12)/2；y3=2*exp(-0.5*x1).*sin(2*pi*x1)；plot(x,y1,'k:',x,y2,'b--',x1,y3,'rp')例29：在0≤x≤2区间内，绘制曲线y1=2e-0.5x和y2=cos(4πx)并给图形添加图形标注。解：x=0:pi/100:2*pi;y1=2*exp(-0.5*x);y2=cos(4*pi*x);plot(x,y1,x,y2)；title('xfrom0to2{\pi}');%加图形标题xlabel('VariableX');%加X轴说明ylabel('VariableY');%加Y轴说明text(0.8,1.5,'曲线y1=2e^{-0.5x}')%在指定位置添加图形说明text(2.5,1.1,'曲线y2=cos(4{\pi}x)');legend(‘y1’,‘y2’)%加图例例30：采用图形保持，在同一坐标内绘制曲线y1=0.2e-0.5xcos(4πx)和y2=2e-0.5xcos(πx)。解：x=0:pi/100:2*pi;y1=0.2*exp(-0.5*x).*cos(4*pi*x);plot(x,y1)；y2=2*exp(-0.5*x).*cos(pi*x);plot(x,y2);例31：分别以条形图、阶梯图、杆图和填充图形式绘制曲线y=2sin(x)。解：x=0:pi/10:2*pi;y=2*sin(x);subplot(2,2,1);bar(x,y,'g');title('bar(x,y,''g'')');axis([0,7,-2,2]);subplot(2,2,2);stairs(x,y,'b');title('stairs(x,y,''b'')');axis([0,7,-2,2]);subplot(2,2,3);stem(x,y,'k');title('stem(x,y,''k'')');axis([0,7,-2,2]);subplot(2,2,4);fill(x,y,'y');title('fill(x,y,''y'')');axis([0,7,-2,2]);例32：绘制r=sin(t)cos(t)的极坐标图，并标记数据点。解：t=0:pi/50:2*pi;r=sin(t).*cos(t);polar(t,r,'-*');例33：绘制y=10x2的对数坐标图，并与直角线性坐标图比较。解：x=0:0.1:10；y=10*x.*x;subplot(2,2,1)；plot(x,y)；title('plot(x,y)')gridon；subplot(2,2,2)；semilogx(x,y)；title('semilogx(x,y)')gridon；subplot(2,2,3)；semilogy(x,y)gridontitle('semilogy(x,y)')；subplot(2,2,4)loglog(x,y)；title('loglog(x,y)')gridon例34：绘制三维曲面图z=sin(x+sin(y))-x/10。解：[x,y]=meshgrid(0:0.25:4*pi);z=sin(x+sin(y))-x/10;mesh(x,y,z);axis([04*pi04*pi-2.51]);例35：绘制三维图形：(1)绘制魔方阵的三维条形图。(2)以三维杆图形式绘制曲线y=2sin(x)。(3)已知x=[2347,1827,2043,3025]，绘制饼图。(4)用随机的顶点坐标值画出五个黄色三角形。解：subplot(2,2,1);bar3(magic(4))；subplot(2,2,2);y=2*sin(0:pi/10:2*pi);stem3(y);subplot(2,2,3);pie3([2347,1827,2043,3025]);subplot(2,2,4);fill3(rand(3,5),rand(3,5),rand(3,5),'y')例36：从不同视点绘制多峰函数曲面。解：subplot(2,2,1);mesh(peaks);view(-37.5,30);%指定子图1的视点title('azimuth=-37.5,elevation=30')；subplot(2,2,2);mesh(peaks);view(0,90);%指定子图2的视点title('azimuth=0,elevation=90')；subplot(2,2,3);mesh(peaks);view(90,0);%指定子图3的视点title('azimuth=90,elevation=0')；subplot(2,2,4);mesh(peaks);view(-7,-10);%指定子图4的视点title('azimuth=-7,elevation=-10')例37：绘制了peaks函数曲面并且将它绕z轴旋转。解：[X,Y,Z]=peaks(30);surf(X,Y,Z)；axis([-3,3,-3,3,-10,10])；axisoff;shadinginterp;olormap(hot);m=moviein(20);%建立一个20列大矩阵fori=1:20view(-37.5+24*(i-1),30)%改变视点m(:,i)=getframe;%将图形保存到m矩阵endmovie(m,2);%播放画面2次例38：分别求3×4矩阵x中各列和各行元素中的最大值，并求整个矩阵的最大值和最小值。解：A=[13,-56,78;25,63,-235;78,25,563;10-1]；a=max(A,[],2)；b=min(A,[],2)；c=max(A)d=min(A)；e=max(max(A))；f=min(min(A))例39：求矩阵A的每行元素的乘积和全部元素的乘积。解：A=[1234;5678;9101112]；a=prod(A)；b=prod(a)例40：某实验对一根长10米的钢轨进行热源的温度传播测试。用x表示测量点0:2.5:10(米)，用h表示测量时间0:30:60(秒)，用T表示测试所得各点的温度(℃)。试用线性插值求出在一分钟内每隔20秒、钢轨每隔1米处的温度TI。解：x=0:2.5:10;h=[0:30:60]';T=[95,14,0,0,0;88,48,32,12,6;67,64,54,48,41];xi=[0:10];hi=[0:20:60]';TI=interp2(x,h,T,xi,hi)例41：用一个三次多项式逼近函数sinx解：x=linspace(0,2*pi,50)；y=sin(x)；p=polyfit(x,y,3)；X=linspace(0,2*pi,20)；Y=sin(X)Y1=polyval(p,X)；plot(X,Y,':O',X,Y1,'-*')例42：已知f(x)(1)计算f(x)=0的全部根。(2)由方程f(x)=0的根构造一个多项式g(x)，并与f(x)进行对比。解：P=[3,0,4,-5,-7.2,5];X=roots(P)%求方程f(x)=0的根G=poly(X)%求多项式g(x)例43：计算二重定积分解：(1)建立一个函数文件fxy.m：functionf=fxy(x,y)；globalki;ki=ki+1;%ki用于统计被积函数的调用次数f=exp(-x.^2/2).*sin(x.^2+y);(2)调用dblquad函数求解。globalki;ki=0;I=dblquad('fxy',-2,2,-1,1)；kiI=1.57449318974494ki=1038例44：用LU分解求解例7-1中的线性方程组。解：A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4];b=[13,-9,6,0]';[L,U]=lu(A);x=U\(L\b)或采用LU分解的第2种格式，命令如下：[L,U,P]=lu(A);x=U\(L\P*b)例45：用QR分解求解例7-1中的线性方程组。解：A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4];b=[13,-9,6,0]';[Q,R]=qr(A);x=R\(Q\b)或采用QR分解的第2种格式，命令如下：[Q,R,E]=qr(A);x=E*(R\(Q\b))例46求f(x)=x-10x+2=0在x0=0.5附近的根。解：步骤如下：(1)建立函数文件funx.m。functionfx=funx(x)fx=x-10.^x+2;(2)调用fzero函数求根。z=fzero('funx',0.5)z=0.3758例47：求下列非线性方程组在(0.5,0.5)附近的数值解。解：(1)建立函数文件myfun.m。functionq=myfun(p)x=p(1);y=p(2);q(1)=x-0.6*sin(x)-0.3*cos(y);q(2)=y-0.6*cos(x)+0.3*sin(y);(2)在给定的初值x0=0.5,y0=0.5下，调用fsolve函数求方程的根。x=fsolve('myfun',[0.5,0.5]',optimset('Display','off'))例48：求f(x)=x3-2x-5在[0,5]内的最小值点与最大值点。解：(1)建立函数文件mymin.m。functionfx=mymin(x）；fx=x.^3-2*x-5;(2)调用fmin函数求最小值点与最大值点。xmin=fmin('mymin',0,5)；xmin=0.8165；xmax=fmin(‘-mymin',0,5)例49：求解下面的有约束优化问题解：x0=[0.5;0.5]；A=[-1,-0.5;-0.5,-1]；b=[-0.4;-0.5]；lb=[0,0]；[x,f]=fmincon('fop',x0,A,b,[],[],lb,[],[],[])；functionf=fop(x)；f=0.4*x(2)+x(1)^2+x(2)^2-x(1)*x(2)+1/30*x(1)^3例50：完成下列操作：(P361,T4)(1)求[100,999]之间能被21整除的数的个数。(2)建立一个字符串向量，删除其中的大写字母。解：clc;clear；A=100:999;length(A(mod(A,21)==0))；ch='SDFHJI23423bfjdbg';Chk=find(ch<='Z'&ch>='A');ch(k)=[]例51：下面是一个线性方程组：(1)求方程的解。(2)将方程右边向量元素b3改为0.53再求解，并比较解的相对变化。(3)计算系数矩阵A的条件数并分析结论。解：clc;clear;A=1./[234;345;456];b1=[0.95;0.67;0.52];x1=inv(A)*b1%x2=A\b1%(2)b2=[0.95;0.67;0.53];x2=inv(A)*b2；x3=x2-x1；B1=norm(x1)%x1的2-范数B2=norm(x2)；B=B2/B1%(3)条件数：D=cond(A)%A的2-范数下的条件数，该数都大于1，但越接近与1越好。if(D>1e3)disp('结论：A的条件数远大于1，为不良矩阵。')例54：求分段函数的值。用if语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y值。解：clc;clear;x=[-5.0-3.01.02.02.53.05.0];if(x<3&x~=-3)y=x^2+x-6;elseif(x>=0&x<5&x~=2&x~=3)；y=x^2-5*x+6;Elsey=x.^2-x-1;end例55.输入一个百分制成绩，要求输出成绩等级A、B、C、D、E。其中90分~100分为A，80分~89分为B，79分~79分为C，60分~69分为D，60分以下为E。要求：(1)分别用if语句和switch语句实现。(2)输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。解：clc;clear;x=input('#if#pleaseinputyourscore:');%useifif(x>=90&x<=100)y='A'elseifx>=80&x<=90；y='B'elseifx>=70&x<=80y='C'elseifx>=60&x<=70y='D'elseifx>=0&x<=60y='E'elsedisp('ERROR!Pleaseiuputaintegernumber.')disp('')end%usesiwtchx=input('#switch#pleaseinputyourscore:');ifx>=0&x<=100switchfix(x/10);case{9,10};y='A'case8;y='B'case7;y='C';case6;y='D';case{5,4,3,2,1,0}y='E';endelsedisp('ERROR!Pleaseiuputaintegernumber.')end例56：硅谷公司员工的工资计算方法如下：(1)工作时数超过120小时者，超过部分加发15%。(2)工作时数低于60小时者，扣发700元。(3)其余按每小时84元计发。试编程按输入的工号和该号员工的工时数，计算应发工资。解：clc;clear;number=input('Pleaseinputyourworknumber:');time=input('Pleaseinoutyourworkhours:');iftime>120;y=120*84+(time-120)*84*1.15;elseiftime>=60&time<=120;y=time*84;elseiftime<60&time>0；y=time*84-700;Elsey='ERROR!'enddisp('')disp(['ThewagesoftheNo.',num2str(number),'employeeis',num2str(y)])例57：设，在同一图形窗口采用子图的形式绘制不同图形：条形图、阶梯图、杆图和全对数坐标图。解：clc;clear;t=-pi:0.25:pi;y=1./(1+exp(-t));figure;subplot(2,2,1);bar(t,y);subplot(2,2,2);stairs(t,y);subplot(2,2,3);stem(t,y);subplot(2,2,4);loglog(t,y);例58：数值与符号计算（1）求极限（2）求不定积分（3）已知线性方程组Ax=b，其中，运用稀疏存储矩阵的方式求其解解：clc;clear;symsabx;%(1)f=exp(x)./(a+b*exp(x));limit(f,x,inf,'left')%(2)symsalphax;f=x*exp(alpha*x);int(f)%(3)B=[-120;-12-1;-12-1;-12-1;02-1];d=[-1;0;1];A=spdiags(B,d,5,5)%稀疏存储矩阵的方式%A=full(A)全矩阵b=[1;0;0;0;0];x=(inv(A)*b)'%转置例58：绘图（1）绘制极坐标图：。（2）绘制曲面图：解：clc;clear;%(1)thrta=-pi/6:0.1:pi/6;rho=3.*sin(thrta).*cos(thrta)./(sin(thrta).^3+cos(thrta).^3);polar(thrta,rho)%(2)x=-3:0.1:3;y=-3:0.1:3;[x,y]=meshgrid(x);z=-5./(1+x.^2+y.^2);figure;surf(x,y,z)figure[x,y,z]=peaks(40);surf(x,y,z)例59：求非线性方程组解：functionF=myfun(X)；x=X(1);y=X(2);F(1)=x^2+y^2-9;F(2)=x+y-1;end命令：clc;clear;fsolve('myfun',[3,0],optimset('Display','off'))例60：已知某精密仪器的某部件轮廓线的数据如下：用三次样条插值法求x每改变0.1时的y值。解：clc;clear;x=[0345791112131415];y=[01.21.72.02.12.01.81.21.01.6];X=0:0.1:15;Y=interp1(x,y,X,'spline')；plot(X,Y)例61：设计一个图形用户界面，其中有一个坐标平面和两个按钮，当单击第一个按钮时，在坐标平面上绘制一副图形，当单击第二个按钮时，可以改变界面的背景颜色。解：figure('name','图形演示系统','numbertitle','off','menubar','none');hplot=uimenu(gcf,'Label','&Plot');uimenu(hplot,'Label','&huatu','Call',['x=0:pi/100:2*pi;','y=sin(x);','plot(x,y);']);hselet=uimenu(gcf,'label','&Selet');uimenu(hselet,'label','&Red','call','set(gcf,"color","r");');uimenu(hselet,'label','&Yellow','call','set(gcf,"color","y");');例62：分别用if语句和switch语句实现以下计算，其中,,abc和x的值从键盘输入。解：clc;clear;a=input('Pleaseinputa=');b=input('Pleaseinputb=');c=input('Pleaseinputc=');x=input('Pleaseinputx=');disp('#if#');ifx>=0.5&x<1.5；y=a*x^2+b*x+c；elseifx<3.5&x>=1.5；y=a*sin(b)^c+xelseifx>=3.5&x<5.5；y=log(abs(b+(c./x)))elsey='ERROR!'enddisp('#switch#')；switch(x*10/5)case{12}y=a*x^2+b*x+ccase{3456}y=a*sin(b)^c+xcase{78910}y=log(abs(b+(c./x)))otherwisey='ERROR!'end例63：求下列矩阵的主对角线元素、上三角阵、下三角阵、秩、范数、条件数和迹。解：clc;clear;A=[1-123;51-42;3052;111509]diag(A)triu(A)tril(A)rank(A)trace(A)norm(A)cond(A)B=[0.43432;-8.9421]diag(B)triu(B)tril(B)rank(B)trace(B)%方阵才有迹。norm(B)cond(B)例64：求微分方程组的数值解，并绘制解的曲线。解：函数文件functiondy=rigid(t,y)dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=-y(1)*y(3);dy(3)=-0.51*y(1)*y(2);end命令：clc;clear[T,Y]=ode45('rigid',[0,20],[0,1,1])plot(T,Y(:,1),'-',T,Y(:,2),'+',T,Y(:,3),'*');例65：已知完成下列运算：(1)B=P1·P2·A；(2)B的逆矩阵并验证结果；(3)包括B矩阵主对角线元素的下三角阵；(4)B的行列式值。解：clc;clear;p1=sym('[010;100;001]');p2=sym('[100;010;101]');A=sym('[abc;def;ghk]');B=p1*p2*A；B1=inv(B)；B2=inv(B1)；tril(B)；det(B)例66：用符号方法求下列极限或导数。（第三问！）解：clc;clear;%(1)symsx;f=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))./sin(x)^3;limit(f)%(2)symsxy;y=(1-cos(2*x))/x;diff(y,x)diff(y,x,2)%(3)symsxyf;f=(x^2-2*x)*exp(-x^2-y^2-x*y);Yx=-diff(f,x)/diff(f,y)dfxy=diff(diff(f,x),y);x=sym('0');y=sym('1');eval(dfxy)例67：(1)最大值、最小值、各数之和。(2)正数、零、负数的个数。解：clc;clear;forn=1:100ifn==1f(n)=1;elseifn==2f(n)=0;elseifn==3f(n)=1;elsef(n)=f(n-1)-2*f(n-2)+f(n-3);endendmax=max(f)min=min(f)sum=sum(f)length(find(f>0))length(find(f==0))length(find(f<0))例68：若两个连续自然数的乘积减1是素数，则称这两个边疆自然数是亲密数对，该素数是亲密素数。例如，2×3-1=5，由于5是素数，所以2和3是亲密数，5是亲密素数。求[2,50]区间内：(1)亲密数对的对数。(2)与上述亲密数对对应的所有亲密素数之和。解：clc;clear;n=0;s=0;fori=2:49b=i*(i+1)-1;m=fix(sqrt(b));forj=2:mifrem(b,j)==0breakendendifj==mn=n+1;s=s+b;end例69：一物理系统可用下列方程组来表示：从键盘输入m1、m2和θ的值，求a1、a2、N1和N2的值。其中g取9.8，输入θ时以角度为单位。要求：定义一个求解线性方程组AX=B的函数文件，然后在命令文件中调用该函数文件。解：函数文件：functionX=mat(m1,m2,thrta)g=9.8;t=thrta*pi/180;A=[m1*cos(t)-m1-sin(t)0m1*sin(t)0cos(t)00m2-sin(t)000-cos(t)1];b=[0;m1*g;0;m2*g];X=inv(A)*b;命令：clc;clear;m1=input('m1=');m2=input('m2=');thrta=input('thrta=');X=mat(m1,m2,thrta)’例70：利用MATLAB提供的rand函数生成30000个符合均匀分布的随机数，然后检验随机数的性质：(1)均值和标准方差。(2)最大元素和最小元素。(3)大于0.5的随机数个数占总数的百分比。解：clc;clear;A=rand(300,100);mean(mean(A))std(std(A))max(max(A))min(min(A))length(find(A>0.5))/30000例71：将100个学生5门功课的成绩存入矩阵P中，进行如下处理：(1)分别求每门课的最高分、最低分及相应学生序号。(2)分别求每门课的平均分和标准方差。(3)5门课总分的最高分、最低分及相应学生序号。(4)将5门课总分按从大到小顺序存入zcj中，相应学生序号存入xsxh。提示：为避免输入学生成绩的麻烦，可用取值范围在[45,95]之间的随机矩阵来表示学生成绩。解：clc;clear;p=45+(95-45)*rand(100,5);[x,l]=max(p)[y,k]=min(p)mean=mean(p)std=std(p)sum=sum(p,2)[X,m]=max(sum)[Y,n]=min(sum)[zcj,xsxh]=sort(p,'descend')例72：某气象观测得某日6:00~18:00之间每隔2h的室内外温度如下表所示。试用三次样条插值分别求出该日室内外6:30~18:30之间每隔2h各点的近似温度（0C）。解：clc;clear;h=6:2:18;t1=[18.0,20.0,22.0,25.0,30.0,28.0,24.0];t2=[15.0,19.0,24.0,28.0,34.0,32.0,30.0];H=6.5:2:18.5;T1=interp1(h,t1,H,'spline');T2=interp1(h,t2,H,'spline');figureplot(H,T1,'-k',H,T2,'--')例73：有3个多项式P1(x)=x4+2x3+4x2+5，P2(x)=x+2，P3(x)=x2+2x+3，试进行下列操作：(1)求P(x)=P1(x)+P2(x)P3(x)。(2)求P(x)的根。(3)当x取矩阵A的每一元素时，求P(x)的值。其中：(4)当以矩阵A为自变量时，求P(x)的值。其中A的值与第(3)题相同。解：clc;clear;p1=[12405];p2=[12];p3=[123];p=p1+[0,conv(p2,p3)]x1=roots(p)A=[-11.2-1.4;0.7523.5;052.5];y1=polyval(p,A)y2=polyvalm(p,A)例74：用数值方法求定积分。解：clc;clear;g=inline('sqrt(cos(t.^2)+4*sin((2*t).^2)+1)');I=quad(g,0,2*pi)g2=inline('log(1+x)./(1+x.^2)');I2=quad(g2,0,2*pi)例75：（1）组合成一个43×的矩阵，第一列为按列顺序排列的a矩阵元素，第二列为按列顺序排列的b矩阵元素，第三列为按列顺序排列的c矩阵元素；（2）按照a、b、c的列顺序组合成一个行向量。clc;clear;a=[42;75];b=[71;83];c=[59;62];A=[a(:)b(:)c(:)]reshape(A,1,12)例76：先求下列表达式的值，然后显示MATLAB工作空间的使用情况。clc;clear;%(1)x=[2,1+2i;-0.45,5];z2=0.5*log(x+sqrt(1+x.^2))%(2)a=-3.0:0.1:3.0;z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)%(3)t=0:0.5:2.5;z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^2-2*t+1)例77：(1)求它们的乘积C。(2)将矩阵C的右下角3×2子矩阵赋给D。(3)查看MATLAB工作空间的使用情况。clc;clear;A=(reshape(1:25,5,5)')B=[3016;17-69;023-4;970;41311];C=A*BD=C(3:5,2:3)例78：(1)y<3时,的最大n值。(2)与(1)的n值对clc;clear;n=0;y=0;while(y<3)y=y+1/(2*n-1);n=n+1;endy=y-1/(2*n-1)n=n-2例79：已知y1=x2，y2=cos(2x)，y3=y1×y2，完成下列操作：(1)在同一坐标系下用不同的颜色和线型绘制三条曲线。(2)以子图形式绘制三条曲线。(3)分别用条形图、阶梯图、杆图和填充图绘制三条曲线。clc;clear;x=-pi:0.01;pi;y1=x.^2;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'-',x,y2,'--',x,y3,'-.');figure;subplot(1,3,1);plot(x,y1,'-');subplot(1,3,2);plot(x,y1,'--');subplot(1,3,3);plot(x,y1,'-.');figuresubplot(2,2,1);bar(x,y1);subplot(2,2,2);stairs(x,y1);subplot(2,2,3);stem(x,y1);subplot(2,2,4);fill(x,y1,'r');例80：绘制极坐标曲线ρ=asin(b+nθ)，并分析参数a、b、n对曲线形状的影响。clc;clear;t=0:pi/100:2*pi;a=input('输入a=');b=input('输入b=');n=input('输入n=')r=a*sin(b+n*t);subplot(4,2,1);polar(t,r,'m')holdon%控制变量法分析三个参数的影响%n的影响input('控制a和b不变，改变n的值')n=input('再输入一个n=')r=a*sin(b+n*t);subplot(4,2,3);polar(t,r,'m');holdoninput('控制a和b不变，再改变n的值')n=input('再输入一个n=')r=a*sin(b+n*t);subplot(4,2,4);polar(t,r,'m');holdon%a的影响input('控制n和b不变，改变a的值')a=input('再输入一个a=')r=a*sin(b+n*t);subplot(4,2,5);polar(t,r,'m');holdoninput('控制n和b不变，再改变a的值')a=input('再输入一个a=')r=a*sin(b+n*t);subplot(4,2,6);polar(t,r,'m');holdon%b的影响input('控制a和n不变，改变b的值')b=input('再输入一个b=')r=a*sin(b+n*t);subplot(4,2,7);polar(t,r,'m');holdoninput('控制a和nb不变，再改变b的值')b=input('再输入一个b=')r=a*sin(b+n*t);subplot(4,2,8);polar(t,r,'m');holdon分析结论：采用控制变量的方法，固定两个参数，然后变一个参数当a,n固定时，图形的形状也就固定了，b只影响图形的旋转角度。当a,b固定是,n只影响图形扇形数，特别的，当n是奇数时，扇形对数就是n，当是偶数时，扇形对数是2n个。当b,n固定时，a影响的是图形大小，特别的，当a是整数时，图形半径就是a。例81：绘制函数的曲线图和等高线。其中x的21个值均匀分布[-5，5]范围，y的31个值均匀分布在[0，10]，要求使用subplot(2,1,1)和subplot(2,1,2)将产生的曲面图和等高线图画在同一个窗口上。clc;clear;x=linspace(-5,5,21);y=linspace(0,10,31);[x,y]=meshgrid(x,y);z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2)/4);subplot(2,1,1);surf(x,y,z);title('曲面图');subplot(2,1,2);surfc(x,y,z);title('等高线图');例82：绘制曲面图形，并进行插值着色处理。clc;clear;s=0:pi/100:pi/2;t=0:pi/100:3*pi/2;[s,t]=meshgrid(s,t);x=cos(s).*cos(t);y=cos(s).*sin(t);z=sin(s);subplot(2,2,1);mesh(x,y,z);title('δ×ÅÉ«μÄͼÐÎ');subplot(2,2,2);surf(x,y,z);title('shadingfaceted(ȱʡ)');subplot(2,2,3);surf(x,y,z);title('shadingflat');shadingflat;subplot(2,2,4);surf(x,y,z);shadinginterp;title('shadinginterp');例83：已知lgx在[1,101]区间10个整数采样点的函数值如下表所示。试求lgx的5次拟合多项式p(x)，并绘制出lgx和p(x)在[1,101]区间的函数曲线。clc;clear;x=1:10:101;y=log10(x);P=polyfit(x,y,5)Y1=polyval(P,x);plot(x,y,':o',x,Y1,'-*')例84：分别用3种不同的数值方法解线性方程组。解：clc;clear;%1invA=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];B=[-4;13;1;11];x1=(inv(A)*B)'%2LU[L,U]=lu(A);LU=L*U;x2=(U\(L\B))'%3QR[Q,R]=qr(A);QR=Q*R;x3=(R\(Q\B))'例85：要求：分别用循环结构和向量运算（使用sum函数）来实现。clc;clear;%Ñ»·y=0;n=input('n=');fori=1:ny=y+1/i/i;endpi=sqrt(6*y)%sumclear;n=input('n=');a=1:n;b=1./a.^2;pi=sqrt(6*sum(b))例86：一个自然数是素数，且它的数字位置经过任意对换后仍为素数。例如13是绝对素数。试求所有两位绝对素数。函数文件：functionf=mat1(x)%是否是素数f=1;m=fix(sqrt(x));fori=2:mifrem(x,i)==0f=0;breakendend命令文件：clc;clear;fori=10:99;j=10*rem(i,10)+fix(i/10);ifmat1(i)==mat1(j)disp(i)endend例87：级数符号求和。clc;clear;%(1)symsn;symsum(1/(2*n-1),1,10)%(2)symsnx;symsum(n^2*x^(n-1),n,1,inf)symsn;symsum(n^2/5^n,1,inf)例88：建立一个图形窗口，使之背景颜色为红色，并在窗口上保留原有的菜单项，而且在按下鼠标器的左键之后显示出LeftButtonPressed字样。hf=figure('Color',[1,0,0],'MenuBar','figure',...'WindowButtonDownFcn','display(''LeftButtonPressed'')')例89：求代数方程的数值解。（有初值）（1）函数：functiong=f(x)g=3*x+sin(x)-exp(x);命令：clc;clear;fzero('f',1.5)（2）函数：FunctionF=fun(X)x=X(1);y=X(2);z=X(3);F(1)=sin(x)+y^2+log(z)-7;F(2)=3*x+2-z^3+1;F(3)=x+y+z-5;命令：clc;clear;X=fsolve(‘fun’,[1,1,1]’,optimset(‘Display’,’off’))例90：求微分方程的数值解。函数文件：functionxdot=sys(x,y)xdot=[y(2);(5*y(2)-y(1))/x];命令：clc;clear;x0=1.0e-9;xf=20;[x,y]=ode45('sys',[x0,xf],[00])例91:函数文件：functiony=fx(x)y=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);命令：clc;clear;y=fx(2)a=[1,2;3,4];y=fx(a)例92：（1）函数文件：functionf1=mat3(n)f1=n+10*log(n*n+5);命令：clc;clear;y=mat3(40)/(mat3(30)+mat3(20))（2）函数文件:functionf2=mat2(n)i=1:n;m=i.*(i+1);f2=sum(m);end命令：clc;clear;y=mat2(40)/(mat2(30)+mat2(20))例93：先利用默认属性绘制曲线y=x2e2x，然后通过图形句柄操作来改变曲线的颜色、线型和线宽，并利用文件对象给曲线添加文字标注。x=-2:0.01:2;y=x.^2.*exp(2*x);h=line(x,y);set(h,'color','r','linestyle',':','linewidth',2)text(1,exp(2),'y=x^2*exp(2*x)')例94：为图形窗口建立快捷菜单用以控制窗口的背景颜色和大小（实验十二第二题）clear;g=figure('Position',[500,330,300,320],'Name','图形','NumberTitle','off',...'menubar','none');h=uicontextmenu;hc=uimenu(h,'Label','颜色');hw=uimenu(h,'Label','大小');uimenu(hc,'Label','红色','Call','set(g，"Color","r");');uimenu(hc,'Label','绿色','Call','set(g，"Color","g");');uimenu(hw,'Label','放大','Call','set(g，"position",[200,130,450,300]);');uimenu(hw,'Label','缩小','Call','set(g，"position",[200,130,200,100]);');set(g,'UIContextMenu',h);例95：求函数在指定区间的极值。函数文件：functionf=g(u)x=u(1);y=u(2);f=2*x.^3+4*x.*y^3-10*x.*y+y.^2;命令：clc;clear;formatlongf=inline('(x^3+cos(x)+x*log(x))/exp(x)');[x,fmin1]=fminbnd(f,0,1)[U,fmin2]=fminsearch('g',[0,0])例96：求微分方程组的通解。clc;clear；[xyz]=dsolve('Dx=2*x-3*y+3*z',...'Dy=4*x-5*y+3*z','Dz=4*x-4*y+2*z','t')例97：设a=[1-23;45-4;5-67]，（1）请编程找出小于1的全部矩阵元素；（2）找到这些矩阵元素的位置，包括单下标和全下标两种形式；（3）并将其单下标转换成全下标。clear,clca=[1-23;45-4;5-67];[x,y]=find(a<0);c=[];fori=1:length(x)c(i,1)=a(x(i),y(i));c(i,2)=x(i);c(i,3)=y(i);c(i,4)=(y(i)-1)*size(a,2)+x(i);End例98：分析下面这组稀疏矩阵，找到其分布规律，使用三元组方法与非三元组法，创建下列矩阵:2080000104006000方法一：clear,clcdata=[28146];ir=[11234];jc=[13421];s=sparse(ir,jc,data,4,4);full(s)方法二：不用三元组法clear,clca=zeros(4,4);a(1,[1,3])=[2,8];a(2,4)=1;a(3,2)=4;a(4,1)=6;例99：设有矩阵A和B，A=[12345;678910;1112131415;1617181920;2122232425]，B=[3016;17-69;023-4;970;41311]（1）求它们的乘积C；（2）将矩阵C的右下角3x2子矩阵赋给D；（3）将矩阵D的左上角1x2子矩阵赋给E；解：A=[12345;678910;1112131415;1617181920;2122232425];B=[3016;17-69;023-4;970;41311];c=A*BD=c([3:5],[2:3])D=c([3,4,5],[2,3])（1）.求[100，999]之间能被61整除的数及其个数;（2）.建立一个字符串向量，删除其中的大写字母；1)clear,clc2)clear,clct=[100:999];b=['asSldSjfkSlGjFhslFf']i=find(rem(t,61)==0);d=find(b<='Z'&b>='A');t(i)b(d)=[]d=size(i,1)*size(i,2)例100：已知A=[2310-780；41-45655；325032；6-549214]，取出其前3行构成矩阵B，其前两列构成矩阵C，其左下角3x2子矩阵构成矩阵D，B与C的乘积构成矩阵E，分别求E<D、E&D、E|D、~E|~D。解：A=[2310-780;41-45655;325032;6-549214];B=A([1:3],:)c=A(:,[1,2])D=A([2:4],[1,2])E=B*cE<D,E&D,E|D,~E|~D例101：编程将这两段文字转换为字符串格式，注意保持两段文字的格式：（1）在英式用法中，引号通常是单引号，如'Fire!'。（2）InGBusagequotationmarksareusuallysingle:'fire!'.解：a=['在英式用法中，引号通常是单引号，如''Fire!''。'];b=['InGBusagequotationmarksareusuallysingle:''fire!''.'];例102：用结构体数组来存储2名学生的基本情况数据，每名学生的数据包括学号、姓名、专业和2门课程的成绩。解：clear,clcxuesheng=struct('num',{222,22},'name',{'mike','hike'},'major','tele','score',{rand(2,1)*100,rand(2,1)*100})xuesheng.score13.求S=20+21+22+23+24+……+210的值。解：a=2*ones(1,10)b=cumprod(a)c=sum(b)+1aa=2*ones(1,11)aa(1)=1bb=cumprod(aa)cc=sum(bb)例103：建立矩阵A并回答有关问题A{1,1}=1;A{1,2}='Brenden';A{2,1}=reshape(1:9,3,3);A{2,2}={12,34,2;54,21,3;4,23,67};使用访问的方式如何将Brenden改写成BRENDEN？分别执行A(3)=[]和A{3}=[]后，A的值各是多少？并说明原因。解：A{1,1}=1;A{1,2}='Brenden';A{2,1}=reshape(1:9,3,3);A{2,2}={12,34,2;54,21,3;4,23,67};A{1,2}=upper(A{1,2})A{3}=[]A(3)=[]例104：利用MATLAB提供的rand函数和圆整函数随机生成4Ｘ3整数矩阵A，进行如下操作1）A各列元素的平均值和中值2）A的最大元素和最小元素及它们的位置3）求A的每行元素的和以及全部元素之和4）分别对A的每列元素按升序、每行元素按降序排列解：A=round(rand(4,3))A_mean=mean(A)A_median=median(A)A_max=max(max(A))A_min=min(min(A))[x,y]=find(A==A_max)[xx,yy]=find(A==A_min)A_maxlocal=[x(1)y(1)]A_minlocal=[xx(1)yy(1)]A_rsum=sum(A,2)A_sum=sum(sum(A))A_csort=sort(A)A_rsort=abs(sort(-A,2))例105.输入一个百分制成绩，要求输出成绩等级A、B、C、D、E。其中90分~100分为A，80分~89分为B，70分~79为C，60分~69分为D，60分以下为E。要求：1）分别用if语句和switch语句实现2）对不合理的成绩应输出出错信息“输入的成绩不合理”（若成绩出现小数，则只能是“.5”）解：score=input('请输入您的分数：');ifrem(score,0.5)==0switchscorecasenum2cell(90:0.5:100)disp('您的成绩等级为A')casenum2cell(80:0.5:89)disp('您的成绩等级为B')casenum2cell(70:0.5:79)disp('您的成绩等级为C')casenum2cell(60:0.5:69)disp('您的成绩等级为D')casenum2cell(0:0.5:59.5)disp('您的成绩等级为E')otherwisedisp('你输入的成绩不合理')endelsedisp('你输入的成绩不合理，若出现小数只能是0.5')end例106：设计程序，完成两位数的加、减、乘、除四则运算，即产生两个两位随机整数，再输入一个运算符号，做相应的运算，显示相应的结果，并要求结果显示类似于“a=x+y=34”。(sy312.m)解：a=input('请输入一个数：');b=input('请再输入一个数：');fuhao=input('请输入一个运算符号（+-*/）：','s');switchfuhaocase{'+'}he=a+b;disp(['输入的两数和=',num2str(a),'+',num2str(b),'=',num2str(he)])case{'-'}he=a-b;disp(['输入的两数差=',num2str(a),'-',num2str(b),'=',num2str(he)])case{'*'}he=a*b;disp(['输入的两数乘积=',num2str(a),'*',num2str(b),'=',num2str(he)])case{'/'}he=a/b;disp(['输入的两数商=',num2str(a),'/',num2str(b),'=',num2str(he)])otherwisedisp('请输入正确的符号')end例107：求下列分段函数的值要求：用if语句实现，分别输出x=-5.0，-3.0，1.0，2.0，2.5，3.0，5.0时的y值。其中x的值以向量的形式从键盘输入。(sy313.m)解：x=[];y=[];x=input('请输入x的值：');fork=1:length(x)ifx(k)<0&x(k)~=-3y(k)=x(k).^2+x(k)-6;elseifx(k)>=0&x(k)<10&x(k)~=2&x(k)~=3y(k)=x(k).^2-5*x(k)+6;elsey(k)=x(k).^2-x(k)-1;endendfori=1:length(y)disp(['y(',num2str(x(i)),')','=',num2str(y(i))]);end例108：使用for循环语句找出最小的n值，使得n!>10100,并求出n!。要求显示的结果为n=70；n!=1.197857e+100>1e100；70!=1.197857e+100>1e100解：用IF语句执行如下：s=1;forn=1:100s=s*n;if(s>10^100)break;endendstr1=['n=',num2str(n)];str2=['n!=',num2str(s)];disp(str1)disp(str2)例109：使用while循环语句找出最小的n值，使得n!>10100,并求出n!。要求显示的结果为n=70；n!=1.197857e+100>1e100；70!=1.197857e+100>1e100用WHILE语句执行如下：n=1;s=1;while(s<=10^100)s=s*n;n=n+1;endstr1=['n=',num2str(n-1)];str2=['n!=',num2str(s)];disp(str1)disp(str2)例110：Fibonacci数组的元素满足Fibonacci规则：ak+2=ak+ak+1,（k=1，2，3，…）；且a1=a2=1。请设计一段程序，求出该数组中第一个大于10000的元素要求显示的结果为：i=21；a(i)=10946；a(21)=10946方法1：a(1)=1;a(2)=1;fork=1:10000%不能是3:10000a(k+2)=a(k)+a(k+1);ifa(k)>10000break;endenddisp(['k=',num2str(k)])disp(['a(k)=',num2str(a(k))])法二：clear,clcs=1;forn=1:infs=s*n;ifs>10^100breakendendfprintf('%d!=%1.6e>1e100\n',n,s)例111：已知f1=1,n=1；f2=0,n=2；f3=1,n=3；fn=fn-1-2fn-2+fn-3,n>3求f1~f50中：1）最大值和最小值及它们的位置，各数之和（最大值给变量MAX，其位置给变量x，最小值给变量MIN，其位置给变量y，各数之和给变量SUM）。2）正数、零、负数的个数（依次赋值给变量positive，zero，negative）。方法一：clear;f(1)=1;f(2)=0;f(3)=1;forn=4:50f(n)=f(n-1)-2*f(n-2)+f(n-3);end[MAX,x]=max(f)%注意不要分号，要执行的语句[MIN,y]=min(f)SUM=sum(f)positive=length(find(f>0))zero=length(find(f==0))negative=length(find(f<0))MAX=406631x=49MIN=-403795y=50SUM=76861positive=25zero=2negative=23方法二：clearall,clcforn=1:50ifn==1f(1)=1;elseifn==2f(2)=0;elseifn==3f(3)=1;elseifn>3f(n)=f(n-1)-2*f(n-2)+f(n-3);endend[MAX,y]=max(f)[Min,x]=min(f)positive=length(find(f>0))zero=length(find(f==0))negative=length(find(f<0))例112：假设已知矩阵，试给出相应的MATLAB命令，将其全部偶数行提取出来，赋给矩阵，用命令生成矩阵，用上述命令检验一下结果是不是正确。>>A=magic(8)A=64236160675795554121351501617474620214342244026273736303133323435292838392541232244451918484915145253111056858595462631>>B=A(2:2:end,:)B=95554121351501640262737363031334123224445191848858595462631例113：用数值方法可以求出，试不采用循环的形式求出和式的数值解。由于数值方法是采用double形式进行计算的，难以保证有效位数字，所以结果不一定精确。试采用运算的方法求该和式的精确值。（1）数值解：i=[0:63];sum(2.^i)ans=1.8447e+019（2）解析解：>>symsk;>>symsum(2^k,0,63)ans=184********709551615例114：根据，求的近似值，当n分别取100,1000,10000时，结果是多少？要求：分别用循环结构和向量运算（使用sum函数）来实现。（1）(循环结构)解：s=0;n=input('n=?');fori=1:ns=s+1/i/i;endPI=sqrt(6*s)Pi°n=?100PI=3.1321ans=3.1416s=0;n=input('n=?');fori=1:ns=s+1/i/i;endPI=sqrt(6*s)pin=?1000PI=3.1406ans=3.1416s=0;n=input('n=?');fori=1:ns=s+1/i/i;endPI=sqrt(6*s)pin=?10000PI=3.1415ans=3.1416（2）向量运算解：n=input('n=?')