Matlab上机实验答案

Matlab上机实验答案实验一MATLAB运算基础1.先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。(1)0122sin851ze>>z1=2*sin(85*pi/180)/(1+exp(2))z1=0.2375(2)221ln(1)2zxx，其中2120.455ix>>x=[21+2i;-0.455];>>z2=1/2*log(x+sqrt(1+x^2))z2=0.7114-0.0253i0.8968+0.3658i0.2139+0.9343i1.1541-0.0044i(3)0.30.330.3sin(0.3)ln,3.0,2.9,,2.9,3.022aaeeazaa>>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)(>>z33=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)可以验证z3==z33，是否都为1）z3=Columns1through50.7388+3.1416i0.7696+3.1416i0.7871+3.1416i0.7913+3.1416i0.7822+3.1416iColumns6through100.7602+3.1416i0.7254+3.1416i0.6784+3.1416i0.6196+3.1416i0.5496+3.1416iColumns11through150.4688+3.1416i0.3780+3.1416i0.2775+3.1416i0.1680+3.1416i0.0497+3.1416iColumns16through20-0.0771+3.1416i-0.2124+3.1416i-0.3566+3.1416i-0.5104+3.1416i-0.6752+3.1416iColumns21through25-0.8536+3.1416i-1.0497+3.1416i-1.2701+3.1416i-1.5271+3.1416i-1.8436+3.1416iColumns26through30-2.2727+3.1416i-2.9837+3.1416i-37.0245+0.0000i-3.0017+0.0000i-2.3085+0.0000iColumns31through35-1.8971+0.0000i-1.5978+0.0000i-1.3575+0.0000i-1.1531+0.0000i-0.9723+0.0000iColumns36through40-0.8083+0.0000i-0.6567+0.0000i-0.5151+0.0000i-0.3819+0.0000i-0.2561+0.0000iColumns41through45-0.1374+0.0000i-0.0255+0.0000i0.0792+0.0000i0.1766+0.0000i0.2663+0.0000iColumns46through500.3478+0.0000i0.4206+0.0000i0.4841+0.0000i0.5379+0.0000i0.5815+0.0000iColumns51through550.6145+0.0000i0.6366+0.0000i0.6474+0.0000i0.6470+0.0000i0.6351+0.0000iColumns56through600.6119+0.0000i0.5777+0.0000i0.5327+0.0000i0.4774+0.0000i0.4126+0.0000iColumn610.3388+0.0000i(4)2242011122123ttzttttt，其中t=0:0.5:2.5>>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)z4=00.250001.25001.00002.25002.已知：1234413134787,2033657327AB求下列表达式的值：(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]>>A=[1234-4;34787;3657];>>B=[13-1;203;3-27];>>A+6*Bans=1852-10467105215349>>I=eye(3);>>A-B+Ians=1231-3328840671(2)>>A*Bans=684462309-72596154-5241>>A.*Bans=1210246802619-13049(3)>>A^3ans=372262338244860424737014918860076678688454142118820>>A.^3ans=172839304-643930434365850327274625343(4)>>A/Bans=16.4000-13.60007.600035.8000-76.200050.200067.0000-134.000068.0000>>B\Aans=109.4000-131.2000322.8000-53.000085.0000-171.0000-61.600089.8000-186.2000(5)>>[A,B]ans=1234-413-13478720336573-27>>[A([1,3],:);B^2]ans=1234-436574511101920-5403.设有矩阵A和B1234530166789101769,111213141502341617181920970212223242541311AB(1)求它们的乘积C。(2)将矩阵C的右下角3×2子矩阵赋给D。(3)查看MATLAB工作空间的使用情况。>>A=(reshape(1:25,5,5))';>>B=[3016;17-69;023-4;970;41311];>>C=A*BC=9315077258335237423520397588705557753890717>>D=C(3:5,2:3)D=520397705557890717>>whosNameSizeBytesClassAttributesA5x5200doubleB5x3120doubleC5x3120doubleD3x248double4.完成下列操作：(1)求[100,999]之间能被21整除的数的个数。(2)建立一个字符串向量，删除其中的大写字母。>>n=100:999;>>l=find(rem(n,21)==0);>>length(l)ans=43>>ch='aegbBOIEG0je23RGnc';>>wz=find(ch>='A'&ch<='Z');>>ch(wz)=[]ch=aegb0je23nc实验二MATLAB矩阵分析与处理1.设有分块矩阵33322322ERAOS，其中E、R、O、S分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22ERRSAOS。>>E=eye(3);>>R=rand(3,2);>>O=zeros(2,3);>>S=diag(1:2);>>A=[E,R;O,S]A=1.0000000.81470.913401.000000.90580.6324001.00000.12700.09750001.0000000002.0000>>A^2ans=1.0000001.62942.740101.000001.81161.8971001.00000.25400.29260001.0000000004.0000>>[E,R+R*S;O,S^2]ans=1.0000001.62942.740101.000001.81161.8971001.00000.25400.29260001.0000000004.0000>>A^2==[E,R+R*S;O,S^2]ans=1111111111111111111111111通过验证，矩阵22ERRSAOS成立。2.产生5阶希尔伯特矩阵H和5阶帕斯卡矩阵P，且求其行列式的值Hh和Hp以及它们的条件数Th和Tp，判断哪个矩阵性能更好。为什么？>>H=hilb(5)H=1.00000.50000.33330.25000.20000.50000.33330.25000.20000.16670.33330.25000.20000.16670.14290.25000.20000.16670.14290.12500.20000.16670.14290.12500.1111>>P=pascal(5)P=111111234513610151410203515153570>>Hh=det(H)Hh=3.7493e-12>>Hp=det(P)Hp=1>>Th=cond(H)Th=4.7661e+05>>Tp=cond(P)Tp=8.5175e+03答：5阶帕斯卡矩阵P的性能好。矩阵的性能是由条件数决定的，条件数越接近于1其性能就越好。由上机操作求得Th=4.7661e+005，Tp=8.5175e+003。Tp的值更接近于1则其性能要好。所以5阶帕斯卡矩阵P的性能好。3.建立一个5×5矩阵，求它的行列式值、迹、秩和范数。>>A=rand(5)A=0.27850.97060.42180.03570.74310.54690.95720.91570.84910.39220.95750.48540.79220.93400.65550.96490.80030.95950.67870.17120.15760.14190.65570.75770.7060>>det(A)ans=-0.1322>>trace(A)ans=3.4127>>rank(A)ans=5>>norm(A)ans=3.27594.已知2961820512885A求A的特征值及特征向量，并分析其数学意义。>>A=[-29,6,18;20,5,12;-8,8,5]A=-2961820512-885>>[V,D]=eig(A)V=0.71300.28030.2733-0.6084-0.78670.87250.34870.55010.4050D=-25.3169000-10.518200016.8351在数学上，特别是线性代数中，对于一个给定的线性变换，它的特征向量（本征向量或称正规正交向量）是这样一个非零的向量v：当v经过这个线性变换的作用之后，得到的新向量（长度也许改变）仍然与原来的v保持在同一条线上。一个特征向量的长度在该线性变换下缩放的比例称为其特征值（本征值）。如果特征值为正，则表示v在经过线性变换的作用后方向也不变；如果特征值为负，说明方向会反转；如果特征值为0，则是表示缩回零点。但无论怎样，仍在同一条直线上。5.下面是一个线性方程组：1231112340.951110.673450.52111456xxx(1)求方程的解。(2)将方程右边向量元素b3改为0.53再求解，并比较b3的变化和解的相对变化。(3)计算系数矩阵A的条件数并分析结论。>>formatrat%用分数格式显示>>A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]A=1/21/31/41/31/41/51/41/51/6>>format%恢复默认格式>>b=[0.95;0.67;0.52]b=0.95000.67000.5200>>X=A\bX=1.20000.60000.6000>>b2=[0.95;0.67;0.53]b2=0.95000.67000.5300>>X2=A\b2X2=3.0000-6.60006.6000>>D=cond(A)D=1.3533e+03矩阵的条件数决定矩阵的性能，条件数越接近于1其性能越好，系数矩阵A的条件数为1.3533e+003，和1相差很大，则其性能不好。因此b矩阵个别元素的微小变动，对方程的解（X值）影响很大。6.建立A矩阵，试比较sqrtm(A)和sqrt(A)，分析它们的区别。实验三选择结构程序设计1.求分段函数的值。2226035605231xxxxyxxxxxxx且且及其他用if语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y值。程序设计：clear;x=[-5.0-3.01.02.02.53.05.0];ifx<0&x~=-3y=x.^2+x-6;elseifx>=0&x<5&x~=2&x~=3y=x.^2-5*x+6;elsey=x.^2-x-1;endy运行结果：y=29.000011.0000-1.00001.00002.75005.000019.00002.输入一个百分制成绩， 要求 对教师党员的评价套管和固井爆破片与爆破装置仓库管理基本要求三甲医院都需要复审吗 输出成绩等级A、B、C、D、E。其中90分~100分为A，80分~89分为B，79分~79分为C，60分~69分为D，60分以下为E。要求：(1)分别用if语句和switch语句实现。程序设计：（一）用if语句clear;score=input('请输入百分制成绩：');a=fix(score/10);ifa==9|a==10grade='A';elseifa==8grade='B';elseifa==7grade='C';elseifa==6grade='D';elsegrade='E';endgrade（二）用if语句clear;score=input('请输入百分制成绩：');ifscore>=90&score<=100grade='A';elseifscore>=80&score<90grade='B';elseifscore>=70&score<80grade='C';elseifscore>=60&score<70grade='D';elsegrade='E';endgrade（三）用switch语句clear;score=input('请输入百分制成绩：');switchfix(score/10)case{9,10}grade='A';case8grade='B';case7grade='C';case6grade='D';otherwisegrade='E';endgrade(2)输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。clear;score=input('请输入百分制成绩：');ifscore>=0&score<=100%判断成绩合理性switchfix(score/10)case{9,10}grade='A';case8grade='B';case7grade='C';case6grade='D';otherwisegrade='E';endgradeelsedisp('error')%成绩不合理时输出出错信息end3.硅谷公司员工的工资计算方法如下：(1)工作时数超过120小时者，超过部分加发15%。(2)工作时数低于60小时者，扣发700元。(3)其余按每小时84元计发。试编程按输入的工号和该号员工的工时数，计算应发工资。程序设计：clear;num=input('请输入员工工号：');time=input('请输入员工工时数：');iftime<60wages=time*84-700;elseiftime>120wages=84*120+(time-120)*84*(1+0.15);elsewages=time*84;enddisp=(['员工工号：'num2str(num)'应发工资：'num2str(wages)])运行结果：请输入员工工号：4请输入员工工时数：35disp=员工工号：4应发工资：22404.设计程序，完成两位数的加、减、乘、除四则运算，即产生两个两位随机整数，再输入一个运算符号，做相应的运算，并显示相应的结果。解：程序设计：clear;a=floor(rand(1)*90+10)b=floor(rand(1)*90+10)c=input('请输入运算符号：','s');ifc=='+'s=a+belseifc=='-'s=a-belseifc=='*'s=a*belseifc=='\'s=a\belsedisp('error')end运行结果：a=30b=77请输入运算符号：-s=-475.建立5×6矩阵，要求输出矩阵第n行元素。当n值超过矩阵的行数时，自动转为输出矩阵最后一行元素，并给出出错信息。程序设计：clear;A=rand(5,6)n=input('n=?')ifn>=1&n<=5B=A(n,1:6);elsedisp('erro!')B=A(5,1:6);endB运行结果：A=0.77920.01190.52850.68920.91330.07820.93400.33710.16560.74820.15240.44270.12990.16220.60200.45050.82580.10670.56880.79430.26300.08380.53830.96190.46940.31120.65410.22900.99610.0046n=?3n=3B=0.12990.16220.60200.45050.82580.1067实验四循环结构程序设计1.根据2222211116123n，求π的近似值。当n分别取100、1000、10000时，结果是多少？要求：分别用循环结构和向量运算（使用sum函数）来实现。程序设计（一）：clear;forn=[100,1000,10000]sum=0;forx=1:nsum=sum+1/(x.^2);endresult=sqrt(6*sum)end运行结果：result=3.1321result=3.1406result=3.1415程序设计（二）：clear;forn=[100,1000,10000]x=1:n;result=sqrt(6*sum(1./x./x))end运行结果：result=3.1321result=3.1406result=3.14152.根据11113521yn，求：(1)y<3时的最大n值。(2)与(1)的n值对应的y值。程序设计：y=0;n=1;whiley<3y=y+1/(2*n-1);n=n+1;endn=n-1;y=y-1/(2*n-1)n=n-1运行结果：y=2.9944n=56【%验证n=56时的y值clearall;clc;n=56;i=1:n;f=1./(2*i-1);y=sum(f)】3.考虑以下迭代 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 ：1nnaxbx其中a、b为正的学数。(1)编写程序求迭代的结果，迭代的终止条件为|xn+1-xn|≤10-5，迭代初值x0=1.0，迭代次数不超过500次。(2)如果迭代过程收敛于r，那么r的准确值是242bba，当(a,b)的值取(1,1)、(8,3)、(10,0.1)时，分别对迭代结果和准确值进行比较。程序设计：clear;a=input('a=?');b=input('b=?');x0=1;x1=a/(b+x0);y=abs(x1-x0);n=1;whiley>10^(-5)&n<=499x0=x1;x1=a/(b+x0);y=abs(x1-x0);n=n+1;endx1%比较收敛准确值和迭代结果r1=(-b+sqrt(b^2+4*a))/2r2=(-b-sqrt(b^2+4*a))/2运行结果：a=?1b=?1x1=0.6180r1=0.6180r2=-1.6180a=?8b=?3x1=1.7016r1=1.7016r2=-4.7016a=?10b=?0.1x1=3.1127r1=3.1127r2=-3.21274.已知12312311021323nnnnfnfnfnffffn求f1~f100中：(1)最大值、最小值、各数之和。(2)正数、零、负数的个数。程序设计：clear;f1=1;f2=0;f3=1;f=[f1,f2,f3];forn=4:100f4=f3-2*f2+f1;f=[ff4];f1=f2;f2=f3;f3=f4;endmaxf=max(f)minf=min(f)sumf=sum(f)positivenum=length(find(f>0))negtivenum=length(find(f<0))zeronum=length(find(f==0))运行结果：maxf=4.3776e+11minf=-8.9941e+11sumf=-7.4275e+11positivenum=49negtivenum=49zeronum=25.若两个连续自然数的乘积减1是素数，则称这两个边疆自然数是亲密数对，该素数是亲密素数。例如，2×3-1=5，由于5是素数，所以2和3是亲密数，5是亲密素数。求[2,50]区间内：(1)亲密数对的对数。(2)与上述亲密数对对应的所有亲密素数之和。程序设计（一）：clearpair=[];form=2:49f=(m+1)*m-1;signf=0;%标记为0表示f为素数fork=2:sqrt(f)if(rem(f,k)==0)signf=1;%标记为1表示f不是素数endendifsignf==0pair=[pair;mm+1];%亲密数对矩阵endendnum=size(pair);susum=sum(pair(:,1).*pair(:,2)-1);disp(['亲密数对的对数',num2str(num(1))])disp(['亲密素数之和',num2str(susum)])运行结果：亲密数对的对数28亲密素数之和21066程序设计：clearsusum=0;num=0;form=2:49f=(m+1)*m-1;n=fix(sqrt(f));fork=2:nif(rem(f,k)==0)break;%f不是素数时跳出循环endendifk==nnum=num+1;%亲密数对个数susum=susum+f;%亲密素数和累加endenddisp(['亲密数对的对数',num2str(num)])disp(['亲密素数之和',num2str(susum)])运行结果：亲密数对的对数28亲密素数之和21066实验五函数文件1.定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用该函数文件。建立函数文件fushu.m：function[e,ln,s,c]=fushu(x)e=exp(x);ln=log(x);s=sin(x);c=cos(x);建立命令文件fushu2.mA=input('请输入任意一个复数A');[a,b,c,d]=fushu(A);%调用fushu.m的函数文件abcdfushu2.m运行结果：>>fushu2请输入任意一个复数A：4ia=-0.6536-0.7568ib=1.3863+1.5708ic=0.0000+27.2899id=27.30822.一物理系统可用下列方程组来表示：11121112220cossin0sin0cos000sin000cos1ammamgmNmNmg从键盘输入m1、m2和θ的值，求a1、a2、N1和N2的值。其中g取9.8，输入θ时以角度为单位。要求：定义一个求解线性方程组AX=B的函数文件，然后在命令文件中调用该函数文件。建立函数文件xxfcz.mfunction[a1,a2,N1,N2]=xxfcz(m1,m2,t)t=t*pi/180;A=[m1*cos(t)-m1-sin(t)0;m1*sin(t)0cos(t)0;0m2-sin(t)0;00-cos(t)1];g=9.8;B=[0;m1*g;0;m2*g];X=inv(A)*B;a1=X(1);a2=X(2);N1=X(3);N2=X(4);建立命令文件xxfcz2.mclearm1=input('请输入m1：');m2=input('请输入m2：');t=input('请输入角度t：');[a1,a2,N1,N2]=xxfcz(m1,m2,t);a1a2N1N2运行结果：>>xxfcz2请输入m1：5请输入m2：3请输入角度t：90a1=9.8000a2=3.7505e-16N1=1.1251e-15N2=29.40003.一个自然数是素数，且它的数字位置经过任意对换后仍为素数。例如13是绝对素数。试求所有两位绝对素数。要求：定义一个判断素数的函数文件。建立判断素数的函数文件sushu.mfunction[x]=sushu(m)%判断数m是否为素数forM=2:(m/2)ifrem(m,M)==0b=0;breakelseb=1;endendifb==1x=m;%判断数m为素数，并返回elsex=[];%数m不为素数endend建立命令文件sushu2.mclear;js=[];%初始绝对素数form=10:99x=sushu(m);ifxx=fix(x/10)+rem(x,10)*10;y=sushu(x);js=[js,y];%绝对素数elsecontinueendendjs执行命令文件sushu2.m>>sushu2js=113171137317379779法二：建立函数文件sushu3.mfunctionp=sushu3(m)b=m(end);p=m;fori=2:sqrt(b)n=find(rem(p,i)==0&p~=i);%将m中能被i整除，且不等于i的元素的下标找出，即np(n)=[];%将下标为n的元素剔除，其余的即为素数end建立命令文件sushu4.mclear;m=10:99;p=sushu3(m);%找出10到99内所有素数p=10*rem(p,10)+(p-rem(p,10))/10;%将p素数矩阵每个元素个位十位调换顺序p=sushu3(p)%再对调换顺序后的素数矩阵找出其中的素数执行命令文件sushu4.m>>sushu4p=1131711373173797794.设2411()(2)0.1(3)0.01fxxx，编写一个MATLAB函数文件fx.m，使得调用f(x)时，x可用矩阵代入，得出的f(x)为同阶矩阵。建立函数文件fx.mfunction[f]=fx(x)f=1./((x-2).^2+0.1)+1./((x-3).^4+0.01);end5.已知(40)(30)(20)fyff(1)当f(n)=n+10ln(n2+5)时，求y的值。(2)当f(n)=1×2+2×3+3×4+...+n×(n+1)时，求y的值。建立函数文件fn1.mfunction[fn]=fn1(n)fn=n+10*log(n.^2+5);end建立函数文件fn2.mfunction[fn]=fn2(n)fn=0;forx=1:nfn=fn+x*(x+1);end求（1）（2）中y的值：y1=fn1(40)/(fn1(30)+fn1(20));y2=fn2(40)/(fn2(30)+fn2(20));运行结果：>>yfny1=0.6390y2=1.7662实验六高层绘图操作1.设23sin0.5cos1xyxx，在x=0~2π区间取101点，绘制函数的曲线。>>x=linspace(0,2*pi,101);>>y=(0.5+3.*sin(x)./(1+x.^2)).*cos(x);>>plot(x,y);2.已知y1=x2，y2=cos(2x)，y3=y1×y2，完成下列操作：(1)在同一坐标系下用不同的颜色和线型绘制三条曲线。>>x=linspace(0,2*pi,101);>>y1=x.^2;>>y2=cos(2*x);>>y3=y1.*y2;>>plot(x,y1,'k:',x,y2,'b--',x,y3,'m-.');(2)以子图形式绘制三条曲线。x=linspace(0,2*pi,101);y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(3,1,1);plot(x,y1);title('y1');subplot(3,1,2);plot(x,y2);title('y2');subplot(3,1,3);plot(x,y3);title('y3');(3)分别用条形图、阶梯图、杆图和填充图绘制三条曲线。subplot(3,4,1);bar(x,y1,'y');title('bar(x,y1,"y")');subplot(3,4,2);stairs(x,y1,'m');title('stairs(x,y1,"m")');subplot(3,4,3);stem(x,y1,'g');title('stem(x,y1,"g")');subplot(3,4,4);fill(x,y1,'c');title('fill(x,y1,"c")');subplot(3,4,5);bar(x,y2,'y');title('bar(x,y2,"y")');subplot(3,4,6);stairs(x,y2,'m');title('stairs(x,y2,"m")');subplot(3,4,7);stem(x,y2,'g');title('stem(x,y2,"g")');subplot(3,4,8);fill(x,y2,'c');title('fill(x,y2,"c")');subplot(3,4,9);bar(x,y3,'y');title('bar(x,y3,"y")');subplot(3,4,10);stairs(x,y3,'m');title('stairs(x,y3,"m")');subplot(3,4,11);stem(x,y3,'g');title('stem(x,y3,"g")');subplot(3,4,12);fill(x,y3,'c');title('fill(x,y3,"c")');3.已知2201ln(1)02xxeyxxx在-5≤x≤5区间绘制函数曲线。x=-5:0.1:5;ifx<=0y=(x+sqrt(pi))./(exp(2));elsey=1/2*(log(x+sqrt(1+x.^2)));endplot(x,y)4.绘制极坐标曲线ρ=asin(b+nθ)，并分析参数a、b、n对曲线形状的影响。a=input('a=');b=input('b=');n=input('n=');theta=0:0.01:2*pi;rho=a*sin(b+n*theta);polar(theta,rho,'k')运行结果：a=1b=1n=15.绘制函数的曲线图和等高线。224coscosxyzxye其中x的21个值均匀分布[-5，5]范围，y的31个值均匀分布在[0，10]，要求使用subplot(2,1,1)和subplot(2,1,2)将产生的曲面图和等高线图画在同一个窗口上。x=-5:0.5:5;y=0:1/3:10;[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);subplot(2,1,2);contour3(x,y,z);6.绘制曲面图形，并进行插值着色处理。coscos3cossin0,022sinxstyststzss=0:0.05:pi/2;t=0:0.05:3*pi/2;[s,t]=meshgrid(s,t);[x,y,z]=peaks(30);x=cos(s).*cos(t);y=cos(s).*sin(t);z=sin(s);surf(x,y,z);shadingflat;