河南城建学院MATLAB上机实验答案.doc

该【河南城建学院MATLAB上机实验答案 】是由【泰山小桥流水】上传分享，文档一共【19】页，该文档可以免费在线阅读，需要了解更多关于【河南城建学院MATLAB上机实验答案 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。百度文库 -让每个人平等地提升自我一 熟悉Matlab工作环境1、熟悉Matlab的5个基本窗口思考题:1)变量如何声明,变量名须遵守什么规则、是否区分大小写。答:变量一般不需事先对变量的数据类型进行声明,系统会依据变量被赋值的类型自动进行类型识别,也就是说变量可以直接赋值而不用提前声明。变量名要遵守以下几条规则:变量名必须以字母开头,只能由字母、数字或下划线组成。变量名区分大小写。变量名不能超过63个字符。关键字不能作为变量名。最好不要用特殊常量作为变量名。2)试说明分号、逗号、冒号的用法。分号:分隔不想显示计算结果的各语句;矩阵行与行的分隔符。逗号:分隔欲显示计算结果的各语句;变量分隔符;矩阵一行中各元素间的分隔符。冒号:用于生成一维数值数组;表示一维数组的全部元素或多维数组某一维的全部元素。3)linspace()称为“线性等分”函数,说明它的用法。LINSPACELinearlyspacedvector. 线性等分函数LINSPACE(X1,X2),X2为末元素平均生成 100个元素的行向量。LINSPACE(X1,X2,N)<2,,X2为末元素平均生成 n个元素的行向量。如果 n<2,返回X2。ClasssupportforinputsX1,X2:float:double,single数据类型:单精度、双精度浮点型。(4)说明函数ones()、zeros()、eye()的用法。ones()生成全1矩阵。zeros()生成全0矩阵。eye()生成单位矩阵。1百度文库 -让每个人平等地提升自我2、Matlab的数值显示格式思考题:(1)3次执行exist(’pi’)的结果一样吗?如果不一样,试解释为什么?>>pi>>pi=0;>>clearans=>>exist('pi')>>exist('pi')ans=ans=>>sin(pi);15>>exist('pi')>>pi>>pians=pi=ans=50答:3次执行的结果不一样。exist()函数是返回变量搜索顺序的一个函数。在第一次执行时返回5代表变量pi是由Matlab构建的变量。在第二次执行时已经通过赋值语句定义了变量pi,返回1代表pi是工作空间变量。第三次执行前清除了工作空间,此时 pi为系统默认常量,和第一次执行时性质一样,所以又返回 5。(2)圆周率pi是系统默认常量,为什么会被改变为 0。pi=0 为赋值语句,此时 pi不再是系统默认常量,而是定义的变量了。MATLAB语言基础1、向量的生成和运算练****使用logspace()创建1~4π的有10个元素的行向量。A=logspace(0,,10)A=2、矩阵的创建、引用和运算(1)矩阵的创建和引用练****创建以下矩阵:A为3×4的全1矩阵、B为3×3的0矩阵、C为3×3的单位矩阵、D为3×3的魔方阵、E由C和D纵向拼接而成、F抽取E的2~5行元素生成、G由F经变形为3×4的矩阵而得、以G为子矩阵用复制函数生成6×8的大矩阵H。>>A=ones(3,4),B=zeros(3,3),C=eye(3,3),D=magic(3)A=11**********B=0**********D=000816C=357100492E=[C;D],F=E(2:5,:),G=reshape(F,3,4)E=100010001816357492F=0**********G=031101568007>>H=repmat(G,2)H=0311031101560156800780070311031130 1 5 6 0 1 5 68 0 0 7 8 0 0 72)矩阵运算练****1)用矩阵除法求下列方程组的解 x=>>A=[634;-257;8-1-3],B=[3;-4;-7]A=B=6343-257-48-1-3-7x=A\Bx=2)求矩阵的秩;>>r=rank(A)r=33)求矩阵的特征值与特征向量>>[X,Lamda]=eig(A)X=Lamda=0000004)矩阵的乘幂(平方)与开方A^2ans=622933434 12 626 22 34A1=sqrtm(A)A1=--+ - -++5)矩阵的指数与对数(以 e为底)Ae=expm(A)Ae=+004*>>Ael=logm(A)Ael=--+ - -++6)矩阵的提取(取右上三角)与翻转(逆时针转 90度)>>a=triu(A)>>a1=rot90(A)a=a1=63447-305735-100-36-283、多维数组的创建及运算练****创建三维数组A,第一页为,第二页为,第三页为。然后用reshape函数重排为数组B,B为3行、2列、2页。>>a=[13;42],b=[12;21],c=[35;71]5>>A=cat(3,a,b,c)A(:,:,1)=A(:,:,2)=A(:,:,3)=131235422171>>B=reshape(A,3,2,2)B(:,:,1)=B(:,:,2)=122741153231三Matlab 数值运算1、多项式运算练****求 的商及余多项式。p1=conv([101],conv([13],[11]))p1=1 4 4 4 3[qr]=deconv(p1,[1021])q=1 4r=0 0 2 -5 -12、多形式插值和拟合有一组实验数据如附表 1-1所示。请分别用拟合(二阶至三阶)和插值(线性和三次样条)的方法来估测X=时Y的值X**********Y163270142260436682101014321960x=1:10;y=[163270142260436682101014321960];p1=polyfit(x,y,1)p1= >>y1=polyval(p1,6y1= p3=+003>> y3=p2=polyfit(x,y,2),y2=polyval(p2, +003p2= >>y4=interp1(x,y,y4=y2= 1696+003 >>y5=spline(x,y,y5=>> 1682p3=polyfit(x,y,3),y3=polyval(p3,3****题(1)用函数roots求方程 的根>>roots([1-1-1])ans=(2) ,在n个节点(n不要太大,如取 5~11)上用分段线性和三次样条插值方法,计算 m个插值点(m可取50~100)的函数值。通过数值和图形输出,将两种插值结果与精度进行比较。适当增加 n,再作比较。x=linspace(0,2*pi,8),y=sin(x)x=0y=0xi=linspace(0,2*pi,100);y0=sin(xi);y1=interp1(x,y,xi);y2=interp1(x,y,xi,'spline');>>plot(xi,y0,'*',xi,y1,'-.',xi,y2)-----10 1 2 3 4 5 6 7e1=y1-y0;e2=y2-y0;plot(xi,e1)----- 1 2 3 4 5 6 78>>plot(xi,e2)--- 1 2 3 4 5 6 7(3)大气压强p随高度x变化的理论公式为 ,为验证这一公式,测得某地大气压强随高度变化的一组数据如表所示。试用插值法和拟合法进行计算并绘图,看那种方法较为合理,且总误差最小。高度/m0300600100015002000压强/Pa插值法:>>x=[0300600100015002000];p=[ ];xi=linspace(0,2000);p0=*exp(-(xi+500)/7756);p1=interp1(x,p,xi,'spline');plot(xi,p0,'*',xi,p1)e1=p1-p0;e=sum(e1.^2)e=拟合法:>>x=[0300600100015002000];p=[ ];9P=log10(p)P=p1=polyfit(x,P,1)p1=b=p1(1)/,a=10.^p1(2)b=a=xi=linspace(0,2000);p0=*exp(-(xi+500)/7756);p2=polyval(p1,xi);P2=10.^p2;e2=P2-p0;e=sum(e2.^2)e=10