资源描述:
《matlab平时作业共享.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、第一章习题13.请指出如下5个变量名中,哪些是合法的?abcd-2,xyz_3,ABCDefghabcd-2xyz_33chana变量ABCDefgh6.指令clear,clf,clc各有什么用处?Clear为清除工作空间变量指令,包括内存中的变量和数据;Clc为清命令窗口指令,而保留工作空间内容;Clf为清除当前图形指令。7.以下两种说法对吗?(1)“MATLAB的数值表达精度与其指令窗中的数据显示精度相同。”(2)MATLAB指令窗中显示的数值有效位数不超过7位。”解:(1)错误;(2)错误,该显示数值有效位可以在菜单中设置或通过vpa(X,n)进行设置。8.
2、想要在MATLAB中产生二维数组,下面哪些指令能实现目的?S=[1,2,3;4,5,6;7,8;9]S=[123;456;789]S=[1,2,3;4,5,6;7,8,9]解:S=[123;456;789];S=[1,2,3;4,5,6;7,8,9]。第三章习题31.要求在闭区间上产生具有10个等距采样点的一维数组。试用两种不同的指令实现。解:a=linspace(0,2*pi,10)a=0:2*pi/9:2*pi2.由指令rand('state',0),A=rand(3,5)生成二维数组A,试求该数组中所有大于0.5的元素的位置,分别求出它们的“全下标”和“单下标”
3、。解:rand('state',0),A=rand(3,5);L(A>0.5)=NAN程序如下:%注释:最终求出的该数组中所有大于0.5的元素的“全下标”信息存储于C中%%%%“单下标”信息存储于D中clear;ii=0;jj=0;x=1;C=zeros(15,2);D=zeros(1,15);rand('state',0),A=rand(3,5);forii=1:3forjj=1:5ifA(ii,jj)>0.5C(x,:)=[ii,jj];D(x)=(ii-1)*5+jj;x=x+1;endendendC,D结果如下:C=1115222425313233340000
4、00000000D=157910111213140000005.在时间区间[0,10]中,绘制曲线。要求分别采取“标量循环运算法”和“数组运算法”编写两段程序绘图。解:(1)标量循环运算法plot(t,y)fort=0:0.01:10y=1-exp(-0.5*t)*cos(2*t);plot(t,y);grid;holdon;end如图:(2)数组运算法t=0:0.01:10;y=1-exp(-0.5*t).*cos(2*t);plot(t,y)如图:6.先运行clear,formatlong,rand('state',1),A=rand(3,3),然后根据A写出两个
5、矩阵:一个对角阵B,其相应元素由A的对角元素构成;另一个矩阵C,其对角元素全为0,而其余元素与对应的A阵元素相同。A={0.952782149656620.598158524172190.836819600676340.704062166775000.840743198113070.518703059724920.953877473592230.442818842235130.02220977857260}B={0.952782149656620000.840743198113070000.02220977857260}C={00.598158524172190.83
6、6819600676340.7040621667750000.518703059724920.953877473592230.442818842235130}第四章习题43.求函数的数值积分,并请采用符号计算尝试复算。解:(1)数值积分:clearsymsxx=0:0.1:pi;y=exp(sin(x).^3);>>trapz(x,y)ans=5.0954(2)符号计算:symsxIsym=vpa(int(exp(sin(x)^3),x,0,pi))Warning:Explicitintegralcouldnotbefound.>Insym.intat58Isym=5
7、.13701355671761418934690897508764.用quad求取的数值积分,并保证积分的绝对精度为。解:Y=quad('exp(-abs(x)).*abs(sin(x))',-5*pi,1.7*pi,1e-9)Y=1.087849437547795.求函数在区间中的最小值点。解:clearsymsxy=((sin(5*x))^2)*(exp(0.06*x^2))-1.5*x*cos(2*x)+1.8*abs(x+0.5);yd=diff(y,x);a=solve(yd)yd_a=vpa(subs(yd,x,xs0),6)y_a=vp