MATLAB习题key

1、练习：（一）：先求下列表达式的值，然后显示在MATLAB工作空间的使用情况并全部保存到变量名为: xitiZ1=sind(85)/(1+exp(2)= 0.1187x=(5+cosd(47)/(1+sqrt(7)-2i)=1.1980 + 0.6572i（3）已知：x=0:0.5:8 ，求y=sinx2x=0:0.5:8y=sin(x.2)plot(x,y)（二）保存到mb_01, mb_02 .1.写出完成下列操作的命令（1）建立3阶单位矩阵A； A=eye(3)（2）建立5x6阶随机矩阵A1; A1=rand(5,6)（3）将矩阵A1第25行中的第1，3，5列元素赋给矩阵B; B=A1(2

2、:5,1),A1(2:5,3),A1(2:5,5)（4）删除矩阵A的第7号元素；A(7)=（三）设：（1） 求:f(x)+g(x) ,f(x)-g(x);（2） f(x)*g(x) ,f(x)/g(x)syms xf1=3*x5-5*x4+2*x3-7*x2+5*x+6f2=3*x2+5*x-3(1) f1+f2= 3*x5-5*x4+2*x3-4*x2+10*x+3 f1-f2= 3*x5-5*x4+2*x3-10*x2+9(2)f1*f2= (3*x5-5*x4+2*x3-7*x2+5*x+6)*(3*x2+5*x-3) f1/f2= (3*x5-5*x4+2*x3-7*x2+5*x+6)

3、/(3*x2+5*x-3)(四)求多项式x4+8x3-10的根clear all;clc;P=1 8 0 0 10 poly2sym(P) X=roots(P) M=poly(X)P = 1 8 0 0 10ans = x4+8*x3+10X = -7.9803 -1.1335 0.5569 + 0.8918i 0.5569 - 0.8918iM = 1.0000 8.0000 -0.0000 0.0000 10.00003. 计算与的数组乘积。a=6,9,3;2,7,5b=2,4,1;4,6,8a.*bans = 12 36 3 8 42 403 对于，如果，求解X。A=4,9,2;7,6,

4、4;3,5,7B=37;26;28X=ABx = -0.51184.04271.33184， 求解多项式x3-7x2+2x+40的根。clear all;clc;P=1 -7 2 40poly2sym(P) X=roots(P) M=poly(X)ans =x3-7*x2+2*x+40X = 5.0000 4.0000 -2.0000M = 1.0000 -7.0000 2.0000 40.00005.求下列极限： ， ,syms a xf=1+a/xlimit(f,x,inf,'left') ans =1syms a x mf=(x(1/m)-a(1/m)/(x-a)limi

5、t(f,x,a)ans =a(1/m)/a/m6.已知：x=0:0.5:8 ，求y=sinx2，并绘制图形图形，要求1）采用红色虚线2）作坐标轴，图形标注，加图例3）加网格线，4）对点（0，0.7）处添加说明曲线，(1)x=0:0.5:8y=sin (x.*x)plot (x,y,'r-')(2)xlabel(x=0:0.5:8)ylabel('y=sin (x.*x)')legend('y=sin (x.*x)') (3)grid on(4)text(0,0.7,'曲线')7.求积分： 和 syms xf=sqrt(x3+x4)

6、int(f,x)ans = -1/48*(x3+x4)(1/2)*(-16*(x2+x)(3/2)+12*(x2+x)(1/2)*x+6*(x2+x)(1/2)-3*log(1/2+x+(x2+x)(1/2)/x/(x+1)*x)(1/2)syms xf=x3/ (x-1)100int(f,x,2,3)ans =8.已知：f(x)=3x5+4x3-5x2-7.2x+5 (1)计算f(x)=0的全部根。 （2）计算f(5), f(7.8), f(9.6), f(12.3)的值。(1)P=3 0 4 -5 -7.2 5poly2sym(P) X=roots(P) X = -0.3046 + 1.6

7、217i -0.3046 - 1.6217i -1.0066 1.0190 0.5967 (2) P=3 0 4 -5 -7.2 5poly2sym(P,5)poly2sym(P,7.8)poly2sym(P,9.6)poly2sym(P,12.3)ans =9719ans =275493997/3125ans =773830729/3125ans =85119516329/1000009.已知：f(x)=3x5-5x4+2x3-7x2-5x+6 g(x)=3x2+5x-3 (1) 求f(x)+g(x)和f(x)-g(x)。P=3 -5 2 -7 -5 6Q=3 5 -3f=poly2sym(

8、P)g=poly2sym(Q)f+gf-gans =3*x5-5*x4+2*x3-4*x2+3ans =3*x5-5*x4+2*x3-10*x2-10*x+9（2）求f(x) g(x)和f(x)/g(x)。simplif(f*g)simplify(f/g)ans =9*x7-28*x5+4*x4-56*x3+14*x2+45*x-18ans =(3*x5-5*x4+2*x3-7*x2-5*x+6)/(3*x2+5*x-3)10. 在一个图形窗口中以子图形式同时绘制正弦、余弦、正切、余切曲线。 x=linspace(0,6)subplot(221)plot(x,sin(x)xlabel('

9、;x')ylabel('sin(x)')title('sin(x)')subplot(222)plot(x,cos(x)xlabel('x')ylabel('cos(x)')title('cos(x)')y1=tan(x);y2=cot(x)subplot(223)plot(x,y1)xlabel('x')title('tan(x)')subplot(224)plot(x,y2)xlabel('x')title('cot(x)')ylabel(

10、'cot(x)')ylabel('tan(x)')11.化简：s=(x2+y2)2+(x2-y2)2syms x ys=(x2+y2)2+(x2-y2)2simplify(s)ans = 2*x4+2*y412求级数之和：syms n xs=symsum(n*xn,n,x,inf)s =1/(x-1)2*(-x(2+x)+2*x*xx)13.已知表达式为： 试编写能对任意一组输入x值求相应y值的程序，令函数名为fun13。function y=fun13(x)x=input('x=?')if x>=1 y=10;elseif x>1&

11、amp;x<1 y=0;else y=-10end14编写一个函数文件，用于求两个矩阵的乘积和点乘，然后在命令文件中调用该函数。假设给定矩阵：1）A为一个三阶单位矩阵，A1为2行3列随机矩阵2）B1为由第一小题的第（3）问求得的结果矩阵B而生成的新矩阵，即抽取B的13行所有列元素得到B1。试计算：A*B1与A.*B1，A1*B1，A1.*B1（可以吗？）并比较结果，体会乘与点乘的区别。建立M文件（1）function z=chengji(x,y)z=x*y（2）function z=diancheng(x,y)z=x.*yA=eye(3)A1=rand(2,3)A1=rand(5,6)B

12、=A1(2:5,1),A1(2:5,3),A1(2:5,5)B1=B(1:3,1),B(1:3,2),B(1:3,3)Q1=chengji(A,B1)Q2=diancheng(A,B1)Q3=chengji(A1,B1)Q4=diancheng(A1,B1)15以子图的形式绘制曲线：（1） （3）（2） （4） x=linspace(0,6)t=linspace(0,6)subplot(221)y1=100./(1+x.2)plot(x,y1)xlabel('x')ylabel('y')title('y=100/(1+x2)')subplot(2

13、22)y2=sqrt(1-x.2)plot(x,y2)xlabel('x')ylabel('y')title('x2+y2=1')subplot(223)y3=1/(2*pi)*exp(-x.2/2)plot(x,y3)xlabel('x')ylabel('y')title(' y=1/(2*pi)*exp(-x2/2)')subplot(224)x=t.2y4=5*t.3plot(x,y4)xlabel('x')ylabel('y')16：已知 ，当n=100时，求 的值。（for循环）17已知某商场对顾客所买的商品实行打折销售，标准如下（商品价格用price来表示）：price<200 没有折扣 3%折扣 5%折扣 8%折扣 10%折扣 14%折扣提示：用switch语句，另外函数fix(x)表示对x向零方向取整数。如fix(2.8)=2; fix(-2.8)=-2。a=fix(price/100)switch a case 1 p=price case 2,3,4 p=price*0.97 case 5,6,