Lab 2 参考答案

合集下载

2第二届--全国--大学生化学实验笔-答案-(真的是全的)

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程序设计与应用第二版习题答案【篇一: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(您输入的成绩不是百分制成绩,请重新输入。

(完整版)MATLAB数学实验第二版答案(胡良剑)

(完整版)MATLAB数学实验第二版答案(胡良剑)

数学实验答案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参考答案

matelab作业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)的图像绘制出来。

操作系统Lab2实验报告

操作系统Lab2实验报告

深入学习 操作系统的原理和实 现细节,包括进程通 信、死锁处理、虚拟 化技术等。
实践应用
希望将所学的操作系 统知识应用到实际项 目中,提高自己的实 践能力。
持续学习
随着技术的不断发展 ,操作系统也在不断 更新和进步,需要保 持持续学习的态度。
探索新领域
对于新兴的操作系统 领域,如云操作系统 、物联网操作系统等 ,也希望有所涉猎和 探索。
学生将了解操作系统内核的组成和功能,以及操作系统如何通过内核 实现各种基本功能。
学生将通过实验深入了解操作系统的内部实现细节,提高自己对操作 系统设计和实现的认知水平。
02
实验内容
进程管理实验
总结词
理解进程状态及其转换、进程控制块PCB的组成、进程调度算法。
详细描述
通过实验,我们深入了解了进程的概念、状态及其转换。进程控制块PCB的组成,包括进程标 识符、进程状态、优先级、CPU寄存器内容等信息。此外,我们还学习了进程调度算法,如先 来先服务、最短作业优先、优先级调度等,并进行了模拟实验。
THANKS
感谢观看
操作系统Lab2实验 报告
汇报人:
202X-01-08
目录
• 实验目的 • 实验内容 • 实验结果与分析 • 问题与改进建议 • 总结与展望
01
实验目的
理解操作系统的基本概念
01
掌握核心概念
02
操作系统是计算机系统的核心软件,负责管理计算机硬件和软件资源 ,为用户提供便利的操作界面和高效的服务。
03
操作系统具有并发性、虚拟性、异步性和随机性的特点,这些特点决 定了操作系统在计算机系统中的重要地位和作用。
04
操作系统的基本功能包括进程管理、内存管理、文件管理和设备管理 ,这些功能是操作系统实现其核心目标的基础。

Matlab程序设计与应用第二版刘卫国课后实验答案

Matlab程序设计与应用第二版刘卫国课后实验答案

实验一: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 课后实验题答案实验一 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 +⎡⎤=⎢⎥⎣⎦。

操作系统lab2实验报告

运行结果四实验体会本实验是关于物理内存管理主要从物理内存和建立页表两个方面设计实验首先了解如何发现系统中的物理内存然后了解如何建立对物理内存的初步管理最后了解页表相关的操作总体来说还是有难度尤其后面需要编程的内容需要很好理解各个文件的代码参考答案进行分析才能更好的理解实验内涵通过实验也能深入理解段页式内存管理机制对上课的内容有了更深理解体会
二、目的
1.理解基于段页式内存地址的转换机制;
2.理解页表的建立和使用方法;
3.理解物理内存的管理方法。
三、实验设计思想和练习题
练习0:填写已有实验
使用eclipse中的diff/merge工具将实验1的代码填入本实验中代码中有“LAB1”的注释相应部分。
练习1:实现first-fit连续物理内存分配算法(需要编程)
list_entry:双向链表指针,指向空闲的物理页;
nr_free:记录当前空闲页的个数的无符号整形变量。
(2)连续物理内存分配思路:
物理内存页管理器顺着双向链表进行搜索空闲内存区域,直到找到一个足够大的空闲区域,这是一种速度很快的算法,因为它尽可能少地搜索链表。如果空闲区域的大小和申请分配的大小正好一样,则把这个空闲区域分配出去,成功返回;否则将该空闲区分为两部分,一部分区域与申请分配的大小相等,把它分配出去,剩下的一部分区域形成新的空闲区。其释放内存的设计思路是把这块区域重新放回双向链表中。
(2)如果找到,获得指向分配的页,重新设置标志位,从空闲链表中删除此页;
(3)判断空闲块大小是否合适;
(4)如果合适,不操作;如果不合适,分割页块;
ห้องสมุดไป่ตู้(5)计算剩余空闲页个数,返回分配的页块地址。
B.具体实现:
注:
算法改进空间:第一页重置标志位操作被注释后,依然可以编译运行。

matlab2022实验2参考答案

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)%%这种方法较第一种解法缺点:在除法运算中,会产生误差,故此题应先将分母的多项式相乘后,再与分子部分的多项式进行运算。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

一般说来,如果你已知某个值是不能改变的(特别是作为函数参数的某个值),你就 应该把它声明为常量(const)。这样,即使在你编程时弄糊涂了甚至试图修改它的值,编 译器也能发现你的错误。 1.1.3 随机数发生函数 C++标准库里包含有一个 rand()函数,用于产生位于 0 至 RAND_MAX(由编译器 定义的整型常数)之间的随机数。但所产生的这些数并非真正意义上的随机数;他们只 是基于一个特定的“种子”数而生成的貌似随实则固定的数字序列。 为了确保我们不会总 是得到相同的随机数序列,我们通常采用当前时间作为种子数。下面举例说明之:
4.1 写一个返回两个整数和的求和函数。 同时写出求两个实数(double 类型)之和的等价 函数,函数命名为 sum. [4 分] 答案 (4 分) 程序如下: 1 int sum( const int x, const int y) { 2 return x + y; 3 } 4 5 double sum( const double x, const double y) { 6 return x + y; 7 }
1 2 3 4 5 6 7 } } return sum ; int sum ( const int numbers [] , const int numbersLen ) { int sum=0; for( int i=0; i< numbersLen ; ++ i) { sum+= numbers [i ];
Lab 2 及参考答案
1.1.1 默认参数 假定你有一个函数,其中有个参数,可能在调用时经常保持不变。例如,我们需要 这样一个函数,它将指定的信息显示 n 次。但大多数我们调用它的时候,只需让它将指 定信息显示出来一次就行了:
1. 2. 3. 4. 5. } } void printNTimes( char *msg,Байду номын сангаасint n) { for ( int i=0; i<n; ++i) { cout << msg;
4.2 请解释,为何在使用 4.1 节的求和函数 sum(1,10.0)时,会出现语法错误。(提示: 考虑在函数调用中参数类型的转换,比如,将第一个参数值 1 看作 double 类型,从而返回 double 类型的和;以及将两个参数均视作 int 型,从而返回 int 型的和)【1 分】
答案(1 分) 混用并匹配 int 和 double 类型,使编译器不能确定到底要调用哪个函数。编译器可 能将 1 强制转换成 double 类型并调用 double 版本的 sum 函数;也可以将 10.0 强制转 换为 int 类型,并调用整型版本的函数。 4.3 再写两个函数,以便完成任何 2 至 4 个整数的求和,函数名及参数写成 sum(num1, num2, ……)的形式。[2+2 分] 答案(2+2 分) 1 int sum( const int x, const int y, const int z) { 2 return x+y+z; 3 } 4 5 int sum( const int a, const int b, const int c, const int d) { 6 return a+b+c+d; 7 } 4.4 现在使用默认参数只写一个函数,便可实现任意 2 至 4 个整数之和的计算及返回。 如果你把这个(函数)定义和 4.3 节中的 3 参数函数放到同一个文件中, 并调用 sum(3,5,7), 会发生什么情况?为什么?【5+1 分】 答案(5+1 分) 1 int sum ( const int a, const int b, const int c=0, const int d=0){
1 2 3 4 5 6 void print( const long &x) { // Pass -by -reference avoids overhead // of copying large number cout << x; } int main() {
7 8 9 10 11 }
long x = 234923592; print(x); // We are guaranteed that x // will not be changed by this return 0;
将函数的一个参数声明为“int n = 1”, 就允许我们用 printNTimes("Some message") 来调用这个函数;编译器会自动插入 1 作为第二个参数。 一个函数中可以有多个默认参数: 1 void printNTimes( char *msg = "\n", int n = 1) { 2 for ( int i=0; i<n; ++i) { 3 cout << msg; 4 } 5} 这样,若需输出一新的行(换行),我们只需写出语句 printNTimes();实现函数调用 就行了。然而,C++不允许跳过参数,所以我们不可用调用语句 printNTimes(k)来实现 k 个换行。若非得这样做,则必须写成 printNTimes("\n", k); 1.1.2 常量参数 有些函数必须对其所带的某些参数视作常量, 这种情况下通常就将它们指定为常量 参数。如同常规的变量定义那样,这时我们就将函数的参数声明为 const: 1. void print( const int n) { 2. cout << n; 3. } 有时我们必须采用传址方法来向函数传递参数但又不希望在函数内部对参数的原 值作任何修改,此时定义常量参数尤其重要:
1 2 3 4 5 6 7 8 9 10 11 #include <iostream > #include <cstdlib > // C standard library // defines rand(), srand(), RAND_MAX #include <ctime> // C time functions -defines time() int main() { srand( time(0) ); // Set the seed; // time(0) returns current time as a number int randNum = rand(); std::cout << "A random number: " << randNum << endl; return 0; }
为了显示一次信息,我们可以将调用语句写成 printNTimes("Some message", 1); 但 C++有更好的办法,那就是允许在函数中定义默认参数,如下:
1. 2. 3. 4. 5. } } void printNTimes( char *msg, int n = 1 ) { for (int i = 0; i < n; ++ i ) { cout << msg;
2 3 } return a+ b+ c+ d;
如果给定的几个定义放在一起,编译器会提示编译错误,因为在缺省参数调用的情 况下,编译器无法确定是调用 3 参数函数还是调用 4 参数函数。 4.5 试写一个 sum 函数,可以实现任意个整数的求和。这个函数必须带有 2 个参数, 包含一个循环,并返回整型和数(提示:要用 2 个参数来表示任意个数的整数集合,你 应该采用哪种数据类型?)【5 分】 答案(5 分)
4.6 现在采用递归代替循环改写 4.5 节的函数。函数的署名不能改变。考虑一下指针的 算术运算会对你有所帮助。【8 分】 答案(8 分)
1 2 3
int sum( const int numbers[], const int numbersLen) { return numbersLen == 0 ? 0 : numbers[0] + sum(numbers + 1, numbersLen -1); }
相关文档
最新文档