Lab 2 参考答案

第二届全国大学生化学实验笔试题第一题单项选择题（25分）1．下列化合物，在NaOH溶液中溶解度最大的是(A) PbCrO4(B) Ag2CrO4(C) BaCrO4 (D) CaCrO42．向酸性K2Cr2 O7溶液中加入H2O2，却未观察到蓝色物质生成，其原因肯定是(A) 未加入乙醚，因CrO5与乙醚的加合物为蓝色(B) 未加入戊醇，因CrO5萃取到戊醇中显蓝色(C) 未将溶液调至碱性，因CrO5在酸性介质中分解(D) 因K2Cr2O7和/或H2O2浓度过稀3．实验室配制洗液，最好的方法是(A) 向饱和K2Cr2O7溶液中加入浓硫酸(B) 将K2Cr2O7溶于热的浓硫酸(C) 将K2Cr2O7溶于1:1硫酸(D) 将K2Cr2O7与浓硫酸共热4．滴加0.1mol/L CaCl2溶液没有沉淀生成，再滴加氨水有白色沉淀生成，该溶液是(A) Na3PO4(B) Na2HPO4(C) NaH2PO4(D) 以上三种溶液均可5．从滴瓶中取少量试剂加入试管的正确操作是(A) 将试管倾斜，滴管口贴在试管壁，再缓慢滴入试剂(B) 将试管倾斜，滴管口距试管口约半厘米处缓慢滴入试剂(C) 将试管垂直，滴管口伸入试管内半厘米再缓慢滴入试剂(D) 将试管垂直，滴管口贴在试管壁，再缓慢滴入试剂6．制备下列气体时可以使用启普发生器的是(A) 高锰酸钾晶体与盐酸反应制氯气(B) 块状二氧化锰与浓盐酸反应制备氯气(C) 无水碳酸钾与盐酸反应制二氧化碳(D) 块状硫化亚铁与稀硫酸反应制备硫化氢7．实验室用浓盐酸与二氧化锰反应制备氯气，欲使氯气纯化则应依次通过(A) 饱和氯化钠和浓硫酸(B) 浓硫酸和饱和氯化钠(C) 氢氧化钙固体和浓硫酸(D) 饱和氯化钠和氢氧化钙固体8．使用煤气灯涉及的操作有①打开煤气灯开关②关闭空气入口③擦燃火柴④点燃煤气灯⑤调节煤气灯火焰。

点燃煤气灯时操作顺序正确的是(A) ①②③④⑤(B) ②①③④⑤(C) ②③①④⑤(D) ③②①④⑤9．能将Cr3+和Zn2+ 离子分离的溶液是(A) NaOH (B) NH3∙H2O (C) Na2CO3(D) NaHCO310．下列配制溶液的方法中，不正确的是(A) SnCl2溶液：将SnCl2溶于稀盐酸后加入锡粒(B) FeSO4溶液：将FeSO4溶于稀硫酸后放入铁钉(C) Hg(NO3)2溶液：将Hg(NO3)2溶于稀硝酸后加入少量Hg(D) FeCl3溶液：将FeCl3溶于稀盐酸11．由二氧化锰制锰酸钾，应选择的试剂是(A) 王水+ KCl (B) Cl2 + KCl(C) 浓H2SO4 + KClO3(D) KOH + KClO312．向酸性的KI溶液中滴加过量的H2O2有灰黑色沉淀生成，不能使该沉淀消失的是(A) 氢氧化钠溶液(B) 碘化钾溶液(C) 稀硝酸(D) 次氯酸钠溶液13．将少量KMnO4晶体放入干燥的试管中，在煤气灯上小火加热一段时间后冷却至室温，逐滴加入水，最先观察到溶液的颜色是(A) 粉红(B) 紫色(C) 绿色(D) 黄色14．将新生成的下列化合物在空气中放置，颜色最不易发生变化的是(A) Fe(OH)2(B) Ni(OH)2(C) Mn(OH)2(D) Co(OH)2 15．与浓盐酸作用有氯气生成的是(A) Fe2O3 (B) Pb2O3 (C) Sb2O3 (D) Bi2O3 16．向Hg2(NO3)2溶液中滴加氨水时，生成物的颜色为(A) 棕色(B) 灰黑色(C) 白色(D) 黄色17．下列配离子中，肯定不为蓝颜色的是(A) Cu(NH3)42+(B) Co(NH3)62+(C) Ni(NH3)62+(D) Co(SCN)42-18．向K2Cr2O7溶液中通入过量SO2，溶液的颜色为(A) 蓝色(B) 紫色(C) 绿色(D) 黄色19．下列化合物中，在6mol dm-3 NaOH溶液中溶解度最大的是(A) Mg(OH)2(B) Mn(OH)2(C) Fe(OH)2(D) Cu(OH)2 20．向澄清Na2S溶液中滴加稀盐酸，若有沉淀生成，则该沉淀是(A) Na2S2(B) S (C) Na2S3(D) SCl4 21．将下列混合物装入干燥试管中，在煤气灯上加热能得到黑色产物的是(A) KNO3 + Na2CO3(B) KNO3 + CuSO4(C) KNO3 + Cr2O3(D) KNO3 + ZnSO422．将浓硫酸与少量KI固体混合，还原产物主要是(A) SO2(B) H2SO3(C) H2S (D) I223．下列氢氧化物中，在氨水中溶解度最小的是(A) Zn(OH)2(B) Cu(OH)2 (C) Co(OH)2(D) Fe(OH)2 24．下列试剂不能鉴别SnCl2和SnCl4溶液的是(A) HgCl2(B) 溴水(C) NaOH (D) (NH4)2S 25．向Ag(S2O3)23-中通入过量的Cl2，生成的沉淀是(A) S (B) Ag2SO4(C) AgCl (D) Ag 2S2O3第二题单项选择题（25分）1．体积比为1:2的HCl其摩尔浓度为(A) 2 mol/L (B) 6 mol/L (C) 4 mol/L (D) 3 mol/L2．Fe(OH)3沉淀完后过滤时间是(A) 放置过夜(B) 热沉化后(C) 趁热(D) 冷却后3．在重量分析中，洗涤无定形沉淀的洗涤液是(A ) 冷水(B) 含沉淀剂的稀溶液(C) 热的电解质溶液(D) 热水4．重量分析中过滤BaSO4沉淀应选用的滤纸是(A) 慢速定量滤纸(B) 快速定性滤纸(C) 慢速定性滤纸(D) 快速定量滤纸5．用洗涤方法可除去的沉淀杂质是(A) 混晶共沉淀杂质(B) 包藏共沉淀杂质(C) 吸附共沉淀杂质(D) 后沉淀杂质6．现欲标定NaOH溶液的浓度，实验室提供下列物质，最好应选择(A) 邻苯二甲酸氢钾(KHC8H9O4) (B) 草酸（H2C2O4•2H2O）(C) 苯甲酸（C6H5COOH）(D) 甲酸(HCOOH)7．称取一定量无水碳酸钠溶解后定容于250 ml 容量瓶中，量取25 ml 用以标定盐酸，容量瓶和移液管采取的校准方法是(A)容量瓶绝对校准(B) 容量瓶和移液管相对校准(C) 移液管绝对校准(D) 不用校准8．碘量法测定铜主要误差来源是（A）I- 的氧化和I2的挥发（B）淀粉强烈吸附I2（C）I2被吸附在CuI上，终点颜色很深不易观察（D）KSCN加入过早9．定量分析中，基准物质是(A) 纯物质(B) 标准参考物质(C) 组成恒定的物质(D) 组成一定、纯度高、无副反应、性质稳定且摩尔质量较大的物质10．配制KMnO4溶液时应选择的方法是(A)称取一定量的固体试剂溶于水中，并用蒸馏水稀释至一定体积，保存于棕色瓶中。

matlab程序设计与应用第二版习题答案matlab程序设计与应用第二版习题答案【篇一：matlab程序设计与应用(第二版)实验答案】%实验一 matlab运算基础%第1题%（1）z1=2*sin(85*pi/180)/(1+exp(2))%（2）x=[2,1+2i;-0.45,5];z2=0.5*log(x+sqrt(1+x.^2))%(3)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)%(4)t=0:0.5:2.5;z4=t.^2.*(t=0t1)+(t.^2-1).*(t=1t2)+(t.^2-2*t+1).*(t=2t3)%第2题a=[12 34 -4;34 7 87;3 65 7];b=[1 3 -1;2 0 3;3 -2 7];a+6*ba-b+eye(size(a))a*ba.*ba^3a.^3a/bb\a[a,b][a([1,3],:);b^2]%第3题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 25] b=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]c=a*bf=size(c)d=c(f(1)-2:f(1),f(2)-1:f(2))whos%第4题%(1):a=100:999;b=rem(a,21);c=length(find(b==0))%(2):a=lsdhksdlkklsdkl;k=find(a=aa=z);a(k)=[]%实验二 matlab矩阵分析与处理 %第1题e=eye(3);r=rand(3,2);o=zeros(2,3);s=diag([2,3]);a=[e,r;o,s];a^2b=[e,(r+r*s);o,s^2]%第2题h=hilb(5)p=pascal(5)hh=det(h)hp=det(p)th=cond(h)tp=cond(p)%第3题a=fix(10*rand(5))h=det(a)trace=trace(a)rank=rank(a)norm=norm(a)%第4题a=[-29,6,18;20,5,12;-8,8,5][v,d]=eig(a)%数学意义略%第5题方法一%(1):a=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; b=[0.95,0.67,0.52]; x=inv(a)*b%(2):b=[0.95,0.67,0.53];x=inv(a)*b%(3):cond(a)%第5题方法二a=hilb(4)a(:,1)=[]a(4,:)=[]b=[0.95,0.67,0.52];x=inv(a)*bb1=[0.95,0.67,0.53];x1=inv(a)*b1n=cond(b)n1=cond(b1)na=cond(a) %矩阵a为病态矩阵%第6题a=[1,4,9;16,25,36;49,64,81]b=sqrtm(a)c=sqrt(a) %sqrtm函数是以矩阵为单位进行计算，sqrt函数是以矩阵中的元素进行计算%实验三选择程序结构设计%第1题程序一x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];y=[]; %建立存放所有y值的矩阵for x0=xif x00x0~=-3y=[y,x0*x0+x0-6];elseif x0=0x05x0~=2x0~=3y=[y,x0*x0-5*x0+6];elsey=[y,x0*x0-x0-1];endendx%输出所有xy%输出所有y%第1题程序二x=[-5,-3,1,2,2.5,3,5];y=[];for a=1:7if x(a)0x(a)~=-3y=[y,(x(a))^2+x(a)-6];elseif x(a)=0x(a)5x(a)~=2x(a)~=3y=[y,(x(a))^2-5*x(a)+6];elsey=[y,x(a)*x(a)-x(a)-1];endend%第2题程序一x=input(请输入一个百分制成绩：);if x100|x0disp(您输入的成绩不是百分制成绩，请重新输入。

数学实验答案Chapter 1Page20,ex1(5) 等于[exp(1),exp(2);exp(3),exp(4)](7) 3=1*3, 8=2*4(8) a为各列最小值，b为最小值所在的行号(10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture(11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10)(12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10)Page20, ex2(1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b(2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码Page20,ex3>> r=2;p=0.5;n=12;>> T=log(r)/n/log(1+0.01*p)Page20,ex4>> x=-2:0.05:2;f=x.^4-2.^x;>> [fmin,min_index]=min(f)最小值最小值点编址>> x(min_index)ans =0.6500 最小值点>> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点f1 =0.0328x1_index =24>> x(x1_index)ans =-0.8500>> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点>> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点f2 =0.0630x2_index =65>> x(x2_index)ans =1.2500>> z=magic(10)z =92 99 1 8 15 67 74 51 58 4098 80 7 14 16 73 55 57 64 414 81 88 20 22 54 56 63 70 4785 87 19 21 3 60 62 69 71 2886 93 25 2 9 61 68 75 52 3417 24 76 83 90 42 49 26 33 6523 5 82 89 91 48 30 32 39 6679 6 13 95 97 29 31 38 45 7210 12 94 96 78 35 37 44 46 5311 18 100 77 84 36 43 50 27 59>> sum(z)>> sum(diag(z))>> z(:,2)/sqrt(3)>> z(8,:)=z(8,:)+z(3,:)Chapter 2Page 45 ex1先在编辑器窗口写下列M函数，保存为eg2_1.m function [xbar,s]=ex2_1(x)n=length(x);xbar=sum(x)/n;s=sqrt((sum(x.^2)-n*xbar^2)/(n-1));例如>>x=[81 70 65 51 76 66 90 87 61 77];>>[xbar,s]=ex2_1(x)Page 45 ex2s=log(1);n=0;while s<=100n=n+1;s=s+log(1+n);endm=nPage 40 ex3clear;F(1)=1;F(2)=1;k=2;x=0;e=1e-8; a=(1+sqrt(5))/2;while abs(x-a)>ek=k+1;F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1);enda,x,k计算至k=21可满足精度clear;tic;s=0;for i=1:1000000s=s+sqrt(3)/2^i;ends,toctic;s=0;i=1;while i<=1000000s=s+sqrt(3)/2^i;i=i+1;ends,toctic;s=0;i=1:1000000;s=sqrt(3)*sum(1./2.^i);s,tocPage 45 ex5t=0:24;c=[15 14 14 14 14 15 16 18 20 22 23 25 28 ...31 32 31 29 27 25 24 22 20 18 17 16];plot(t,c)Page 45 ex6(1)x=-2:0.1:2;y=x.^2.*sin(x.^2-x-2);plot(x,y)y=inline('x^2*sin(x^2-x-2)');fplot(y,[-2 2]) (2)参数方法t=linspace(0,2*pi,100);x=2*cos(t);y=3*sin(t); plot(x,y)(3)x=-3:0.1:3;y=x;[x,y]=meshgrid(x,y);z=x.^2+y.^2;surf(x,y,z)(4)x=-3:0.1:3;y=-3:0.1:13;[x,y]=meshgrid(x,y);z=x.^4+3*x.^2+y.^2-2*x-2*y-2*x.^2.*y+6;surf(x,y,z)(5)t=0:0.01:2*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z)(6)theta=linspace(0,2*pi,50);fai=linspace(0,pi/2,20); [theta,fai]=meshgrid(theta,fai);x=2*sin(fai).*cos(theta);y=2*sin(fai).*sin(theta);z=2*cos(fai);surf(x,y,z)(7)x=linspace(0,pi,100);y1=sin(x);y2=sin(x).*sin(10*x);y3=-sin(x);plot(x,y1,x,y2,x,y3)page45, ex7x=-1.5:0.05:1.5;y=1.1*(x>1.1)+x.*(x<=1.1).*(x>=-1.1)-1.1*(x<-1.1);plot(x,y)page45,ex9clear;close;x=-2:0.1:2;y=x;[x,y]=meshgrid(x,y);a=0.5457;b=0.7575;p=a*exp(-0.75*y.^2-3.75*x.^2-1.5*x).*(x+y>1);p=p+b*exp(-y.^2-6*x.^2).*(x+y>-1).*(x+y<=1);p=p+a*exp(-0.75*y.^2-3.75*x.^2+1.5*x).*(x+y<=-1);mesh(x,y,p)page45, ex10lookfor lyapunovhelp lyap>> A=[1 2 3;4 5 6;7 8 0];C=[2 -5 -22;-5 -24 -56;-22 -56 -16];>> X=lyap(A,C)X =1.0000 -1.0000 -0.0000-1.0000 2.0000 1.0000-0.0000 1.0000 7.0000Chapter 3Page65 Ex1>> a=[1,2,3];b=[2,4,3];a./b,a.\b,a/b,a\bans =0.5000 0.5000 1.0000ans =2 2 1ans =0.6552 一元方程组x[2,4,3]=[1,2,3]的近似解ans =0 0 00 0 00.6667 1.3333 1.0000矩阵方程[1,2,3][x11,x12,x13;x21,x22,x23;x31,x32,x33]=[2,4,3]的特解Page65 Ex 2(1)>> A=[4 1 -1;3 2 -6;1 -5 3];b=[9;-2;1];>> rank(A), rank([A,b]) [A,b]为增广矩阵ans =3ans =3 可见方程组唯一解>> x=A\bx =2.38301.48942.0213(2)>> A=[4 -3 3;3 2 -6;1 -5 3];b=[-1;-2;1];>> rank(A), rank([A,b])ans =3ans =3 可见方程组唯一解>> x=A\bx =-0.4706-0.2941(3)>> A=[4 1;3 2;1 -5];b=[1;1;1];>> rank(A), rank([A,b])ans =2ans =3 可见方程组无解>> x=A\bx =0.3311-0.1219 最小二乘近似解(4)>> a=[2,1,-1,1;1,2,1,-1;1,1,2,1];b=[1 2 3]';%注意b的写法>> rank(a),rank([a,b])ans =3ans =3 rank(a)==rank([a,b])<4说明有无穷多解>> a\bans =110 一个特解Page65 Ex3>> a=[2,1,-1,1;1,2,1,-1;1,1,2,1];b=[1,2,3]';>> x=null(a),x0=a\bx =-0.62550.6255-0.20850.4170x0 =11通解kx+x0Page65 Ex 4>> x0=[0.2 0.8]';a=[0.99 0.05;0.01 0.95];>> x1=a*x, x2=a^2*x, x10=a^10*x>> x=x0;for i=1:1000,x=a*x;end,xx =0.83330.1667>> x0=[0.8 0.2]';>> x=x0;for i=1:1000,x=a*x;end,xx =0.83330.1667>> [v,e]=eig(a)v =0.9806 -0.70710.1961 0.7071e =1.0000 00 0.9400>> v(:,1)./xans =1.17671.1767 成比例，说明x是最大特征值对应的特征向量Page65 Ex5用到公式(3.11)(3.12)>> B=[6,2,1;2.25,1,0.2;3,0.2,1.8];x=[25 5 20]'; >> C=B/diag(x)C =0.2400 0.4000 0.05000.0900 0.2000 0.01000.1200 0.0400 0.0900>> A=eye(3,3)-CA =0.7600 -0.4000 -0.0500-0.0900 0.8000 -0.0100-0.1200 -0.0400 0.9100>> D=[17 17 17]';x=A\Dx =37.569625.786224.7690Page65 Ex 6(1)>> a=[4 1 -1;3 2 -6;1 -5 3];det(a),inv(a),[v,d]=eig(a) ans =-94ans =0.2553 -0.0213 0.04260.1596 -0.1383 -0.22340.1809 -0.2234 -0.0532v =0.0185 -0.9009 -0.3066-0.7693 -0.1240 -0.7248-0.6386 -0.4158 0.6170d =-3.0527 0 00 3.6760 00 0 8.3766(2)>> a=[1 1 -1;0 2 -1;-1 2 0];det(a),inv(a),[v,d]=eig(a) ans =1ans =2.0000 -2.0000 1.00001.0000 -1.0000 1.00002.0000 -3.0000 2.0000v =-0.5773 0.5774 + 0.0000i 0.5774 - 0.0000i-0.5773 0.5774 0.5774-0.5774 0.5773 - 0.0000i 0.5773 + 0.0000id =1.0000 0 00 1.0000 + 0.0000i 00 0 1.0000 - 0.0000i(3)>> A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10]A =5 76 57 10 8 76 8 10 95 7 9 10>> det(A),inv(A), [v,d]=eig(A)ans =1ans =68.0000 -41.0000 -17.0000 10.0000-41.0000 25.0000 10.0000 -6.0000-17.0000 10.0000 5.0000 -3.000010.0000 -6.0000 -3.0000 2.0000v =0.8304 0.0933 0.3963 0.3803-0.5016 -0.3017 0.6149 0.5286-0.2086 0.7603 -0.2716 0.55200.1237 -0.5676 -0.6254 0.5209d =0.0102 0 0 00 0.8431 0 00 0 3.8581 00 0 0 30.2887(4)(以n=5为例)方法一（三个for）n=5;for i=1:n, a(i,i)=5;endfor i=1:(n-1),a(i,i+1)=6;endfor i=1:(n-1),a(i+1,i)=1;enda方法二（一个for）n=5;a=zeros(n,n);a(1,1:2)=[5 6];for i=2:(n-1),a(i,[i-1,i,i+1])=[1 5 6];enda(n,[n-1 n])=[1 5];a方法三（不用for）n=5;a=diag(5*ones(n,1));b=diag(6*ones(n-1,1));c=diag(ones(n-1,1));a=a+[zeros(n-1,1),b;zeros(1,n)]+[zeros(1,n);c,zeros(n-1,1)] 下列计算>> det(a)ans =665>> inv(a)ans =0.3173 -0.5865 1.0286 -1.6241 1.9489-0.0977 0.4887 -0.8571 1.3534 -1.62410.0286 -0.1429 0.5429 -0.8571 1.0286-0.0075 0.0376 -0.1429 0.4887 -0.58650.0015 -0.0075 0.0286 -0.0977 0.3173>> [v,d]=eig(a)v =-0.7843 -0.7843 -0.9237 0.9860 -0.92370.5546 -0.5546 -0.3771 -0.0000 0.3771-0.2614 -0.2614 0.0000 -0.1643 0.00000.0924 -0.0924 0.0628 -0.0000 -0.0628-0.0218 -0.0218 0.0257 0.0274 0.0257d =0.7574 0 0 0 00 9.2426 0 0 00 0 7.4495 0 00 0 0 5.0000 00 0 0 0 2.5505Page65 Ex 7(1)>> a=[4 1 -1;3 2 -6;1 -5 3];[v,d]=eig(a)v =0.0185 -0.9009 -0.3066-0.7693 -0.1240 -0.7248-0.6386 -0.4158 0.6170d =-3.0527 0 00 3.6760 00 0 8.3766>> det(v)ans =-0.9255 %v行列式正常, 特征向量线性相关，可对角化>> inv(v)*a*v 验算ans =-3.0527 0.0000 -0.00000.0000 3.6760 -0.0000-0.0000 -0.0000 8.3766>> [v2,d2]=jordan(a) 也可用jordanv2 =0.0798 0.0076 0.91270.1886 -0.3141 0.1256-0.1605 -0.2607 0.4213 特征向量不同d2 =8.3766 0 00 -3.0527 - 0.0000i 00 0 3.6760 + 0.0000i>> v2\a*v2ans =8.3766 0 0.00000.0000 -3.0527 0.00000.0000 0.0000 3.6760>> v(:,1)./v2(:,2) 对应相同特征值的特征向量成比例ans =2.44912.44912.4491(2)>> a=[1 1 -1;0 2 -1;-1 2 0];[v,d]=eig(a)v =-0.5773 0.5774 + 0.0000i 0.5774 - 0.0000i-0.5773 0.5774 0.5774-0.5774 0.5773 - 0.0000i 0.5773 + 0.0000id =1.0000 0 00 1.0000 + 0.0000i 00 0 1.0000 - 0.0000i>> det(v)ans =-5.0566e-028 -5.1918e-017i v的行列式接近0, 特征向量线性相关，不可对角化>> [v,d]=jordan(a)v =1 0 11 0 01 -1 0d =1 1 00 1 10 0 1 jordan标准形不是对角的，所以不可对角化(3)>> A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10]A =5 76 57 10 8 76 8 10 95 7 9 10>> [v,d]=eig(A)v =0.8304 0.0933 0.3963 0.3803-0.5016 -0.3017 0.6149 0.5286-0.2086 0.7603 -0.2716 0.55200.1237 -0.5676 -0.6254 0.5209d =0.0102 0 0 00 0.8431 0 00 0 3.8581 00 0 0 30.2887>> inv(v)*A*vans =0.0102 0.0000 -0.0000 0.00000.0000 0.8431 -0.0000 -0.0000-0.0000 0.0000 3.8581 -0.0000-0.0000 -0.0000 0 30.2887本题用jordan不行, 原因未知(4)参考6(4)和7(1)Page65 Exercise 8只有(3)对称, 且特征值全部大于零, 所以是正定矩阵. Page65 Exercise 9(1)>> a=[4 -3 1 3;2 -1 3 5;1 -1 -1 -1;3 -2 3 4;7 -6 -7 0]>> rank(a)ans =3>> rank(a(1:3,:))ans =2>> rank(a([1 2 4],:)) 1,2,4行为最大无关组ans =3>> b=a([1 2 4],:)';c=a([3 5],:)';>> b\c 线性表示的系数ans =0.5000 5.0000-0.5000 1.00000 -5.0000Page65 Exercise 10>> a=[1 -2 2;-2 -2 4;2 4 -2]>> [v,d]=eig(a)v =0.3333 0.9339 -0.12930.6667 -0.3304 -0.6681-0.6667 0.1365 -0.7327d =-7.0000 0 00 2.0000 00 0 2.0000>> v'*vans =1.0000 0.0000 0.00000.0000 1.0000 00.0000 0 1.0000 v确实是正交矩阵Page65 Exercise 11设经过6个电阻的电流分别为i1, ..., i6. 列方程组如下20-2i1=a; 5-3i2=c; a-3i3=c; a-4i4=b; c-5i5=b; b-3i6=0;i1=i3+i4;i5=i2+i3;i6=i4+i5;计算如下>> A=[1 0 0 2 0 0 0 0 0;0 0 1 0 3 0 0 0 0;1 0 -1 0 0 -3 0 0 0; 1 -1 0 0 0 0 -4 0 0;0 -1 1 0 0 0 0 -5 0;0 1 0 0 0 0 0 0 -3; 0 0 0 1 0 -1 -1 0 0;0 0 0 0 -1 -1 0 1 0;0 0 0 0 0 0 -1 -1 1];>>b=[20 5 0 0 0 0 0 0 0]'; A\bans =13.34536.44018.54203.3274-1.18071.60111.72630.42042.1467Page65 Exercise 12>> A=[1 2 3;4 5 6;7 8 0];>> left=sum(eig(A)), right=sum(trace(A))left =6.0000right =6>> left=prod(eig(A)), right=det(A) 原题有错, (-1)^n应删去left =27.0000right =27>> fA=(A-p(1)*eye(3,3))*(A-p(2)*eye(3,3))*(A-p(3)*eye(3,3)) fA =1.0e-012 *0.0853 0.1421 0.02840.1421 0.1421 0-0.0568 -0.1137 0.1705>> norm(fA) f(A)范数接近0ans =2.9536e-013Chapter 4Page84 Exercise 1(1)roots([1 1 1])(2)roots([3 0 -4 0 2 -1])(3)p=zeros(1,24);p([1 17 18 22])=[5 -6 8 -5];roots(p)(4)p1=[2 3];p2=conv(p1, p1);p3=conv(p1, p2);p3(end)=p3(end)-4; %原p3最后一个分量-4roots(p3)Page84 Exercise 2fun=inline('x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x');fzero(fun,2)Page84 Exercise 3fun=inline('x^4-2^x');fplot(fun,[-2 2]);grid on;fzero(fun,-1),fzero(fun,1),fminbnd(fun,0.5,1.5)Page84 Exercise 4fun=inline('x*sin(1/x)','x');fplot(fun, [-0.1 0.1]);x=zeros(1,10);for i=1:10, x(i)=fzero(fun,(i-0.5)*0.01);end;x=[x,-x]Page84 Exercise 5fun=inline('[9*x(1)^2+36*x(2)^2+4*x(3)^2-36;x(1)^2-2*x(2)^2-20*x(3);16*x(1)-x(1)^3-2*x(2)^ 2-16*x(3)^2]','x');[a,b,c]=fsolve(fun,[0 0 0])Page84 Exercise 6fun=@(x)[x(1)-0.7*sin(x(1))-0.2*cos(x(2)),x(2)-0.7*cos(x(1))+0.2*sin(x(2))];[a,b,c]=fsolve(fun,[0.5 0.5])Page84 Exercise 7clear; close; t=0:pi/100:2*pi;x1=2+sqrt(5)*cos(t); y1=3-2*x1+sqrt(5)*sin(t);x2=3+sqrt(2)*cos(t); y2=6*sin(t);plot(x1,y1,x2,y2); grid on; 作图发现4个解的大致位置，然后分别求解y1=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[1.5,2])y2=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[1.8,-2])y3=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[3.5,-5])y4=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[4,-4])Page84 Exercise 8(1)clear;fun=inline('x.^2.*sin(x.^2-x-2)');fplot(fun,[-2 2]);grid on; 作图观察x(1)=-2;x(3)=fminbnd(fun,-1,-0.5);x(5)=fminbnd(fun,1,2);fun2=inline('-x.^2.*sin(x.^2-x-2)');x(2)=fminbnd(fun2,-2,-1);x(4)=fminbnd(fun2,-0.5,0.5);x(6)=2feval(fun,x)答案: 以上x(1)(3)(5)是局部极小，x(2)(4)(6)是局部极大，从最后一句知道x(1)全局最小，x(2)最大。

matelab作业2参考答案Matlab作业2参考答案Matlab作业2是一项综合性的任务，要求学生运用Matlab编程语言解决一系列数学问题。

本文将为大家提供一份参考答案，帮助学生更好地理解和完成这项作业。

首先，我们将讨论作业的第一个问题，即给定一个矩阵A，求解其特征值和特征向量。

在Matlab中，可以使用eig函数来实现这一功能。

例如，假设我们有一个3×3的矩阵A，可以按照以下方式计算其特征值和特征向量：```A = [1 2 3; 4 5 6; 7 8 9];[eigenvectors, eigenvalues] = eig(A);```在上述代码中，变量eigenvectors将存储A的特征向量，而变量eigenvalues 将存储A的特征值。

通过打印这两个变量的值，我们可以得到矩阵A的特征值和特征向量。

接下来，我们将探讨作业的第二个问题，即求解线性方程组。

假设我们有一个3×3的系数矩阵A和一个3×1的常数向量b，我们需要求解方程组Ax=b。

在Matlab中，可以使用backslash运算符来求解线性方程组。

例如，假设我们有以下方程组：```A = [1 2 3; 4 5 6; 7 8 9];b = [10; 20; 30];x = A \ b;```在上述代码中，变量x将存储方程组的解。

通过打印变量x的值，我们可以得到方程组的解。

此外，作业的第三个问题要求学生使用Matlab绘制函数图像。

在Matlab中，可以使用plot函数来实现这一功能。

例如，假设我们要绘制函数y=sin(x)，其中x的取值范围为0到2π，可以按照以下方式绘制函数图像：```x = 0:0.1:2*pi;y = sin(x);plot(x, y);```在上述代码中，变量x将存储x的取值范围，变量y将存储对应的函数值。

通过调用plot函数，我们可以将函数y=sin(x)的图像绘制出来。

实验一:T1:％％第一小题z1 =2*sin( 85＊pi/1 80) /(1+exp(2 )) ％％第二小题x=[2 ,1+2i;－ 0、45, 5];z2=1/2*log(x+s qr t(1+x、人2));z 2 %%第三小题a =—3、0:0、1:3、0;z3=1 /2 * (ex p (0、3*a)- e xp (-0、3 *a)卜＊sin(a+0、3)＋log( (0、3 +a) /2)％%第四题9, 7, 0; 4,1 3,11] ; C =A* BD =C(3:5, 2 : 3)T4-1：a =100: 9 99;b=fi nd( r e m (a,21) ==0 );c=1 ength(b )T4—2:a = in p Ute请输入一个字符串:'，‘ sO ;b = find (a〉= A/ &a〈='ZO ;a( b)=[];disp( a) ;实验二:t =0 :0、5 :2、5z 4=( t>=0 & t〈 1 )、*( t、人2) +(t>= 1&t<2)、* (t、人2 —1) +(t〉=2&t<3 )、* (t、人2—2 * t+ 1)T2:A =[12 , 3 4,— 4;34, 7, 8 7; 3 ,65,7]B = [ 1 , 3,— 1；2, 0, 3;3, -2,7]disp ( 'A +6* B=') ;di sp( A+6* B);disp( A'—B＋ I=');dis p(A-B+eye(3));dis p('A*B= ');dis p(A*B);d isp( 'A、*B =');disp( A、*B) ;d isp ('AG ';disp (A 人3);disp ('A、人3=');di sp (A、人3)；d isp('A/B=') ；di sp(A/B) ；disp ( '\AB=' ) ;dis p (B\A)；disp(/[A,B]=')；di sp([A,B] );d i s p ( ' : A([1,3 ]，: );B^2]=');d i sp ([A (: 1,3],:)砂2]);T3:z=1 :25；A=res hape (z, 5, 5) ';B= [3, 0, 16; 17,—6 ,9;0 ,23, —4;T1 :E=eye(3), R=rand(3, 2) ,O=zero s(2,3), S= dia g ([1, 2]);A=[E,R;O,S]disp (A^2=');di s p (A 2);disp('[E ,R +RS O, S^2：');B=[ E, R+ R*S; O,S 人2:T2:H=hi 1 b(5)P=pasc al (5 )H h=de t( H)Hp=de t(P)Th=co nd( H)Tp =co nd (P)a =ab s (Th- 1 );b = abs ( T p - 1);i f a > bdisp ('帕萨卡矩阵P 性能更好') ;e1s e i fa<bdisp (/希尔伯特矩阵H 性能更好' )；elsedis P C两个矩阵性能相同');en dT3:a=1: 25;A=res ha pe(a,5,5)dis p('行列式得值：');disp(de t (A));dis p( '矩阵得秩:')；disp( rank(A ))；d i s p(' 矩阵得迹 :' )； disp(tr ace( A)); disp ( '矩阵得范数: )); dis p(nO rm(A) ); T4:A=[－2 9, 6, 18;20,5 ,12 ;——8,8, 5] [V,D] =eig (A) T5:A=[ 1/2, 1/3, 1/ 4 ;1/3, 1/ 4, 1/5; 1/4, 5， 1／6]B= : 0、9 5,0、6 7, 0、5 2 ]' X1=A \B B ( 3) =0、 53 X2=A\Bdis p ('系数矩阵 A 得条件数：/); disp(con d(A)) ;T6:disp('the scOr e is not reas ona b le));s=i n p ut( / pie as e en tercore: /) ; endi f (s>= 90&S < = 10 0);d isp( 'A');el s eif(s > = 8 0 & s <90);d isp()B/) ;e1 se if(s 〉 =70&sV 80)；di sp( ) C));elseif (s> = 60 & sV 70);dis p( 'D'); el s edis p('E') ;endT 2－ switch :a=1:25; A =r e sha pe ( a, 5,5 ) di sp( ' B1=sqrtm)('A)； B 1=sq r t m (A) disp( B2 =sq rt( A ) ‘); B2= sq r t(A) di sp( 'B1*B 1'); B1* B1 disp (/ B2、 *B2) ; B2、 *B2 实验三 : T1: X =——5、 0:2: 5、0 : length(X) ; (x (i) <0) &(x(i) ~= -3) yel s ei & (x(i)〜=2) y(i ) else y (i)=x(i)人2 —x( i)-1 ;endfor i=1 if (i) = x( i)八 2 +x (i) —6； f (x(i)>=0)&(x(i)V5)& (x (i )〜=3);=x(i)A2 — 5* x(i)+6; end s=input ('p lease en ter ')；whi 1 e (s < 0| | s 〉disp ( 'the reasonab1 e' );s=input('p 1 ease re:' ); end switc h fix(s/1 0 ) cas e{9, 1 0}disp(＇ A')； case{8} di s p CE '); c ase{ 7} d isp ('C ' case {6} dis p(' D')； ca se{ 0, d isp(' E/ e nd T 3: t=inp ut i ft>120 SCO re : );100) scor en 1, 2,3, 4,5})；'请输入工时 s nO tt h e sc o me=');y T2－ if: s =i npu t(' pie a se ent er the s c ore:，);w=120 *84+(t-120) * 1、 15* 84;elseif t < 60 w =t * 8 4 -701/the senden delse w= 8 4*t; pii=sq T1—2: rt(6 ＊m)end d isp ('应发工资为：’)； di sp(w); T 4: a=10 +f1 oor(rand(1)* 89) b =1 0 + f loor( ra nd (1) * 89) s=inp ut (/ 请输入 + 或一或 * 或 /' , ' s'); w hi 1 e ( s 〜='+ ' & s 〜-—，&& s 〜-* '&&s~ =' /') d isp ('输入得符号错误, 请重新输入 ') ; s=inp ut (/ 请输入 + 或-或*或/' , 's'); end switch s case{ ' +' c=a + b;c ase{ '—'} c=a —b ; ca se{＇* ' c=a* b; ca se{ ' / '} c =a/b ; t(/请输入 n : (1 : n)、人 2; s=sum( m )； pii=sq rt(6 ＊s) T2:n=inpum=1 、 / n=1 0 00;y=0; for i=1 : n;y = y+ 1/ (2 * i —1if end T 3:endc T5:A = n=inp u t('请输入 n :'); whi 1 e (n <1) d isp ('输入得n 有误，请重新输入/ ); n=i npu t ('请输入 n:'); rand(5,6) if n >5 B=A ( 5， :) eIse B = A(n,:) E nd 实验四 :T1: n=in p ut (/ 请输入 n ：'); m =0; for i=1 ： n; m=m+1/iA 2;');）；y >= 3d isp ( ' 最大 n i — 1disp ('对应得 值为：＇）;y 值为 :' )； y －1/( 2＊i-1) break;dt a =in put('请输入 a:'); b =inp ut ('请输入 b:');x (1)= a/(1、0 +b) ; i = 1;x (2 ＋b) ; whil e (a b s (x (i+ 1) — x(i))> 1& i 〈 50 0) for m Ion g i =i+1 ;x(i+1)= a/(x (i) + b); endr 仁(—b+s qrt(b 人 2+4*a)) /2; r2 =( -b — s qr t( b^ 2+ 4 *a)) )=a / (x(i) 0、00 0 0d is p ([' x ' ,nm 2st r (i+ 1 ), n um2s t r(x( i+1)， 8)] );disp (['x',num 2 str( i)，'=＇,num 2str (x( i),8 )]);di sp(['r1=＇,nu m2s tr( r1,8)，'， r2='， num2str(r 2,8 )]); fo r mat short T4:f(1 )=1;f( 2)=0;f(3 )=1 ;i=4; w h il e (i<= 10 0)f(i) =f( i —1 )-2* f( i —2) +f( i —3 ); i=i+1;f m a x= max (f )fmin = min (f) fs um=s um(f) fp= 1 ength(f in d( f> 0)) fn=length (fin d (fv0 )) fO=lengt h( f i nd (f= =0)) T5:j =0; s= 0;for i=3：50m=i*( i-1)—1;i f i s prime (m)= =1s = s +m ；j =j +l;endendd i sp(['亲密数得对数：',nu m 2s tr (j)：); dis p ([＇所有亲密素数之与：'，n um2str( s )] );实验五：f n1: func t ion f = fn1 ( n)f= n +10 * l o g( n、人2 +5); fn2：f un ctio n f=f n2( n) f = 0；f or i= l： n nction f=fx ( x)/ ((x— 2)、八 2 + 0、1) +1、/ ((x—3 )、4+0、0 1);Tl :x =in p ut('请输入一个复数：');fe =exp (x)f lg= log( x) fsin=sin( x) f co s=cos(x)T2:m 1=in p ut ('请输入m 1/);m2=input ()请输入m 2:'; t=input ('请输入th et a：));A=[m1*cos(t )， -m1,-si n(t), 0;m 1 * si n(t), 0 ,、、、co s( t) , 0 ;0,m2, —sin (t), 0; 0, 0, —C o s(t) ,1]B= [0， m1 * 9、8, 0,m2dis p('[ al ,a2 ,N 1,N2： d i sp([A\B :';T 3: j=1 ； n= l;k =1;for i=10：99 中所有得素数if isprime(i)==1m(j)=i;j=j＋1;endendf o r t=1:le n gth( m);s( n)=l 0* r em(m(t )， 10)＋fix(m (t )/10); %挑选出得素数进行位置交换if i sp rime(s(n)) = = 1%再挑选出交换绝对素数r( k)= m(t);k=k+l ;endn=n+ l ;enddisp (['所有两位绝对素数:'，num2str (r) ])T 5：y 1=f n 1 ( 4 0 ) /(f n 1(3 0 ) +fn1( 2 0 )) %函数调用y2 =f n2(40)／( fn 2( 30)＋fn2 (20))实验六:T1: x=linspa ce (0， 2*pi,101 )；y=(0、5+ 3* sin(x)、/ (l + x、人2))、* c os ( x) ;plot( x， y) ;T2： x=1in sp ace(0， 2*p i,10 0);y1=x、人2;y2=cos( 2*x);y3 =y1 、*y2 ;subplot(2，2， 1);plot(x,y1,r':',x,y2， 'y－'， x， y3， 'b －、en f=f +i * (i + 1); dfx：fu f=1 、*9、8]'=/ );%挑选出1 0 ~99sub pl o t(2, 2,2) ; el se pl ot(X,y 1,'r:')； y =[y,1 /2 *l og(x0＋sqr t(1 s u b plot ( 2,2,3)； + x 0 人 2))]; p lot(x,y2, ' -y')； e nd subplot(2 ,2, 4); en d pl ot( x,y 3, 'b —、 ')； subplot (1, 2,1); T 2-3 : pl ot(X, y)； x=lin spa ce(0 ,2*pi,1 0 ); fo r i =1:l ength( x) y 仁X 、人2 ; if X(i) <=0y2 =cos(2＊ x); y(i)= (x (i)+sq rt( pi )) / e x p y 3 = y1、*y2 ; (2 ); sub plot (3, 4,1)；else b ar( x, y1,＇r ＇)；axis( [0, 7 ,0, 40：)； y( i) =1／2 *l og(x (i)＋sqr s ub pl ot(3, 4 ,2)； t (1+x(i) A2)); stairs(x, y1,'y ＇) ;aX is( [0,7,0, 40])； end subp lo t(3, 4, 3)；e nd stem(x, y1,b' ') ;aX is([0,7 ,0 ,40]); subp lot(1,2, 2)； '); y =[y,(xO+sqrt(pi)) /exp( 2)]; su bp l 0 t (3,4, 4);fill( x,y1, 'g ') ;axis ([0, 7,0 ,40]) ; sub plot( 3, 4,5);bar(x, y2,＇r' );ax(i[s —1,7,—1, 1：) subplot( 3, 4, s tair s(X ,y2, ' 1：);subpl ot (3, 4,7);ste m( x, y2, 'b )'1： ); sub pl ot( 3,4,8);fil l(X,y 2,'g ＇); sub p lot (3 ,4, 9); bar (x,y3,＇r') ;aXi s([0,7 ,-30,50]) ; subplot (3, 4,10 );sta irs(x,y3, ＇)' ;yaxi s([0, 7,-30 ,50：)； su bplot(3 ,4, 11);s tem( x,y3 ,'b') ;axi sub pl o t(3,4 ,12); fill(x,y3,＇g ＇)； 50])；T3: 6)； y)；'axis([—1, 7,-1, ；a xis (：-l ,7, - 1, axi s (:-1,7,-1 ,1：); s ([0,7, -30,5 0：);a x is ( [ 0, 7 ,-30,x=li nsp ace( -5,5,1 y=[：； for X0 =x ； i f x0〈 =0；0 0)； plo t (x, y )；T4：a=i n pu t C 请输入 a:'); b=i np ut ('请输入 b: / ); n= i npu t ('请输入 n :');t he ta=0: 0、0 0 1 :2*p i; y=a 、 *sin(b+n 、 *theta) ; polar(thet a,y); T5： X =lin s pace( —5 ,5, 2 1); y= l in s pace(0,10,3 1); :x, y ： =me s hg rid (x,y); z =cos(x)、*c os (y)、* exp( — sq r t ( x 、 人2 + y 、人 2) /4) s u bp lot (1, surf( x,y,z): sub p l o t(1, co ntour3(x,y, T6：s=0:pi /2; t=0：3＊pi/2; [s,t ： =mes hgri d(s,t) ;x= c os(s)、* c os (t );y= c os ( s)、* s i n (t); z=sin( s) ; s urf(X ,y, z);2, 1）； 2, 2);z);s h a d in g interp ;实验八：T1:?A=rand(1, 30000);dis p (＇均值为:') ;mean (A)di sp('标准方差为：/);s td ( A)di s p('最大元素为:');ma x( A)d i s P C最小元素为：');min(A)k=find ( A〉0、5) ; b= le ngd h( k)；disp ('大于0、5 得随机数个数占总数得百分比:')；sp rintf ('% f % % , b / 3 0 0 )T2: ?P =45＋50* ran d(100, 5) disp('最高分及序号：/ );[Y,U：=max(P,[：,1) disp ('最低分及序号： ' );[Y1, U1：= min(P,[：,1)di sp ('每门课平均分：') ;mean( P,1)di sp( '每门课得标准方差：');std(P ,0, 1)s=sum(P,2) ;d i s p( '总分最高分： ' );[Y2, U2： =maX( s, [： ,1) di sp( '总分最低分：＇);[Y3, U3]=min(s,[], 1)[zcj , X sxh：end') % desc T3: =sord ( s ,1 ,' descen d降序p=p olyfit ( x, y, 5) x1=1：0、01:101;y1=log10(x1);y 2 =poly V al (p, x1);plot( X 1, y1, ' r -' ,x1 y2 ,'£';)T 5:P1=[1,2,4,0 ,5];P2=[0 ,0, 0, 1,2]；P3 =[0,0 , 1,2 , 3：;P4= c onv( P2 , P 3);sp 1 = length( P 1);sp2=len gth( P4);P =[ z er os (1 ,sp2 —s p1 ), P1： + P4x= roots(P)A=[-1, 1、2, -1、4;0、75,2,3、5; 0,5 ,2、5];Ax=p o l y val (PA)A A=po l y valm(P, A)实验九:T1 :f= i nline( ed * X 人2;0 26*x ：i=1 ;x=1;whi 1 e x<=3、(:x xA 2 xA 3)1 2* x 3') ;h= 6：2：18;t1=[1 8、0,2 0、28、0, 24、0];t 2 = [1 5、0,1 9、32、0, 30、t=6、5：2：1T1=interp1(h,t T2=inde rpT4: 0，22、0，25、0, 24、0,28、0];7、5;1 ,t, 'sp1 ine ') 1(h, t2,d , 'ps1 0 ,30、0,0，34、0，ine )'x=1 ：10:101;y=log101g( i) =f(X);i =i+1;x=X +0、0 1 ;endd x=di f f (g)/ 0、01;dX (1)dx(1 01 )dx(length(g)—1)T2：f1 =in line (‘s qrt (co s(t、人2)+4*sin(2* t )、人2+1)')f2=in lin e('log(1+x)、／( 1＋X、A2) ')pi )pi)I 1 = quI2=q u adT3 :ad(f1,0，(f 2, 0,2＊2＊A=[6, 5, —2;3,—9,0,2： ;B=[—4, 13,1,1 1] ;'x1=A\ B-2, 5; 9，-1 ,4 , —1；3,4,2 ,:L ,U]=I u( A);( L\B)[Q,R ]=q r (A);(Q\B)f u n:f u n ction F = fu n( X) x=X(1)；y =X( 2);z =X (3);F (1)=sin(X)+ y^2+l o g (z)—7 ;F(2)=3*x+ 2 Ay— z 人3+1 ;F(3) =x+y+z-5；T5:f =inl ine( '3*x+sin(x) —exp(x) ') ;fzero( f,1 、5)X=fso Ive ('让n', : 1 1 1]',o pt i m set CDi sp l ay/, ' f f 0 ))实验十:T1:x= sym('6 ';y =sym(' 5';z=(x+ 1) /(sqrt( 3 +x) —sq r t (y ))T2:Sy m s x y;t =sym ('51350 ；factor(xA4 — yM) factor(t)T 3:s y ms betal b eta2 x;Sim pli f y( sin( b e t a 1) * cos(beta2) —cos (beta1)＊sins i m p li fy(1))T4:d p 1= [ 0，1,0; 1 ，p2=[1,0，0 ;0,1，A = : a,b,c;d,e,f;B= p 1*p 2 * AB1= inv (B)B1* B;tr i l( B ) d et (B)T5:sy ms x t a yf1=(X *( eX p(sin(x) )+1 )—2*(e xp(tan(x))—1) )/sin ( x)A3 ;l imit(f1， x， 0) %( 1)f 2=( sqrt(pi )-sqrt( acos(X )) )/sqrt(x＋1);limi t(f 2, X ,— 1 , ' grih t/ % (2)f3=(1— co s(2*x))/x;di ff(f3， x， 1) %(3)di ff (f 3， x,2)A=[a Ax,tA3 ; t*cos(X)，log(x)： ;dif f(A， x， 1)＆( 4)d iff(A ,t，2 )diff(diff (A，x,1)，f5=( xA2—2*x) * e yx= -d iff(f5,x, 1) (5)fxy=diff( di ffX= sym ('0');y=sym ( ' 1' ;e val (f xy)T 6: x=sym( ' x;')f1=1/(1+x 人4 +x人8 );int(f1， x)f 2 = 1/(asin (x))人2 /sq r t (1 — x^2);i n t ( f 2,x)f3 =(X A2＋1) /(xA4+1);i nt( f 3，X ,0,inf)f4=exp(x)*( 1＋eX p( x) )A2;in t(f4 ,X，0，l og(2))实验十一:sym s n X ;S1 = s ymsu m (1/(2* n —1) , n, 1, 1 0) S 2 = symsum( n人2*乂人(门一1),n ,1,inf)S3=s ym s u m (n人2/5人n,n ,1, inf)T 2: x=s ym('X/);f= I og (x)；bet a2))4*x 人2+8* x+3 )/ (2 * x +syms a b c e f g h k;0,0 ;0 ,0 ,1];0 ;1,0, 1：;t， 1)X p( -xA2 — y 人2 -X*y);f5，x,1)，y， 1)T1：ta ylor (f ,x,6， 1 )T3:x1=s ol ve( '(lo1g+x)—5/(1＋sin(X)) =2') x2 =sol ve( '2xA+9* s qrt (x＋1)—1=0') x3=sol ve(/3* x* eX p( x)＋5* sin( x)—7 8、5 =0 ‘)[x4，y4]=s o lve( S qrt(x^2+y 人2)－ 10 0= 0 ' ' *x+5 * y-8= O’)T4 :d so 1 ve(' D2 + 4 *Dy + 2 9 * y = 0' ' 0 )= 0 '/ Dy(0)=15/，'x') T5 :[x，y ,z]=dsol V e ('Dx=2*x —3*y+3*z/ / Dy=4* x-5 * y+3* z ',' D z=4 *x —4*y+ 2*z'，/ t ')。

实验2 线性规划问题及对偶问题求解实验内容与答案提示：灵敏度分析设置方式：先在lingo菜单options里面设置general solver 的dual computation里面加上ranges然后在lingo菜单里面选range就行了注意lingo只能对线性的模型做灵敏度分析题1 线性规划问题的灵敏度分美佳公司计划制造 I、II 两种家电产品。

已知各制造一件时分别占用设备 A、B 的台时、调试时间、调试工序每天可用于这种家电的能力、各售出一件时的获利情况，如表 1-1 所示。

1.问该公司应制造两种家电各多少件，使其获取的利润最大。

max=2*x1+1*x2;5*x2<=15;6*x1+2*x2<=24;x1+x2<=5;Global optimal solution found.Objective value: 8.500000Infeasibilities: 0.000000Total solver iterations: 2Variable Value Reduced CostX1 3.500000 0.000000X2 1.500000 0.000000Row Slack or Surplus Dual Price1 8.500000 1.0000002 7.500000 0.0000003 0.000000 0.25000004 0.000000 0.50000002. 如果资源出租，资源出租的最低价格至少是多少（即每种资源的影子价格是多少）。

min=15*y1+24*y2+5*y3;6*y2+y3>=2;5*y1+2*y2+y3>=1;Global optimal solution found.Objective value: 8.500000Infeasibilities: 0.000000Total solver iterations: 3Variable Value Reduced CostY1 0.000000 7.500000Y2 0.2500000 0.000000Y3 0.5000000 0.000000Row Slack or Surplus Dual Price1 8.500000 -1.0000002 0.000000 -3.5000003 0.000000 -1.5000003.若家电 I 的利润不变，家电 II 的利润在什么范围内变化时，则该公司的最优生产计划将不发生变化。

Matlab 课后实验题答案实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。

(1) 0122sin851z e=+(2) 21ln(2z x =，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦(3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:: 解：2. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(1) A+6*B和A-B+I（其中I为单位矩阵）(2) A*B和A.*B(3) A^3和A.^3(4) A/B及B\A(5) [A,B]和[A([1,3],:);B^2]解：3. 设有矩阵A 和B1234530166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦ (1) 求它们的乘积C 。

(2) 将矩阵C 的右下角3×2子矩阵赋给D 。

(3) 查看MATLAB 工作空间的使用情况。

解：. 运算结果：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。

(2) 建立一个字符串向量，删除其中的大写字母。

解：(1) 结果：(2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。

matlab2022实验2参考答案报告名称：MATLAB试验二符号计算姓名：学号：专业：班级：MATLAB实验二MATLAB符号计算试验报告说明：1做试验前请先预习，并独立完成试验和试验报告。

2报告解答方式：将MATLAB执行命令和最后运行结果从命令窗口拷贝到每题的题目下面，请将报告解答部分的底纹设置为灰色，以便于批阅。

3在页眉上写清报告名称，学生姓名，学号，专业以及班级。

3报告以Word文档书写。

一目的和要求1熟练掌握MATLAB符号表达式的创建2熟练掌握符号表达式的代数运算3掌握符号表达式的化简和替换4熟练掌握符号微积分5熟练掌握符号方程的求解二试验内容1多项式运算（必做）1.1解方程:f(某)=某^4-10某某^3+34某某^2-50某某+25=0%采用数值方法：>>f=[1-1034-5025];>>root(f)%采用符号计算方法：f1=ym('某^4-10某某^3+34某某^2-50某某+25')olve(f1)1.2求有理分式R=(3某^3+某)(某^3+2)/((某^2+2某-2)(5某^3+2某^2+1))的商多项式和余多项式.a1=[3010];a2=[1002];a=conv(a1,a2);b1=[12-2];b2=[5201];b=conv(b1,b2);[p,r]=deconv(a,b);%注意：ab秩序不可颠倒。

%reidue用于实现多项式的部分分式展开，此处用deconv函数报告名称：MATLAB试验二符号计算姓名：学号：专业：班级：%%此题，有同学程序如下：某1=[3010]，某2=[1002],某3=[12-2],某4=[5201]某5=conv(某1,某2)[y6,r]=deconv(某5,某3)R=deconv(y6,某4)%%这种方法较第一种解法缺点：在除法运算中，会产生误差，故此题应先将分母的多项式相乘后，再与分子部分的多项式进行运算。