资源描述
精选优质文档-倾情为你奉上1.完成下列操作：(1) 求100,999之间能被21整除的数的个数。(2) 建立一个字符串向量，删除其中的大写字母。(1)m=100:999;n=find(mod(m,21)=0);length(n)(2)ch=SVHBFguisdyihnFNKSDFBNJDFNnjkdf;k=find(ch=A&ch=Z);ch(k)=2. 设有分块矩阵，其中E、R、O、S分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证。E=eye(3,3);R=rand(3,2);O=zeros(2,3);S=diag(1,2);A=E,R;O,S;a=E,R+R*S;O,S2A2A2-a3.下面是一个线性方程组：(1) 求方程的解。(2) 将方程右边向量元素b3改为0.53再求解，并比较解的相对变化。(3) 计算系数矩阵A的条件数并分析结论。 (1)a=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6;b1=0.95,0.67,0.52;x1=inv(a)*b1 (2)b2=0.95,0.67,0.53;x2=inv(a)*b2 (3)c=cond(a)4. 求分段函数的值。用if语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y值。x=input(please input x:);if x=0&x5&x=2&x=3 y=x*x-5.*x+6;else y=x*x-x-1;endy5. 输入一个百分制成绩，要求输出成绩等级A、B、C、D、E。其中90分100分为A，80分89分为B，79分79分为C，60分69分为D，60分以下为E。要求：(1) 分别用if语句和switch语句实现。(2) 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。score=input(please input score:);if score=0 switch fix(score/10) case 9,10; grade=A; case 8 grade=B; case 7 grade=C; case 6 grade=D; otherwise grade=E; end gradeelse disp(错误);end6. 硅谷公司员工的工资计算方法如下：(1) 工作时数超过120小时者，超过部分加发15%。(2) 工作时数低于60小时者，扣发700元。(3) 其余按每小时84元计发。试编程按输入的工号和该号员工的工时数，计算应发工资。n=input(please input number:);t=input(please input time:);switch t case t120 w=t*84+(t-120)*15/100*84; case t=0.5&x1.5 y=a*x2+b*x+c;elseif x=1.5 y=a.*sin(b)c+x;elseif x=3.5&x0)length(find(f=0)length(find(f0.5n=n+1;endendp=n/3000023. 将100个学生5门功课的成绩存入矩阵P中，进行如下处理：(1) 分别求每门课的最高分、最低分及相应学生序号。(2) 分别求每门课的平均分和标准方差。(3) 5门课总分的最高分、最低分及相应学生序号。(4) 将5门课总分按从大到小顺序存入zcj中，相应学生序号存入xsxh。提示：为避免输入学生成绩的麻烦，可用取值范围在45,95之间的随机矩阵来表示学生成绩。%(1)A=45+51*rand(100,5);Y,U=max(A)a,b=min(A)%(2)m=mean(A)s=std(A)%(3)sum(A,2)Y,U=max(ans)a,b=min(ans)%(4)zcj,xsxh=sort(ans)24. 某气象观测得某日6:0018:00之间每隔2h的室内外温度如下表所示。时间h 6 8 10 12 14 16 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试用三次样条插值分别求出该日室内外6:3018:30之间每隔2h各点的近似温度（0C）。h=6:2:18;x=6.5:2:18.5;t1=18,20,22,25,30,28,24;t2=15,19,24,28,34,32,30;T1=spline(h,t1,x)T2=spline(h,t2,x)25. 有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)题相同。%(1)p1=1,2,4,0,5;p2=1,2;p3=1,2,3;p=p1+0,conv(p2,p3) %为使两向量大小相同，所以补0%(2)A=roots(p)%(3)A=-1,1.2,-1.4;0.75,2,3.5;0,5,2.5;polyval(p,A)%(4)polyvalm(p,A)26. 用数值方法求定积分。(1) 的近似值。(2) clc;clear;f=inline(sqrt(cos(t.2)+4*sin(2*t).2+1);I1=quad(f,0,2*pi) g=inline(log(1+x)./(1+x.2);I2=quad(g,0,2*pi) 27、将矩阵、和组合成两个新矩阵：（1）组合成一个的矩阵，第一列为按列顺序排列的矩阵元素，第二列为按列顺序排列的矩阵元素，第三列为按列顺序排列的矩阵元素；（2）按照、的列顺序组合成一个行向量。a=4 2;5 7; b=7 1;8 3; c=5 9;6 2; d=a(:) b(:) c(:) e=a(:);b(:);c(:)28、有一组测量数据满足，t的变化范围为010，用不同的线型和标记点画出、和三种情况下的曲线，并在图中添加标题、用箭头线标识出各曲线的取值和图例框。 t=linspace(0,10)y1=exp(-0.1*t);y2=exp(-0.2*t);y3=exp(-0.5*t);plot(t,y1,s-,t,y2,d,t,y3,v:)title( y=exp(-a*t);text(t(30),y1(30),leftarrow a=0.1); text(t(40),y2(40),leftarrow a=0.2); text(t(30),y3(30),leftarrow a=0.5);legend(a=0.1,a=0.2,a=0.5)29. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况。(1) ，其中(2) (3) ，其中t=0:0.5:2.5%(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.2.*(t=0&t=1&t=2&t3)30. 设有矩阵A和B：(1) 求它们的乘积C。(2) 将矩阵C的右下角32子矩阵赋给D。(3) 查看MATLAB工作空间的使用情况。A=1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25B=3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11C=A*BF=size(C)D=C(F(1)-2:F(1),F(2)-1:F(2)whos31. 根据，求：(1) y3时的最大n值。(2) 与(1)的n值对应的y值。y=0;n=1;while(y1e-5 Xn=Xn1; Xn1=a/(b+Xn); n=n+1; if n=500 break; endendnXn1r1=(-b+sqrt(b*b+4*a)/2r2=(-b-sqrt(b*b+4*a)/233. 已知y1=x2，y2=cos(2x)，y3=y1y2，完成下列操作：(1) 在同一坐标系下用不同的颜色和线型绘制三条曲线。(2) 以子图形式绘制三条曲线。(3) 分别用条形图、阶梯图、杆图和填充图绘制三条曲线。%（1）x=linspace(-2*pi,2*pi,100);y1=x.2;y2=cos(2*x);y3=y1.*y2;plot(x,y1,b-,x,y2,r:,x,y3,y-); text(4,16,leftarrow y1=x2);text(6*pi/4,-1,downarrow y2=cos(2*x);text(-1.5*pi,-2.25*pi*pi,uparrow y3=y1*y2);%（2）x=linspace(-2*pi,2*pi,100);y1=x.2;y2=cos(2*x);y3=y1.*y2;subplot(1,3,1);%分区plot(x,y1);title(y1=x2);%设置标题subplot(1,3,2);plot(x,y2);title(y2=cos(2*x);subplot(1,3,3);plot(x,y3);title(y3=x2*cos(2*x);%（3）x=linspace(-2*pi,2*pi,20);y1=x.2;subplot(2,2,1);%分区bar(x,y1);title(y1=x2的条形图);%设置标题subplot(2,2,2);stairs(x,y1);title(y1=x2的阶梯图);subplot(2,2,3);stem(x,y1);title(y1=x2的杆图);subplot(2,2,4);fill(x,y1,r);%如果少了r则会出错title(y1=x2的填充图);%其他的函数照样做。34. 绘制极坐标曲线=asin(b+n)，并分析参数a、b、n对曲线形状的影响。a=input(a=);b=input(b=);n=input(n=);t=-2*pi:0.01:2*pi;r=a*sin(b+n*t);polar(t,r)35. 绘制函数的曲线图和等高线。其中x的21个值均匀分布-5，5范围，y的31个值均匀分布在0，10，要求使用subplot(2,1,1)和subplot(2,1,2)将产生的曲面图和等高线图画在同一个窗口上。x=linspace(-5,5,21);y=linspace(0,10,31); x,y=meshgrid(x,y);%在-5,5*0,10的范围内生成网格坐标z=cos(x).*cos(y).*exp(-sqrt(x.2+y.2)/4);subplot(2,1,1);surf(x,y,z);subplot(2,1,2);contour3(x,y,z,50);%其中50为高度的等级数，越大越密36. 绘制曲面图形，并进行插值着色处理。ezsurf(cos(s)*cos(t),cos(s)*sin(t),sin(s),0,0.5*pi,0,1.5*pi); %利用ezsurf隐函数shading interp %进行插值着色处理37. 已知lgx在1,101区间10个整数采样点的函数值如下表所示。 x 1 11 21 31 41 51 61 71 81 91 101lgx 0 1.0414 1.3222 1.4914 1.6128 1.7076 1.7853 1.8513 1.9085 1.9510 2.0043试求lgx的5次拟合多项式p(x)，并绘制出lgx和p(x)在1,101区间的函数曲线。x=1:0.1:101;y1=log10(x);p=polyfit(x,y1,5)y2=polyval(p,x);plot(x,y1,:,x,y2,-)38. 分别用3种不同的数值方法解线性方程组。A=6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2;b=-4,13,1,11;L,U=lu(A);x=U(Lb)%直接方法x=Ab%QR分解法Q,R=qr(A);x=R(Qb)%Cholesky分解(必须是正定矩阵)（此题为非正定矩阵）%R=chol(A);39. 用符号方法求下列积分。%(1)x=sym(x);f=1/(1+x4+x8);int(f)%(2) x=sym(x);f=1/(asin(x)2*sqrt(1-x2)1/2);int(f)%(3)x=sym(x);f=(x2+1)/(x4+1);int(f,0,inf)%(4)x=sym(x);f=exp(x).*(1+exp(x)2;int(f,0,log(2)40. 根据，求的近似值。当n分别取100、1000、10000时，结果是多少？要求：分别用循环结构和向量运算（使用sum函数）来实现。%循环实现clearn=input(place enter n: );sum=0;for i=1:n sum=sum+1/(i*i);pi=sqrt(6*sum);endpi%向量实现clearn=input(pi=place enter n: );i=1:n;f=1./i.2;y=sum(f);pi=sqrt(6*y);pi41. 一个自然数是素数，且它的数字位置经过任意对换后仍为素数。例如13是绝对素数。试求所有两位绝对素数。 function flag=sushu(n)for k=2:sqrt(n) flag=1; if rem(n,k)=0 flag=0; break; endend %在命令窗口输入clearfor i=10:99 j=10*rem(i,10)+fix(i/10); if sushu(i)&sushu(j) disp(i) endend42. 级数符号求和。(1) 计算。(2) 求级数的和函数，并求之和。%（1）n=sym(n);s=symsum(1/(2*n-1),n,1,10)%（2）n=sym(n);x=sym(x);s=symsum(n2*x(n-1),n,1,inf)s=symsum(n2/5n,n,1,inf)43. 建立一个图形窗口，使之背景颜色为红色，并在窗口上保留原有的菜单项，而且在按下鼠标器的左键之后显示出Left Button Pressed字样。44. 求代数方程的数值解。(1) 3x+sinx-ex=0在x0=1.5附近的根。(2) 在给定的初值x0=1，y0=1，z0=1下，求方程组的数值解。%(1)function f=fz(x)f=3*x+sin(x)-exp(x); %在命令窗口输入fzero(fz,1.5) %(2)function F=myfun(X)x=X(1);y=X(2);z=X(3);F(1)=sin(x)+y2+log(z)-7;F(2)=3*x+2y-z3+1;F(3)=x+y+z-5; %在命令窗口输入X=fsolve(myfun,1,1,1,optimset(Display,off)45. 求微分方程的数值解。function ydot=sys(x,y)ydot=(5*y(1)-y(2)/x;y(1); %在命令窗口输入x0=0;xf=20;x,y=ode45(sys,x0,xf,0,0);x,y46. 设，编写一个MATLAB函数文件fx.m，使得调用f(x)时，x可用矩阵代入，得出的f(x)为同阶矩阵。function y=fx(x)y=1./(x-2).2+0.1)+1./(x-3).4+0.01); %在命令窗口输入y=fx(2)a=1,2,3,4;y=fx(a)47.已知(1) 当f(n)=n+10ln(n2+5)时，求y的值。(2) 当f(n)=12+23+34+.+n(n+1)时，求y的值。function f1=m(n)f1=n+10*log(n*n+5) %在命令窗口输入y=m(40)/(m(30)+m(20)function f2=m(n)i=1:n;f2=sum(i.*(i+1); %在命令窗口输入y=m(40)/(m(30)+m(20)48. 先利用默认属性绘制曲线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=x2*exp(2*x)49. 用符号方法求下列极限或导数。%(1)syms x;f=(pi1/2-acos(x)1/2)/(x+1)1/2;limit(f,x,-1,right)%(2)syms a x t;A=ax t3;t*cos(x) log(x);diff(A,x)diff(A,t,2)diff(A,x)*diff(A,t)%(3)syms x zf=(x2-2*x)*exp(-x2-z2-x*z); Zx=-diff(f,x)/diff(f,z)dfxz=diff(diff(f,x),z);x=sym(0);z=sym(1);eval(dfxz) 50. 求微分方程初值问题的符号解，并与数值解进行比较。clear all;clc;dsolve(D2y+4*Dy+29*y,y(0)=0,Dy(0)=15,x)51. 为图形窗口建立快捷菜单，用以控制窗口的背景颜色和大小。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);52. 设计一个对话框，其中有一个编辑框和按钮，当单击按钮时，使编辑框的内容加5。 53. 求函数在指定区间的极值。(1) 在(0,1)内的最小值。(2) 在0,0附近的最小值点和最小值。 %1function f=g(x)f=(x.3+cos(x)+x*log(x)/exp(x); %在命令窗口输入x,fmin=fminbnd(g,0,1)%2function f=f(x,y)f=2*x3+4*x*y3-10*x*y+y2; %在命令窗口输入x,y,fmin=fminbnd(fx,0,0)54. 求下列方程的符号解。clear all;clc;x1=solve(log(x+1)-5/(1+sin(x)=2)x2=solve(x2+9*sqrt(x+1)-1)x3=solve(3*x*exp(x)+5*sin(x)-78.5)x4 y4=solve(sqrt(x2+y2)-100,3*x+5*y-8)55. 求微分方程组的通解。clear all;clc;x y z=dsolve(Dx=2*x-3*y+3*z,. Dy=4*x-5*y+3*z,Dz=4*x-4*y+2*z,t)56. 在图形窗口默认菜单条上增加一个Plot菜单项，利用该菜单项可以在本窗口绘制三维曲面图形。screen=get(0,ScreenSize);figure(Position,200,130,300,120,Name,图形,NumberTitle,off,.menubar,none);hplot=uimenu(gcf,Label,&plot);uimenu(hplot,Label,&ssk,Call,. x=-pi:pi/20:pi;,x,y=meshgrid(x);,. z=sin(x).*cos(y);,surf(x,y,z););专心-专注-专业
展开阅读全文